Cortical Regions Encoding Hardness Perception Modulated by Visual Information Identified by Functional Magnetic Resonance Imaging With Multivoxel Pattern Analysis

Recent studies have revealed that hardness perception is determined by visual information along with the haptic input. This study investigated the cortical regions involved in hardness perception modulated by visual information using functional magnetic resonance imaging (fMRI) and multivoxel pattern analysis (MVPA). Twenty-two healthy participants were enrolled. They were required to place their left and right hands at the front and back, respectively, of a mirror attached to a platform placed above them while lying in a magnetic resonance scanner. In conditions SFT, MED, and HRD, one of three polyurethane foam pads of varying hardness (soft, medium, and hard, respectively) was presented to the left hand in a given trial, while only the medium pad was presented to the right hand in all trials. MED was defined as the control condition, because the visual and haptic information was congruent. During the scan, the participants were required to push the pad with the both hands while observing the reflection of the left hand and estimate the hardness of the pad perceived by the right (hidden) hand based on magnitude estimation. Behavioral results showed that the perceived hardness was significantly biased toward softer or harder in >73% of the trials in conditions SFT and HRD; we designated these trials as visually modulated (SFTvm and HRDvm, respectively). The accuracy map was calculated individually for each of the pair-wise comparisons of (SFTvm vs. MED), (HRDvm vs. MED), and (SFTvm vs. HRDvm) by a searchlight MVPA, and the cortical regions encoding the perceived hardness with visual modulation were identified by conjunction of the three accuracy maps in group analysis. The cluster was observed in the right sensory motor cortex, left anterior intraparietal sulcus (aIPS), bilateral parietal operculum (PO), and occipito-temporal cortex (OTC). Together with previous findings on such cortical regions, we conclude that the visual information of finger movements processed in the OTC may be integrated with haptic input in the left aIPS, and the subjective hardness perceived by the right hand with visual modulation may be processed in the cortical network between the left PO and aIPS.

Recent studies have revealed that hardness perception is determined by visual information along with the haptic input. This study investigated the cortical regions involved in hardness perception modulated by visual information using functional magnetic resonance imaging (fMRI) and multivoxel pattern analysis (MVPA). Twenty-two healthy participants were enrolled. They were required to place their left and right hands at the front and back, respectively, of a mirror attached to a platform placed above them while lying in a magnetic resonance scanner. In conditions SFT, MED, and HRD, one of three polyurethane foam pads of varying hardness (soft, medium, and hard, respectively) was presented to the left hand in a given trial, while only the medium pad was presented to the right hand in all trials. MED was defined as the control condition, because the visual and haptic information was congruent. During the scan, the participants were required to push the pad with the both hands while observing the reflection of the left hand and estimate the hardness of the pad perceived by the right (hidden) hand based on magnitude estimation. Behavioral results showed that the perceived hardness was significantly biased toward softer or harder in >73% of the trials in conditions SFT and HRD; we designated these trials as visually modulated (SFTvm and HRDvm, respectively). The accuracy map was calculated individually for each of the pair-wise comparisons of (SFTvm vs. MED), (HRDvm vs. MED), and (SFTvm vs. HRDvm) by a searchlight MVPA, and the cortical regions encoding the perceived hardness with visual modulation were identified by conjunction of the three accuracy maps in group analysis. The cluster was observed in the right sensory motor cortex, left anterior intraparietal sulcus (aIPS), bilateral parietal operculum (PO), and occipito-temporal cortex (OTC).

INTRODUCTION
Tactile texture perception provides essential information for not only recognition of objects (Klatzky et al., 1987) but also manipulation of objects (Johansson and Flanagan, 2009). Hardness is one of the fundamental dimensions that determine the perception of the tactile texture properties, along with roughness, stickiness, and temperature (Okamoto et al., 2013;Bensmaia, 2016). Previous studies have revealed that when one actively touches a deformable object with the hand, the compressional force to the fingers provides important information regarding the perceived hardness (softness) of the objects (Srinivasan and LaMotte, 1995;Friedman et al., 2008). This information is likely conveyed by the activity of slowly adapting type 1 (SA1) fibers, which innervate mechanoreceptors in superficial layers of the skin and respond to steady skin indentation with a sustained discharge (Srinivasan and LaMotte, 1996). Thus, it is clear that hardness perception is determined primarily by haptic signals from the periphery.
However, recent studies have revealed that hardness perception by active touch is also affected by visual information along with haptic inputs. Using virtual reality techniques, Kuschel et al. (2010) and Cellini et al. (2013) investigated the effect of vision on hardness (softness) perception. In these studies, visual information was presented via movies created by computer graphics, in which a deformable surface was pushed by a non-corporeal object, such as a ball or cylinder. Haptic information was provided to participants by the squeezing of two levers with simulated repulsive force (Kuschel et al., 2010) or by the touching of silicone specimens with the hand (Cellini et al., 2013). The indentation of the virtual surface on the display was controlled such that it was congruent or incongruent with that simulated by the haptic deformation. The results indicated that the perceived hardness in the incongruent condition was biased toward harder or softer when the indentation of the virtual surface on the display was smaller or larger, respectively. Furthermore, Punpongsanon et al. (2015) showed that the perceived softness is enhanced when participants push an object with their hand while observing an augmented (more indented) surface induced by a projection mapping technique. Together, these studies demonstrate how incongruent visual information can distort our perception of hardness.
These studies employed virtual reality techniques to experimentally manipulate the visual and haptic stimuli. In our previous study (Katsuyama et al., 2018), we attempted to examine the effect of visual information on hardness perception by active touch under a more natural condition. For this purpose, we used a mirror visual feedback (MVF) paradigm, in which the reflection of the hand appears to be the hand hidden behind the mirror. Participants were required to touch a polyurethane foam pad with both hands placed at the front and back of a mirror at the midline. To the hand at the front side of the mirror, one pad of three different levels of hardness (hard, medium, and soft) was presented in a given trial, while only the medium pad was presented to the hidden hand in all trials. The participants were required to estimate the hardness of the pad perceived by the hidden hand while observing the reflection. The result indicated that the perceived hardness assessed by magnitude estimation was significantly larger or smaller when participants observed the reflection of the hand pushing the harder or softer pad, respectively, at the same time as the hidden hand was pushing a pad of constant hardness. The modulation of the perceived hardness disappeared when participants touched the pads with both hands at different times or with their eyes closed, indicating that the visual modulation was not induced by the interaction of the motor commands for or haptic inputs from both hands in the brain. Further examination revealed that most participants utilized the visual cues regarding both the hand and the pad, such as the finger displacement and indentation of the surface of the pad, to infer the hardness. These results indicate that the modulation of hardness perception is induced by observation of the reflection of a hand pushing a pad of different hardness at the same time as the hidden hand was presented with a pad of constant hardness. Collectively, these findings indicate that hardness perception by active touch is affected by visual information along with haptic inputs.
This study aimed to investigate the cortical regions involved in hardness perception by active touch induced by multimodal integration. For this purpose, we conducted a functional magnetic resonance imaging (fMRI) experiment with a wholebrain searchlight multivoxel pattern analysis (MVPA) along with a conventional univariate analysis. However, it is difficult to distinguish the cortical activity involved in the perception of hardness from that involved in the motor commands for the finger movements. When one touches a deformable object with the hand, the compressional force to the fingers varies with the finger movements, particularly the finger displacement into the object. Since perceived hardness (softness) relies on the peak compressional force to the fingers (Srinivasan and LaMotte, 1995;Friedman et al., 2008), it is plausible that the cortical activity related to hardness perception and motor commands for finger movements may also co-vary, which allows them to be indistinguishable. This is particularly valid when applying the conventional subtraction technique and correlation analysis because these methods are based on the assumption that the intensity of cortical activity differs among cognitive conditions (Petersen et al., 1989;Friston et al., 1996;Amaro and Barker, 2006). To overcome this problem, we employed the hardness estimation task under the MVF that we used in our previous study. Our previous study showed that finger displacement of the hidden hand was constant across all conditions, which suggests that cortical activity related to the motor command and haptic input from the hidden hand should be constant as well while those involved in the perceived hardness cannot. This may allow us to successfully separate the cortical regions specific to the perceived hardness from those specific to the motor command.
Many studies have demonstrated that tactile illusions, including the rubber hand illusion, referred sensation, and MVF, are accompanied by illusory ownership toward the fake hands (Botvinick and Cohen, 1998;Ehrsson et al., 2004;Tsakiris et al., 2006;Bertamini et al., 2011;Bekrater-Bodmann et al., 2014). In particular, when the dummy and hidden hands are moving simultaneously, participants also experience the sensation that the dummy hand is moved by their own will; this is known as the sense of agency (Sanchez-Vives et al., 2010;Kalckert and Ehrsson, 2012;Jenkinson and Preston, 2015;Katsuyama et al., 2018). Therefore, we also investigated whether participants experienced the sense of hand ownership of the image of the hand in the mirror during the hardness estimation task with MVF.

