Activation of Bilateral Secondary Somatosensory Cortex With Right Hand Touch Stimulation: A Meta-Analysis of Functional Neuroimaging Studies

Background: Brain regions involved in processing somatosensory information have been well documented through lesion, post-mortem, animal, and more recently, structural and functional neuroimaging studies. Functional neuroimaging studies characterize brain activation related to somatosensory processing; yet a meta-analysis synthesis of these findings is currently lacking and in-depth knowledge of the regions involved in somatosensory-related tasks may also be confounded by motor influences. Objectives: Our Activation Likelihood Estimate (ALE) meta-analysis sought to quantify brain regions that are involved in the tactile processing of the right (RH) and left hands (LH) separately, with the exclusion of motor related activity. Methods: The majority of studies (n = 41) measured activation associated with RH tactile stimulation. RH activation studies were grouped into those which conducted whole-brain analyses (n = 29) and those which examined specific regions of interest (ROI; n = 12). Few studies examined LH activation, though all were whole-brain studies (N = 7). Results: Meta-analysis of brain activation associated with RH tactile stimulation (whole-brain studies) revealed large clusters of activation in the left primary somatosensory cortex (S1) and bilaterally in the secondary somatosensory cortex (S2; including parietal operculum) and supramarginal gyrus (SMG), as well as the left anterior cingulate. Comparison between findings from RH whole-brain and ROI studies revealed activation as expected, but restricted primarily to S1 and S2 regions. Further, preliminary analyses of LH stimulation studies only, revealed two small clusters within the right S1 and S2 regions, likely limited due to the small number of studies. Contrast analyses revealed the one area of overlap for RH and LH, was right secondary somatosensory region. Conclusions: Findings from the whole-brain meta-analysis of right hand tactile stimulation emphasize the importance of taking into consideration bilateral activation, particularly in secondary somatosensory cortex. Further, the right parietal operculum/S2 region was commonly activated for right and left hand tactile stimulation, suggesting a lateralized pattern of somatosensory activation in right secondary somatosensory region. Implications for further research and for possible differences in right and left hemispheric stroke lesions are discussed.


INTRODUCTION
Somatosensory function is crucial for daily life, guiding our interactions with the world around us through the detection, discrimination and recognition of body sensations (1). Somatosensation is important not only for perception, but also for goal-directed action (2,3). For example, somatosensation contributes to the fundamental pinch grip-lift-and hold task (4) and is important in dexterous movement of the hand (5). Following stroke, reduced functional arm use is contributed to by motor and somatosensory deficits. Somatosensory impairment has a negative impact on grasp and manipulation of objects (6) and is associated with reduced arm use (7). Further, somatosensory brain regions have been implicated in motor recovery (8). It has been suggested that somatosensory processing for the guidance of action can be dissociated from the processing that leads to perception (2). Here we focus on brain regions involved in somatosensation, specifically tactile stimulation of the hand, without motor confounds.
The neuroanatomy of somatosensory processing is well established through a large body of lesion, post-mortem, animal and structural neuroimaging studies (9)(10)(11)(12). Reproducible functional activation in the contralateral primary somatosensory cortex (S1) has been demonstrated in healthy controls when asked to perceive a touch stimulus to their fingertips (13). Technological advances in recent years have even allowed mapping of individual fingers to corresponding areas of S1 (14) and the temporal acuity of anticipation of a tactile stimulus originating in the ipsilateral S1 (15). Different patterns of activation and lateralization emerge when examining somatosensory processing in the secondary somatosensory cortex (S2). Median nerve stimulation has been shown to activate bilateral S2 regions, including parietal operculum, regardless of the hand being stimulated, but only the contralateral S1 (16). This has also been seen in other stimulation studies. Lee et al. (17) recently examined the differential neural activations associated with vibrotactile, pressure and temperature stimulation of right palm, showing common activation in the contralateral S1 and bilateral S2/insula regardless of stimulation type. Bilateral S2 region activation has also been seen with vibrotactile stimulation irrespective of other cognitive demands (18). It has been suggested that serial somatosensory processing occurs from contralateral S1 to contralateral S2 in response to electrical stimulation, but when stimulation becomes more intense or painful there is an increase in hemispheric integration (19). A meta-analysis of studies examining the functional role of S2 in somatosensory processing divided the area into OP1 (parietal operculum 1), OP2, OP3, and OP4 (10). While OP1 is reported to represent the human homolog of macaque area S2 and was generally more responsive to pure somatosensory (tactile) stimuli (10), overall the areas were all implicated in different somatosensory processes (20). A thorough review of the functional role of S2, from the bi-laterality of activation with unilateral stimulation, to the mapping of the hand area spread of OP1-OP4, has been provided by Eickhoff et al. (10).
When examining the literature it becomes clear that the functional activation of somatosensory processing in the brain is still a developing area. There are various stimulation techniques to investigate reflexive neural activity, for example vibrotactile stimulation (18,21) as opposed to MNS median nerve stimulation (16,22,23), that can yield different results. Somatosensory stimuli are applied to various body parts, including the face, upper limb, and lower limb (10, 24), but may not be performed on each hand separately (25,26). Finally, studies have often been confounded by motor contributions to the task, e.g., involving movement intention and/or execution (27)(28)(29).
Our aim was to characterize and synthesize the somatosensory brain activation network during touch sensation, with potential influence of motor contributions eliminated. We employed the ALE meta-analytic technique to provide a statistically-based likelihood estimation of the brain regions that are consistently activated during tactile stimulation of the hands. Firstly, studies were limited to those that involved only tactile stimulation of the right (RH) or left hand (LH) separately in order to allow interpretation of networks that account for hemispheric dominance. Following this, studies which incorporated any motor movements during the stimulation task were excluded, to address confounding motor influence during somatosensory task performance. Lastly, to characterize neural correlates specific to touch sensation, studies involving other somatosensory modalities, such as pain or proprioception, were excluded.

