- Research article
Spatiotemporal integration of tactile information in human somatosensory cortex
BMC Neurosciencevolume 8, Article number: 21 (2007)
Our goal was to examine the spatiotemporal integration of tactile information in the hand representation of human primary somatosensory cortex (anterior parietal somatosensory areas 3b and 1), secondary somatosensory cortex (S2), and the parietal ventral area (PV), using high-resolution whole-head magnetoencephalography (MEG). To examine representational overlap and adaptation in bilateral somatosensory cortices, we used an oddball paradigm to characterize the representation of the index finger (D2; deviant stimulus) as a function of the location of the standard stimulus in both right- and left-handed subjects.
We found that responses to deviant stimuli presented in the context of standard stimuli with an interstimulus interval (ISI) of 0.33s were significantly and bilaterally attenuated compared to deviant stimulation alone in S2/PV, but not in anterior parietal cortex. This attenuation was dependent upon the distance between the deviant and standard stimuli: greater attenuation was found when the standard was immediately adjacent to the deviant (D3 and D2 respectively), with attenuation decreasing for non-adjacent fingers (D4 and opposite D2). We also found that cutaneous mechanical stimulation consistently elicited not only a strong early contralateral cortical response but also a weak ipsilateral response in anterior parietal cortex. This ipsilateral response appeared an average of 10.7 ± 6.1 ms later than the early contralateral response. In addition, no hemispheric differences either in response amplitude, response latencies or oddball responses were found, independent of handedness.
Our findings are consistent with the large receptive fields and long neuronal recovery cycles that have been described in S2/PV, and suggest that this expression of spatiotemporal integration underlies the complex functions associated with this region. The early ipsilateral response suggests that anterior parietal fields also receive tactile input from the ipsilateral hand. The lack of a hemispheric difference in responses to digit stimulation supports a lack of any functional asymmetry in human somatosensory cortex.
The spatiotemporal integration of tactile inputs from different skin regions and across body parts is an important function of human somatosensory cortex. For example, the integration of inputs across the digits is vital for the successful manual manipulation and identification of objects. However, the mechanism of this integration is not well understood. Good candidates for the performance of this function are the second somatosensory area, S2, and the parietal ventral area, PV. These two fields are mirror symmetric representations of the body's surface [1–4], joined at the representation of the hand. Because it is difficult to distinguish between the hand representations of S2 and PV using functional imaging techniques we refer to this region as S2/PV.
Electrophysiological recording studies in monkeys indicate that neuronal receptive fields in S2 [2, 5, 6] and PV [2, 4, 7] are large, encompassing multiple digits or even the entire hand. Further, these receptive fields are often bilateral, including, for example, both the contra- and ipsilateral hand. Studies of neuroanatomical connections also show that, beside the local homotopic connections in the ipsilateral hemisphere, both S2 [1, 7] and PV [1, 4, 7] have dense bilateral connections.
There is also evidence from human imaging studies that inputs from different skin regions interact in S2/PV. For example, bilateral integration of inputs from the hands takes place in human S2/PV [8–12]. In addition, the activation in human S2 evoked by stimulation of one finger can be modulated by simultaneous  or non-simultaneous [8, 9, 12] stimulation of other fingers of the same hand.
Another possibility is that integration of inputs from different skin regions takes place, at least to some extent, in anterior parietal areas 3a, 3b, 1 and 2 [9, 11–15]. Previous MEG studies have confirmed that tactile stimulation to the human finger evokes responses in both contralateral anterior parietal fields as well as bilateral S2/PV [16, 17]. Each finger has a distinguishable somatotopic representation in contralateral anterior parietal cortex [18, 19]. However, there is some evidence that interaction among digit representations may also occur in anterior parietal cortex. For example, the strength of the early response to stimulation of a single finger was attenuated by simultaneous [11, 14, 15] or non-simultaneous stimulation of another finger at short (<100 ms) ISIs . Further, the magnitude of spatial integration decreased with the distance of separation of the digit representations in anterior parietal fields [11, 14]. There is also evidence from electrophysiological recording  and optical imaging  studies in monkeys that simultaneous stimuli from different skin regions could be merged together into a single activation zone in anterior parietal cortex.
However, in contrast to S2/PV, human MEG studies suggest that responses in anterior parietal cortex to stimulation of one hand are not affected by stimulation of the opposite hand [10–12]. This finding is supported by neuroanatomical results in monkeys indicating that the hand representation in 3b is largely acallosal [1, 22–25]. Previous MEG work also suggests that the spatial integration of inputs across the digits is much stronger in S2/PV than in anterior parietal fields [9, 11, 12]. However, the extent of spatial integration in S2/PV, and the difference in extent and timing of integration between anterior parietal fields and S2/PV has not been quantified. These differences may represent different steps in the complex process of spatiotemporal integration.
In this study, we measured somatosensory evoked fields (SEFs) during non-simultaneous tactile stimulation of digits of both hands with an oddball paradigm. We compared the extent of spatial integration in the hand representations in human S2/PV and anterior parietal fields, and discuss the difference in the timing of integration in these two areas.
Contra- and ipsilateral responses to index finger stimulation
1). Dipole orientation and localization
In all subjects, dipoles were fit for both the early (30–70 ms) and late (70–130 ms) response windows. Figure 1 contains the actual contour plots of MEG sensor data and the plots of the averaged magnetic fields measured using 275 axial gradiometers from a single right-handed subject showing responses to index finger stimulation (average of 'deviants' alone at low rate, ISI = 2s). The left and right columns show responses to right index finger (RD2) and left index finger (LD2) stimulation, respectively. Identified dipole sources of responses in the right hemisphere were superimposed on this subject's MRI and are shown in Figure 2.
Interestingly, bilateral activation appeared during both early and late time periods in this subject for both LD2 and RD2 stimulation, though the early responses were of larger amplitude and shorter latency in the contralateral (e.g. the first peak at 45 ms in Figure 1e) versus the ipsilateral (e.g. the first peak at 51 ms in Figure 1f) hemisphere. The dipole direction and location of the ipsilateral (left hemisphere) response to LD2 stimulation were similar to those of the contralateral (left hemisphere) response to RD2 stimulation and likewise for the right hemisphere responses. The dipole directions tended to be mirror images in the two hemispheres. Dipole analysis and co-registration of the MEG localization information with the MRI scans verified that the anterior parietal area (the posterior wall of the central sulcus) was the primary contributor to both contralateral and ipsilateral early (30–70 ms) responses and that the S2 region (the upper bank of the Sylvian fissure) was the primary contributor to the late (70–130 ms) responses in this subject (Figure 2).
The average current dipole positions were specified in Talairach coordinates (x, y, z in mm): the early contralateral response dipoles in left and right hemispheres were (-49.2 ± 1.6, -19.2 ± 2.5, 45.8 ± 2.4) and (48.0 ± 1.4, -14.5 ± 1.4, 51.3 ± 3.7); the late contralateral response dipoles in left and right hemispheres were (-61.9 ± 3.7, -16.6 ± 3.4, 6.9 ± 3.8) and (56.2 ± 2.4, -5.4 ± 4.7, 13.9 ± 4.8), respectively. These localizations in Talairach space were consistent with the locations of anterior parietal fields and S2/PV in both hemispheres.
Early ipsilateral responses for LD2 stimulation were found in 12/15 right- and 2/6 left-handed subjects, while 7/15 right- and 2/6 left-handed subjects showed early ipsilateral responses to stimulation of RD2. Figure 3 shows two more examples of early ipsilateral responses (one right-handed and one left-handed). The source of early ipsilateral responses was identified in anterior parietal fields. Late ipsilateral responses were recorded from all 21 subjects. Like the late contralateral responses, the late ipsilateral responses also originated from the upper banks of the Sylvian fissure, in the S2 region.
The peak latencies of responses from a single subject (contralateral early and late responses, ipsilateral early and late responses) are shown in Figure 1. The early contralateral responses from both hemispheres of the same subject peaked at 45 ms, while late contralateral responses peaked at about 100 ms (Figure 1a, b and 1d, e). Early and late ipsilateral responses in the right hemisphere peaked at 60 and 120 ms (Figure 1a and 1c). The early ipsilateral responses in the left hemisphere peaked at 51 ms. Later responses, which looked more complex with multiple peaks, peaked around 110 ms (Figure 1d and 1f). Thus, both early and late ipsilateral response peaks appeared later than contralateral peaks in both hemispheres in this subject.
The average early and late contralateral response latencies for all subjects (right- and left-handed) were not significantly different across hemispheres (p > 0.05). The ipsilateral early (56.2 ± 6.6 ms) and late (116.6 ± 13.1 ms) responses appeared significantly later than contralateral early (46.8 ± 6.9 ms) and late (108.8 ± 12.5 ms) responses in all subjects. The mean delay of the early ipsilateral responses relative to the contralateral responses was 10.7 ± 6.1 ms (range: 3.5–24 ms), while the delay of the late ipsilateral responses relative to the contralateral responses was 12.7 ± 13.1 ms (range 3.0–30.8 ms). The mean difference between delays for the early versus late responses was not significant (p > 0.05).
The subject's data from Figure 1 clearly show the early ipsilateral responses above the noise level (Figure 1c and 1f). The amplitudes of this subject's early ipsilateral responses were 47.3 and 35.8 fT (RMS), while the noise levels were 13.6 and 9.8 fT (RMS), in the left and right hemispheres, respectively. The average amplitude of the early ipsilateral response (37.9 ± 2.0 fT) across all subjects was significantly higher (p < 0.01) than the average noise level (14.5 ± 0.9 fT) based on the pre-stimulation period of 50 ms.
The average amplitude of early and late responses for all subjects did not show a significant difference between hemispheres (p > 0.05), nor was there a significant difference in amplitude of response for right- versus left-handers (p > 0.05). As the strength of early and late responses did not show a significant difference between hemispheres, the responses to left and right hand stimulation were combined for further analysis.
The mean amplitude of the early response was significantly smaller (p < 0.01) for the early ipsilateral response (37.9 ± 2.0 fT) than for the early contralateral response (66.9 ± 3.3 fT). Similarly, the late ipsilateral response (50.4 ± 4.8 fT) was significantly weaker (p < 0.01) than the late contralateral response (69.3 ± 5.6 fT). Further, the mean amplitude of the early contralateral response (66.9 ± 3.3 fT) was not significantly different from the late contralateral response (69.3 ± 5.6 fT; p > 0.05), while the late ipsilateral response (50.4 ± 4.8 fT) was significantly stronger than the early ipsilateral response (37.9 ± 2.0 fT; p < 0.01; see Figure 4).
Spatial integration between digit representations
Figure 5 is an example of the spatial interaction from a typical right-handed subject. This figure shows the averaged magnetic field responses recorded from all 275 sensors under the different stimulus conditions. Each line is the average over all trials for a given sensor. The left column in Figure 5 shows contralateral (right hemisphere) sensor waveforms for LD2 stimulation, and the right column shows the contralateral (left hemisphere) sensor waveforms for RD2 stimulation. The latencies of both early and late responses did not change under different stimulation conditions.
In an example subject, low rate index finger stimulation (mean ISI: 2s) clearly elicited both early and later responses (Figure 5a and 5f). The response amplitude was markedly decreased for high rate stimulation (ISI: 0.33 s), and the late response (S2/PV) almost disappeared in both hemispheres (Figure 5b and 5g). The intervened standard stimuli (D3, D4 and opposite D2) suppressed the late response for D2 (deviant) stimulation in both hemispheres. There appeared to be the greatest attenuation for the same (D2) digit high rate stimulation compared to alternate standard digit (D3, D4 and opposite D2) stimulation, and greater attenuation for adjacent (D3) compared to non-adjacent (D4 and opposite D2) digit stimulation (Figure 5b–e and 5g–j). The standard stimuli reduced the amplitude of the early response for D2 (deviant) stimulation in some conditions for this subject (Figure 5c and 5j).
Grand averaged sensor amplitude (RMS; contralateral hemisphere only) and dipole moment (nAm) response values to deviant stimuli are displayed in Figure 4. Both sensor and moment values of early responses were statistically significantly smaller for high rate index finger stimulation (ISI: 0.33s) compared to low rate index finger stimulation (mean ISI: 2s). Both contralateral and ipsilateral early amplitude of responses to single finger stimulation (D2) were not significantly affected by stimulation of the other fingers (D3, D4 and opposite D2).
As with the early responses, sensor and moment values of the late responses were statistically significantly smaller for high rate index finger stimulation (ISI: 0.33s) compared to low rate index finger stimulation (mean ISI: 2s). The late response decreased more than the early response regardless of handedness (p < 0.01; also see Figure 5).
Further, the intervened standard stimuli (D3, D4, and opposite D2) significantly reduced the late contralateral and ipsilateral responses to deviant (D2) stimuli. In the contralateral hemisphere this effect was greater for the adjacent finger (D3) compared to the non-adjacent fingers (D4 or opposite D2). The attenuation of the late contralateral response was significantly greater for the same finger high rate stimulation (D2 alone at high rate) compared to non-adjacent finger (D4 or opposite D2) stimulation. For the ipsilateral response, there were no differences in amplitude reduction related to the location of intervened standard stimulation. Only opposite D2 stimulation differed from high rate D2 stimulation in the late ipsilateral response. The suppressive effect of opposite D2 stimulation was significantly stronger for the contralateral than the ipsilateral late response (Figure 4).
In summary, we examined the spatiotemporal integration across digit representations within anterior parietal somatosensory areas and S2/PV, as well as interhemispheric integration, using non-simultaneous tactile stimulation. We have shown 1) spatial integration of responses in S2/PV, with a decrease in interaction with the increase in separation of digit representations; 2) greater temporal integration of inputs within one digit representation in S2/PV versus anterior parietal fields; 3) both contra- and ipsilateral early responses in anterior parietal fields to cutaneous mechanical stimulation. The early ipsilateral response has a longer latency and lower amplitude than the early contralateral response; and 4) processing of cutaneous inputs is symmetrical across hemispheres. What follows is a discussion of these findings in light of previous work from human and non-human primates.
Spatiotemporal integration between digit representations in S2/PV
The distance-dependence of spatial integration between digit representations in the S2/PV area was examined in this study. We found that the amplitude of the contralateral late responses to deviant stimuli (D2) decreased when standard stimuli (D3, D4 and opposite D2) were intervened. Like the distance-dependence of spatial integration in anterior parietal areas [11, 14], the attenuation decreased with the increasing distance of separation between receptive fields, toward greater attenuation for adjacent fingers (D3) compared to non-adjacent fingers (D4 and opposite D2).
In agreement with previous MEG work on digit response interaction in S2/PV [8, 9, 11, 12, 26], our data show that the standard stimuli significantly attenuate responses in bilateral S2/PV to deviant (D2) stimulation. The pattern of spatial integration in S2/PV is consistent with the fact that neurons in S2 and PV have large receptive fields, frequently encompassing multiple digits [2, 4–7, 27], while cells in anterior parietal areas do not [28, 29].
We also found significant temporal integration in the S2/PV region. The amplitude of the late response was smaller for high rate versus low rate index finger stimulation, and the late responses decreased more than the early responses with high rate stimulation, as previous SEP  and SEF [26, 31–34] investigations have shown. This result is consistent with previous studies revealing the shorter recovery cycle or refractory period for neurons in anterior parietal areas than in S2 [35, 36].
Some non-human primate studies have suggested serial information processing from anterior parietal areas to S2/PV [37, 38], with tactile inputs from anterior parietal fields converging onto S2 [27, 39]. However, other studies suggest hierarchical equivalence of anterior parietal areas and S2 for tactile processing in cats and new-world primates [40–42]. To date, the debate about serial versus parallel organization of somatosensory cortex in humans is unresolved. Both early  and late [44, 45] concurrent activities in anterior parietal and lateral sulcal areas have been reported in humans. In contrast, a recent MEG study revealed that the onset latency of response following transcutaneous electrical stimulation was longer in S2/PV than in anterior parietal areas . However, we cannot rule out the contribution of activity in S2/PV to early responses and anterior parietal areas to late responses. It is possible that both serial and parallel processing play a role in tactile perception.
Serial convergence may account for the spatial and temporal specialties in S2/PV. Unlike anterior parietal areas with restricted and defined neural receptive fields, the large, complex receptive fields within S2/PV may play an important role in integrating tactile inputs over space and time. For example, the digit representations in S2 and surrounding fields lack fine somatotopy, and some cells responded better to proprioceptive rather than cutaneous stimuli . In addition, studies have found that 23% of cells in the S2 region were orientation tuned, often across multiple digit pads . These receptive field properties, tuned for relatively long-range spatial and temporal integration, go some way towards explaining S2's purported involvement in functions like tactile object exploration and identification [47–50], bimanual coordination [2, 4, 6, 7], and tactile learning and memory [49, 50]. Consistent with this idea is evidence in humans that damage to parietotemporal cortices (including S2) causes impairment of tactile object recognition in the absence of basic somesthetic dysfunction [51, 52].
Spatiotemporal integration between digit representations in anterior parietal somatosensory areas
We found that the amplitude of responses in anterior parietal areas to deviant alone stimuli applied to a single digit was significantly lower at high rate (the ISI is also 0.33s) versus low rate (ISI: 2s) stimulation, as previous studies [26, 30–34] have shown. This decrease of amplitude resulting from increased rate of stimulation was much smaller in anterior parietal areas than in S2/PV. However, the average amplitude of early responses in anterior parietal areas for deviant stimulation did not show a significant decrease when the standard stimuli were interleaved between deviants (ISI: 0.33s), though the amplitude decrease was seen in some cases. Previous work using non-simultaneous stimuli with long ISIs (1–4 s) also showed similar results [26, 53]. However, an interaction between digit representations has been reported in anterior parietal areas when using simultaneous [11, 14, 15] or non-simultaneous stimuli with short (<100 ms) ISIs . These findings are in line with the short recovery cycle time constant (~110 ms) for neurons in anterior parietal areas, which was estimated by fitting an exponential curve to SEF intensities for responses to electrical stimulation of the median nerve with different ISIs (100–500 ms) . With a long ISI, the intervened standard inputs may arrive out of the recovery period of neurons in anterior parietal areas, with little effect on responses to deviant stimuli.
It is known that there is an inhibitory surround structure of receptive fields in monkey area 3b [54–58]. A three-component model has been proposed in which a lagged inhibitory surround overlaps the excitatory center and its one or more fixed inhibitory flankers. The fixed (as a spatial filter) and lagged (as a temporal filter) surround inhibitory components confer the spatial and temporal selectivity of neurons in anterior parietal areas . The spatiotemporal property of inhibitory surround structure of neural receptive fields may account for the restricted neural receptive field and rapid recovery cycle in anterior parietal areas observed in the present study.
Ipsilateral responses in anterior parietal areas
In the present investigation we clearly recorded early ipsilateral responses with high goodness-of-fit to single finger tactile stimulation (both right- and left-handers). Though the amplitude of the ipsilateral response was very weak in some cases, on average it was 2.5 times greater than the noise level. An ipsilateral response from anterior parietal field has also been observed in a small number of subjects in several previous studies [59–63] using MEG or EEG to measure responses to median nerve stimulations. A growing body of evidence from monkey studies [64, 65] and human fMRI studies [66, 67] also report ipsilateral responses in anterior parietal fields, suggesting that anterior parietal fields receive tactile input from the ipsilateral hand. Ipsilateral responses were relatively robust in our study, likely due to the cutaneous mechanical stimulation used in our experiment. Recent fMRI studies on both monkeys and humans have shown that responses from ipsilateral anterior parietal areas are more robust for mechanical versus median nerve electrical stimulation .
Single cell recordings from monkeys indicate that some neurons representing multiple digits in areas 2 and 5 have bilateral or ipsilateral receptive fields , and a few nociceptive neurons at the junction of areas 3b and 1 have large, bilateral receptive fields . Thus the source of the early ipsilateral response is likely not only area 3b, but may also include input from cells in area 1, and possibly areas 2 and 5 as well. Recent evidence from electrophysiological recording in monkeys suggests that responses from ipsilateral anterior parietal fields are elicited by feedback rather than feedforward afferents , which may be related to the latency delay of the ipsilateral response observed in previous SEP  and the current MEG study. It has been suggested that transcallosal pathways may mediate responses from ipsilateral anterior parietal areas both in humans [59–61, 67] and in monkeys [64, 65]. However, anatomical evidence from macaque monkeys indicates that anterior parietal areas, in particular area 3b, have sparse or absent callosal connections for the hand representation [1, 22, 25]. Callosal connections become denser in the more caudal fields, from areas 3b to 1 and 2 . Thus, transcallosal pathways mediating the early ipsilateral responses would be through these caudal areas, or S2/PV.
In addition, some previous functional imaging studies suggest functional dominance, or an increased source amplitude, in the left versus the right hemisphere in human primary [70–72] and secondary [73–77] somatosensory cortex. An expanded hand representation in left anterior parietal fields has also been proposed [72, 78]. However, morphometric and cytoarchitectonic measurements show no lateralized differences in somatosensory cortices . Psychophysical tests show no differences in spatial acuity between hands [80, 81], although there does appear to be a left-hemisphere advantage for tactile processing of simultaneity judgment . Furthermore, as in our study, previous somatosensory evoked potential studies show no difference in topography or response amplitude between the two hemispheres [83, 84]. Thus, it is necessary to reassess the hemispheric asymmetry of somatosensory cortex.
The differences between studies may be due to the different types of stimuli used. The electrical stimulation used in previous studies showing asymmetry is not receptor specific, and elicits a more complex cortical activation. A recent MEG study showed that movement might be involved in the lateral asymmetry of somatosensory cortex, and the intensity of electrical stimulation often exceeded the motor threshold in previous studies [70–75, 77]. Passive finger movement was found to evoke laterally asymmetrical responses in somatosensory cortex . Thus, the increased response in left somatosensory cortex evoked by electric stimulation may reflect the lateral asymmetry of movement rather than tactile information processing.
This human MEG study revealed that the bilateral spatiotemporal integration in S2/PV takes place over a large cortical area and over a long time period. Further, the strength of integration in this region is distance-dependent. The wide overlap of digit receptive fields in S2/PV might account for complex functions such as manual exploration and bimanual coordination, thought to be crucially dependent on S2/PV. In contrast, the properties of spatiotemporal integration in anterior parietal areas were different from those in S2/PV with significant temporal integration while spatial integration was reduced, which is consistent with the critical role of anterior parietal fields in spatial discrimination. The early ipsilateral response observed in over half of subjects suggests that anterior parietal fields do accept tactile inputs from the ipsilateral side of body. In addition, no tactile response difference between hemispheres indicated functional symmetry in human somatosensory cortex.
SEFs were recorded in fifteen right-handed and six left-handed healthy adult subjects (10 male and 5 female right-handers; 2 male and 4 female left-handers; age range 22–40 years), using an Omega 2000 Whole-Cortex MEG System (CTF Systems Inc. Port Coquitlam, Canada; 275 DC SQUID first-order axial gradiometers). All subjects signed an informed consent as approved by the Committee on Human Research of the University of California, San Francisco. The Edinburgh inventory  was used to determine the direction and degree of handedness for each subject. The subjects were seated comfortably and maintained their head position during MEG testing in a magnetically shielded room. Subjects' hands were placed palm up on opposite armrests. The subjects closed their eyes, wore earplugs and were asked not to pay attention to the tactile stimulation.
Tactile stimuli were pneumatically driven pulses (~140 ms duration) applied to the tips of the digits with balloon diaphragms. A digit oddball paradigm [9, 16, 30] was used to examine the representations of the index fingers (D2, infrequent deviant stimulus) while varying the location of the frequent standard stimuli across adjacent (D3), non-adjacent (D4) and contralateral (D2) digits. In this oddball paradigm, infrequent deviant stimuli interspersed with frequent standard stimuli were presented. The intensity of all stimuli was well above detection threshold at 17 PSI (pounds per square inch). The location of the standard was changed between blocks and the order (including deviant alone block) was randomized across subjects. In each block of 600–900 trials with an ISI of 0.33s and a jitter of 10 ms, standard and deviant stimuli were presented at probabilities of 0.83 and 0.17 respectively. Deviants were always followed by 3–7 standards to allow for an adaptation to the standard. To compare the tactile response to single finger stimulation in different hemispheres, deviants were also presented alone at two rates: one equal to that for deviants (mean ISI: 2 s) and the other equal to that for the compound of deviants and standards (mean ISI: 0.33s) in the mismatch paradigm (Figure 6).
Data recording and analysis
Data were collected at a sample rate of 1200 Hz. The filtering passband for data analysis was 2–40 Hz. About 100 artifact-free trials for deviant stimuli (D2) were averaged in each test block. Head position relative to the MEG sensors was determined before and after each test block by means of three small coils placed at landmark sites (nasion, left and right preauricular points).
MEG data was analyzed using an equivalent current dipole (ECD) embedded in a spherical conducting medium. The locations of index finger representations in somatosensory cortex were determined using a dipole fit with contralateral index finger stimulation alone. SEFs in somatosensory cortex peaking in the time window up to 200 ms following stimulus onset were analyzed. The early (30–70 ms) response was analyzed for activation in anterior parietal fields; and the late response (70–130 ms) was analyzed for activation of S2/PV [11, 26]. Sensors recording from the hemisphere contralateral to the stimulated index finger were chosen to determine the ECD of the most dominant source. The position and orientation of the ECD corresponding to the early response were first found and fixed; then another dipole corresponding to the late response was added with the early one fixed, then was fitted and fixed successively. Only sources with high goodness of fit (> 85%) were accepted. Dipoles matching anterior parietal fields and S2/PV in each hemisphere were identified and fixed based on contralateral index finger stimulation, and they were used for analysis of data from other stimulation conditions. The response latencies, amplitudes (root-mean-square value, RMS) of each hemisphere and dipole moments (Q value) for all of four dipole locations were estimated for the different stimulation conditions based on peaks within the early and late time periods.
Eight of the subjects were also scanned using a 1.5T MRI scanner (GE Medical System, Milwaukee, WI) to acquire a 3D structural image (flip angle = 40°, TR = 27 ms, TE = 6 ms, FOV = 240 × 240 mm, 1.5 mm slice thickness, 256 × 256 × 124 pixels). Three fiducials were placed on the subject at the same three locations as the localizing coils in MEG. This information was used to coregister the MEG data to the MRI image. Each subject's structural MRI was normalized to MNI space using Neurodynamic Utility Toolbox for MEG . Then the MNI coordinates of dipole positions were converted to Talairach coordinates, using a non-linear transform .
Paired Student t-tests were used to assess significant difference between condition pairs. In cases where the number of conditions exceeds two, repeated measures ANOVAs were used for assessing statistical significance between responses across different stimulation conditions. In these cases, Bonferroni post-hoc tests were used to assess significance between specific condition pairs, whereby after correction for multiple comparisons a significance threshold of either p < 0.01 or p < 0.05 was used. Data are presented as mean values ± standard error of the mean throughout.
parietal ventral area
somatosensory evoked fields
equivalent current dipole
root mean square
left index finger
right index finger
Krubitzer LA, Kaas JH: The organization and connections of somatosensory cortex in marmosets. J Neurosci. 1990, 10 (3): 952-74.
Krubitzer L, Clarey J, Tweedale R, Elston G, Calford M: A redefinition of somatosensory areas in the lateral sulcus of macaque monkeys. J Neurosci. 1995, 15 (5 Pt 2): 3821-39.
Disbrow E, Roberts T, Krubitzer L: Somatotopic organization of cortical fields in the lateral sulcus of Homo sapiens: evidence for SII and PV. J Comp Neurol. 2000, 418 (1): 1-21. 10.1002/(SICI)1096-9861(20000228)418:1<1::AID-CNE1>3.0.CO;2-P.
Qi HX, Lyon DC, Kaas JH: Cortical and thalamic connections of the parietal ventral somatosensory area in marmoset monkeys (Callithrix jacchus). J Comp Neurol. 2002, 443 (2): 168-82. 10.1002/cne.10113.
Whitsel BL, Petrucelli LM, Werner G: Symmetry and connectivity in the map of the body surface in somatosensory area II of primates. J Neurophysiol. 1969, 32 (2): 170-83.
Burton H, Sinclair R: Somatosensory cortex and tactile perceptions. Touch and Pain. 1996, London: Academic Press
Disbrow E, Litinas E, Recanzone GH, Padberg J, Krubitzer L: Cortical connections of the second somatosensory area and the parietal ventral area in macaque monkeys. J Comp Neurol. 2003, 462 (4): 382-99. 10.1002/cne.10731.
Hari R, Hamalainen H, Hamalainen M, Kekoni J, Sams M, Tiihonen J: Separate finger representations at the human second somatosensory cortex. Neuroscience. 1990, 37 (1): 245-9. 10.1016/0306-4522(90)90210-U.
Forss N, Jousmaki V, Hari R: Interaction between afferent input from fingers in human somatosensory cortex. Brain Res. 1995, 685 (1–2): 68-76. 10.1016/0006-8993(95)00424-O.
Disbrow E, Roberts T, Poeppel D, Krubitzer L: Evidence for interhemispheric processing of inputs from the hands in human S2 and PV. J Neurophysiol. 2001, 85 (5): 2236-44.
Hoechstetter K, Rupp A, Stanèák A, Meinck HM, Stippich C, Berg P, Scherg M: Interaction of tactile input in the human primary and secondary somatosensory cortex – a magnetoencephalographic study. Neuroimage. 2001, 14 (3): 759-67. 10.1006/nimg.2001.0855.
Simões C, Mertens M, Forss N, Jousmäki V, Lütkenhöner B, Hari R: Functional overlap of finger representations in human SI and SII cortices. J Neurophysiol. 2001, 86 (4): 1661-5.
Huttunen J, Ahlfors S, Hari R: Interaction of afferent impulses in the human primary sensorimotor cortex. Electroencephalogr Clin Neurophysiol. 1992, 82 (3): 176-81. 10.1016/0013-4694(92)90165-E.
Biermann K, Schmitz F, Witte OW, Konczak J, Freund HJ, Schnitzler A: Interaction of finger representation in the human first somatosensory cortex: a neuromagnetic study. Neurosci Lett. 1998, 251 (1): 13-6. 10.1016/S0304-3940(98)00480-7.
Tanosaki M, Suzuki A, Takino R, Kimura T, Iguchi Y, Kurobo Y, Haruta Y, Hoshi Y, Hashimoto I: Neural mechanisms for generation of tactile interference effects on somatosensory evoked magnetic fields in humans. Clin Neurophysiol. 2002, 113 (5): 672-80. 10.1016/S1388-2457(02)00052-4.
Hari R, Forss N: Magnetoencephalography in the study of human somatosensorycortical processing. Philos Trans R Soc Lond B Biol Sci. 1999, 354 (1387): 1145-54. 10.1098/rstb.1999.0470.
Kakigi R, Hoshiyama M, Shimojo M, Naka D, Yamasaki H, Watanabe S, Xiang J, Maeda K, Lam K, Itomi K, Nakamura A: The somatosensory evoked magnetic fields. Prog Neurobiol. 2000, 61 (5): 495-523. 10.1016/S0301-0082(99)00063-5.
Nakamura A, Yamada T, Goto A, Kato T, Ito K, Abe Y, Kachi T, Kakigi R: Somatosensory homunculus as drawn by MEG. Neuroimage. 1998, 7 (4 Pt 1): 377-86. 10.1006/nimg.1998.0332.
Maeda K, Kakigi R, Hoshiyama M, Koyama S: Topography of the secondary somatosensory cortex in humans: a magnetoencephalo-graphic study. Neuroreport. 1999, 10 (2): 301-6. 10.1097/00001756-199902050-00018.
Gardner EP, Costanzo RM: Spatial integration of multiple-point stimuli in primary somatosensory cortical receptive fields of alert monkeys. J Neurophysiol. 1980, 43 (2): 420-43.
Chen LM, Friedman RM, Roe AW: Optical imaging of a tactile illusion in area 3b of the primary somatosensory cortex. Science. 2003, 302 (5646): 881-5. 10.1126/science.1087846.
Jones EG, Powell TP: Connexions of the somatic sensory cortex of the rhesus monkey. II. Contralateral cortical connexions. Brain. 1969, 92 (4): 717-30. 10.1093/brain/92.4.717.
Karol EA, Pandya DN: The distribution of the corpus callosum in the rhesus monkey. Brain. 1971, 94 (3): 471-786. 10.1093/brain/94.3.471.
Killackey HP, Gould HJ, Cusick CG, Pons TP, Kaas JH: The relation of corpus callosum connections to architectonic fields and body surface maps in sensorimotor cortex of new and old world monkeys. J Comp Neurol. 1983, 219 (4): 384-419. 10.1002/cne.902190403.
Padberg J, Disbrow E, Krubitzer L: The organization and connections of anterior and posterior parietal cortex in titi monkeys: do New World monkeys have an area 2?. Cereb Cortex. 2005, 15 (12): 1938-63. 10.1093/cercor/bhi071.
Mertens M, Lütkenhöner B: Efficient neuromagnetic determination of landmarks in the somatosensory cortex. Clin Neurophysiol. 2000, 111 (8): 1478-87. 10.1016/S1388-2457(00)00349-7.
Fitzgerald PJ, Lane JW, Thakur PH, Hsiao SS: Receptive field properties of the macaque second somatosensory cortex: evidence for multiple functional representations. J Neurosci. 2004, 24 (49): 11193-204. 10.1523/JNEUROSCI.3481-04.2004.
Nelson RJ, Sur M, Felleman DJ, Kaas JH: Representations of the body surface in postcentral parietal cortex of Macaca fascicularis. J Comp Neurol. 1980, 192 (4): 611-43. 10.1002/cne.901920402.
Sripati AP, Yoshioka T, Denchev P, Hsiao SS, Johnson KO: Spatiotemporal receptive fields of peripheral afferents and cortical area 3b and 1 neurons in the primate somatosensory system. J Neurosci. 2006, 26 (7): 2101-14. 10.1523/JNEUROSCI.3720-05.2006.
Kekoni J, Hämäläinen H, Saarinen M, Gröhn J, Reinikainen K, Lehtokoski A, Näätänen R: Rate effect and mismatch responses in the somatosensory system: ERP-recordings in humans. Biol Psychol. 1997, 46 (2): 125-42. 10.1016/S0301-0511(97)05249-6.
Tiihonen J, Hari R, Hamalainen M: Early deflections of cerebral magnetic responses to median nerve stimulation. Electroencephalogr Clin Neurophysiol. 1989, 74 (4): 290-6. 10.1016/0168-5597(89)90059-2.
Wikstrom H, Huttunen J, Korvenoja A, Virtanen J, Salonen O, Aronen H, Ilmoniemi RJ: Effects of interstimulus interval on somatosensory evoked magnetic fields (SEFs): a hypothesis concerning SEF generation at the primary sensorimotor cortex. Electroencephalogr Clin Neurophysiol. 1996, 100 (6): 479-87.
Mauguiere F, Merlet I, Forss N, Vanni S, Jousmaki V, Adeleine P, Hari R: Activation of a distributed somatosensory cortical network in the human brain: a dipole modelling study of magnetic fields evoked by median nerve stimulation. Part II: Effects of stimulus rate, attention and stimulus detection. Electroencephalogr Clin Neurophysiol. 1997, 104 (4): 290-5. 10.1016/S0013-4694(97)00018-7.
Nagamine T, Makela J, Mima T, Mikuni N, Nishitani N, Satoh T, Ikeda A, Shibasaki H: Serial processing of the somesthetic information revealed by different effects of stimulus rate on the somatosensory-evoked potentials and magnetic fields. Brain Res. 1998, 791 (1–2): 200-8. 10.1016/S0006-8993(98)00095-X.
Schnitzler A, Volkmann J, Enck P, Frieling T, Witte OW, Freund HJ: Different cortical organization of visceral and somatic sensation in humans. Eur J Neurosci. 1999, 11 (1): 305-15. 10.1046/j.1460-9568.1999.00429.x.
Hamada Y, Otsuka S, Okamoto T, Suzuki R: The profile of the recovery cycle in human primary and secondary somatosensory cortex: a magnetoencephalography study. Clin Neurophysiol. 2002, 113 (11): 1787-93. 10.1016/S1388-2457(02)00258-4.
Pons TP, Garraghty PE, Friedman DP, Mishkin M: Physiological evidence for serial processing in somatosensory cortex. Science. 1987, 237 (4813): 417-20. 10.1126/science.3603028.
Pons TP, Garraghty PE, Mishkin M: Serial and parallel processing of tactual information in somatosensory cortex of rhesus monkeys. J Neurophysiol. 1992, 68 (2): 518-27.
Friedman DP, Murray EA, O'Neill JB, Mishkin M: Cortical connections of the somatosensory fields of the lateral sulcus of macaques: evidence for a corticolimbic pathway for touch. J Comp Neurol. 1986, 252 (3): 323-47. 10.1002/cne.902520304.
Zhang HQ, Murray GM, Turman AB, Mackie PD, Coleman GT, Rowe MJ: Parallel processing in cerebral cortex of the marmoset monkey: effect of reversible SI inactivation on tactile responses in SII. J Neurophysiol. 1996, 76 (6): 3633-55.
Zhang HQ, Murray GM, Coleman GT, Turman AB, Zhang SP, Rowe MJ: Functional characteristics of the parallel SI- and SII-projecting neurons of the thalamic ventral posterior nucleus in the marmoset. J Neurophysiol. 2001, 85 (5): 1805-22.
Zhang HQ, Zachariah MK, Coleman GT, Rowe MJ: Hierarchical equivalence of somatosensory areas I and II for tactile processing in the cerebral cortex of the marmoset monkey. J Neurophysiol. 2001, 85 (5): 1823-35.
Karhu J, Tesche CD: Simultaneous early processing of sensory input in human primary (SI) and secondary (SII) somatosensory cortices. J Neurophysiol. 1999, 81 (5): 2017-25.
Allison T, McCarthy G, Wood CC, Williamson PD, Spencer DD: Human cortical potentials evoked by stimulation of the median nerve. II. Cytoarchitectonic areas generating long-latency activity. J Neurophysiol. 1989, 62 (3): 711-22.
Stancak A, Hoechstetter K, Tintera J, Vrana J, Rachmanova R, Kralik J, Scherg M: Source activity in the human secondary somatosensory cortex depends on the size of corpus callosum. Brain Res. 2002, 936 (1–2): 47-57. 10.1016/S0006-8993(02)02502-7.
Inui K, Wang X, Tamura Y, Kaneoke Y, Kakigi R: Serial processing in the human somatosensory system. Cereb Cortex. 2004, 14 (8): 851-7. 10.1093/cercor/bhh043.
Garcha HS, Ettlinger G: The effects of unilateral or bilateral removals of the second somatosensory cortex (area SII): a profound tactile disorder in monkeys. Cortex. 1978, 14 (3): 319-26.
Garcha HS, Ettlinger G: Tactile discrimination learning in the monkey: the effects of unilateral or bilateral removals of the second somatosensory cortex (area SII). Cortex. 1980, 16 (3): 397-412.
Mishkin M: Analogous neural models for tactual and visual learning. Neuropsychologia. 1979, 17 (2): 139-51. 10.1016/0028-3932(79)90005-8.
Murray EA, Mishkin M: Relative contributions of SII and area 5 to tactile discrimination in monkeys. Behav Brain Res. 1984, 11 (1): 67-83. 10.1016/0166-4328(84)90009-3.
Caselli RJ: Rediscovering tactile agnosia. Mayo Clin Proc. 1991, 66 (2): 129-42.
Caselli RJ: Ventrolateral and dorsomedial somatosensory association cortex damage produces distinct somesthetic syndromes in humans. Neurology. 1993, 43 (4): 762-71.
Huttunen J, Hari R, Leinonen L: Cerebral magnetic responses to stimulation of ulnar and median nerves. Electroencephalogr Clin Neurophysiol. 1987, 66 (4): 391-400. 10.1016/0013-4694(87)90208-2.
Mountcastle VB, Powell TP: Neural mechanisms subserving cutaneous sensibility, with special reference to the role of afferent inhibition in sensory perception and discrimination. Bull Johns Hopkins Hosp. 1959, 105: 201-32.
Iwamura Y, Tanaka M, Sakamoto M, Hikosaka O: Functional subdivisions representing different finger regions in area 3 of the first somatosensory cortex of the conscious monkey. Exp Brain Res. 1983, 51: 315-326.
DiCarlo JJ, Johnson KO, Hsiao SS: Structure of receptive fields in area 3b of primary somatosensory cortex in the alert monkey. J Neurosci. 1998, 18 (7): 2626-45.
DiCarlo JJ, Johnson KO: Velocity invariance of receptive field structure in somatosensory cortical area 3b of the alert monkey. J Neurosci. 1999, 19 (1): 401-19.
DiCarlo JJ, Johnson KO: Spatial and temporal structure of receptive fields in primate somatosensory area 3b: effects of stimulus scanning direction and orientation. J Neurosci. 2000, 20 (1): 495-510.
Schnitzler A, Salmelin R, Salenius S, Jousmaki V, Hari R: Tactile information from the human hand reaches the ipsilateral primary somatosensory cortex. Neurosci Lett. 1995, 200 (1): 25-8. 10.1016/0304-3940(95)12065-C.
Korvenoja A, Wikstrom H, Huttunen J, Virtanan J, Laine P, Aronen HJ, Seppalainen AM, Ilmoniemi RJ: Activation of ipsilateral primary sensorimotor cortex by median nerve stimulation. Neuroreport. 1995, 6 (18): 2589-93. 10.1097/00001756-199512150-00033.
Korvenoja A, Huttunen J, Salli E, Pohjonen H, Martinkauppi S, Palva JM, Lauronen L, Virtanen J, Ilmoniemi RJ, Aronen HJ: Activation of multiple cortical areas in response to somatosensory stimulation: combined magnetoencephalographic and functional magnetic resonance imaging. Hum Brain Mapp. 1999, 8 (1): 13-27. 10.1002/(SICI)1097-0193(1999)8:1<13::AID-HBM2>3.0.CO;2-B.
Noachtar S, Lüders HO, Dinner DS, Klem G: Ipsilateral median somatosensory evoked potentials recorded from human somatosensory cortex. Electroencephalogr Clin Neurophysiol. 1997, 104 (3): 189-98. 10.1016/S0168-5597(97)00013-0.
Kanno A, Nakasato N, Hatanaka K, Yoshimoto T: Ipsilateral area 3b responses to median nerve somatosensory stimulation. Neuroimage. 2003, 18 (1): 169-77. 10.1006/nimg.2002.1283.
Lipton ML, Fu K, Branch CA, Schroeder CE: Ipsilateral hand input to area 3b revealed by converging hemodynamic and electrophysiological analyses in macaque monkeys. J Neurosci. 2006, 26 (1): 180-185. 10.1523/JNEUROSCI.1073-05.2006.
Tommerdahl M, Simons SB, Chiu JS, Favorov O, Whitsel BL: Ipsilateral input modifies the primary somatosensory cortex response to contralateral skin flutter. J Neurosci. 2006, 26 (22): 5970-7. 10.1523/JNEUROSCI.5270-05.2006.
Nihashi T, Naganawa S, Sato C, Kawai H, Nakamura T, Fukatsu H, Ishigaki T, Aoki I: Contralateral and ipsilateral responses in primary somatosensory cortex following electrical median nerve stimulation – an fMRI study. Clin Neurophysiol. 2005, 116 (4): 842-8. 10.1016/j.clinph.2004.10.011.
Hlushchuk Y, Hari R: Transient suppression of ipsilateral primary somatosensory cortex during tactile finger stimulation. J Neurosci. 2006, 26 (21): 5819-24. 10.1523/JNEUROSCI.5536-05.2006.
Iwamura Y, Iriki A, Tanaka M: Bilateral hand representation in the postcentral somatosensory cortex. Nature. 1994, 369 (6481): 554-6. 10.1038/369554a0.
Kenshalo DR, Isensee O: Responses of primate SI cortical neurons to noxious stimuli. J Neurophysiol. 1983, 50 (6): 1479-96.
Rossini PM, Narici L, Martino G, Pasquarelli A, Peresson M, Pizzella V, Tecchio F, Romani GL: Analysis of interhemispheric asymmetries of somatosensory evoked magnetic fields to right and left median nerve stimulation. Electroencephalogr Clin Neurophysiol. 1994, 91 (6): 476-82. 10.1016/0013-4694(94)90168-6.
Buchner H, Ludwig I, Waberski T, Wilmes K, Ferbert A: Hemispheric asymmetries of early cortical somatosensory evoked potentials revealed by topographic analysis. Electromyogr Clin Neurophysiol. 1995, 35 (4): 207-15.
Jung P, Baumgärtner U, Bauermann T, Magerl W, Gawehn J, Stoeter P, Treede RD: Asymmetry in the human primary somatosensory cortex and handedness. Neuroimage. 2003, 19 (3): 913-23. 10.1016/S1053-8119(03)00164-2.
Forss N, Hari R, Salmelin R, Ahonen A, Hämäläinen M, Kajola M, Knuutila J, Simola J: Activation of the human posterior parietal cortex by median nerve stimulation. Exp Brain Res. 1994, 99 (2): 309-15. 10.1007/BF00239597.
Kany C, Treede RD: Median and tibial nerve somatosensory evoked potentials: middle-latency components from the vicinity of the secondary somatosensory cortex in humans. Electroencephalogr Clin Neurophysiol. 1997, 104 (5): 402-10. 10.1016/S0168-5597(97)00045-2.
Simões C, Alary F, Forss N, Hari R: Left-hemisphere-dominant SII activation after bilateral median nerve stimulation. Neuroimage. 2002, 15 (3): 686-90. 10.1006/nimg.2001.1007.
Alary F, Simões C, Jousmäki V, Forss N, Hari R: Cortical activation associated with passive movements of the human index finger: an MEG study. Neuroimage. 2002, 15 (3): 691-6. 10.1006/nimg.2001.1010.
Mäkelä JP, Illman M, Jousmäki V, Numminen J, Lehecka M, Salenius S, Forss N, Hari R: Dorsal penile nerve stimulation elicits left-hemisphere dominant activation in the second somatosensory cortex. Hum Brain Mapp. 2003, 18 (2): 90-9. 10.1002/hbm.10078.
Sörös P, Knecht S, Imai T, Gürtler S, Lütkenhöner B, Ringelstein EB, Henningsen H: Cortical asymmetries of the human somatosensory hand representation in right- and left-handers. Neurosci Lett. 1999, 271 (2): 89-92. 10.1016/S0304-3940(99)00528-5.
White LE, Andrews TJ, Hulette C, Richards A, Groelle M, Paydarfar J, Purves D: Structure of the human sensorimotor system. II: Lateral symmetry. Cereb Cortex. 1997, 7 (1): 31-47. 10.1093/cercor/7.1.31.
Sathian K, Zangaladze A: Tactile spatial acuity at the human fingertip and lip: bilateral symmetry and interdigit variability. Neurology. 1996, 46 (5): 1464-6.
Vega-Bermudez F, Johnson KO: Differences in spatial acuity between digits. Neurology. 2001, 56 (10): 1389-91.
Nicholls ME, Lindell AK: A left hemisphere, but not right hemispace, advantage for tactual simultaneity judgments. Percept Psychophys. 2000, 62 (4): 717-25.
Kakigi R, Shibasaki H: Effects of age, gender, and stimulus side on scalp topography of somatosensory evoked potentials following median nerve stimulatiom. J Clin Neurophysiol. 1991, 8 (3): 320-330.
Kakigi R, Shibasaki H: Effects of age, gender, and stimulus side on scalp topography of somatosensory evoked potentials following posterior tibial nerve stimulatiom. J Clin Neurophysiol. 1992, 9 (3): 431-440.
Oldfield RC: The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia. 1971, 9 (1): 97-113. 10.1016/0028-3932(71)90067-4.
Dalal SS, Zumer JM, Agrawal V, Hild KE, Sekihara K, Nagarajan SS: NUTMEG: a neuromagnetic source reconstruction toolbox. Neurol Clin Neurophysiol. 2004, 2004: 52.
Brett M: The MNI brain and the Talairach atlas. [http://imaging.mrc-cbu.cam.ac.uk/imaging/MniTalairach]
This work was supported by NIH grant NS044590 to EAD and from the Dana Foundation to SSN. We thank Anne Findlay and Susanne Honma for their excellent technical help with data collection and analysis. We also thank Sarang S. Dalal for help with data analysis using the Neurodynamic Utility Toolbox for MEG.
ZZ, EAD, DJM and SSN designed the study. ZZ and EAD drafted the manuscript. JMZ, DJM and SSN helped to draft the manuscript. ZZ and JMZ recorded and analyzed MEG and MRI data. All authors read and approved the final manuscript.