Participants
Twenty-two healthy individuals (mean age 21.0 ± 0.64 years, eight females) participated in this study. They were all righthanded (handedness score = 89.6 ± 11.8, measured with the Edinburgh Handedness Inventory, Oldfield, 1971). They had normal or corrected to-normal vision. All participants provided written informed consent before the experiment. This study was approved by the local Ethical Committee for the Faculty of Dentistry of the Tokyo Medical and Dental University (D2016-051) Primate Research Institute of Kyoto University (2016-13), and Tamagawa University (N27-42), and corresponded to the Human Subjects Guidelines of the Declaration of Helsinki.

Apparatus
Three polyurethane foam pads (9 × 7 × 4 cm, length × width × height) of varying degrees of hardness were used (Inoac Corporation, Nagoya, Japan). They were labeled S (soft), M (medium), and H (hard) in increasing order of hardness. The physical hardness of each pad was calculated on the basis of the elastic modulus (S: 1871.7, M: 3434.6, H: 9803.9 N/cm 2 ). They were wrapped in a green cloth such that the participants were unable to judge the hardness of the pads based on the appearance.
To perform MVF lying in the scanner, an acrylic platform (Figures 1A,Ba) with a mirror (Figures 1A,Bb) on the top table was prepared. The platform was placed over the body of the participants while they were lying in the MR scanner. The left and right hands were placed at the front and back of the mirror, respectively. They could observe the reflection of the left hand on the platform in another mirror (Figure 1Ac) attached on the top of the head coil, while the right hand was hidden behind the mirror on the platform ( Figure 1B). The position, height, and angle of the mirror on the platform were adjusted individually such that the reflection of the left hand appeared to be the right hand. The position of both hands was also adjusted at the beginning of every run. Visual cues were presented on a flat monitor placed 5 m from the entrance of the bore of the scanner. Participants could observe the cues in the mirror on the head coil. The horizontal and vertical visual angle of the mirror was 36.2 and 25.0 degrees, respectively. All participants could observe the visual cues in their entirety and the mirror reflection of the left hand on the platform through the mirror. The oral responses of participants were recorded by a voice-capture system for MRI experiments (FOMRI TM II, Optoacoustics Ltd., Or Yehuda, Israel) via a microphone attached to the head coil (Figure 1Ad).

Pilot Study
We conducted a pilot study before the main experiment, to confirm whether the visual modulation of hardness perception can be induced in a supine position in the MR scanner (see Supplementary Material for details of the study). All apparatus, including the scanner, pads, and platform, were the same as those used in the main experiment. The participants performed the hardness estimation task while lying in the scanner. The procedure for the task was basically the same as that in the main experiment described below. Along with the three pad conditions (SFT, MED, and HRD), two timing conditions (SYNC and ASYNC) were employed. Participants touched the pad with both hands synchronously and asynchronously in the SYNC and ASYNC conditions, respectively, while observing the mirror reflection of the left hand. The results were consistent with our previous behavioral study (Katsuyama et al., 2018): the perceived hardness significantly increased as the hardness of the pad in the mirror reflection increased in the SYNC condition and remained constant in the ASYNC condition (Supplementary Figure 1). This result indicated that the visual modulation of hardness perception was induced even while the participant was lying in the MR scanner. To acquire enough number of trials with visual modulation for multivariate analysis, we used only the SYNC condition in the main experiment (see the Discussion section for details).

Main Experiment
During scanning runs, participants were required to perform a hardness estimation task. Figure 2 illustrates the schematic diagram of the task. The task was performed in a block design, and the scanner and a computer for the behavioral task were synchronized by a trigger signal from the scanner. Each run began with an initial rest (14 s). In the initial 10 s of the rest, participants were instructed to read Japanese characters presented on the monitor aloud for the sensitivity adjustment of the voice-capture system in the scanner noise. After the initial rest, an open square was presented on the monitor for 8 s (cue period). During the cue period, an experimenter beside the scanner presented a pad to each of the participants' hands. To the left hand, one of three pads (S, M, or H) was presented in a given trial, while only M was presented to the Participants could observe the reflection of the left hand through another mirror on the head coil (not shown), whereas the right hand was hidden behind the mirror on the platform.
FIGURE 2 | Schematic diagram of the hardness estimation task. The time sequence of one trial is depicted. An experimenter beside the scanner gave a pad to each of the participants' hands during the cue period. They were required to touch the pad with both hands simultaneously while watching the mirror reflection of the left hand during the touch period. In the following rest period, the experimenter removed the pads, and participants answered the hardness of the pad perceived by the right (hidden) hand orally. During the rest period, they were required to fixate the fixation target (a white cross).
right hand in all trials. Participants were instructed to place the middle three digits on the top surface of the pad and hold them. When the open square was replaced with a solid square on the monitor, the participants touched the pad with both hands while observing the reflection of the left hand for 6 s (touch period). The task was to estimate the hardness of the pad perceived with the right (hidden) hand during the period. Participants were told to investigate how hard (soft) the pad felt, not to estimate the physical hardness of the pad. At the end of the touch period, the experimenter removed the pads from the participants' hands, and they reported the estimated hardness based on magnitude estimation (Gescheider, 1997), such that the impression of the hardness of the pad matched the size of the positive number. Any positive numbers, including decimals and fractions, could be used for the estimation. Assuming participants would make an absolute ratio judgment of sensory magnitude (Gescheider, 1997), we did not use a standard stimulus (modulus) in the present study. Furthermore, no limitation was employed for the maxima and minima of the numbers. The responses were recorded and served as the offline analysis. After another rest period (6 s), the subsequent cue period followed. In all rest periods, the participants were asked to look at a white cross presented on the monitor in the mirror on the head coil. They were also instructed to put their hands on the platform in a rest position during the period. Before data collection began, participants received enough training to be able to make a stable hardness estimation while lying in the scanner.

