Decomposing Neural Representational Patterns of Discriminatory and Hedonic Information during Somatosensory Stimulation

Abstract The ability to interrogate specific representations in the brain, determining how, and where, difference sources of information are instantiated can provide invaluable insight into neural functioning. Pattern component modeling (PCM) is a recent analytic technique for human neuroimaging that allows the decomposition of representational patterns in brain into contributing subcomponents. In the current study, we present a novel PCM variant that tracks the contribution of prespecified representational patterns to brain representation across areas, thus allowing hypothesis-guided employment of the technique. We apply this technique to investigate the contributions of hedonic and nonhedonic information to the neural representation of tactile experience. We applied aversive pressure (AP) and appetitive brush (AB) to stimulate distinct peripheral nerve pathways for tactile information (C-/CT-fibers, respectively) while patients underwent functional magnetic resonance imaging (fMRI) scanning. We performed representational similarity analyses (RSAs) with pattern component modeling to dissociate how discriminatory versus hedonic tactile information contributes to population code representations in the human brain. Results demonstrated that information about appetitive and aversive tactile sensation is represented separately from nonhedonic tactile information across cortical structures. This also demonstrates the potential of new hypothesis-guided PCM variants to help delineate how information is instantiated in the brain.


Introduction
One of the holy grails of cognitive neuroscience is to be able to "read out" the content of cognition from nothing but patterns of brain activation. Whereas multivariate approaches that emphasize "decoding" are effective at classification or mapping distance in representational space, current efforts are focusing on how to model the specific content being instantiated by patterns of BOLD response across voxels (Kriegeskorte et al., 2008;Kriegeskorte and Kievit, 2013;Visser et al., 2013Visser et al., , 2015Visser et al., , 2016. Here, we present one such approach, focusing on using a variant of pattern component modeling (PCM; Diedrichsen et al., 2011Diedrichsen et al., , 2018Kriegeskorte and Kievit, 2013) that fits combinations of prespecified conceptual similarity patterns to observed neural representational patterns to investigate instantiations of sensory and hedonic tactile information in the CNS.
Sensory experiences, such as the embrace of a loved one or the pain of a stubbed toe, can be broken down into two central components: discrimination of the sensory information and the associated hedonic response. The information processed by sensory systems is typically viewed as objective, forming representations of a tangible external environment that are not yet colored by emotional evaluations. Yet, in the somatosensory system, there is strong evidence that hedonic information (good vs bad) is coded by peripheral afferents before any cortical processing (Iggo, 1959(Iggo, , 1960Vallbo et al., 1999). This indicates that there are pathways of tactile sensation that signal pleasant and painful emotional content, rather than "pure" sensory experience (Craig, 2015;Kryklywy et al., 2020). Yet while discriminatory and hedonic information are clearly dissociable in the peripheral nervous system (PNS), investigating the dissociability versus integration of these signals as instantiated in the CNS remains a challenge, in part because of the limited analytic techniques for parcellating overlapping patterns of representation in the brain. The current work presents a novel approach to the analysis of human neuroimaging data that decomposes multivariate representational patterns (described via representational similarity analyses; RSAs) into predefined representational patterns of interest (POIs) that model the information content instantiated by the BOLD response. Specifically, in the current study, we interrogated BOLD patterns evoked by aversive pressure (AP) and pleasant caress to examine instantiation of hedonic and discriminative somatosensory information in the CNS.
The somatosensory system contains multiple functional subsystems, with specific peripheral nerves serving as labeled lines for information traveling into the CNS (McGlone and Reilly, 2010;McGlone et al., 2014). Fast large-diameter myelinated afferent fibers (A-fibers) carry signals pertaining to sensory discrimination, while slower unmyelinated fiber pathways [C-fibers, Qiu et al., 2006; C-tactile (CT) fibers, Löken et al., 2009;Marshall et al., 2019] carry signals of the hedonic value of the sensation (pain and pleasure, respectively). A-fibers predominantly convey information about the timing and location of cutaneous sensory stimulation (McGlone and Reilly, 2010), with some fibers specialized for nociception (Nagi et al., 2019). Additional small-diameter unmyelinated afferent pathways support the hedonic response to touch, conveying information about affective aspects of aversive touch and nociception. At the cortical level, primary somatosensory cortex (S1) is the dominant entry point for information carried along myelinated cutaneous pathways. The extent to which hedonic and nonhedonic peripheral pathways remain segregated on their entry into cortical areas, including S1, remains unclear, in part because of an inability to differentiate between multiple distinct representations patterns versus a single integrated representational pattern instantiated in a single brain region.
There is evidence that integration of hedonic information into discriminatory touch representations in these early sensory structures occurs through centrally-mediated appraisal of pain and pleasure (Bushnell et al., 1999;Gazzola et al., 2012). This observation is consistent with conventional views positing that modulation of sensory information by emotion relies on re-entrant projections from higher order structures assessing hedonic value [e.g., prefrontal cortices (PFCs), insula (Ins), amygdala] to the sensory cortices (Pessoa and Adolphs, 2010;Rolls, 2019). Yet, evidence for peripheral labeling of affective information suggests that not all affective modulation of sensory signals is the result of central feedback (Qiu et al., 2006). Considerable evidence exists to support a neural bases of pain-coding in the periphery (McGlone and Reilly, 2010;Nagi et al., 2019). Anatomical projection studies in primates indicate that information carried along unmyelinated C-fiber pathways does not project to the entirety of S1, as observed for A-fibers. Rather, it projects to an anterior region of S1 (insula-adjacent area 3a; Whitsel et al., 2009;Vierck et al., 2013), with additional direct projections to the insula, anterior cingulate cortex (ACC), and PFC (Baumgärtner et al., 2006;Qiu et al., 2006). Similarly, recently identified pathways, C-tactile (CT) fibers, have been shown to carry information about caress or pleasant touch (Löken et al., 2009;Marshall et al., 2019). These fibers originate from mechanoreceptors located in hairy skin rather than the glabrous (i.e., hairless) skin of the palms, where previous research had focused (McGlone and Reilly, 2010;McGlone et al., 2014;Marshall et al., 2019). These CT-fiber afferents respond preferentially to touch that is subjectively perceived as pleasant caress (Olausson et al., 2002;Löken et al., 2009;Croy et al., 2016). Thus, in the cutaneous system there may be distinct parallel representations for tactile stimulation beginning from the point of contact, carried through hedonic labeled lines that independently inform the experience of hedonic value.
Previous functional magnetic resonance imaging (fMRI) studies examining neural substrates of affective-tactile processing have investigated either C-or CT-fiber pathways, but not both. The independent examination of Cand CT-fibers does not allow for the dissociation between these two distinct systems and is unable to discriminate valence-specific hedonic information from general tactile salience and arousal. Moreover, these studies have also relied predominantly on univariate statistical approaches (Olausson et al., 2002;Löken et al., 2009;McGlone et al., 2014). Univariate approaches have limited ability to discriminate specific information that is represented within a region, particularly perceptual and hedonic information represented by different sensory systems (Chikazoe et al., 2014;. By contrast, multivariate analyses, including representational similarity analyses (Kriegeskorte et al., 2008) allow the examination of population-based neuronal coding in a multidimensional representational space. When further analyzed through pattern component modeling (PCM; Kriegeskorte and Kievit, 2013;Diedrichsen et al., 2018), this representational space can be decomposed into weighted subcomponents of experience (Diedrichsen et al., 2018), thus representing neural activity as an integration of multiple heterogeneous sets of overlapping neural representations. In the present study, we implemented an innovative analytic approach that uses theory-guided components to perform PCM. We generated combinations of predefined similarity matrices in an approach that is a variation of previously described applications of a priori conceptual models of representational dissimilarity matrices (mRDM; Kriegeskorte et al., 2008;Kriegeskorte and Kievit, 2013;Popal et al., 2019). We refer to our hypothesized similarity models as patterns of interest (POIs), each of which characterizes a correlation pattern that would be expected given perfect representation of a single vector of information, e.g., painful touch. Whereas previous approaches fit single hypothesized conceptual similarity/dissimilarity matrices to the observed data, the current POI  looks for the simultaneous best fit of combinations of multiple models to observed representational patterns extracted from predefined regions of interest (ROIs).
Functional neuroimaging data were collected while participants received aversive pressure on the thumb or appetitive caress stimulation on the forearm (Fig. 1A) and viewed images of faces with neutral expressions. RSA conducted on the BOLD signal identified similarity/dissimilarity between tactile conditions that not only represent distinct discriminatory patterns but also putatively stimulate distinct fiber pathways for aversive and appetitive tactile Figure 1. A, Experimental time course. Participants completed tactile visual conditioning tasks. Only data collected during CS-US paired blocked will be presented. B, Representational similarity analyses (RSA). RSA was conducted correlating all experimental trials independently. Resultant Pearson correlation coefficients were averaged across conditions (removing autocorrelation) to create a 6 Â 6 condition similarity matrix comparing all conditions of interest. C, Pattern of interest (POI) identification. A greedy best-first search algorithm (GBFS) with Bayesian information criterion (BIC) was used to identify the combination of POIs that best predicted experimental data. Level-1 search independently fit 13 unique POIs to observed representation similarity. The POI with optimal fit was combined with all remaining models (Level 2) to determine the combination of POI represented in the ROI. This process was done iteratively until the addition of none of the POIs improved model fit. D, POI weighting. Multiple regression allowed extraction of b value indicating the representational strength for each contributing POI. experience (Fig. 1B). POIs were created for task-relevant information content and included specific aspects of tactile and emotional experience. Bayesian information criterion (BIC) analyses were then performed to characterize the combination of POIs that best predicted observed similarity patterns of neural activity in each ROI. This allowed us to identify and weigh dissociable representations of discriminative and hedonic tactile signals, revealing potential C-and CT-fiber pathways projection. Critically, these analyses allow us to fit multiple overlapping pattern components simultaneously and identify the best fitting combination, rather than fitting each component alone and comparing the representational strength of each component in isolation. We believe that this feature affords a nuanced and faithful glimpse into how a brain region may hold and process multiple sources of information.
We predicted that modeled representations of the hedonic components of the tactile stimulation in frontotemporal cortices, including vmPFC, ACC, and insula would be distinct from representations in primary somatosensory cortices. This would demonstrate the dual coding of somatosensation. We expected that dissociable representational patterns for appetitive versus aversive tactile stimulation would be identified in the insula and vmPFC, as these regions may receive direct unprocessed information from hedonic-labeled lines. Patterns observed in the ACC were predicted to be most heavily weighted toward representation of aversive touch, consistent with this region's preferential activation in response to modality-general pain.

Participants
A total of 121 participants (x age = 21.1, SD = 2.8; 41F) were recruited from a pool of 488 participants at Cornell University who had completed a large-scale behavioral study examining individual differences in responsiveness to reward and punishment (Moore, 2017). Of the data collected for the present study, we were unable to complete preprocessing for 54 participants. Multiecho (ME) preprocessing (described below) failed to generate the required denoised datasets (i.e., medn files) for 25 participants, while we were unable to obtain convergence in the independent component analyses (ICA) for an additional 14 participants. Data files (e.g., imaging run, stimulus onset timings, motion correction files) were missing for 12 participants. Two participants were excluded because of motion artifacts and a final participant was excluded to nonstandardized data collection. Results from the remaining 67 participants are reported. All participants gave written, informed consent and had normal or corrected-to-normal vision. Participants were prescreened for a history of anxiety and depression as well as other psychopathology, epilepsy, and brain surgery. Prescreening was followed up in person by an additional interview to ensure inclusion criteria were met. As this study was conducted as part of larger research program, all participants provided saliva samples for genotyping, and fecal sample for microbiome analyses. The experiment was performed in accordance with the Institutional Review Board for Human Participants at Cornell University.

Stimuli and apparatus
Photographs of three male and three female faces with neutral expressions were chosen from the Karolinska directed emotional faces picture set (Goeleven et al., 2008). These face images were used as conditioned stimuli (CS) in two classical conditioning paradigms, one aversive and one appetitive. Analysis here focuses only on the trials that included tactile stimulation = paired with the faces, as effects of conditioning are reported elsewhere (Ehlers et al., 2021). Unconditioned stimuli (US) consisted of either aversive pressure delivered to the right thumb, or appetitive caress to the participant's left forearm. These tactile manipulations were aimed to maximally activate C-and C-tactile somatosensory afferents, respectively. Aversive pressure stimuli were delivered using a custom designed hydraulic device (Giesecke et al., 2004;López-Solà et al., 2010) capable of transmitting controlled pressure to 1-cm 2 surface placed on the subjects' right thumbnail. Applied pressure levels were individually calibrated for each participant before the experiment to ensure that the pressure intensity was experienced as aversive but not excessively painful. Light appetitive caresses lasting ;4 s were manually applied to the left forearm with a brush by a trained experimenter to maximally activate CT-fiber pathways (McGlone et al., 2014). Only participants who had previously demonstrated reliable affective responses to the tactile manipulations in the larger behavioral study (e.g., positive responses to caress; see below for details) were invited to participant in the scanning session, as this indicated likely recruitment of C-and CT-fiber pathways of the hedonic tactile manipulations (which can be subject to interparticipant variability).

Procedure
While undergoing functional MR scanning, participants completed two separate conditioning tasks (appetitive conditioning and aversive conditioning), each involving a series of tactile and visual pairings ( Fig. 1A; Visser et al., 2015). In each task, participants completed seven CSonly blocks interleaved with six CS-US paired blocks. Single blocks of either the CS-only or the CS-US pairing contained one presentation of each facial stimulus (i.e., three face stimuli, two CS1 and one CSÀ, per block of each conditioning task). Individual trials consisted of an initial fixation period (19,500 ms) followed by the presentation of a face (4000 ms). A fixed and long interstimulus interval (19,500 ms) was included in the experimental design to reduce intrinsic noise correlations and enable trial by trial analyses by means of RSA (Visser et al., 2013(Visser et al., , 2016. During CS-only trials, all faces were presented without tactile stimulation. During CS-US paired trials, two of three facial stimuli presentations overlapped with tactile stimulation, thus creating two CS1 and one CSÀ. The US was delivered from the midpoint of the face presentation (2000 ms postonset), remained for the rest of the time the face was visible (2000 ms) and persisted following the offset (2000 ms; total US = 4000 ms). The order of face presentation was randomized within each CS-US paired block. Participants completed two experimental tasks (one for each US, order counterbalanced across participants), totaling 26 blocks (six CS-US paired and seven CS only blocks for each US type).
Note that this paradigm was designed to target two distinct questions: (1) How and where are discriminatory and hedonic tactile signals represented in the brain? and (2) How do novel affective associations change the neural representation of conditioned stimuli? The current work focuses on the former question. Importantly, trials described as CSÀ during CS-US pairings are not independent of tactile stimulation. Rather, these trials lacked pain or brush stimulation, but participants still experienced tactile stimulation from the scanner, their clothes etc.

MRI acquisition and preprocessing
MR scanning was conducted on a 3 Tesla GE Discovery MR scanner using a 32-channel head coil. For each subject, a T1-weighted MPRAGE sequence was used to obtain high-resolution anatomic images [TR = 7 ms, TE = 3.42 ms, field of view (FOV) 256 Â 256 mm, slice thickness 1 mm, 176 slices]. Functional tasks were acquired with the following multiecho (ME) EPI sequence: TR = 2000 ms, TE1 = 11.7 ms, TE2 = 24.2 ms and TE3 = 37.1 ms, flip angle 77°; FOV 240 Â 240 mm. These parameters are consistent with recent work demonstrating improved effect-size estimation and statistical power for multiecho acquisition parameters (Lombardo et al., 2016). Specifically, the multiecho sequence was chosen because of its enhanced capacity for differentiating BOLD and non-BOLD signal (Kundu et al., 2012(Kundu et al., , 2014, as well as its sensitivity for discrimination of small nuclei in areas susceptible to high signal dropout (Markello et al., 2018). A total of 468 volumes (102 slices, thickness 3.5 mm; 72 Â 72 acquisition matrix, 3.33 Â 3.33 mm) was acquired for each functional run. Pulse and respiration data were acquired with scanner-integrated devices.
Preprocessing and analysis of the fMRI data were conducted using Analysis of Functional NeuroImages software (AFNI; Cox, 1996) and the associated toolbox meica.py (Kundu et al., 2014(Kundu et al., , 2017. For maximal sensitivity during multivariate pattern detection, no spatial smoothing was performed on the data (Haynes, 2015). Preprocessing of multiecho imaging data followed the procedural steps outlined by Kundu and colleagues (Kundu et al., 2012(Kundu et al., , 2013, as described below, and used independent component analyses to define a set of components using their TE dependence to be individually classified as BOLD or non-BOLD (e.g., motion artifact). An optimally combined (OC) dataset was generated from the functional multiecho data by taking a weighted summation of the three echoes, using an exponential T2* weighting approach (Posse et al., 1999). Multiecho principal components analysis (PCA) was first applied to the OC dataset to reduce the data dimensionality. Spatial independent components analysis (ICA) was then applied and the independent component time-series were fit to the preprocessed time-series from each of the three echoes to generate ICA weights for each echo. These weights were subsequently fitted to the linear TE-dependence and TE-independence models to generate F-statistics and component-level k and r values, which, respectively, indicate BOLD and non-BOLD weightings. The k and r metrics were then used to identify non-BOLD-like components to be regressed out of the OC dataset as noise regressors. Regressor files of interest were generated for all individual trials across the experiment, modeling the time course of each stimulus presentation during each run (36 total events: two tasks Â six CS-US blocks Â three CS, with each event beginning at the face presentation onset). The relevant hemodynamic response function was fit to each regressor for linear regression modeling. This resulted in a b coefficient and t value for each voxel and regressor. To facilitate group analysis, each individual's data were transformed into the standard brain space of the Montreal Neurologic Institute (MNI).

fMRI analyses: structural regions of interest
To assess tactile [pressure and caress and nonspecific touch (nST)] and hedonic (pressure vs caress) representations in neural patterns, nine bilateral regions of interest (ROIs) were generated from the standard anatomic atlas (MNIa_caez_ml_18) implemented with AFNI. Sensory ROIs included primary somatosensory cortex (S1), secondary somatosensory cortex (S2), primary visual cortex (V1) and ventral visual structures (VVS) while integrative ROIs included amygdalae, ventromedial prefrontal cortex (vmPFC; Posse et al., 1999), anterior cingulate cortex (ACC), and both the anterior and posterior divisions or insular cortex (Ins; for review of the functional and histologic divisions of this region, see Nieuwenhuys, 2012). S1 and V1 were selected as the primary sites of tactile and visual information, respectively. VVS were chosen because of their role in visual classification (Kanwisher et al., 1997;Kravitz et al., 2013). Amygdala, vmPFC, ACC and posterior/anterior Ins divisions were selected for their hypothesized roles in affect and pain representations subdivisions (Anderson and Phelps, 2002; for rationale behind multiple insular ROIs, see Cauda et al., 2012;Chikazoe et al., 2014;Orenius et al., 2017;Kragel et al., 2018). For extended details on defining our ROI, see Table 1.

fMRI analyses: RSA
In order to identify and compare representational patterns elicited by the experimental conditions, representational similarity analysis (RSA; Mur et al., 2009;Kriegeskorte and Kievit, 2013;Fig. 1B) was performed using the PyMVPA Python package (Hanke et al., 2009). For each participant, a vector was created containing the spatial patterns derived from b coefficients from each voxel related to each face-touch pairing in each ROI. Pairwise Pearson coefficients were calculated between all vectors of a single ROI, thus resulting in a similarity matrix containing correlations between patterns of b weights for all trials for each participant (i.e., how closely the pattern of voxel activation elicited in one trial resembles the patterns of voxel activation observed in all other trials). Fisher transformations were performed on all similarity matrices to allow comparisons between participants. Following the identification of trial-by-trial correlations, each matrix weas down-sampled to represent the condition-by-condition similarity. Thus, each cell in the presented similarity matrices is the average value of all correlations between b -values elicited during the relevant trial types. The average correlation within each condition (i.e., correlation between all repetitions of a single face-touch pairing) is the mean of 30 individual correlations (6 Â 6 trial matrix trials with auto-correlations removed), while the average correlation between conditions (i.e., correlations between two different face-touch pairings) is the mean of 36 individual correlations (complete 6 Â 6 trial matrix). Correlation matrix transformations were performed using MATLAB (The MathWorks).

fMRI analyses: patterns of interest
Pattern component modeling is a method to decompose multivariate similarity patterns into distinct subcomponents (Kriegeskorte et al., 2008;Diedrichsen et al., 2011Diedrichsen et al., , 2018Kriegeskorte and Kievit, 2013). Here, we developed a novel theory-guided implementation of PCM that determines the contribution of known sources of information to overall representational patterns observed in the brain. Patterns of interest (POIs) were created illustrating the similarity matrices that would be observed if the data were to ideally represent a single type of information perfectly. POIs are functionally similar to conceptual dissimilarity models (mRDM) for modeling content of RSA patterns (Kriegeskorte et al., 2008;Kriegeskorte and Kievit, 2013;Popal et al., 2019) in that they act to fit observed data to theory-based prespecified representational patterns. Critically, however, the current approach is distinct from these in that it acts to fit combinations of conceptual patterns of interest to a given brain region, rather than individual models.
Thirteen POIs were designed to represent unique taskrelevant representational patterns could be expected to be elicited during the current task. For example, one POI would illustrate what the data pattern would look like if the voxels across a given ROI were to represent only the positive or negative valence of the touch as opposite poles of a bipolar scale. Note that the POIs were designed to capture information about the specific manipulations of the current experiment as well as cognitive processes associated with the ROIs defined in this study. The POIs we described have a particular focus on potential representations of tactile experience. Thus, they may not accurately reflect information represented in brain regions processing predominantly nontactile information and are not intended to be an exhaustive model of all potential representational states. The following is an extended description of each POI, with additional description and visualization presented in Table 2.

Experimental task (ET)
Experimental task was conceptualized as representational overlap for all trials contained in one conditioning task (i.e., the aversive vs appetitive task). This included both trials with and without tactile stimulation. Ideal representation of Experimental task was defined as perfect correlation (r = 1) between all trials in each specific conditioning task. All correlations for trial across conditioning tasks were set to r = 0.  • i.e., within separate aversive and appetitive tasks, all the trials (in the presence and absence of touch) are represented as similar.
Nonspecific touch nST r = 1 for all comparisons between all trials where a tactile manipulation occurred.
• i.e., shared representation both within and between aversive pressure and appetitive caress trials.
Specific touch ST r = 1 for all comparisons between all trials where an identical tactile experience occurred.
• i.e., shared representation within, but not between trials with aversive pressure, trials with appetitive caress and unpaired trials, respectively.
Appetitive brush AC r = 1 for all comparisons between all trials that involved the delivery of an appetitive caress to the participant's arm.
Aversive pressure AP r = 1 for all comparisons between all trials that involved the delivery of aversive pressure to the participant's thumb.
Touch valence TV r = 1 for all comparisons between trials involving tactile stimulation within each conditioning task. r = À1 for all comparisons between trials involving tactile stimulation between conditioning tasks. NOTE: This reflects a linear representation of hedonic tactile information.
Positive events PE r = 1 for all comparisons between all trials experienced as positively valenced relative to its experimental task.
• i.e., tactile stimulation trials in the appetitive task are represented as similar to trials with no stimulation in the aversive task.
Negative events NE r = 1 for all comparisons between all trials experienced as negatively valenced relative to its experimental task.
• i.e., tactile stimulation trials in the aversive task are represented as similar to trials with no stimulation in the appetitive task.
All valence AV r = 1 for all comparisons between trials containing positive events and between trials containing negative events. r = À1 for all comparisons between trials of positive events and negative events. *NOTE: This reflects a linear representation of all hedonic information.

Salience
Sa r = 1 for all comparisons between all trial with highly tactile salience.
• shared representation both within and between trials with appetitive and aversive stimulation. r = 1 for all comparisons between all trial with minimal tactile salience.
• shared representation both within and between trials with no appetitive/aversive stimulation. r = À1 for all correlation between highly and minimally salient trials.
Research Article: Methods/New Tools Nonspecific touch (nST) Nonspecific touch was conceptualized as identical representations for all trials where a tactile manipulation occurred. Importantly, this POI does not represent any discriminable information between the two tactile manipulations (aversive vs appetitive), but rather defines them as sharing representational space. This POI may be derived from information carried by nonhedonic (A-fiber) or hedonic (C/C-tactile) peripheral channels; nonhedonic signals may indicate quality/strength of tactile information while hedonic signals may indicate emotional saliency of tactile experience. Note, however, that this pattern is unlikely to occur as a result of first order hedonic projection (which are separated into unique positive and negative valence signals), as it includes, by definition, a shared representational space across conditions. Ideal representation of nonspecific touch was modeled as perfect representational overlap (r = 1) between each appetitive brush (AB)-paired trial and all other appetitive brush (AB-AB) trials, each aversive pressure-paired trial and other aversive-pressure (AP-AP) trials, as well as aversive pressure and appetitive brush (AP-AB) trials. All correlations involving no tactile manipulations were defined as r = 0.

Specific touch (ST)
Specific touch was conceptualized as representing each tactile experience as occupying a unique representation space. Critically, in addition to a unique representational space for appetitive and aversive stimulation, this POI also defines the tactile experience when there was active manipulation as occupying its own representation space. Thus, a pattern similar to this POI cannot be derived from hedonic tactile channels, but rather must rely on alternate signals of discriminative touch (i.e., information carried by A-fibers). Ideal representation of specific touch entailed perfect correlation (r = 1) between all trials in which similar tactile experience occurred. This included AP-AP trials, AB-AB trials, and no manipulation (i.e., generic scanner sensation) trials. All additional correlations crossing tactile experience types were set to r = 0.

Appetitive brush (AB)
Appetitive brush was conceptualized as representing information carried by C-tactile fiber pathways and the associated experience of pleasant touch. Ideal representation of appetitive brush was defined as perfect correlation (r = 1) between all trials that involved the delivery of an appetitive brush stroke to the participant's arm (i.e., AB-AB trials). It should be noted that similarity between these trials could be ascribed to either the tactile sensation of brushing or the hedonic nature of the experience. All additional correlations were set to r = 0.

Aversive pressure (AP)
Aversive pressure was conceptualized as representing information carried by C-fiber nociceptive pathways and the associated experience of painful stimulation. Ideal representation of aversive pressure was defined as perfect correlation (r = 1) between all trials that involved the delivery of aversive pressure to the thumbnail (i.e., AP-AP trials). As with Appetitive brush, it should be noted that similarity between these trials could ascribed to either the tactile sensation of pressure or the hedonic nature of the experience. All additional correlations were set to r = 0.

Touch valence (TV)
Touch valence represents painful and pleasant tactile experience as opposite ends of a shared linear representational space. Ideal representation of touch valence was defined as a linear representation of the hedonic aspect of the tactile manipulation. For this, all correlations between all AB-AB trials and all AP-AP trials were set to r = 1. To represent the contrasting hedonic experience of AB-AP trials, these correlations were set to r = À1. All additional correlations were set to r = 0.

Positive events (PE)
Positive events were conceptualized as a general representation for individual trials that are positively valenced • i.e., distinct representation for each of the 6 faces (3 faces Â 2 tasks).
Violation of expectation VE r = 1 for all comparisons correlation between all trials involving the less probable outcome. NOTE: As there were two faces paired with tactile stimulation compared with one face never paired with tactile stimulation for each experimental task, the less probable outcome was always the unpaired trials.
Temporal adjacency TA r = 1 for all comparisons between all comparisons that included trials that were temporally contained within the same block (i.e., temporally adjacent exposures). NOTE: Because of the removal of autocorrelation from within condition averaging, these did not contain temporally adjacent trials.
*General matrix structure can be found in Figure 1B; red represents r = 1, blue represents r = À1 and yellow represents r = 0.
relative to the context in which they are situated. This includes both the positive valence related to the presence of pleasant brushing during the appetitive conditioning task and the absence of painful pressure during the aversive conditioning task. Ideal representation of positive events was defined as perfect correlation (r = 1) between all trials where the event was positive relative to its context. Positive trials (PT) included the delivery of a brush stroke during the appetitive conditioning task as well as the absence of pressure in the aversive conditioning task. All additional correlations were set to r = 0.
Negative events (NE) Negative events were conceptualized as a general representation for individual trials that negatively valenced relative to the context in which they are situated. This includes both the feeling of painful pressure during the aversive conditioning task, as well as the absence of pleasurable brushing during the appetitive conditioning task. Ideal representation of negative events was defined as perfect correlation (r = 1) between all trials where the event was negative relative to the rest of the conditioning task. Negative trials (NT) included the delivery of pressure during the aversive conditioning task as well as the absence of a brush stroke in the appetitive conditioning task. All additional correlations were set to r = 0.

All valence (AV)
All valence was conceptualized as a linear representation of the hedonic aspects of the experimental procedure relative to the task in which they were situated. Here, "positive events" and "negative events" are independent POIs that are considered to be opposite ends of a shared linear representational space. To represent this, all PT-PT and NT-NT correlations were set to r = 1. All PT-NT correlations were set to r = À1.

Salience (Sa)
Salience was conceptualized as a linear representation of the salience of individual trials, with "high salience" and "low salience" considered as opposite ends of a shared linear representational space. For this, correlations between all highly salient trials (i.e., AB and AP trials) were set to r = 1. In addition, correlations between all minimally salient trials (i.e., CSÀ trials) were also set to r = 1. All correlations between highly and minimally salient trials were set to r = À1.

Face stimulus (FS)
Face stimulus was conceptualized as a unique representational space for each independent facial identity. To model this, FS defines correlations for all trials contained an identical visual CS (six in total, three/experimental task) as r = 1. All additional correlations between different CS were set to r = 0.

Violation of expectation (VE)
Violation of expectation was conceptualized as a shared representation space for the less likely tactile outcome in each conditioning task As there were two CS-US paired trials for each CS-only trials for each experimental task, the less probable outcome was always the CS-only trials. As such, the VE POI modeled all CS-minus-CSminus correlations (within and across experimental task) as r = 1. All additional correlations were set to r = 0.

Temporal adjacency (TA)
Temporal adjacency was conceptualized as representing any task-irrelevant cognitive processes that may extend beyond single modeled events, and thus inform the representational pattern of temporally adjacent trials. Ideal representation of temporal adjacency was defined as perfect correlation (r = 1) between all comparisons that included trials that were temporally contained within the same block (i.e., within each group of three CS presentations). Because of the removal of autocorrelation from within-condition averaging, these did not contain temporally adjacent trials, nor did any correlation of trial between experimental tasks. Thus, correlations for both of these comparisons were set to r = 0.
In consideration of the dominant contralateral input to S1 and the lateralized tactile stimulation across tasks (aversive pressure applied to the RIGHT thumbnail, and appetitive caress with a brush applied to the LEFT forearm), two additional POIs were described to facilitate unilateral investigation of this region. There were as follows.

Right specific touch (rST)
This POI represents the discriminable tactile experience expected elicited from the right side of the body. It represents both aversive pressure pain (applied to the right thumb) and generic scanner touch (including that occurring on the right side during appetitive brush trials) as dissociable tactile experiences. This POI is predicted to manifest in representation in left S1.

Left specific touch (lST)
This POI represents the discriminable tactile experience expected elicited from the left side of the body. It represents both appetitive brushing (applied to the left forearm) and generic scanner touch (including that occurring on the left side during the application of aversive pressure) as dissociable tactile experiences. This POI is predicted to manifest in representation in right S1.

fMRI analyses: patterns component modeling
To determine the POI combinations that best explained the observed data patterns for each ROI, we conducted Monte Carlo cross-validated (CV) PCM using Bayesian information criterion (BIC) to fit our pattern component models. An uninformed greedy best-first search (GBFS) algorithm (Doran and Michie, 1966) was implemented to identify the best fitting POI combination in a step wise manner (see Fig.  1C). This allows us to distinguish between subtle nuances in representational patterns by iteratively layering POIs in a step-wise manner rather than all at once, thus comparing the fit of overtly similar component patterns both independent from, and in combination with, each other. Initial model testing was conducted by fitting each independent POI to the observed similarity between conditions for a given ROI (Level 1). Upon identification of the best fitting POI (POI B1 ), model fitting was conducted on each remaining POI in combination with POI B (Level 2) . The POI combination (i.e., POI B1 1 POI B2 ) that provided the best fit to the ROI data would be held as a constant for model fitting in Level 3. This process was repeated iteratively until no addition of remaining POIs led to an improved fit to the ROI. A DBIC . 2 was defined as indicative of an improved fit (Fabozzi, 2014). Following similar equivalency criteria, all POI combinations at a given search level with DBIC scores , 2 to the best fitting combination were also extended to path completion (Fabozzi, 2014). This approach allowed for the decomposition of observed representational patterns into multiple unique contributing sources of information. The independent pattern components identified as contributing to the representational space were subsequently fit using linear regression modeling to the observed similarity in the original ROI (see Fig. 1D) to determine the weight of each contributing POI to the observed similarity from the full ROI.
To ensure that regression fits were not a product of overfitting, the following cross-validation procedures were performed: Initial model fitting was performed on a randomly selected sample of participants ("randomsample": RS = 60), with the identified components fit as a predictor to data from the remaining participants held-out of this initial sample (the "hold-out": HO = 7). Monte Carlo cross-validation (CV; Picard and Cook, 1984) parameters were chosen to maximize cross-validation performance by minimizing cross-validation variance while maximizing model selection accuracy (Arlot and Celisse, 2010). These analyses identified POIs contributing to representational patterns for each ROI in the RS. Beta coefficients and intercepts, determined by fitting these POIs as predictors to the experimental data, were used to create a reconstructed and averaged dataset. The reconstructed dataset was then fitted as a predictor to the hold-out, with each iteration approximating a single fold of a 10-fold validation.

Results
Results described here are generated from a Monte Carlo cross-validation procedure (1000 iterations) conducted on a random sample of 60 participants (RS = 60). A weighted re-combination of these POIs was then validated against the remaining seven participants (HO = 7; see

POI identification and weighting
To determine the combinations of POIs that best explained the observed patterns of data for each ROI, PCM was conducted through Bayesian information criterion (BIC) analyses and an uninformed greedy best-first search (GBFS) algorithm. This strategy of iterative model fitting with progressively increasingly complex predictor combinations contrasts previous RSA approaches which typically do not combine potential predictors (Kragel et al., 2018;Dobs et al., 2019). Through the current analyses, we are able to determine the degree information was integrated versus separate in any specific region. To ensure that regression fits were not a product of overfitting, we performed Monte Carlo cross-validation. Specific outputs of interest included the proportion of cross-validation iterations (i = 1000) in which a POI was identified as a contributing component of the experimental data in the random samples, the representational weight of those components identified at a rate significantly greater than chance, and the model fit of the reconstructed RS components to the holdout (Fig. 2; for complete summary of the crossvalidation results, see Table 3). The proportion of iterations for POI identification was compared with chance identification for each region of interest (ROIs; i.e., number of POIs identified/total number of POIs). For each iteration of the MCCV procedure, the total number of paths required for a given search is defined as the n-path.

Sensory regions of interest
Primary somatosensory cortex (S1). In S1, a POI modeling nonspecific aspects of tactile experience (nST) was the strongest individual predictor of S1 representational patterns in the random samples (b nST = 0.129). Additional POIs identified contributing to overall representational patterns in S1 modeled experimental task and nonhedonic tactile experience (b ET = 0.055 and b ST = 0.025, respectively). In the held-out sample, POIs identified in the random sample explained 24.8% of the variance (R 2 = 0.248, F (1,145) = 52.13, p , 0.001). This pattern indicates a representation of distinct discriminative tactile experiences rather than hedonic value in S1 (Fig. 3A).
Two additional PCM analyses were conducted in S1 ROIs that interrogated right and left S1 ROIs independently. Importantly, these analyses added two additional POIs, modeling the left and right lateralized components of nonhedonic tactile experience. In a pattern similar to that observed across the bilateral ROI, data from the random sample were predicted most strongly by POIs modeling nonspecific aspects of tactile experience (Left S1: b nST = 0.140; Right S1: b nST = 0.129), with secondary contributions from the POI modeling global similarity within, but not between, experimental tasks (Left S1: b ET = 0.052; Right S1: b ET = 0.066). This similarity, however, was not observed for representations of nonspecific tactile experience observed in the bilateral ROI. Left S1 did not represent aversive pressure as isolated from other forms of nonhedonic tactile states (b AP = 0.032) and was lacking general representation for right lateralized nonhedonic touch. By contrast, right S1 represented nonhedonic touch experience (i.e., appetitive caress and scanner-generic touch as distinct states b lST = 0.045). In both unilateral ROIs, POIs contributing to representational patterns in the random sample also significantly predicted the representational patterns of the held-out participants (Left S1: R 2 = 0.237, F (1,145) = 49.68, p , 0.001; Right S1: R 2 = 0.240, F (1,145) = 49.97, p , 0.001).
Secondary somatosensory cortex (S2). For S2, the strongest predictor of representational patterns was nonspecific touch (b nST = 0.185). Additional POIs identified as contributing to overall representational patterns were aversive touch (b AP = 0.094), experimental task (b ET = 0.041), and the task-specific positive experience (i.e., caress, or safety; b PE = 0.002). Combined, weighted POIs identified as being represented in the random sample Figure 2. Cross-validation across brains. A 1000 iteration Monte Carlo cross-validation determined (1) that identified POIs from the whole sample data (n = 67) were reliably identified when the procedure was replicated on subsets of the sample (n = 60) and (2) that reconstructed data generated through POI identification and weighting accurately predicted activational similarity pattern in the held-out participants (n = 7). Results from each MCCV iteration are represented as a row of data, with the identified POI noted and the dark blue, and the fit to the HO shown in the center-right column for each ROI. Data summaries collapsed across all MCCV iterations are shown in the red box for each ROI.   accounted for an average of 41.1% of the variance in the held-out participants (R 2 = 0.411, F (1,145) = 113.07, p , 0.001). This suggests that S2 may receive hedonic signals that are not represented in S1 (Fig. 3B). Visual cortices. The most predictive POI in both V1 and ventral visual stream modeled task-related changes in experience. No other POIs were identified at a rate significantly greater than chance. In V1, this component explained an average of 5.9% of the variance in the initial random sample, while in VVS, it explained an average of 6.1% of the variance. This weighted component, however, failed to significantly predict patterns observed in the held-out participants in cross-validation procedures (V1: R 2 = 0.021, F (1,145) = 4.27, p = 0.10; VVS: R 2 = 0.028, F (1,145) = 5.27, p = 0.07). This suggests that representational patterns in visual cortices may reflect visual attentional demands of the experimental task and are relatively uninformative of nonvisual information, regardless of its hedonic value (Fig. 3C,D). The inability to replicate significant findings in the held-out participants, however, suggests that activational patterns in these regions are likely not driven by information modeled in the current set of POIs, but may instead be more accurately represented by untested representational patterns related to specific aspects of visual (rather than tactile) experience

Integrative regions of interest
Amygdalae. A combination of three POIs were identified as contributing to the representational patterns of bilateral amygdalae in the initial random sample. These POIs modeled a linear spectrum of tactile valence (b TV = 0.014), nonspecific tactile experiences (b nST = 0.019), and global differences in experimental tasks (b ET = 0.011), respectively. POIs identified in the random sample accounted for 11.6% of the observed variance in the held-out participants (R 2 = 0.116, F (1,145) = 20.86, p = 0.0049). Notably, activation patterns in the amygdalae represented valence on a linear spectrum, where appetitive and aversive touch were most dissimilar, polar opposites of a shared representational space (Fig. 4A).
Ventromedial prefrontal cortex. In the vmPFC, representational patterns were predicted by multiple POIs, most notably the two distinct hedonic touch POIs: Aversive pressure (b AP = 0.050) and Appetitive caress (b AB = 0.033). Additional POIs identified as contributing to the representational patterns in the random sample modeled differences in experimental task (b ET = 0.011), nonspecific tactile representations (b nST = 0.015), and temporal order of trial events (i.e., whether or not the similarity matrix reflected temporally adjacent events; b TA = 0.005). Total variance accounted for in the held-out participant by the models identified in the random samples was on average 7.4% (R 2 = 0.074, F (1,145) = 12.98, p = 0.023). This demonstrates that vmPFC activity contains information about the hedonic value of the tactile stimulation, representing positive and negative values as distinctly independent and nonopposing signals (Fig. 4B). Furthermore, the heterogeneity of nonhedonic POI identification in this region this suggests that while vmPFC does consistently represents aversive pressure and appetitive caress, there may be extensive interparticipant variability regarding what information can be in this region.
Anterior cingulate cortex. The aversive pressure (AP) POI was the strongest predictor of ACC activation patterns (b AP = 0.067). Additional POIs identified as contributing to the representational patterns included those modeling nonspecific tactile experience (b nST = 0.055) and experimental task (b ET = 0.034). A weighted combination of the POIs identified in the random samples predicted an average of 19.1% of the variance in the held-out participants (R 2 = 0.191, F (1,145) = 37.84, p , 0.001). This suggests that the ACC represents general tactile information but is particularly sensitive to tactile information associated with pain (Fig. 4C).
Insula. The insula was anatomically subdivided at the anterior commissure into distinct nonoverlapping anterior/posterior regions. Representational patterns in the anterior insula (aIns) were significantly predicted by four POIs. In order of representational strength, these POIs modeled aversive tactile experience (b AP = 0.107), nonspecific tactile experience (b nST = 0.050), experimental task (b ET = 0.029), and task-contextualized negative events (i.e., aversive touch OR lack of appetitive touch; b NE = À0.016). POIs identified in the random sample predicted an average of 15.3% of the variance in the held-out participants (R 2 = 0.153, F (1,145) = 28.74, p , 0.0023).
As in the aIns, in pIns, activity was significantly predicted by POIs modeling aversive tactile experience (b AP = 0.059), nonspecific tactile experience (b nST = 0.102), and experimental task (b ET = 0.043). An additional POI modeling tactile valence (b TV = 0.009) was also identified in this region. In pIns, combinations of POIs identified in the random sample predicted an average of 36.3% of the variance in the held-out participants (R 2 = 0.363, F (1,145) = 90.82, p , 0.001). This demonstrates that whereas the general type of information processed across the insula may be similar for the anterior and posterior sections, each region Provided are the most likely combination of pattern component models (POIs) to explain observed representational patterns for each ROI. Model fitting was performed by identifying the best individually fitting POI, then iteratively adding remaining POIs. POIs held from prior levels of analyses are indicated in under "Included1" (forth column from left), with the BIC score resulting from combining the held POIs independently with all remaining POIs displayed in the right columns. Significantly improved fit because of the addition of an additional POI is indicated by DBIC . 2. Note, identified final POI combinations predicted the observed data significantly better than a combination of all POIs ("All"), and a single uniform predictor ("None").
sensitive to both hedonic and nonhedonic signals, the precise nature and dominance of these representations differs (Fig. 4D,E).

n-Path analyses
To assess the robustness of POI contributions, a oneway ANOVA was conducted on the average number of search paths required to find the best fitting component combination for each iteration (i.e., n-path data). This identified a significant main effect of region (F (8,7992) = 195.507, p , 0.001). A follow-up series of independent sample t tests (all reported p values are Bonferroni-corrected) identified four distinct clusters of ROIs characterized by their n-path. A lower search path likely indicates either a poor fit of the POI models (if only a single model is frequently identified; e.g., V1/VVS), or robust representations for a specific subset of models (if identified POI . 1; e.g., vmPFC). By contrast, a higher n-path likely indicates more overlapping representational space (e.g., insular subdivisions). Visual areas required expansion of fewer paths than any other area (all p , 0.001), yet they did not differ significantly from each other (p = 1.0). vmPFC had more branching than either V1 or VVS but less than all other ROIs (all p , 0.001). ACC, amygdala, and S2 did not significantly differ from each other, yet required less expansion of search paths than S1 or either insular ROI (all ps , 0.001). Finally, while the anterior insula did not significantly differ from either the posterior insula or S1 (both ps = 1), the posterior insula displayed greater branching than S1 (p , 0.001). The greater n-paths in these regions suggests that there is likely greater overlap of representational space in these regions between the modeled components compared with representations in other regions.

Discussion
In this study, we applied a novel form of pattern component modeling with representational similarity analysis to demonstrate that by fitting predefined representational patterns (POIs) to observed patterns of similarity between voxel activation patterns in the brain, we are able to model the specific nature of information represented in a brain region, teasing apart isolated or integrated signals of valence and experience. Specifically, in the current work, we were able to dissociate how discriminatory versus hedonic tactile information, coded at the somatosensory receptors and carried by A-and C-/CT-fibers, respectively, contribute to population coded representations in the human brain. Distinct representations of hedonic information were observed in frontal and temporal structures, including ventromedial prefrontal cortex (vmPFC), insula (Ins) and anterior   With no clear POI providing optimum fit at the 2nd and 3rd level of GBFS-BIC analyses (i.e., POI 1 -POI 2 , 2), search paths for each statistically equivalent option were completed. The POI combination providing the lowest BIC at completion of the path was defined as the best fitting combination.
cingulate cortex (ACC), as well as in secondary somatosensory cortex (S2). Primary somatosensory cortex (S1) did not represent all tactile information coded by peripheral receptors, with limited representation of some hedonic signaling. Specifically, we did not observe any representation of positive hedonic touch signals and only limited representation of negative hedonic signals (carried by CT-fiber afferents and C-fiber afferents, respectively. Visual areas, including primary/secondary visual cortex (V1) and ventral visual structures (VVS), displayed no representation of either affective or discriminative touch information. Together, the findings support the hypothesis that processing of sensation carried by hedonic-labeled tactile signals from C-and CT-fiber pathways, despite their salience and homeostatic significance, does not depend on representation of these signals in S1. Rather, this information is represented predominantly in frontotemporal structures more typically implicated in interoception (Craig, 2011;Pollatos et al., 2016;Strigo and Craig, 2016) and the central mediation of emotional relevance (McFarland and Sibly, 1975;Rolls, 2000;Todd et al., 2020). Of note, however, negative hedonic information made a minor contribution to representational patterns in S1 contralateral to the tactile stimulation, indicating that some nociceptive information does reach this area independent of frontotemporal processing. Thus, some propagation of tactile hedonic information is distinct from traditional exteroceptive signals highlighting nontraditional mechanisms (Kryklywy et al., 2020) by which prioritized information may be incorporated into emotionally-guided cognitive processes.

Cortical representations for nonhedonic touch
In S1, neural activity displayed representational patterns that discriminated tactile experiences, as well as nonvalence specific components, of tactile manipulations. The finding of strong representation of specific touch experiences in this area are consistent with S1's traditional role in exteroception, processing discriminatory tactile information as carried by out A-fiber afferents (McGlone and Reilly, 2010;McGlone et al., 2014). The POI "specific touch" is defined such that representational patterns for trials with no hedonic manipulation (i.e., the somatosensory experience generated by lying in a scanner) have equivalent contribution to the global representation of touch experience as trials with hedonic manipulation. Thus, its manifestation is unlikely to be generated by peripheral hedonic signaling alone. In addition to specific touch, S1 strongly represented the experience of "nonspecific touch" indicating activity in this region was also influenced by salient tactile experiences and contained shared representational space for both appetitive and aversive tactile manipulations. This suggests that these representations are likely not shaped by information carried by C-and CT-fiber in isolation, as the two distinct peripheral signals are represented by overlapping activation representational patterns. Nonspecific touch representation is likely to be either a discriminatory representation of body location (i.e., arm; not dependent on C-or CT-fiber activation) or general tactile salience (may or may not integrate information from C-and CT-fiber activation; i.e., hedonic salience). In support of the latter interpretation, there is evidence that S1 likely integrates reentrant hedonic signals from multisensory emotion-related regions (Vuilleumier, 2005;Pessoa and Adolphs, 2010;Orenius et al., 2017).
Components indexing unprocessed projections of hedonic-labeled afferent pathways (i.e., "aversive pressure" and "appetitive caress") were absent in bilateral representational patterns observed in S1, which may on first inspection be interpreted as hedonic signals not being instantiated in traditional somatosensory processing structures. Yet, unilateral investigation of left S1 identified representation patterns associated with aversive pressure in absence of specific touch experience, indicating that discriminative and nociceptive information may be integrated before reaching S1 (for candidate regions, see Abraira et al., 2017;Marshall and McGlone, 2020;Neubarth et al., 2020). Alternatively, it may be that sustained changes in tonic firing of rates of slow adapting mechanoreceptors (for review, see Knibestol, 1975;Abraham and Mathew, 2019) in response to the strong pressure manipulation (right hand), result in a distinct tactile representation in left S1 during CSÀ trials (pressure task). This representation would be distinct from the representation of scanner-generic sensation (right hand) experienced during the caress (which was applied to the left hand). Given substantive evidence to support both interpretations, it is likely that the observed POI contributions reflect a combination of these processes. No evidence of unique representational patterns in S1 for signals of appetitive hedonic information was observed.
Taken together, nonhedonic tactile representation, likely of signals carried along A-fiber pathways (McGlone and Reilly, 2010), appear to dominate activity in early somatosensory cortices. While some evidence for hedonic representations in these regions exists, it appears that precortical integration of A-and C-/CT-fiber pathways (Abraira et al., 2017;Marshall and McGlone, 2020;Neubarth et al., 2020), or re-entrant feedback from higher order integrative structures (Vuilleumier, 2005;Pessoa and Adolphs, 2010) are the most probable source of these representations.

Cortical representations for hedonic touch
Amongst all regions investigated, only the vmPFC displayed independent representation of both appetitive and aversive touch (POIs: AC/AP). This suggests that this region either (1) receives information carried along C-and C-tactile fiber afferents as distinct signals before their integration with each other or other tactile information, or (2) has decomposed an integrated hedonic representation back into distinct signals of positive and negative value to inform situation specific behaviors and decision. The potential of first order representation of peripherally labeled hedonic signals in vmPFC is particularly intriguing considering the critical role these ventral medial structures play in appraising emotional salience to guide value-based decision-making (Euston et al., 2012;Dixon et al., 2017;Hiser and Koenigs, 2018). Propagation of hedonic-labeled tactile signals to these regions independent of any prior cortical processing would act as a mechanism to facilitate the prioritization of evolutionarily relevant sensation (Kryklywy et al., 2020), and allow for expedited integration of action-outcomes into value appraisal to guide decision-making processes.
The absence of distinct representations for pleasurable tactile signals in both the anterior and posterior insula is notable, as these regions have been highlighted as potential cortical recipients of C-tactile fiber signaling (Olausson et al., 2002;Rolls et al., 2003). While this may be because of variability in response to the appetitive touch manipulation used in the current design, the identification of a clear appetitive caress contribution to representational patterns in vmPFC indicated that this is unlikely. It is important to note much of the prior work identifying modulation of insula activity by pleasurable touch has been performed either independent of aversive touch (Olausson et al., 2002), or examined by treating the two signals orthogonally without direct comparison (Rolls et al., 2003). This leaves open the possibility that previous results were driven by general affective salience of the tactile cue (as observed in the current work), rather than pleasurable sensation alone. Related to this, the potential for sensory adaptation or habituation to hedonic tactile experience should be considered (McBurney and Balaban, 2009;Morrison, 2016). The current design relies on the averaged representational similarity across six distinct tactile exposures for each condition. Thus, it is possible that differences in the time course of downregulating tactile signaling between painful and pleasurable stimulation following multiple exposures may underlie some of the differences in representational strength of hedonic information. This question could be targeted by future studies that dissociate early and late exposure to hedonic touch.
Although distinct representation of appetitive touch was identified only in the vmPFC, distinct representations of aversive pressure were identified within the ACC as well as the anterior and posterior insula. Notably, both of these regions are heavily implicated in the representation of painful experience (Corradi-Dell'Acqua et al., 2016;Kragel et al., 2018) and are postulated to underlie awareness of one's own internal homeostatic balance (Craig, 2011(Craig, , 2015Pollatos et al., 2016;Strigo and Craig, 2016). One potential explanation for this pattern of results is that hedonic-labeled peripheral afferents are not processed as tactile signals in the traditional view of sensation (Pinel and Barnes, 2018;Gazzaniga et al., 2019). That is, they may not be instantiated in neocortex as representing the experience of contact with external objects in the environment. Rather, information carried along these pathways indicates internal concerns about homeostatic threat or social safety (Craig, 2011(Craig, , 2015 and manifest cognitively as emotional feelings congruent with these states. Information about the internal state acts can then act as an immediate mechanism for motivating response, independent of its representation as an exteroceptive tactile experience or any other form of cognitive processing.

Integrated representation of tactile experience
Multiple regions displayed patterns of activity that integrated signals from both C-and CT-fiber afferents. These representational patterns showed distinct similarity/dissimilarity between hedonic conditions, rather than maintaining the independence of their sources, indicative of prior processing and representational integration of this information. Specifically, in the ACC as well as the anterior and posterior insula, patterns of neural activity were found to represent nonspecific touch beyond those identified for aversive pressure. These integrated representations likely indicate that processing of sensory information has occurred before their affective representation the insula and ACC, in a manner consistent with traditional models for emotional prioritization (Rolls, 2000(Rolls, , 2019Vuilleumier, 2005;Pessoa and Adolphs, 2010).
While none of the anterior or posterior insula, ACC, or vmPFC were found to display opposing representations of hedonic valence, in the amygdala, a distinct representation of the hedonic experience of touch valence was observed. Here, signals of tactile valence were represented as a single linear vector, with positive and negative hedonic conditions represented as polar ends of a single valence spectrum. Thus, before its representation, or as part of its processing in the amygdala, tactile information must be integrated into the same representational space. These amygdalar bi-polar valence representations are consistent with those identified in the olfactory domain (Jin et al., 2015), but have not been reported for either gustatory or visual hedonic information (Chikazoe et al., 2014(Chikazoe et al., , 2019. This divergence indicates a probable modal-specificity of hedonic processing in the amygdala rather than a centralized a-modal representation of emotional information (Miskovic and Anderson, 2018). In the current study, this unidimensional hedonic vector may be related to the association of the tactile sensation with the concurrent visual stimuli rather than the tactile signals in isolation, as multiple studies in both humans and nonhuman primates have implicated this region in guiding affect-biased attention  and emotional learning in vision (Morris et al., 1998;Everitt et al., 2003). Of note, however, there is substantial evidence to suggest that both hedonic responding and attentional biases may be heavily influenced by individual differences between subjects (Nielsen et al., 2009;Todd et al., 2012;Harjunen et al., 2017). Future work conducted in a larger sample population could provide additional insight into this possibility.

Analytic considerations
By iteratively fitting combinations of predefined models of representational patterns to observed data, rather than fitting each pattern in isolation, we were able to interrogate the representation of information with more conceptual resolution than previous work. Specifically, this approach allows researchers to determine whether individual brain regions represented information as either isolated or integrated patterns, as well as the strength of each representational pattern in the context of others. Applying this technique to targeted hypothesis-specific representations, we were able to identify meaningful differences in the cortical representation of discriminative versus hedonic somatosensory signals.
It is notable that with the current implementation, all patterns of interest were defined as matrices of either 0 or 1 (i.e., perfect or no correlation). While the reliability of results we obtained (see the cross-validation procedures), and their consistency to hypothesized functioning in the tested regions (Craig, 2011(Craig, , 2015Cauda et al., 2012;Jin et al., 2015;Dixon et al., 2017;Miskovic and Anderson, 2018), support its implementation in this form, future work using similar technique may look to add additional refinement to these models. In many cases, there may be a well-supported reason to expect partial, or scaled correlation between conditions, rather than absolutes. While the coarse approach to modeling implemented in the current work does provide meaningful insight into the nature of tactile representation, further refinement of the representational pattern investigated may provide additional details to this picture.
An additional limitation to the current work that by constraining our search to a greedy best-first approach, we have by algorithmic definition limited the combinations of POIs tested. We believe that our current approach provides greater refinement toward understanding how multiple sources of information contribute in parallel to observed neural representational patterns than previous approaches, and does so in a computationally efficient manner. That said, we also acknowledge the importance of future work to explore and compare alternate search algorithms as well as computationally intense brute-force approaches. We believe that POI-guided pattern component modeling represents an exciting new approach to multivariate analyses, and that its application to human neuroimaging data can allow for increasingly detailed "read out" into how and where information in processed within the human brain.
In conclusion, in this work, we have outlined an exciting extension of traditional representational similarity analyses. Using this approach, we demonstrated that hedonic tactile information is not processed in the same fashion as nonhedonic tactile information. The full spectrum of hedonic tactile information, signals carried in the periphery by C-and CT-fiber pathways, is not represented in primary sensory cortices but is represented in some frontotemporal structures typically associate affective processing. Additionally, we believe that the POI-based approach to pattern component modeling outlined here reflects an exciting new avenue for future multivariate neuroimaging analyses. It provides a tool to decompose observed patterns of representational similarity into intuitive, theory-guided representational subcomponents (i.e., POIs), thus allowing researcher to probe deeper into the content of information instantiated by the brain.