High-resolution fMRI investigations of the fingertip somatotopy and variability in BA3b and BA1 of the primary somatosensory cortex
Introduction
S1 consists of four BAs which differ in their cytoarchitectonic properties (Brodmann, 1909). The cortical sensory representation areas of the stimulated body surface were defined as the regions with the largest electrical response amplitude, based on the investigation of single neurons with electrophysiological methods. A separate somatotopic arrangement of the contralateral representation areas was found in BA1 and BA3b, the fingertips lined up sequentially along the postcentral gyrus, with D1 (thumb) located most inferior (Kaas et al., 1979). Optical BOLD imaging was used to reveal the large-scale structure of the activation pattern in monkey cortices invasively based on the hemodynamic coupling to neural activity, exhibiting a modular structure of the representation areas in BA3b (Chen et al., 2001); a structure similar to the continuous V1 visual orientation maps in BA1 (Friedman et al., 2004), where the frequency of the stimulation is the analogous property to the orientation. This corresponds to larger fingertip representations areas in BA1 with a higher overlap compared to BA3b and a more frequent response to multiple fingers of the fingertip areas in BA1, as found in electrophysiological studies (Iwamura et al., 1980, Iwamura et al., 1983b).
fMRI investigations of the human S1 fingertip somatotopy are also based on the BOLD effect and the hemodynamic coupling. The distances between the fingertip positions revealed important associations to behavioral properties in healthy participants and patients. Examples are the associations with the tactile hyper acuity thresholds in healthy participants (Duncan and Boynton, 2007) and CRPS in patients (Di Pietro et al., 2015). In those studies S1 was often not separated into BAs. Recent 7-Tesla fMRI examinations (Besle et al., 2014, Martuzzi et al., 2014) found separate fingertip representations across the BAs in S1 and raised the question whether analogous examinations are possible with 3-Tesla MRI. The cortex in S1 possesses a complex folding pattern (Hopkins et al., 2014). However, the inter finger distances were in most examinations determined using 3D Euclidean distances, which neglect the cortical folding and may obscure fine-scale details. Virtually all fMRI studies on the fingertip somatotopy observed a large inter individual variability. However, whether the observed variability is due to anatomical or functional differences between the participants has not been clarified.
Our examinations address those issues using an automated approach for investigations of the human fingertip somatotopy in BA3b (Pfannmöller et al., 2016) and extended the approach to include BA1. We investigate differences between the fingertip representation in BA1 and BA3b, using cytoarchitectonic probability maps for separation of S1 into BAs (Fischl et al., 2008). In order to maximize the mapping of fine-scale details and to include the cortical anatomy into the distance determination, we compare three different methods for the determination of cortical distances. Exact shortest distances are applied as a gold standard method for the determination of distances along the cortical surface (Balasubramanian et al., 2009). We assess whether Dijkstra (Dijkstra, 1959) or 3D Euclidean distances are valid approximations. Furthermore, we use the nonlinear surface-based FreeSurfer normalization (Fischl et al., 1999) to investigate the relative contribution of the anatomical variability to the total variability of the fingertip distances. This allows us to make inferences about the functional origin of the variability in the somatotopy between participants and on its left- and right-hand symmetry, a central quantity in examinations involving the fingertip somatotopy.
Section snippets
Participants and tactile stimulation
We investigated 18 healthy participants [mean age: 30 ± 11 years; 11 women; average handedness score of 92; range 77–100 (Oldfield, 1971)]. All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional research committee, and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. The authors declare that they have no actual or potential conflicts of interest with this work. Informed consent
Fingertip representation in BA3b and BA1
Concerning the peak value approach: As shown in Table 1, at least 78% of the participants exhibited artifact-free cortical maps in BA3b, while only 56% of the participants revealed artifact-free maps in BA1. In more than 80% of the valid BA3b maps the fingertips revealed a sequential arrangement along the postcentral gyrus, while there were less than 30% for BA1. The fingertip positions in BA3b were confined to a specific sub region after nonlinear normalization to FSaverage, as depicted in Fig.
Fingertip representation in BA3b and BA1
Previous fMRI investigations focused on the BA-specific somatotopies (Besle et al., 2014, Martuzzi et al., 2014) and noticed a difference between the maps for the fingertip representation across the BAs in S1. We investigated the human fingertip somatotopy of D1, D2 and D3 in BA3b and BA1. An excellent agreement with the expected pattern for the fingertip representations was found in BA3b applying the peak value approach and a good agreement in BA1 in the center of mass approach. Although the
Conclusion
Our results provide support for a modular functional structure in BA3b and an orientation column-like structure analogous to V1 in BA1, as found in monkey cortices. Valid surface distances were found to improve the fine-scale feature mapping in the somatotopy considerably. Application of exact shortest distances and nonlinear normalization allowed disentangling anatomical and functional contributions to the variability in the somatotopy. We found a dominant contribution of functional
Acknowledgements
M. Lotze was funded by a grant from the German Research Community (DFG Lo795/7-1). Markus Oelschläger is acknowledged for assistance during data acquisition. J. P. Pfannmöller participated in data acquisition, data evaluation and preparation of the examinations. M. Greiner and M. Balasubramanian participated in data evaluation. M. Lotze participated in the preparation of the examinations. All authors contributed to the preparation of the manuscript.
References (43)
- et al.
What aspect of the fMRI BOLD signal best reflects the underlying electrophysiology in human somatosensory cortex?
Clin Neurophysiol
(2003) - et al.
A new SPM toolbox for combining probabilistic cytoarchitectonic maps and functional imaging data
NeuroImage
(2005) - et al.
Areas 3a, 3b, and 1 of human primary somatosensory cortex 1. Microstructural organization and interindividual variability
NeuroImage
(1999) - et al.
Accurate and robust brain image alignment using boundary-based registration
NeuroImage
(2009) - et al.
Overlapping representation of fingers in the somatosensory cortex (area 2) of the conscious monkey
Brain Res
(1980) - et al.
Evaluation of 14 nonlinear deformation algorithms applied to human brain MRI registration
NeuroImage
(2009) - et al.
The cutaneous sensory system
Neurosci Biobehav R
(2010) The great brain versus vein debate
NeuroImage
(2012)The assessment and analysis of handedness: the Edinburgh inventory
Neuropsychologia
(1971)- et al.
Functional MRI indicates consistent intra-digit topographic maps in the little but not the index finger within the human primary somatosensory cortex
NeuroImage
(2011)
Finger representations in human primary somatosensory cortex as revealed by high-resolution functional MRI of tactile stimulation
NeuroImage
Evaluation of 2D multiband EPI imaging for high-resolution, whole-brain, task-based fMRI studies at 3T: sensitivity and slice leakage artifacts
NeuroImage
How much cortex can a vein drain? Downstream dilution of activation-related cerebral blood oxygenation changes
NeuroImage
Brain morphometry with multiecho MPRAGE
NeuroImage
Optimizing the mapping of finger areas in primary somatosensory cortex using functional MRI
Magn Reson Imaging
Intra- and inter- subject variability of high field fmri digit maps in somatosensory area 3b of new world monkeys
Neuroscience
Exact geodesics and shortest paths on polyhedral surfaces
IEEE T Pattern Anal
Controlling the false discovery rate – a practical and powerful approach to multiple testing
J Roy Stat Soc B Met
Event-related fMRI at 7T reveals overlapping cortical representations for adjacent fingertips in S1 of individual subjects
Hum. Brain Mapp.
Vergleichende Lokalisationslehre der Grosshirnrinde in ihren Prinzipien dargestellt auf Grund des Zellenbaues
Fine-scale organization of SI (Area 3b) in the squirrel monkey revealed with intrinsic optical imaging
J Neurophysiol
Cited by (20)
Primary somatosensory cortex organization for engineering artificial somatosensation
2024, Neuroscience ResearchFunctional characteristics of the human primary somatosensory cortex: An electrostimulation study
2021, Epilepsy and BehaviorCitation Excerpt :Thus, the S1 is characterized by heterogeneous functional organizations across BAs. Recent advances in neuroimaging (PET, fMRI, etc.) [11-17] and electrophysiology (EEG, MEG) [18-21] techniques increasingly applied in lots of studies, to some extent, have opened a window to the functional characteristics of the deep surface of S1, but do not identify directly the necessity of a given cortex area in a particular function. Electrical stimulation still has unique advantages and critical contribution to producing highly reliable outcomes at well-defined anatomical sites [1,22-26].
Functional organization of the human primary somatosensory cortex: A stereo-electroencephalography study
2021, Clinical NeurophysiologyCitation Excerpt :However, the depth of the central sulcus is approximately twice the width of the postcentral gyrus (Woolsey, 1979), which means that the major component of S1 is the sulcal cortex, including BA3 and BA2. Although recent functional imaging studies (Akselrod et al., 2017; Kurth et al., 2000, 1998; Martuzzi et al., 2014; Pfannmoller et al., 2016) have made relevant contributions to the functional investigation of the folded cortex, it is yet less superior to arbitrarily explored a specific area of the brain than DES (Desmurget et al., 2013). Different from previous studies, we first utilized SEEG technology and DES to explore the functional organization of the entire S1.
Cerebral plasticity as the basis for upper limb recovery following brain damage
2019, Neuroscience and Biobehavioral ReviewsAnatomical and functional properties of the foot and leg representation in areas 3b, 1 and 2 of primary somatosensory cortex in humans: A 7T fMRI study
2017, NeuroImageCitation Excerpt :Ultra-high field fMRI has already been applied to the study of the retinotopic organization of the human primary visual cortex (Olman et al., 2010), the tonotopic organization of the human primary auditory cortex (Da Costa et al., 2011), as well as the somatotopic organization in primary somatosensory cortex (S1) (Besle et al., 2014; Martuzzi et al., 2014). Considering the large inter-subject variability of S1 somatotopy, it has been proposed that, in particular, S1 responses should be studied at the single subject level (Pfannmöller et al., 2016), which can be achieved with high resolution at 7 T. However, previous somatotopic descriptions of human S1 were so far limited to the representations of single fingers with only few neuroimaging studies targeting important S1 representations of body parts other than the hand, such as the foot and the leg (Bao et al., 2012; Nakagoshi et al., 2005; Saadon-Grosman et al., 2015). S1 is somatotopically organized to form a cortical map of the body, meaning that there is a point-to-point correspondence between the skin and S1 representations (Penfield and Boldrey, 1937; Rasmussen and Penfield, 1947).