Identification of Studies for Meta-Analysis
The meta-analysis of neuroimaging studies was conducted according to the PRISMA statement and recorded using the suggested checklist (30). A thorough literature search was conducted using Web of Science database (conducted December 12, 2017) and the following search terms: (fMRI OR MRI OR PET OR "functional magnetic resonance imaging" OR "positron emission tomography" OR neuroimaging OR "brain imaging" OR "neural activation") AND (somatosen * OR sens * OR tactile) AND (hand OR "upper limb" OR finger) AND (health * OR control). These papers were then crosschecked with papers identified in the Sleuth functional database (31)(32)(33). The Sleuth database was searched for "somesthesis perception" in the behavioral domain and for "activation only" studies. These were reviewed using the strict inclusion criteria (see Figure 1).

Activation Likelihood Estimation Meta-Analysis
The meta-analysis was performed using Activation Likelihood Estimation (ALE) on the activation voxel coordinates reported by the selected study (34)(35)(36). Analyses were conducted using GingerALE (version 2.3.6) (37) software (downloaded from http://brainmap.org/ale), with coordinates in Montreal Neurological Institute (MNI) space (38,39). Coordinates reported in Talairach space (40) were converted to MNI space using the "icbm2tal transform" (41,42). To minimize withinexperiment and within-group effects we utilized the modified algorithm described in Turkeltaub,Eickhoff (36) and, thus, were able to include multiple contrasts from within the one study. The calculated ALE map had a cluster forming threshold of p < 0.001 with 1000 permutations, corrected for multiple comparisons using the Family Wise Error (FWE) p < 0.05 (20,37,43,44). Contrast and conjunction analyses were calculated to compare activation associated with task type, first by creating an image of two tasks pooled together (e.g., RH and LH) with an uncorrected threshold of p < 0.01 at 10,000 permutations, and then subtracting each original task analysis   from the pooled image in an iterative process (45,46). Contrast analyses permitted identification of regions of difference between groups while conjunction analyses quantify regions of overlap.
To maximize accurate localization and interpretation, images created in GingerALE were also imported into the SPM Anatomy Toolbox (47)(48)(49) to permit localization of the ALE images with 3-dimensional probabilistic cytoarchitectonic mapping (50,51). This regional cytoarchitectonic classification of ALE maps complements the GingerALE localization that uses peak MNI co-ordinates.

RESULTS
A total of n = 45 studies were determined to be suitable for inclusion (see Table 1). Of the 45 studies, 29 were used to perform the RH whole-brain meta-analysis, seven were used for the LH whole-brain meta-analysis (three studies involved stimulation of both LH and RH independently), and 12 studies examined RH stimulation in a ROI analysis. As can be seen in Table 1, for the 29 RH whole-brain studies, a total of n = 375 participants were included (n = 173 males, however n = 3 studies did not report sex) aged 18-76 years. The RH ROI studies included n = 123 participants (n = 63 males) aged 18-66 years. The seven LH studies included n = 85 participants (n = 52 males) aged 18-80 years. The most common form of stimulation was vibration (n = 12 studies), followed by compressed air (n = 8), textures (n = 7), brush stroke (n = 7), Von Frey filaments (n = 4), and pressure (n = 4).
The RH whole-brain studies, RH ROI studies and the LH whole-brain studies were analyzed separately, as presented in Table 2 and Figure 2. For the RH whole-brain stimulation studies, the contralateral (left) primary and secondary somatosensory areas were significant, with a large cluster containing the parietal operculum (92), somatosensory (93), and motor (94) cortices. The ipsilateral (right) secondary somatosensory cortex, S2, was also significant, largely comprising the parietal operculum (92) and inferior parietal cortex (95,96), in addition to a small cluster in the anterior cingulate. The RH ROI studies revealed visually smaller contralateral (left) clusters in the primary and secondary somatosensory regions, with a smaller ipsilateral (right) cluster within S2. The contralateral (left) clusters were separated into a large superior cluster containing the primary somatosensory (93) and motor (94) cortices, and a smaller inferior cluster containing primarily the parietal operculum (92). The ipsilateral (right) cluster contained similar areas to RH whole brain, namely the parietal operculum (92) and inferior parietal cortex (95,96). With the small number of LH stimulation studies, only two clusters were significant in the contralateral (right) primary (containing somatosensory (93) and motor (94) cortices) and secondary somatosensory regions [primarily parietal operculum (92)] and primary auditory cortex (97) ( Table 2).
Contrast analyses were then performed, as presented in Table 3 and Figure 3. When contrasted with the LH whole-brain studies, RH whole-brain studies revealed two clusters in the contralateral (left) primary and secondary somatosensory areas. The largest cluster contained primary somatosensory (93,98) and motor (94) cortices, while the smaller cluster contained primarily the auditory cortex (97), insula (99), and parietal operculum (92). When contrasted with RH whole-brain studies, the LH whole-brain studies activated three small clusters in S1 quite similar to those found in the standalone LH whole-brain analysis, all containing primary somatosensory areas (93,98). In our analysis of conjoined areas (i.e. areas of overlap) for RH and LH whole-brain studies, only one significant cluster was present, in the right secondary somatosensory region, primarily parietal operculum areas OP1, OP3 and OP4 (92). There were no significant differences in the contrast analysis between RH whole-brain studies and RH ROI studies. However, when the two groups were conjoined, significant common regions of activation were identified, with clusters revealed in the left primary (93) and secondary somatosensory areas (92), and the right secondary somatosensory area (92), including OP1.

DISCUSSION
In two important ways our ALE meta-analysis allowed us to examine the brain regions consistently activated during tactile stimulation of the hands in order to characterize functional somatosensory regions and networks, without the influence of motor function. Firstly, the meta-analysis allowed us to characterize and compare areas involved in right hand and left hand tactile stimulation studies separately. Secondly, it revealed the similarities and differences between functional activation studies that focus on specific brain regions (RH ROI studies) and what is actually occurring throughout the brain (RH wholebrain studies). Unfortunately very few studies (n = 7) examined LH stimulation separate to the RH and without the influence of motor activity, making a statistical comparison between the hands difficult and exploratory. For the RH whole-brain stimulation studies (n = 29) not only did we find two large clusters in the contralateral (left) primary (93) and secondary [specifically within parietal operculum areas OP1, OP3 and OP4 (92)] somatosensory cortices as expected, but activation was also revealed in the ipsilateral (right) secondary somatosensory region involving OP1 and OP4 (92) in addition to the anterior cingulate. Bilateral activation of secondary somatosensory S2 region, involving parietal operculum (92) to unilateral stimulation of the right hand is consistent with previous reports (100). From the few LH studies included, two small but significant clusters were revealed in the contralateral (right) S1 and S2. While each hand had significantly greater activation in the contralateral S1 and S2 in comparison to the other hand, the only significant area of overlap was in the right S2, specifically OP1, OP3 and OP4 (92). Lateralized differences have been reported for different sensory modalities, with right hemisphere being more spatially oriented toward the dorsal perceptual/sensory systems (101). Overlap in right S2 is consistent with hemispheric asymmetry involving righthemisphere-based bilateral representation of the body (101), right-sided asymmetry for tactile processing (102) and robust bilateral responses to unilateral stimulation in S2 (100). Due to the difference in numbers of studies included for each hand, this comparison is considered exploratory and highlights the need for more studies to examine LH tactile stimulation separately. Nevertheless, it is an interesting trend and could have significant implications for better understanding somatosensory function and dysfunction.
Activation in the contralateral S1, when using a tactile stimulus on the hand, is quite consistent with previous research (13,14). The pattern of activation shown in the RH wholebrain studies is consistent with research showing contralateral S1 activation only, and studies that have shown bilateral activation in S2 regardless of the hand being stimulated (16,17,103). It is surprising that bilateral S2 activation was not seen for the meta-analysis of LH studies also. However, this may have been attributable to the low number of studies stimulating the LH alone.
The role of S2 both contralateral and ipsilateral to the hand being stimulated is particularly interesting and may have important implications. The secondary somatosensory cortex of nonhuman primates is located on the parietal operculum, and the anatomical cytoarchitectonic maps of OP 1-4 of the human parietal operculum correlate with the functionally defined human somatosensory cortex (92), with OP 1 constituting the putative human homologue of area S2 (92). Further, OP1 is closely connected with the parietal networks for higher order somatosensory processing, while OP 4 is more closely integrated 3 | Anatomical location, summary statistics and MNI coordinates of ALE identified areas for contrast analyses: RH Whole-Brain greater than LH Whole-Brain, LH Whole-Brain greater than RH Whole-Brain, RH Whole-Brain conjoined with LH Whole-Brain, and RH Whole-Brain conjoined with RH ROI studies (p < 0.01, 10,000 p-value permutations, 100 mm cluster threshold).  with areas responsible for basic sensorimotor processing and action control (104). Bilateral secondary somatosensory cortex, in particular, has demonstrated a role in complex integrative processes of stimulus elaboration and attention following stimulation of right hand (105). Tame, Braun (103) have demonstrated bilateral activation in both S1 and S2 regardless of which hand was stimulated, suggesting that these areas may be involved in integrating somatosensory input from both sides of the body. Some may attribute the involvement of ipsilateral S2 to a more cognitive role in sensory processing, and while it is important to consider the cognitive aspects of sensorimotor control, such as planning and strategy (106), bilateral S2 activation has been demonstrated in somatosensory studies regardless of the level of cognitive demand (18).

SPM Anatomy Toolbox region location MNI GingerALE peak location
The involvement of S2 is particularly interesting in the context of aging, somatosensory dysfunction, and sensory rehabilitation. Age-related changes in activation have been seen, with decreased activation in S2 with tactile stimulation evident in elderly participants who are known to experience behavioral decline in somatosensory thresholds (54). The relationship of bilateral S2 with tactile sensation must also be considered in fields such as stroke research, where the location of the lesion has been demonstrated to impact both the type of somatosensory dysfunction (107), and also the ability to recover after stroke (108). Our finding of overlap in activation of right secondary somatosensory region for RH and LH tactile stimulation, may have particular relevance after stroke. For example, FIGURE 3 | ALE images displayed in neurological convention. (A) RH whole-brain activation greater than LH whole-brain activation; (B) LH whole-brain activation greater than RH whole-brain activation; (C) RH whole-brain conjoined with LH whole brain activation; (D) RH whole-brain conjoined with RH ROI activation. a stroke survivor with an infarct in the right hemisphere affecting S2 might not only experience the typically expected impairment of sensation in the contralateral hand (i.e., LH), but also impairment in the ipsilateral right hand; as has been described clinically (2). Further, recent evidence of altered functional connectivity in stroke survivors with impaired touch sensation following left or right hemisphere lesions, highlighted increased laterality indices in ipsilateral (contralesional) S2 relative to healthy controls following lesion of either hemisphere (109). Further, functional connectivity research has demonstrated that an increase in connectivity from contralesional S2 to contralesional thalamus correlates with better somatosensory function 6-months post-stroke (110).
Evaluation of the RH ROI studies (n = 12) revealed that only contralateral (left) S1 and bilateral S2 were examined by studies which predefined the areas thought to be involved in somatosensory processing of the hand. In comparison, the RH whole-brain studies also revealed anterior cingulate activation, and much larger clusters were involved with tactile stimulation. This suggests that when researchers set out to examine the functional activation of a tactile stimulus, if they limit the focus to a-priori areas, this may not capture the entire neural functional process related to the sensation. Anterior cingulate activation may play a significant role in sensory processing. For example, pleasant human touch is represented in anterior cingulate cortex (111). In addition, while attention differentially modulates signal amplitudes in the human somatosensory cortex, at higher intensities activation is also seen in the anterior cingulate cortex, consistent with attention to tactile stimuli in the current studies (112). It has been suggested that Von Economo (spindle) neurons found in cingulate cortex (113), and linked with insula, may have a role as part of a salience network (114). Network analyses identify anterior cingulate as a hub region and common co-activation of anterior cingulate and insula support the interpretation of a saliency network devoted to the integration of information from internal and external sensory environments (115). Further, interhemispheric connections between bilateral thalami occur via the anterior cingulate (113) and healthy controls show interhemispheric functional connectivity between a number of regions associated with somatosensory processing, including anterior cingulate (107), highlighting the contribution of both hemispheres and the broader somatosensory system. Interestingly, cingulate cortex has also been implicated in rats sensory recovery after lesions (116).
Other areas identified in this meta-analysis included inferior parietal lobe, insula, supramarginal gyrus and temporal lobe. Inferior parietal lobe (IPL) of the right hemisphere was identified for both RH whole-brain and ROI analyses. The location included OP1 and OP4. IPL has been associated with multimodal sensory information integration (117,118) and is reported to be part of the larger somatosensory network (119). The insula was also identified using the GingerALE peak maps, although this region was frequently assigned to the parietal operculum using the Anatomical toolbox. The insula has been identified as having a role in recognition, perception and learning in functional models of somatosensory processing (2). S2 is reciprocally connected with granular fields of the insula, reported to be devoted to somatic processing in monkeys (120). The close proximity of locations highlight the importance of the combined parietal opercular-insula region. Supramarginal gyrus is similarly located close to the parietal operculum/S2 region. The SMG is part of the somatosensory association cortex which has a role in interpretation of tactile sensory information as well as in perception of space and limbs location (121). Right SMG was found for RH whole-brain, RH ROI and LH whole-brain, and for the conjoined analyses. Right SMG is associated with spatial processing (121), consistent with tasks requiring localization of stimuli and/or involving spatial features of textures. Activation of left temporal gyrus, including auditory cortex and granular insula area 2, was greater in RH than LH whole-brain studies. Left temporal cortex has been linked with structural and semantic knowledge of body representation (122).
Each of the regions identified above have been implicated in stroke tactile impairment and recovery, potentially highlighting their broader importance. For example, change in functional connectivity from ipsilesional right S1 to right inferior parietal lobe was found in stroke survivors with impaired touch sensation compared to healthy controls (109). In addition, increased interhemispheric connectivity between the S2 region of interest and somatosensory association cortex (involving insula, parietal operculum and SMG) and temporal gyrus was found in healthy age-matched controls compared to stroke survivors with tactile deficits (109). Further, following tactile training, patients with lesions of sensory thalamus and/or internal capsule demonstrated activation in ipsiliesional insula, extending to the temporal pole, and supramarginal gyrus post-intervention (108). Interestingly, the regions identified have a role in the broader interpretation of tactile stimuli, including multi-modal integration, perception and learning, spatial processing and semantic knowledge and appear to be accessed as part of a wider somatosensory network.
There are limitations to this meta-analysis when examining the demographic information regarding the participants (see Table 1). Most of the LH studies (with the exception of one) included young participants (18-43 years). Aside from this, the cohorts were fairly well controlled, with the majority being RH dominant, and with tasks controlled for motor and other influences. Variable naming across studies can also contribute to confusion with interpretation. For example, terms such as secondary somatosensory cortex, secondary somatosensory region and secondary somatosensory area are often used interchangeably, although differences have been defined (10). To maximize accuracy and comparison across studies and broader literature in the field, we have reported on the MNI co-ordinates and peak location ALE results as well as the Anatomy Toolbox regional activation results.
The aim of this meta-analysis was to determine the convergence of foci reported from functional neuroimaging studies of touch sensation, separate to motor contributions and/or confounds. The findings advance our understanding of the separate, but potentially complementary, contributions of brain regions involved in processing touch sensation. Given the role of somatosensation and the somatosensory system in goaldirected actions of the upper limb and recovery after stroke, in depth knowledge of the role of key regions in the network is critical. The importance of bilateral S2 activation with right hand touch stimulation is highlighted, with a potential lateralization of activation in right S2 for right and left hand stimulation. This has implication for possible differences in unilateral vs. bilateral patterns of somatosensory impairment following right or left hemisphere lesion stroke. It may also identify a region with scope to contribute to recovery.
In conclusion, while research has established a role for S1 and S2 contralateral to the hand being stimulated (13,14), this metaanalysis has demonstrated the need to also examine the bilateral activation in S2 with right hand stimulation in order to further delineate the role of this area in tactile processing. Additional studies examining LH tactile processing separate to the RH would be beneficial to further examine whether this same pattern of activation is seen. These two advances in understanding would in turn further research into somatosensory dysfunction and rehabilitation.

AUTHOR CONTRIBUTIONS
LC conceived the study. GL, LC, EL, and AB contributed to the design of the study. GL, SP, EL, and AB conducted the search and extraction of data. GL and PG conducted the meta-analysis. GL and LC interpreted the findings and drafted the manuscript. All authors critically reviewed and revised the manuscript for important intellectual content, provided approval for publication and agree to be accountable for all aspects of the work.