Triplet-Based Presentation
In the following, the condition (trial) in which pads S, M, and H were presented to the participants' left hands is referred to as condition SFT, MED, and HRD, respectively. The order of the three conditions was arranged in a tripletbased manner: any one of conditions SFT, MED, and HRD was assigned to three successive trials in a pseudorandomized order (Figure 3). One run consisted of eight triplets, which resulted in 24 trials in a run (eight trials/condition). Six triplets 1 were arranged into overall runs and assigned to participants in a pseudorandomized order. Each participant completed eight runs (64 trials/condition/participant). FIGURE 3 | Schematic description of the triplet-based presentation of the stimulus. The tactile stimulation was presented such that any one of conditions SFT, MED, and HRD was assigned to three successive trials in a pseudorandomized order. The hardness ratio was calculated by dividing the perceived hardness in a trial of conditions SFT and HRD by that in condition MED (control) within a given triplet. For example, hardness ratio of the first (condition HRD) and third trial (condition SFT) in the triplet 1 is 2.0 (10/5) and 0.6 (3/5), respectively. Trials with a hardness ratio of <1 in condition SFT and >1 in condition HRD were defined as visually-modulated (colored in pink). In contrast, trials with a hardness ratio of ≥1 in condition SFT and ≤1 in condition HRD were defined as not-visually-modulated (colored in green). HRD, condition HRD; MED, condition MED; SFT, condition SFT; H, hard pad; M, medium pad; S, soft pad.

Behavioral Data
Because the visual and haptic information was incongruent in conditions SFT and HRD, but congruent in condition MED, we expected visual information to modulate hardness perception in conditions SFT and HRD, but not in condition MED. Therefore, condition MED was used as the control, and visual modulation was quantified by dividing the perceived hardness in a trial of conditions SFT and HRD by that in condition MED within a given triplet (the hardness ratio). A hardness ratio of <1 in condition SFT and >1 in condition HRD indicated that participants considered the pad to be softer and harder than that of the control (condition MED), respectively. We defined these trials as visually-modulated, and trials in condition SFT with a ratio of ≥1 and those in condition HRD with a ratio of ≤1 as not-visually-modulated (Figure 3). Finally, we also calculated the incidence of visuallymodulated trials in conditions SFT and HRD in every run (modulation ratio).

Multivariate analysis
Preprocessing was performed using the SPM8 software package (Wellcome Trust Centre for Neuroimaging) 2 . The volumes of the initial rest period were discarded. Motion correction was performed by realigning all volumes to the mean of the functional images for each participant followed by slice timing correction. Parameters for co-registration of anatomical images to the mean functional image were calculated and saved for each run. After preprocessing, we ran individual general linear model (GLM) analyses to estimate the cortical activity to the following six conditions per run: (1) SFTvm: visually-modulated trials in condition SFT, (2) SFTnm: not-visually-modulated trials in condition SFT, (3) HRDvm: visually-modulated trials in condition HRD, (4) HRDnm: not-visually-modulated trials in condition HRD, (5) MED: trials in condition MED, and (6) NULL: missed trials. The head motions were additionally modeled as "of no interest" conditions along with condition NULL.
An MVPA was carried out using The Decoding Toolbox (TDT; Hebart et al., 2015) and custom-written programs written in MATLAB (The MathWorks, Natick, MA, United States). A classification accuracy map was calculated for the pairwise comparisons of (SFTvm vs. MED), (HRDvm vs. MED), and (SFTvm vs. HRDvm) using the GLM parameter estimates for each of the three conditions SFTvm, HRDvm, and MED. A whole-brain searchlight analysis (Kriegeskorte et al., 2006) was performed using 9 mm-radius spheres centered around a given voxel for all voxels. A support vector machine (LIVSVM; Chang and Lin, 2011) was trained on seven out of the eight runs and tested on the untrained run, yielding a classification accuracy map. This process was performed iteratively until all runs had been tested (a leave-one-run-out cross-validation scheme), and a mean accuracy map was calculated by averaging the eight accuracy maps for each of the three comparisons.
Group analysis was performed using SPM8. The mean accuracy maps for the three comparisons were individually normalized to the standard Montreal Neurological Institute (MNI) brain template using the parameters obtained by the co-registration, and smoothed with an 8 mm full-width-halfmaximum (FWHM) isotropic Gaussian kernel. The normalized and smoothed maps were entered into a random-effect analysis to investigate the voxels with significantly higher accuracy than chance using a one-way within-subject ANOVA implemented in the second-level routine of SPM8. The resultant statistical map of T-value for the comparisons of (SFTvm vs. MED), (HRDvm vs. MED), and (SFTvm vs. HRDvm) was pseudocolored and overlaid on the structural image. To isolate brain regions encoding the perceived hardness with visual modulation, we performed conjunction analysis of the three T-value maps against the conjunction null hypothesis (Nichols et al., 2005). This analysis seeks the significant voxels present in all three maps. For the group analysis, a statistical threshold of p < 0.05 corrected for familywise error rate was applied. No extent threshold was used. The peak coordinates of the resultant clusters were represented in the MNI coordinate system. The locations of the clusters were identified using an anatomical toolbox for SPM8 (SPM Anatomy Toolbox v. 2.2b, Eickhoff et al., 2005) and a human brain atlas (Talairach and Tournoux, 1988).

Univariate analysis
All analyses were performed using SPM8. The preprocessing was conducted as follows: Motion correction by realigning all volumes to the mean functional image, slice timing correction, co-registering of the anatomical image to the mean functional image, normalization of the all functional images to the MNI template, and spatial smoothing with an 8 mm FWHM isotropic Gaussian kernel.
After preprocessing, we ran individual GLM analyses using the same conditions as in the multivariate analysis (SFTvm, SFTnm, HRDvm, HRDnm, MED, NULL, and the realigned parameters) to estimate the cortical activity per condition. For the group analysis, contrast images of conditions SFTvm, HRDvm, and MED were entered into second-level random-effect analysis by a one-sample t-test on a voxel-tovoxel basis, and the cortical regions activated by contrasts (SFTvm > MED), (HRDvm > MED), (SFTvm > HRDvm), and (HRDvm > SFTvm) were investigated. We also performed a correlation analysis. We then ran further individual GLM analyses using the following parameters as the parametric modulation: (1) The physical hardness of the pad. This analysis was conducted to identify the cortical regions that exhibited activity changes with varying hardness of the pad regardless of visual modulation. Two conditions, PAD and NULL, were employed. The condition PAD indicates the sum of conditions SFTvm, SFTnm, HRDvm, HRDnm, and MED. The physical hardness of the three pads described earlier was entered as the parametric modulation of all trials of condition PAD, and the contrast images were used for the population analysis; (2) the raw values of the perceived hardness; and (3) the hardness ratio. These analyses were performed to identify the cortical regions that exhibited activity that correlated with the perceived hardness with visual modulation. Four conditions, VM, NM, MED, and NULL, were employed.
Conditions VM and NM indicated the sum of conditions SFTvm and HRDvm, and SFTnm and HRDnm, respectively. The behavioral parameters were entered as the parametric modulation. A population analysis was conducted using the contrast images of condition VM. A statistical threshold of p < 1 × 10 −5 , uncorrected for multiple comparisons, was applied for the group analysis. An extent threshold >50 voxels was used.

Hand Ownership and Agency Test
After MRI, a questionnaire regarding the ownership and agency of the mirror image of the hand was administered. The items were invented with reference to previous studies (Botvinick and Cohen, 1998;Longo et al., 2009;Kalckert and Ehrsson, 2014;Kokkinara and Slater, 2014;Jenkinson and Preston, 2015;Ma and Hommel, 2015).
The items were as follows: The questionnaire items were translated into Japanese. The participants were required to rate their agreement with each item on a seven-point scale ranging from −3 (strongly disagree) through 0 (neither agree nor disagree) to +3 (strongly agree). For items 3 and 6, a different scale ranging from −100 (strongly left hand) through 0 (equally left and right hands) to +100 (strongly right hand) was used. Any intermediate value could be used. Agreement or disagreement was tested by comparing the mean score to 0 using a one sample t-test at a significance level of p < 0.05. Furthermore, the participants answered the following questions in a postexperimental interview: 1. Did you know how many different pads were presented to each of your hands in all runs? 2. Which hand was presented with a wider range of pads of varying hardness; your right hand or your left hand? 3. Did you find that the physical hardness of the pad presented to your right hand was constant throughout all runs?

Results of the Hardness Estimation Test During Scanning
In the present study, condition MED was used as a control to test for visual modulation in conditions SFT and HRD. Thus, the ratios of perceived hardness SFT/MED and HRD/MED were calculated for each triplet (hardness ratio) and used to determine whether the perceived hardness was modulated by visual information in a given trial. Figure 4 illustrates the distribution of the hardness ratio in conditions SFT and HRD across all participants and runs. Trials of condition SFT with hardness ratios of <1 indicated that the pad was perceived to be softer than the control. Likewise, trials of condition HRD with ratios >1 indicated that the pad was perceived to be harder than the control. We then defined these trials as visuallymodulated (SFTvm and HRDvm). The ratio of the visuallymodulated trials in conditions SFT and HRD (modulation ratio) averaged over all eight runs for each participant is provided in Table 1. The mean modulation ratio was 73.4% ± 17.7% and 74.0% ± 18.2% for conditions SFT and HRD, respectively. There was no significant difference between these ratios [Welch's t-test, t (42) = 0.11, p = 0.91]. Furthermore, we calculated the modulation ratio of the SFT and HRD conditions for each of the six triplets overall runs and the participants with two-way repeated measures ANOVA, for computing the main effects of the condition (SFT vs. HRD) and the triplet (6 triplets), and found no significant difference between the ratios [F (1, 21) = 0.08, p = 0.78 for the main effect of the condition; F (5, 105) = 1.73, p = 0.13 for the main effect of the triplet; and F (5, 105) = 0.83, p = 0.53 for the interaction between the main effects]. The modulation ratios for each participant, condition, and run are provided in Supplementary Table 1.
The post-experimental interview revealed that most participants did not find that the physical hardness of the pad presented to the hidden hand was constant during scan, until they were informed of the "gimmick" of the task after the MRI experiment. However, two participants were vaguely aware of the possibility. The number of pads of different degrees of hardness presented to the left and right hand was 6.3 ± 2.4 and 6.0 ± 3.5, respectively. Although there was

Results of the fMRI Multivariate Analysis
Cortical regions encoding the perception that the pad is softer than the control were investigated by multivariate analysis using the pair-wise comparison of (SFTvm vs. MED) ( Figure 5A). A large cluster was observed in the sensory motor cortex (SMC) extending from the precentral gyrus to the intraparietal sulcus FIGURE 4 | Distribution of the hardness ratio across all participants and runs. The visually-modulated trials in conditions SFT and HRD were defined as conditions SFTvm and HRDvm, respectively (colored in pink). In contrast, the not-visually-modulated trials in conditions SFT and HRD were defined as conditions SFTnm and HRDnm, respectively (colored in green). The surface view of the cortical regions encoding the perception that the pad is softer than the control condition. The cortical regions whose activity patterns were different between the conditions SFTvm and MED were investigated by multivariate analysis. (B) Brain regions encoding the perception that the pad is harder than the control condition. The areas whose activity patterns were different between the conditions HRDvm and MED were investigated. (C) Cortical regions encoding the perceived hardness. Brain areas whose activity patterns were different between the conditions SFTvm and HRDvm were investigated. The lateral view of the left and right hemisphere is shown in the left and right column, respectively, in all panels. The significance level was set at p < 0.05 corrected for familywise error rate. No extent threshold was applied. The color bar indicates the T-value. A, anterior; P, posterior.
(IPS), parietal operculum (PO), insula, and occipito-temporal cortex (OTC) of the right hemisphere. In the left hemisphere, a cluster was observed in the anterior IPS (aIPS), PO, and OTC. There were no clusters found in the left SMC. Similarly, the brain regions encoding the perception that the pad is harder than the control were investigated by the comparison of (HRDvm vs. MED) (Figure 5B). A large cluster was observed in the SMC, PO and visual areas in the occipital cortex of the right hemisphere. In the left hemisphere, a cluster was found in the anterior insula, aIPS, PO, and OTC. There were no clusters obtained in the left SMC either.
Finally, differences in perceived hardness between SFT and HRD conditions were investigated by the comparison of (SFTvm vs. HRDvm) (Figure 5C). In this analysis, a large cluster covering the motor, sensory, and visual cortices was observed in the right hemisphere. Similarly, a large cluster was obtained in the left hemisphere as well, however, only a few clusters were found in the left motor cortex.
Next, we investigated the cortical regions encoding the perceived hardness with visual modulation by conjunction analysis of the three comparisons of (SFTvm vs. MED), (HRDvm vs. MED), and (SFTvm vs. HRDvm). Figure 6 shows the result. A large cluster extending from the posterior portion of the precentral gyrus to the aIPS through the postcentral gyrus and sulcus was observed in SMC of the right hemisphere (Figures 6B,C). The local maxima of the cluster (x, y, z = 54, -24, 44) was located in the postcentral gyrus ( Figure 6B). In the left hemisphere, no cluster was observed in the SMC. However, a cluster was identified in the aIPS (Figure 6C). Another cluster was obtained in the posterior PO ( Figure 6D) and OTC (Figure 6E) of both hemispheres. The cluster in the OTC was located in the posterior portion of the inferior temporal sulcus. It is noteworthy that the cluster in the aIPS, PO, and OTC was symmetrical in both hemispheres. The location, size, Brodmann's area (BA), MNI coordinates, T-value, and mean accuracy over all participants of the clusters are indicated in Table 2.

Results of the fMRI Univariate Analysis
Univariate analysis was performed using the same data set as was used in the multivariate analysis. First, we analyzed the contrast images of (SFTvm > MED) and (HRDvm > MED) to identify the cortical regions that were more strongly activated in the visually-modulated conditions than in the control. However, we did not observe any activation at the threshold in the population analysis. No activation was observed in the contrast images of (SFTvm > HRDvm) and (HRDvm > SFTvm) either. We then applied a correlation analysis to investigate the cortical regions that exhibited altered activity with varying physical hardness of the pad irrespective of visual modulation. Some participants exhibited significant activation with a BOLD signal that was negatively correlated with the hardness of the pad in the right SMC. However, no significant cluster was obtained in the population analysis of all participants. Finally, we attempted to identify the cortical regions wherein the activity correlated with the perceived hardness modulated by visual information, by entering the raw values of the perceived hardness or hardness ratio as the parametric modulation, but found no significant activation. Thus, we did not obtain any significant activation related to the physical or perceived hardness of the pad in the univariate analysis in the present study.

Results of the Hand Ownership and Agency Test
After the fMRI experiment, a questionnaire regarding illusory hand ownership and sense of agency toward the hand reflection was administered (Figure 7). Items 1-3 were related to illusory ownership. The results revealed that the vast majority of the participants felt that the reflection of the left hand was their own right hand. Items 4-6 were related to the sense of agency. Participants reported that the motion of the hand in the reflection was the motion intended for their own hand. However, the score of item 6, which asked whether the motion of the reflected hand was the one they intended to make with their right or left hand, was biased toward the right hand, but not significantly. The other items (7-11) were the control questions to items 1-6. However, the participants awarded significantly negative scores to all items with the exception of item 10. These results indicate that, although the participants felt that the hand in the reflection was their own right hand and the motion was caused by their own will, it was not necessarily the motion they intended to make with the right hand.

Hardness Perception Modulated by Visual Information
The aim of the present study was to identify the cortical regions involved in hardness perception modulated by visual information. The participants performed the hardness estimation task using MVF while lying in an MR scanner. They were required to estimate the hardness of a pad presented to the right hand that was placed behind the mirror while observing the reflection of the left hand pushing a pad simultaneously. To the left hand, one of three pads (soft, medium, and hard) was presented in a given trial (condition SFT, MED, and HRD, respectively), whereas only the medium pad was presented to the right hand across all trials. The behavioral result indicated that the participants felt that the pad presented to the hidden hand was softer or harder than the control condition (MED) in >73 and 74% of the trials in conditions SFT and HRD, respectively. This result is consistent with that of our previous study demonstrating that the modulation of hardness perception was induced by observation of the reflection of a hand pushing a pad of different hardness. In the study, we excluded the possibility that the modulation was induced by sensory assimilation and bimanual coordination by demonstrating that the modulation diminished when participants touched the pad with both hands at different times or with their eyes closed. Furthermore, no significant difference was observed in the finger displacement of the hidden hand during touch (Katsuyama et al., 2018). In this study, we excluded control conditions to validate the effect of sensory assimilation and bimanual coordination on the behavioral task, to acquire a sufficient number of trials with visual modulation for multivariate analysis, which requires multiple runs (in the present study, participants underwent 192 trials, over eight runs in the supine position). The introduction of control conditions simply doubles the number of total trials. This might cause deterioration of the quality of data due to fatigue. Furthermore, it was impossible to measure the finger movements in the MRI scanner using the electromagnetic tracker used in the previous study. However, we considered the effects of sensory assimilation and bimanual coordination negligible, based on the results of the pilot study (Supplementary Material). The study was conducted to examine whether visual modulation can be induced in a supine position in an MR scanner. Based on our previous study, two touch timing conditions (ASYNC and SYNC) were employed, along with the pad-hardness condition. Since earlier studies have shown that sensory assimilation can be induced when two stimuli are presented to participants sequentially as well as concurrently (Gescheider, 1997;Roberts, 2013), if the modulation of hardness perception was induced by the assimilation of tactile input from both hands, modulation should have been observed in the ASYNC condition also, in which participants touched the pads of varying hardness with both hands in an alternating fashion.
However, similar to our previous study, modulation was induced in the SYNC condition, but not in the ASYNC condition, in the pilot study. Although we lack explicit evidence, we believe that it is unlikely that participants used different strategies and/or moved the fingers in a different manner to reproduce the same behavioral result in the upright and supine positions. Therefore, we concluded that the modulation of perceived hardness in the present study was induced by observation of the reflection, even in the supine position.

Result of the Univariate Analysis of the fMRI Data
In the present study, we analyzed the same dataset with a uniand multivariate approach. First, we investigated whether there were any cortical regions where the activity was higher in the visually-modulated conditions than in the control condition by subtraction analysis. However, we did not identify any cortical regions. Next, we tried to investigate cortical regions that exhibit activity that correlates with the physical hardness of the pad by correlation analysis. However, no such region was found, either.
In the individual analysis, we observed significant activity in the right SMC that was negatively correlated with the hardness of the pad in some participants. Since the finger displacement to push the pad increased as the physical hardness of the pad decreased, it is plausible that the activity reflected the difference in the motor commands for the finger movements and the sensory inputs to the left hand when touching pads of varying hardness during the scan. However, the activity in the right SMC did not exceed the threshold in the group analysis. To the contrary, as detailed below, the informative clusters were obtained in several cortical regions including the right SMC by multivariate analysis. The reason why no activity was found by univariate analysis is not straightforward, because little is known about how hardness (softness) perception is encoded by neurons in the somatosensory cortex. The firing rate of SA1 fibers, which convey the information of the compressional force to the fingers, increases as the indentation of the skin increases (Mountcastle et al., 1966;Harrington and Merzenich, 1970;Muniak et al., 2007). This response property is kept in some S1 neurons (Bensmaia et al., 2008), and may be detectable by univariate analysis. However, cumulative evidence reveals that information processing in S1 is more complicated than thought. Convergence of the information conveyed by different mechanoreceptor afferents including SA1, rapidly adapting (RA) 3 , and The inconsistency in the activation patterns identified by uniand multivariate approaches was also reported for stickiness perception Yeon et al., 2017). As described earlier, hardness perception is determined by the compressional force to the skin, and the information is conveyed by SA1 fibers from periphery (Srinivasan and LaMotte, 1995;Friedman et al., 2008). Likewise, the sense of stickiness is evoked when the skin is stretched by an adhesive object, and the information about the skin deformation is conveyed by several mechanoreceptor afferents including SA1, slowly adapting type 2 (SA2) 5 , and RA fibers . How sensory perception, such as hardness and stickiness, is encoded by cortical neurons and how the neuronal representation is converted to a BOLD signal utilized by fMRI are important problems in the study of the sensory system.

Result of the Multivariate Analysis of the fMRI Data: The Right Hemisphere
In the present study, we investigated cortical regions involved in the perceived hardness modulated by visual information by using a whole-brain searchlight MVPA (Kriegeskorte et al., 2006). For this purpose, the accuracy map of pair-wise comparisons of (SFTvm vs. MED), (HRDvm vs. MED), and (SFTvm vs. HRDvm) were calculated individually, and the maps were used for conjunction analysis in the group analysis. This analysis seeks the significant voxels present in all three maps. As described earlier, the condition MED was designated as the control to conditions SFTvm and HRDvm for the visual modulation. Therefore, the T-value map of the pair-wise comparison of (SFTvm vs. MED) -showing brain regions whose activity patterns are significantly different between conditions SFTvm and MED -indicates cortical regions encoding the perception that the pad is softer than the control in the group analysis. Likewise, the map of (HRDvm vs. MED) illustrates the brain regions encoding the perception that the pad is harder than the control. On the other hand, the T-value map of (SFTvm vs. HRDvm) shows the brain regions sensitive to the difference in the perceived hardness. In the present study, we defined such cortical regions as areas encoding the perceived hardness. However, the map also includes those cortical regions used in perceiving hardness without visual input. Therefore, to isolate brain regions involved in integrating visual inputs to perceive hardness from those that do not, we performed conjunction analysis of the three T-value maps of (SFTvm vs. MED), (HRDvm vs. MED), and (SFTvm vs. HRDvm).
One may argue that the individual accuracy map should be calculated for the pair-wise comparisons of (SFTvm vs. SFTnm) and (HRDvm vs. HRDnm), because the perceived hardness was different between the conditions, whereas the physical stimuli were the same. However, this analysis proved difficult. The modulation ratio was more than 73% in both conditions SFT and HRD, and the modulated trials of both conditions (SFTvm and HRDvm) were contained in all runs and participants. On the other hand, the trials of conditions SFTnm and HRDnm were not in all runs: the mean number of runs with trials of conditions SFTnm and HRDnm was 6.2 and 6.0 out of 8 runs (range: 0-8 for both conditions) per participant, respectively. If the SVM is trained on a condition pair with an unbalanced number of runs, the labeling by the SVM in the untrained run will be biased toward the condition with the most runs (SFTvm and HRDvm in this case), resulting in deterioration of classification accuracy. Therefore, we restricted our analysis to the comparisons of (SFTvm vs. MED), (HRDvm vs. MED), and (SFTvm vs. HRDvm) in the present study.
The searchlight MVPA revealed a large cluster in the SMC of the right hemisphere. The peak voxel of the region was in the postcentral gyrus (S1), and it extended into the precentral gyrus (M1) rostrally and to IPS caudally through the postcentral sulcus, and to the PO ventrally. In the hardness estimation task, the participants touched pads of varying physical hardness with the left hand. Therefore, it is clear that the large cluster in the right hemisphere reflected cortical activity related to the motor command of the finger movements and sensory response to the tactile and proprioceptive inputs of the left hand. The cluster extended from the postcentral gyrus to the IPS via the postcentral sulcus, covering Brodmann's area 1 and 2. Neurons in Brodmann's area 1 and 2 in monkeys represent several fingers, and the receptive fields were overlapped (Iwamura et al., 1980(Iwamura et al., , 1985a. This observation was also confirmed in the human somatosensory cortex (Kurth et al., 2000). These characteristics of the receptive fields of neurons in Brodmann's area 1 and 2 form a functional surface suitable for various aspects of active touch, such as touching, pinching, and grasping (Iwamura et al., 1985b). As described earlier, when one pushes objects of varying hardness, finger displacement decreases as the physical hardness of the object increases, and the peak compressional force to fingers and perceived hardness would also co-vary with the displacement. Thus, the brain activity pattern obtained in the right hemisphere was depictive of the cortical activity related to the motor commands and somatosensory responses of the left hand in the hardness estimation task.
Contrary to the findings in the right hemisphere, no clusters were found on M1 and S1 of the left hemisphere. Considering that the finger displacement of the right hand was constant across all trials, this result was not unexpected. One possible disadvantage of the task design of the present study is that it may have underestimated the role of S1 in hardness perception. Although no cluster was obtained on the left S1, we cannot exclude the possibility that S1 may play a role in the perception of hardness with visual modulation, because recent studies have suggested that S1 is involved in illusory hand ownership and tactile perception in humans and animals (Chen et al., 2003;Blankenburg et al., 2006;Schaefer et al., 2006;Hötting et al., 2009;Shokur et al., 2013). In the present study, participants actually touched a pad with their right hand. The continuous tactile input from the periphery might mask the illusion-related signal in S1 and make it difficult to detect by either univariate or multivariate analysis. The absence of activity in S1 in illusory ownership and tactile perception has been observed in previous fMRI studies, in which participants received tactile stimulation in both test and control conditions (for example, Ehrsson et al., 2004;Matthys et al., 2009;Chambon et al., 2013;Kitada et al., 2014a;Limanowski et al., 2014). However, when a tactile illusion is induced on a part of the body, by tactile stimulation delivered to another part of the body and at different times, the illusionrelated signal is detected on the somatotopically corresponding site at S1 by fMRI (Blankenburg et al., 2006).
Although sensitivity to S1 might be low, considering the difference in activation patterns between the right and left SMC, we concluded that the cortical regions encoding the perception of hardness with visual modulation can be successfully isolated from those involved in the motor commands, using the hardness-estimation task and multivariate analysis as described in the Introduction.

Result of the Multivariate Analysis of the fMRI Data: The Bilateral aIPS
A cluster was found in the aIPS of the left hemisphere. It is well known that the IPS plays a crucial role in multimodal integration. Recent neuroimaging studies have demonstrated that visual, auditory, somatosensory, and vestibular inputs converge in the IPS of both humans and macaque monkeys (Guipponi et al., 2013;Huang and Sereno, 2018). In monkeys, it has been revealed that neurons in the bottom of the intraparietal sulcus (ventral intraparietal area, area VIP) respond to both visual and tactile stimuli (Colby and Duhamel, 1991;Duhamel et al., 1998), and that the activity of some bimodal neurons in VIP is higher when the visual and tactile stimulation was presented at the same time than when each stimulus was presented individually (Avillac et al., 2007). It is noteworthy that the cluster in the left aIPS obtained in the present study was located near the bottom of the sulcus. In humans, many studies have demonstrated that the aIPS is commonly activated by visual and haptic cues in shape and texture recognition tasks (Macaluso and Driver, 2001;Kitada et al., 2006;Stilla and Sathian, 2008;Gentile et al., 2011;Joanne Jao et al., 2014). The activity of the aIPS is enhanced by visuo-tactile stimulation that is more than the sum of the visual and tactile stimulation alone (Gentile et al., 2011). The activity of the aIPS is also observed in object discrimination in a delayed-matching-to-sample task in which the sample and test stimuli were presented visually and haptically, respectively, and vice versa (Grefkes et al., 2002;Kassuba et al., 2013). Pasalar et al. (2010) reported that application of transcranial magnetic stimulation (TMS) to the aIPS disrupts visuo-tactile integration. They found that the detection of tactile stimulation in a finger is enhanced by observation of a visual cue pointing to the same finger. However, the enhancement is eliminated by the application of TMS to the junction of the postcentral sulcus and the aIPS.
The aIPS and posterior parietal cortex (PPC) is also known as a higher order somatosensory cortex. Activation of the aIPS is observed in haptic object (shape) recognition tasks (Pietrini et al., 2004;Hinkley et al., 2009;Podrebarac et al., 2014;Hernández-Pérez et al., 2017) and vibration stimulation to a finger (Schmidt et al., 2014). With regard to the present study, Leib et al. (2016) suggested that stiffness perception may be represented in the PPC. In their study, participants estimated the stiffness of a surface projected on a display by "pushing" the virtual surface with the hand via a haptic device. They applied thetaburst TMS to the PPC of participants and demonstrated that the stiffness perception was disrupted by TMS but the hand movements were not. Moreover, Newman et al. (2005) reported that imagery of the material dimension of objects (roughness, hardness, and temperature) activated the posterior IPS. Finally, it is noteworthy that the cluster was observed in the corresponding site of the aIPS in both hemispheres in the present study. In the hardness estimation task, participants actually touched the pads of varying physical hardness with the left hand, and touched the pad of constant hardness with the right hand. Nevertheless, the perceived hardness of the pad presented to the right hand was modulated by visual information. Thus, the present result that the cluster was located in the corresponding site of the aIPS suggests that this region may be involved in the perception of hardness. Taken together, these results indicate that the aIPS plays a crucial role in not only the integration of visual and haptic information, but also in hardness perception in the present study.

Result of the Multivariate Analysis of the fMRI Data: The Bilateral PO
A cluster was observed in the PO of both hemispheres. This region is known to be the anatomical site of somatosensory area 2 (S2). Human PO is divided into four cytoarchitectonically distinct sub-regions termed OP1 -4. On the basis of the somatotopic arrangement in the lateral fissure, Eickhoff et al. (2006) suggested that areas OP1, OP3, and OP4 may be the human homologs to non-human primate areas S2, ventral somatosensory (VS; Cusick et al., 1989), and parietal ventral (PV; Krubitzer et al., 1995), respectively. The cluster obtained in the present study was located in the posterior and lateral portion of PO in both hemispheres. On the basis of the topography of the four sub-regions in PO (Eickhoff et al., 2006), the cluster may be in OP1 (S2).
Several studies have shown that human PO is involved in somatosensory processing. Electrical stimulation to the Sylvian fissure evokes somatosensory sensations (Penfield and Jasper, 1954;Lüders et al., 1985). Clinical studies have reported that a lesion in PO impairs somatic sensation including pain (Caselli, 1991;Bowsher et al., 2004). By using an fMRI technique, Servos et al. (2001) also found that PO is activated by the hardness identification task, although the location was more anterior to the cluster obtained in the present study. Furthermore, activation of PO has been observed by tactile or electrical stimulation to the hand via fMRI (Disbrow et al., 2000;Ruben et al., 2001;Eickhoff et al., 2007;Hinkley et al., 2007;Burton et al., 2008;Mazzola et al., 2012;Tamè et al., 2012;Björnsdotter et al., 2014;Morrison, 2016), EEG (Stančák et al., 2005), and MEG (Disbrow et al., 2001;Lin and Forss, 2002) studies.
It has been suggested that PO (S2) is a higher order area than S1. Physiological studies with monkeys have shown that neurons in S2 have more complex response properties than S1: they have a larger receptive field and respond to both contra-and ipsilateral stimulation; some S2 neurons respond to stimulation to multiple digits; about one-fourth of S2 neurons have their receptive field in more than one of the four body parts (head, trunk, forelimb, and hindlimb); finally, some S2 neurons are silent to passive stimulation but active in active touch task (Sinclair and Burton, 1993;Fitzgerald et al., 2006a,b;Taoka et al., 2016). Bilateral response has been also confirmed in human PO by fMRI and EEG/MEG studies (Disbrow et al., 2001;Eickhoff et al., 2007;Burton et al., 2008;Mazzola et al., 2012;Tamè et al., 2012). Anatomical studies have shown that S2 has connections with the insula, PPC (area 7), motor cortex, and prefrontal cortex in monkeys (Pandya and Kuypers, 1969;Mufson et al., 1981;Friedman et al., 1986;Matelli et al., 1986;Carmichael and Price, 1995). Functional imaging studies have shown that PO has anatomical and functional connectivity with these cortical regions in humans as well (Eickhoff et al., 2010;Mȃlîia et al., 2018). Binkofski et al. (1999) revealed that PO is involved in object recognition by hand manipulation. In the study, activation of PO was observed without activation of S1 when participants manipulated complex objects. It is noteworthy that a cluster was observed in the left PO while no clusters were found in S1 during the hardness estimation task by active touch in the present study. Activation of PO was also observed in the haptic texture identification task (Stilla and Sathian, 2008). PO is also engaged during a tactile illusion: when the tendon of the wrist extensor muscle is stimulated by vibration, an object touched by the stimulated hand is perceived to move with the hand without any real motion of the hand (hand-object illusion) (Naito and Ehrsson, 2006). Furthermore, PO, but not S1, is activated while watching a movie of someone being touched on the leg. Interestingly, the same region in PO is also activated by tactile stimulation to the same leg of participants (Keysers et al., 2004). The activation observed in the study by Keysers et al. was located close to the cluster obtained in the present study. All these observations indicate that human PO is involved in somatosensory perception, and this region may play a key role in perceiving hardness with visual modulation in the present study.
Considering that PO is sensitive to not only the contralateral but also ipsilateral tactile stimulation, one may argue that the cluster in the left PO may be related to the haptic input from the left hand. This is a valid point, and we do not exclude this possibility. However, we suggest that the region in the left PO also conveys substantial information on the subjective hardness perceived by the right hand for the following reasons. First, although PO is sensitive to tactile stimulation to either side, the response is stronger to contralateral than ipsilateral stimulation in most cases (Burton et al., 1999(Burton et al., , 2008Disbrow et al., 2000Disbrow et al., , 2001Ruben et al., 2001;Servos et al., 2001;Stančák et al., 2005;Eickhoff et al., 2007;Mazzola et al., 2012). Second, a large cluster was observed in the right SMC, reflecting that participants actually touched the pads of varying physical hardness with the left hand throughout all trials. To the contrary, no cluster was observed in the left S1, suggesting that very little information is sent from S1 to PO in the left hemisphere. These results predict that more information about the perceived hardness should be processed in the right PO than the left. However, there was no difference in the size and mean accuracy of the cluster in PO between the left and right hemispheres. The most plausible explanation for this discrepancy is that a certain amount of information about the perceived hardness by the right hand is processed in the left PO. It has been shown that human PO has both anatomical and functional connectivity with aIPS (Eickhoff et al., 2010). Taken together with the earlier discussion about the aIPS, it seems reasonable to suppose that the haptic input from the right hand and visual information are integrated in the left aIPS, and the perceived hardness modulated by the visual information is processed in the cortical network between aIPS and PO.

Result of the Multivariate Analysis of the fMRI Data: The Bilateral OTC
We also observed a cluster in the bilateral OTC. It was located in the inferior temporal sulcus, and the location was symmetrical in both hemispheres. By comparison of the coordinates with those from previous studies, we propose that the region is in the extrastriate body area (EBA). EBA was initially described as one of the extrastriate visual areas that selectively responds to images of body parts except the face (Downing et al., 2001). The mean coordinates of the local maxima of EBA averaged over the previous 11 studies were (x, y, z) = (−47.1 ± 2.6, −72.2 ± 2.7, −3.4 ± 6.0), and (47.1 ± 3.1, −67.9 ± 3.0, −3.1 ± 5.1) in left and right hemisphere, respectively (based on data in Ferri et al., 2013). In comparison, our cluster was located more anteriorly to the classical EBA at the local maxima (left: −50, −62, −4; right: 50, −60, −4). However, the extent of EBA identified by the localizer task (based on the method by Downing et al., 2001) is sometimes larger, covering the anterior occipital to posterior temporal cortices (for example, Spiridon et al., 2006;Orlov et al., 2010;Costantini et al., 2011;Kitada et al., 2014b). Moreover, studies on the topographical organization of EBA showed that the hand is represented anterior to the other body parts, such as the foot and trunk (Astafiev et al., 2004;Orlov et al., 2010). Therefore, the cluster in the OTC may be included in, or overlap with, the anterior portion of EBA that represents the hands somatotopically.
Previous studies have shown that EBA is one of the cortical regions related to illusory hand (body) ownership (Matthys et al., 2009;Limanowski et al., 2014;Limanowski and Blankenburg, 2015), along with the ventral premotor cortex (vPM), right temporo-parietal junction (rTPJ), and PPC (Ehrsson et al., 2004(Ehrsson et al., , 2005Arzy et al., 2006;Tsakiris et al., 2008;Gentile et al., 2011;Ionta et al., 2011;Brozzoli et al., 2012). Using MVF, Matthys et al. (2009) found activation in the OTC when participants observed a mirror reflection of one hand while the other hand was hidden behind the mirror. Because participants felt the mirror reflection was their right hand in the present study, the cluster in the OTC may be related to the illusory hand ownership.
The EBA is also positioned in the cortical network for visuomotor interactions; it processes visual information of one's own body, such as the position, posture, and movements, for actions (Astafiev et al., 2004;Bracci et al., 2012;Orlov et al., 2014;Zimmermann et al., 2018). In the hardness estimation task of the present study, the displacement, speed, and appearance of the strain of the fingers to touch the pad were different between all three conditions, which was detectable by multivariate analysis in the pair-wise comparisons. Furthermore, our previous study showed that participants actually used the visual cues on the finger movements for the hardness estimation along with those on the surface deformation of the pad (Katsuyama et al., 2018). Therefore, the cluster in OTC may also be engaged in processing visual information about finger movements during the hardness estimation task.

Illusory Ownership and Agency of the Reflection Hand
The result of the hand ownership test indicated that the participants felt as if the reflection of the left hand was the right hand behind the mirror. Illusory ownership toward a reflection of a hand has been previously reported (Longo et al., 2009;Mancini et al., 2011;Takasugi et al., 2011;Hoermann et al., 2012;Jenkinson and Preston, 2015). The present result is consistent with such findings. However, the result of the sense of agency test was somewhat different to that of the illusory hand ownership test. The result indicated that the participants felt that the motion of the hand in the reflection was caused by their own will. However, the answer to the question of whether the hand they intended to move was the left or right hand was biased toward the right, but not significantly. These results indicate that participants felt as if the hand in the reflection was their right hand, but the motion was not necessarily the same as that intended for the right hand. Using the moving rubber hand illusion, a previous study revealed that ownership and agency toward a dummy hand is dissociable (Kalckert and Ehrsson, 2012). In the study, participants observed a dummy hand tapping in synchrony with their hidden hand (the dummy and hidden hands were mechanically connected). When the hidden and dummy hands were moved passively in synchrony by an experimenter, the participants experienced ownership but not agency of the dummy hand. However, when the participants observed the dummy hand moving in synchrony with their hidden hand, the participants experienced agency but not ownership of the dummy hand, if the spatial position of the dummy hand was anatomically incongruent. In the present study, some participants reported that they experienced a peculiar feeling during the scan when they realized that the finger displacements observed in the mirror reflection appeared to differ from those perceived by the proprioceptive sensation. The dissociation of the observed and experienced finger displacements may have been due to the diminished agency of the reflection. Previous studies have revealed that cortical regions in the prefrontal and parietal cortices related to intention and monitoring of actions, such as pre-supplementary motor area (pre-SMA), dorsolateral prefrontal cortex (DLPFC), TPJ, and PPC are also involved in the sense of agency (Tsakiris et al., 2010;Yomogida et al., 2010;Sperduti et al., 2011;Chambon et al., 2013;Kühn et al., 2013). The absence of clusters in these cortical regions in the present study may be due to the imperfections in the agency illusion.
As described earlier, the illusory ownership of dummy hands is encoded in several cortical regions, such as vPM, rTPJ, PPC, and EBA. In the present study, although the participants experienced significant ownership of the hand reflection, no clusters were found in these cortical areas, except for the EBA. One possible reason for this is that the object to which participants experienced the ownership of the right hand was the mirror reflection of the left hand in all three conditions (SFTvm, HRDvm, and MED), making it difficult to detect differences in the pattern of cortical activation between the conditions with respect to the illusory hand ownership. Another reason would be that the illusory hand ownership observed in the present study was not accompanied by a change in the proprioceptive sensation of the hand position. In previous studies, the dummy hand was placed at an anatomically plausible, but apparently different position from that of participant's real (hidden) hand. After application of synchronous tactile stimulation to the real and dummy hands, the participants felt as though their real hand was in the same position as the dummy hand. The relative displacement between the observed and experienced position of the hands is defined as proprioceptive drift, and has been used as an objective measure of ownership of dummy hands (Botvinick and Cohen, 1998;Holmes and Spence, 2005;Tsakiris and Haggard, 2005). The previous functional imaging studies showing the involvement of vPM, PPC, rTPJ, and EBA in illusory hand ownership employed experimental paradigms including proprioceptive drift (Ehrsson et al., 2004(Ehrsson et al., , 2005Tsakiris et al., 2008;Brozzoli et al., 2012;Limanowski et al., 2014;Limanowski and Blankenburg, 2015). Thus, both illusory ownership and proprioceptive drift may be necessary for activation of specific cortical regions. In contrast, in the present study, the participants were requested to adjust their hands at each side of the mirror such that the reflection of the left hand appeared to be the right hand at the beginning of each run, and the hands were placed at symmetrical positions to the mirror. In such conditions, proprioceptive drift could not be induced. Recent studies have revealed that subjective hand ownership assessed by questionnaires may be dissociable from proprioceptive sensation of the hand position (Holmes et al., 2006;Longo et al., 2008;Rohde et al., 2011;Romano et al., 2015). The exact relationship between illusory hand ownership and the proprioceptive sense of the hand position, and the underlying cortical mechanisms are interesting points for future studies.

CONCLUSION
In the present study, we investigated the cortical regions involved in the subjective perception of hardness achieved by integration of visual and haptic information. Participants were required to touch a polyurethane foam pad of constant physical hardness with their right (hidden) hand while observing the mirror reflection of their left hand touching a pad of different hardness, and to estimate the hardness of the pad perceived by the right hand. The behavioral results showed that when they observed the mirror reflection touching the softer or harder pad, they perceived the pad in their right hand as softer or harder, respectively. An fMRI experiment using multivariate analysis suggested that the visual information of finger movements processed in the bilateral OTC may be integrated with haptic input in the left anterior intraparietal sulcus (aIPS), and that the subjective hardness with visual modulation perceived by the right hand may be processed in the cortical network between the left aIPS and parietal operculum (PO).

ETHICS STATEMENT
This study was approved by the local Ethical Committee for the Faculty of Dentistry of the Tokyo Medical and Dental University (D2016-051), Primate Research Institute of Kyoto University (2016-13), Tamagawa University (N27-42), and corresponded to the Human Subjects Guidelines of the Declaration of Helsinki.

AUTHOR CONTRIBUTIONS
NK, MT, and KN planned the research project and acquired funding. YK, NU, MT, and NK designed the study. YK, AM, TH, KM, and NK collected the experimental data. YK and NK analyzed the experimental data and prepared the manuscript with supports from NU and KN.

FUNDING
This study was supported by a Grant-in-Aid for Scientific Research on Innovative Areas (No. 25135711) from the Ministry of Education, Science, Culture, Sports, and Science, Japan (NK), and a Primate Research Institute Project Grant to Kyoto University.

ACKNOWLEDGMENTS
MT (deceased as of July 8, 2017) read and approved the early version of the manuscript in July 2017. We sincerely wish to dedicate this work to the memory of MT. We would also like to thank Editage (www.editage.jp) for English language editing.