Regional structural differences across functionally parcellated Brodmann areas of human primary somatosensory cortex

Cortical field maps across human sensory cortex

Cortical processing pathways for sensory information in the mammalian brain tend to be organized into topographical representations that encode various fundamental sensory dimensions. Numerous laboratories have now shown how these representations are organized into numerous cortical field maps (CMFs) across visual and auditory cortex, with each CFM supporting a specialized computation or set of computations that underlie the associated perceptual behaviors. An individual CFM is defined by two orthogonal topographical gradients that reflect two essential aspects of feature space for that sense. Multiple adjacent CFMs are then organized across visual and auditory cortex into macrostructural patterns termed cloverleaf clusters. CFMs within cloverleaf clusters are thought to share properties such as receptive field distribution, cortical magnification, and processing specialization. Recent measurements point to the likely existence of CFMs in the other senses, as well, with topographical representations of at least one sensory dimension demonstrated in somatosensory, gustatory, and possibly olfactory cortical pathways. Here we discuss the evidence for CFM and cloverleaf cluster organization across human sensory cortex as well as approaches used to identify such organizational patterns. Knowledge of how these topographical representations are organized across cortex provides us with insight into how our conscious perceptions are created from our basic sensory inputs. In addition, studying how these representations change during development, trauma, and disease serves as an important tool for developing improvements in clinical therapies and rehabilitation for sensory deficits.

Two-Dimensional Population Receptive Field Mapping of Human Primary Somatosensory Cortex

Functional magnetic resonance imaging can provide detailed maps of how sensory space is mapped in the human brain. Here, we use a novel 16 stimulator setup (a 4 × 4 grid) to measure two-dimensional sensory maps of between and within-digit (D2-D4) space using high spatial-resolution (1.25 mm isotropic) imaging at 7 Tesla together with population receptive field (pRF) mapping in 10 participants. Using a 2D Gaussian pRF model, we capture maps of the coverage of digits D2-D5 across Brodmann areas and estimate pRF size and shape. In addition, we compare results to previous studies that used fewer stimulators by constraining pRF models to a 1D Gaussian Between Digit or 1D Gaussian Within Digit model. We show that pRFs across somatosensory areas tend to have a strong preference to cover the within-digit axis. We show an increase in pRF size moving from D2-D5. We quantify pRF shapes in Brodmann area (BA) 3b, 3a, 1, 2 and show differences in pRF size in Brodmann areas 3a-2, with larger estimates for BA2. Generally, the 2D Gaussian pRF model better represents pRF coverage maps generated by our data, which itself is produced from a 2D stimulation grid.

The 3D Structural Architecture of the Human Hand Area Is Nontopographic

The functional topography of the human primary somatosensory cortex hand area is a widely studied model system to understand sensory organization and plasticity. It is so far unclear whether the underlying 3D structural architecture also shows a topographic organization. We used 7 Tesla (7T) magnetic resonance imaging (MRI) data to quantify layer-specific myelin, iron, and mineralization in relation to population receptive field maps of individual finger representations in Brodman area 3b (BA 3b) of human S1 in female and male younger adults. This 3D description allowed us to identify a characteristic profile of layer-specific myelin and iron deposition in the BA 3b hand area, but revealed an absence of structural differences, an absence of low-myelin borders, and high similarity of 3D microstructure profiles between individual fingers. However, structural differences and borders were detected between the hand and face areas. We conclude that the 3D structural architecture of the human hand area is nontopographic, unlike in some monkey species, which suggests a high degree of flexibility for functional finger organization and a new perspective on human topographic plasticity.SIGNIFICANCE STATEMENT Using ultra-high-field MRI, we provide the first comprehensive in vivo description of the 3D structural architecture of the human BA 3b hand area in relation to functional population receptive field maps. High similarity of precise finger-specific 3D profiles, together with an absence of structural differences and an absence of low-myelin borders between individual fingers, reveals the 3D structural architecture of the human hand area to be nontopographic. This suggests reduced structural limitations to cortical plasticity and reorganization and allows for shared representational features across fingers.

Empirical transmit field bias correction of T1w/T2w myelin maps

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

Metabotropic Glutamate Receptor 5 in the Dysgranular Zone of Primary Somatosensory Cortex Mediates Neuropathic Pain in Rats

The primary somatosensory cortex (S1) plays a key role in the discrimination of somatic sensations. Among subdivisions in S1, the dysgranular zone of rodent S1 (S1DZ) is homologous to Brodmann's area 3a of primate S1, which is involved in the processing of noxious signals from the body. However, molecular changes in this region and their role in the pathological pain state have never been studied. In this study, we identified molecular alteration of the S1DZ in a rat model of neuropathic pain induced by right L5 spinal nerve ligation (SNL) surgery and investigated its functional role in pain symptoms. Brain images acquired from SNL group and control group in our previous study were analyzed, and behaviors were measured using the von Frey test, acetone test, and conditioned place preference test. We found that metabotropic glutamate receptor 5 (mGluR5) levels were significantly upregulated in the S1DZ contralateral to the nerve injury in the SNL group compared to the sham group. Pharmacological deactivation of mGluR5 in S1DZ ameliorated symptoms of neuropathic allodynia, which was shown by a significant increase in the mechanical paw withdrawal threshold and a decrease in the behavioral response to cold stimuli. We further confirmed that this treatment induced relief from the tonic-aversive state of chronic neuropathic pain, as a place preference memory associated with the treatment-paired chamber was formed in rats with neuropathic pain. Our data provide evidence that mGluR5 in the S1DZ is involved in the manifestation of abnormal pain sensations in the neuropathic pain state.

Somatotopic Mapping of the Fingers in the Somatosensory Cortex Using Functional Magnetic Resonance Imaging: A Review of Literature

Multiple studies have demonstrated finger somatotopy in humans and other primates using a variety of brain mapping techniques including functional magnetic resonance imaging (fMRI). Here, we review the literature to better understand the reliability of fMRI for mapping the somatosensory cortex. We have chosen to focus on the hand and fingers as these areas have the largest representation and have been the subject of the largest number of somatotopic mapping experiments. Regardless of the methods used, individual finger somatosensory maps were found to be organized across Brodmann areas (BAs) 3b, 1, and 2 in lateral-to-medial and inferior-to-superior fashion moving from the thumb to the pinky. However, some consistent discrepancies are found that depend principally on the method used to stimulate the hand and fingers. Therefore, we suggest that a comparative analysis of different types of stimulation be performed to address the differences described in this review.

Supraspinal Mechanisms Underlying Ocular Pain

Supraspinal mechanisms of pain are increasingly understood to underlie neuropathic ocular conditions previously thought to be exclusively peripheral in nature. Isolating individual causes of centralized chronic conditions and differentiating them is critical to understanding the mechanisms underlying neuropathic eye pain and ultimately its treatment. Though few functional imaging studies have focused on the eye as an end-organ for the transduction of noxious stimuli, the brain networks related to pain processing have been extensively studied with functional neuroimaging over the past 20 years. This article will review the supraspinal mechanisms that underlie pain as they relate to the eye.

Analysis of Intracerebral Activity during Reflex Locomotion Stimulation According to Vojta’s Principle

Vojta’s therapy is a widely used approach in both the prevention and therapy of musculoskeletal disorders. Changes in the musculoskeletal system have been described repeatedly, but the principles of the approach have not yet been clarified. The objective of our study was to evaluate changes of intracerebral activity using electromagnetic tomography (sLORETA) that arise during reflex locomotion stimulation of the breast trigger zone according to Vojta’s therapy. Seventeen healthy women took part in the experiment (aged 20–30 years old). EEG activity was recorded 5 min prior to the reflex locomotion stimulation, during stimulation, and 5 min after the stimulation. The obtained data were subsequently processed in the sLORETA program and statistically evaluated at the significance level p ≤ 0.05. The analysis found statistically significant differences in the frequency bands alpha-2, beta-1, and beta-2 between the condition prior to stimulation and the actual stimulation in BAs 6, 7, 23, 24, and 31 and between the resting condition prior to stimulation, and the condition after the stimulation was terminated in the frequency bands alpha-1, alpha-2, beta-1, and beta-2 in BAs 3, 4, 6, and 24. The results showed that reflex locomotion stimulation according to Vojta’s therapy modulates electrical activity in the brain areas responsible for movement planning and regulating and performing the movement.

Population Receptive Field Characteristics in the between- and Within-Digit Dimensions of the Undominant Hand in the Primary Somatosensory Cortex

Somatotopy is an important guiding principle for sensory fiber organization in the primary somatosensory cortex (S1), which reflects tactile information processing and is associated with disease-related reorganization. However, it is difficult to measure the neuronal encoding scheme in S1 in vivo in normal participants. Here, we investigated the somatotopic map of the undominant hand using a Bayesian population receptive field (pRF) model. The model was established in hand space with between- and within-digit dimensions. In the between-digit dimension, orderly representation was found, which had low variability across participants. The pRF shape tended to be elliptical for digits with high spatial acuity, for which the long axis was along the within-digit dimension. In addition, the pRF width showed different change trends in the 2 dimensions across digits. These results provide new insights into the neural mechanisms in S1, allowing for in-depth investigation of somatosensory information processing and disease-related reorganization.

Human Somatosensory Processing and Artificial Somatosensation

In the past few years, we have gained a better understanding of the information processing mechanism in the human brain, which has led to advances in artificial intelligence and humanoid robots. However, among the various sensory systems, studying the somatosensory system presents the greatest challenge. Here, we provide a comprehensive review of the human somatosensory system and its corresponding applications in artificial systems. Due to the uniqueness of the human hand in integrating receptor and actuator functions, we focused on the role of the somatosensory system in object recognition and action guidance. First, the low-threshold mechanoreceptors in the human skin and somatotopic organization principles along the ascending pathway, which are fundamental to artificial skin, were summarized. Second, we discuss high-level brain areas, which interacted with each other in the haptic object recognition. Based on this close-loop route, we used prosthetic upper limbs as an example to highlight the importance of somatosensory information. Finally, we present prospective research directions for human haptic perception, which could guide the development of artificial somatosensory systems.

A touch of hierarchy: population receptive fields reveal fingertip integration in Brodmann areas in human primary somatosensory cortex

Several neuroimaging studies have shown the somatotopy of body part representations in primary somatosensory cortex (S1), but the functional hierarchy of distinct subregions in human S1 has not been adequately addressed. The current study investigates the functional hierarchy of cyto-architectonically distinct regions, Brodmann areas BA3, BA1, and BA2, in human S1. During functional MRI experiments, we presented participants with vibrotactile stimulation of the fingertips at three different vibration frequencies. Using population Receptive Field (pRF) modeling of the fMRI BOLD activity, we identified the hand region in S1 and the somatotopy of the fingertips. For each voxel, the pRF center indicates the finger that most effectively drives the BOLD signal, and the pRF size measures the spatial somatic pooling of fingertips. We find a systematic relationship of pRF sizes from lower-order areas to higher-order areas. Specifically, we found that pRF sizes are smallest in BA3, increase slightly towards BA1, and are largest in BA2, paralleling the increase in visual receptive field size as one ascends the visual hierarchy. Additionally, we find that the time-to-peak of the hemodynamic response in BA3 is roughly 0.5 s earlier compared to BA1 and BA2, further supporting the notion of a functional hierarchy of subregions in S1. These results were obtained during stimulation of different mechanoreceptors, suggesting that different afferent fibers leading up to S1 feed into the same cortical hierarchy.

Tactile acuity of fingertips and hand representation size in human Area 3b and Area 1 of the primary somatosensory cortex

Intracortical mapping in monkeys revealed a full body map in all four cytoarchitectonic subdivisions of the contralateral primary somatosensory cortex (S1), as well as positive associations between spatio-tactile acuity performance of the fingers and their representation field size especially within cytoarchitectonic Area 3b and Area 1. Previous non-invasive investigations on these associations in humans assumed a monotonous decrease of representation field size from index finger to little finger although the field sizes are known to change in response to training or in disease. Recent developments improved noninvasive functional mapping of S1 by a) adding a cognitive task during repetitive stimulation to decrease habituation to the stimuli, b) smaller voxel size of fMRI-sequences, c) surface-based analysis accounting for cortical curvature, and d) increase of spatial specificity for fMRI data analysis by avoidance of smoothing, partial volume effects, and pial vein signals. We here applied repetitive pneumatic stimulation of digit 1 (D1; thumb) and digit 5 (D5; little finger) on both hands to investigate finger/hand representation maps in the complete S1, but also in cytoarchitectonic Areas 1, 2, 3a, and 3b separately, in 21 healthy volunteers using 3T fMRI. The distances between activation maxima of D1 and D5 were evaluated by two independent raters, blinded for performance parameters. The fingertip representations showed a somatotopy and were localized in the transition region between the crown and the anterior wall of the post central gyrus agreeing with Area 1 and 3b. Participants were comprehensively tested for tactile performance using von Freyhair filaments to determine cutaneous sensory thresholds (CST) as well as grating orientation thresholds (GOT) and two-point resolution (TPD) for spatio-tactile acuity testing. Motor performance was evaluated with pinch grip performance (Roeder test). We found bilateral associations of D1-D5 distance for GOT thresholds and partially also for TPD in Area 3b and in Area 1, but not if using the complete S1 mask. In conclusion, we here demonstrate that 3T fMRI is capable to map associations between spatio-tactile acuity and the fingertip representation in Area 3b and Area 1 in healthy participants.

Relation between palm and finger cortical representations in primary somatosensory cortex: A 7T fMRI study

Many studies focused on the cortical representations of fingers, while the palm is relatively neglected despite its importance for hand function. Here, we investigated palm representation (PR) and its relationship with finger representations (FRs) in primary somatosensory cortex (S1). Few studies in humans suggested that PR is located medially with respect to FRs in S1, yet to date, no study directly quantified the somatotopic organization of PR and the five FRs. Importantly, the link between the somatotopic organization of PR and FRs and their activation properties remains largely unexplored. Using 7T fMRI, we mapped PR and the five FRs at the single subject level. First, we analyzed the cortical distance between PR and FRs to determine their somatotopic organization. Results show that PR was located medially with respect to D5. Second, we tested whether the observed cortical distances would predict the relationship between PR and FRs activations. Using three complementary measures (cross-activations, pattern similarity and resting-state connectivity), we show that the relationship between PR and FRs activations were not determined by their somatotopic organization, that is, there was no gradient moving from D5 to D1, except for resting-state connectivity, which was predicted by the somatotopy. Instead, we show that the representational geometry of PR and FRs activations reflected the physical structure of the hand. Collectively, our findings suggest that the spatial proximity between topographically organized neuronal populations do not necessarily predicts their functional properties, rather the structure of the sensory space (e.g., the hand shape) better describes the observed results.

Somatotopic Arrangement of the Human Primary Somatosensory Cortex Derived From Functional Magnetic Resonance Imaging

Functional magnetic resonance imaging (fMRI) was used to estimate neuronal activity in the primary somatosensory cortex of six participants undergoing cutaneous tactile stimulation on skin areas spread across the entire body. Differences between the accepted somatotopic maps derived from Penfield's work and those generated by this fMRI study were sought, including representational transpositions or replications across the cortex. MR-safe pneumatic devices mimicking the action of a Wartenberg wheel supplied touch stimuli in eight areas. Seven were on the left side of the body: foot, lower, and upper leg, trunk beneath ribcage, anterior forearm, middle fingertip, and neck above the collarbone. The eighth area was the glabella. Activation magnitude was estimated as the maximum cross-correlation coefficient at a certain phase shift between ideal time series and measured blood oxygen level dependent (BOLD) time courses on the cortical surface. Maximally correlated clusters associated with each cutaneous area were calculated, and cortical magnification factors were estimated. Activity correlated to lower limb stimulation was observed in the paracentral lobule and superomedial postcentral region. Correlations to upper extremity stimulation were observed in the postcentral area adjacent to the motor hand knob. Activity correlated to trunk, face and neck stimulation was localized in the superomedial one-third of the postcentral region, which differed from Penfield's cortical homunculus.

Hemodynamic response varies across tactile stimuli with different temporal structures

Tactile stimuli can be distinguished based on their temporal features (e.g., duration, local frequency, and number of pulses), which are fundamental for vibrotactile frequency perception. Characterizing how the hemodynamic response changes in shape across experimental conditions is important for designing and interpreting fMRI studies on tactile information processing. In this study, we focused on periodic tactile stimuli with different temporal structures and explored the hemodynamic response function (HRF) induced by these stimuli. We found that HRFs were stimulus‐dependent in tactile‐related brain areas. Continuous stimuli induced a greater area of activation and a stronger and narrower hemodynamic response than intermittent stimuli with the same duration. The magnitude of the HRF increased with increasing stimulus duration. By normalizing the characteristics into topographic matrix, nonlinearity was obvious. These results suggested that stimulation patterns and duration within a cycle may be key characters for distinguishing different stimuli. We conclude that different temporal structures of tactile stimuli induced different HRFs, which are essential for vibrotactile perception and should be considered in fMRI experimental designs and analyses.

A probabilistic atlas of finger dominance in the primary somatosensory cortex

With the advent of ultra-high field (7T), high spatial resolution functional MRI (fMRI) has allowed the differentiation of the cortical representations of each of the digits at an individual-subject level in human primary somatosensory cortex (S1). Here we generate a probabilistic atlas of the contralateral SI representations of the digits of both the left and right hand in a group of 22 right-handed individuals. The atlas is generated in both volume and surface standardised spaces from somatotopic maps obtained by delivering vibrotactile stimulation to each distal phalangeal digit using a travelling wave paradigm. Metrics quantify the likelihood of a given position being assigned to a digit (full probability map) and the most probable digit for a given spatial location (maximum probability map). The atlas is validated using a leave-one-out cross validation procedure. Anatomical variance across the somatotopic map is also assessed to investigate whether the functional variability across subjects is coupled to structural differences. This probabilistic atlas quantifies the variability in digit representations in healthy subjects, finding some quantifiable separability between digits 2, 3 and 4, a complex overlapping relationship between digits 1 and 2, and little agreement of digit 5 across subjects. The atlas and constituent subject maps are available online for use as a reference in future neuroimaging studies.

A Data-Driven Multi-scale Technique for fMRI Mapping of the Human Somatosensory Cortex

A previously introduced Bayesian non-parametric multi-scale technique, called iterated Multigrid Priors (iMGP) method, is used to map the topographic organization of human primary somatosensory cortex (S1). We analyze high spatial resolution fMRI data acquired at ultra-high field (UHF, 7T) in individual subjects during vibrotactile stimulation applied to each distal phalange of the left hand digits using both a travelling-wave (TW) and event-related (ER) paradigm design. We compare the somatotopic digit representations generated in S1 using the iMGP method with those obtained using established fMRI paradigms and analysis techniques: Fourier-based analysis of travelling-wave data and General Linear Model (GLM) analysis of event-related data. Maps derived with the iMGP method are similar to those derived with the standard analysis, but in contrast to the Fourier-based analysis, the iMGP method reveals overlap of activity from adjacent digit representations in S1. These findings validate the use of the iMGP method as an alternative to study digit representations in S1, particularly with the TW design as an attractive means to study cortical reorganization in patient populations such dystonia and carpal tunnel syndrome, where the degree of spatial overlap of cortical finger representations is of interest.

A putative human homologue of the macaque area PEc

The cortical area PEc is anatomically and functionally well-defined in macaque, but it is unknown whether it has a counterpart in human. Since we know that macaque PEc, but not the nearby posterior regions, hosts a lower limb representation, in an attempt to recognize a possible human PEc we looked for the existence of leg representations in the human parietal cortex using individual cortical surface-based analysis, task-evoked paradigms and resting-state functional connectivity. fMRI images were acquired while thirty-one participants performed long-range leg movements through an in-house MRI-compatible set-up. We revealed the existence of multiple leg representations in the human dorsomedial parietal cortex, here defined as S-I (somatosensory-I), hPE (human PE, in the postcentral sulcus), and hPEc (human PEc, in the anterior precuneus). Among the three "leg" regions, hPEc had a unique functional profile, in that it was the only one responding to both arm and leg movements, to both hand-pointing and foot pointing movements, and to flow field visual stimulation, very similar to macaque area PEc. In addition, hPEc showed functional connections with the somatomotor regions hosting a lower limb representation, again as in macaque area PEc. Therefore, based on similarity in brain position, functional organization, cortical connections, and relationship with the neighboring areas, we propose that this cortical region is the human homologue of macaque area PEc.

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.

Addressing challenges of high spatial resolution UHF fMRI for group analysis of higher-order cognitive tasks: An inter-sensory task directing attention between visual and somatosensory domains

Functional MRI at ultra‐high field (UHF, ≥7 T) provides significant increases in BOLD contrast‐to‐noise ratio (CNR) compared with conventional field strength (3 T), and has been exploited for reduced field‐of‐view, high spatial resolution mapping of primary sensory areas. Applying these high spatial resolution methods to investigate whole brain functional responses to higher‐order cognitive tasks leads to a number of challenges, in particular how to perform robust group‐level statistical analyses. This study addresses these challenges using an inter‐sensory cognitive task which modulates top‐down attention at graded levels between the visual and somatosensory domains. At the individual level, highly focal functional activation to the task and task difficulty (modulated by attention levels) were detectable due to the high CNR at UHF. However, to assess group level effects, both anatomical and functional variability must be considered during analysis. We demonstrate the importance of surface over volume normalisation and the requirement of no spatial smoothing when assessing highly focal activity. Using novel group analysis on anatomically parcellated brain regions, we show that in higher cognitive areas (parietal and dorsal‐lateral‐prefrontal cortex) fMRI responses to graded attention levels were modulated quadratically, whilst in visual cortex and VIP, responses were modulated linearly. These group fMRI responses were not seen clearly using conventional second‐level GLM analyses, illustrating the limitations of a conventional approach when investigating such focal responses in higher cognitive regions which are more anatomically variable. The approaches demonstrated here complement other advanced analysis methods such as multivariate pattern analysis, allowing UHF to be fully exploited in cognitive neuroscience.

Neural Basis of Touch and Proprioception in Primate Cortex

The sense of proprioception allows us to keep track of our limb posture and movements and the sense of touch provides us with information about objects with which we come into contact. In both senses, mechanoreceptors convert the deformation of tissues-skin, muscles, tendons, ligaments, or joints-into neural signals. Tactile and proprioceptive signals are then relayed by the peripheral nerves to the central nervous system, where they are processed to give rise to percepts of objects and of the state of our body. In this review, we first examine briefly the receptors that mediate touch and proprioception, their associated nerve fibers, and pathways they follow to the cerebral cortex. We then provide an overview of the different cortical areas that process tactile and proprioceptive information. Next, we discuss how various features of objects-their shape, motion, and texture, for example-are encoded in the various cortical fields, and the susceptibility of these neural codes to attention and other forms of higher-order modulation. Finally, we summarize recent efforts to restore the senses of touch and proprioception by electrically stimulating somatosensory cortex. © 2018 American Physiological Society. Compr Physiol 8:1575-1602, 2018.

Parcellating Cerebral Cortex: How Invasive Animal Studies Inform Noninvasive Mapmaking in Humans

The cerebral cortex in mammals contains a mosaic of cortical areas that differ in function, architecture, connectivity, and/or topographic organization. A combination of local connectivity (within-area microcircuitry) and long-distance (between-area) connectivity enables each area to perform a unique set of computations. Some areas also have characteristic within-area mesoscale organization, reflecting specialized representations of distinct types of information. Cortical areas interact with one another to form functional networks that mediate behavior, and each area may be a part of multiple, partially overlapping networks. Given their importance to the understanding of brain organization, mapping cortical areas across species is a major objective of systems neuroscience and has been a century-long challenge. Here, we review recent progress in multi-modal mapping of mouse and nonhuman primate cortex, mainly using invasive experimental methods. These studies also provide a neuroanatomical foundation for mapping human cerebral cortex using noninvasive neuroimaging, including a new map of human cortical areas that we generated using a semiautomated analysis of high-quality, multimodal neuroimaging data. We contrast our semiautomated approach to human multimodal cortical mapping with various extant fully automated human brain parcellations that are based on only a single imaging modality and offer suggestions on how to best advance the noninvasive brain parcellation field. We discuss the limitations as well as the strengths of current noninvasive methods of mapping brain function, architecture, connectivity, and topography and of current approaches to mapping the brain's functional networks.

Somatotopy in the Human Somatosensory System

Previous functional magnetic resonance imaging (fMRI) studies have demonstrated digit somatotopy in primary somatosensory cortex (SI), and even shown that at high spatial resolution it is possible to resolve within-digit somatotopy. However, fMRI studies have failed to resolve the spatial organisation of digit representations in secondary somatosensory cortex (SII). One of the major limitations of high spatial resolution fMRI studies of the somatosensory system has been the long acquisition time needed to acquire slices spanning both SI and SII. Here, we exploit the increased blood oxygenation level dependent contrast of ultra-high-field (7 Tesla) fMRI and the use of multiband imaging to study the topographic organisation in SI and SII with high spatial resolution at the individual subject level. A total of n = 6 subjects underwent vibrotactile stimulation of their face, hand digits and foot (body imaging) and their individual hand digits (digit mapping) for each left and right sides of the body. In addition, n = 2 subjects participated only in the body imaging experiment on both their left and right sides. We show an orderly representation of the face, hand digits and foot in contralateral primary cortex in each individual subject. In SII, there is clear separation of the body areas of the face, hand and foot but the spatial organisation varies across individual subjects. However, separate representation of the individual digits of the hand in SII could not be resolved, even at the spatial resolution of 1.5 mm due to largely overlapping representations.

Detailed T1-Weighted Profiles from the Human Cortex Measured in Vivo at 3 Tesla MRI

Studies into cortical thickness in psychiatric diseases based on T1-weighted MRI frequently report on aberrations in the cerebral cortex. Due to limitations in image resolution for studies conducted at conventional MRI field strengths (e.g. 3 Tesla (T)) this information cannot be used to establish which of the cortical layers may be implicated. Here we propose a new analysis method that computes one high-resolution average cortical profile per brain region extracting myeloarchitectural information from T1-weighted MRI scans that are routinely acquired at a conventional field strength. To assess this new method, we acquired standard T1-weighted scans at 3 T and compared them with state-of-the-art ultra-high resolution T1-weighted scans optimised for intracortical myelin contrast acquired at 7 T. Average cortical profiles were computed for seven different brain regions. Besides a qualitative comparison between the 3 T scans, 7 T scans, and results from literature, we tested if the results from dynamic time warping-based clustering are similar for the cortical profiles computed from 7 T and 3 T data. In addition, we quantitatively compared cortical profiles computed for V1, V2 and V7 for both 7 T and 3 T data using a priori information on their relative myelin concentration. Although qualitative comparisons show that at an individual level average profiles computed for 7 T have more pronounced features than 3 T profiles the results from the quantitative analyses suggest that average cortical profiles computed from T1-weighted scans acquired at 3 T indeed contain myeloarchitectural information similar to profiles computed from the scans acquired at 7 T. The proposed method therefore provides a step forward to study cortical myeloarchitecture in vivo at conventional magnetic field strength both in health and disease.

Somatosensory maps

Somatosensory areas containing topographic maps of the body surface are a major feature of parietal cortex. In primates, parietal cortex contains four somatosensory areas, each with its own map, with the primary cutaneous map in area 3b. Rodents have at least three parietal somatosensory areas. Maps are not isomorphic to the body surface, but magnify behaviorally important skin regions, which include the hands and face in primates, and the whiskers in rodents. Within each map, intracortical circuits process tactile information, mediate spatial integration, and support active sensation. Maps may also contain fine-scale representations of touch submodalities, or direction of tactile motion. Functional representations are more overlapping than suggested by textbook depictions of map topography. The whisker map in rodent somatosensory cortex is a canonic system for studying cortical microcircuits, sensory coding, and map plasticity. Somatosensory maps are plastic throughout life in response to altered use or injury. This chapter reviews basic principles and recent findings in primate, human, and rodent somatosensory maps.

Using diffusion MRI to discriminate areas of cortical grey matter

Cortical area parcellation is a challenging problem that is often approached by combining structural imaging (e.g., quantitative T1, diffusion-based connectivity) with functional imaging (e.g., task activations, topological mapping, resting state correlations). Diffusion MRI (dMRI) has been widely adopted to analyse white matter microstructure, but scarcely used to distinguish grey matter regions because of the reduced anisotropy there. Nevertheless, differences in the texture of the cortical 'fabric' have long been mapped by histologists to distinguish cortical areas. Reliable area-specific contrast in the dMRI signal has previously been demonstrated in selected occipital and sensorimotor areas. We expand upon these findings by testing several diffusion-based feature sets in a series of classification tasks. Using Human Connectome Project (HCP) 3T datasets and a supervised learning approach, we demonstrate that diffusion MRI is sensitive to architectonic differences between a large number of different cortical areas defined in the HCP parcellation. By employing a surface-based cortical imaging pipeline, which defines diffusion features relative to local cortical surface orientation, we show that we can differentiate areas from their neighbours with higher accuracy than when using only fractional anisotropy or mean diffusivity. The results suggest that grey matter diffusion may provide a new, independent source of information for dividing up the cortex.

•

Diffusion MRI provides architectonic contrast between a variety of cortical areas.

•

dMRI provides complementary information in areas of low myelin content.

•

High order rotational invariants provide more reliable parcellation results.

Measuring the effects of attention to individual fingertips in somatosensory cortex using ultra-high field (7T) fMRI

Attention to sensory information has been shown to modulate the neuronal processing of that information. For example, visuospatial attention acts by modulating responses at retinotopically appropriate regions of visual cortex (Puckett and DeYoe, 2015; Tootell et al. 1998). Much less, however, is known about the neuronal processing associated with attending to other modalities of sensory information. One reason for this is that visual cortex is relatively large, and therefore easier to access non-invasively in humans using tools such as functional magnetic resonance imaging (fMRI). With high-resolution fMRI, however, it is now possible to access smaller cortical areas such as primary somatosensory cortex (Martuzzi et al., 2014; Sanchez-Panchuelo et al., 2010; Schweisfurth et al. 2014; Schweizer et al. 2008). Here, we combined a novel experimental design and high-resolution fMRI at ultra-high field (7T) to measure the effects of attention to tactile stimulation in primary somatosensory cortex, S1. We find that attention modulates somatotopically appropriate regions of S1, and importantly, that this modulation can be measured at the level of the cortical representation of individual fingertips.

The impact of ultra-high field MRI on cognitive and computational neuroimaging

The ability to measure functional brain responses non-invasively with ultra high field MRI (7 T and above) represents a unique opportunity in advancing our understanding of the human brain. Compared to lower fields (3 T and below), ultra high field MRI has an increased sensitivity, which can be used to acquire functional images with greater spatial resolution, and greater specificity of the blood oxygen level dependent (BOLD) signal to the underlying neuronal responses. Together, increased resolution and specificity enable investigating brain functions at a submillimeter scale, which so far could only be done with invasive techniques. At this mesoscopic spatial scale, perception, cognition and behavior can be probed at the level of fundamental units of neural computations, such as cortical columns, cortical layers, and subcortical nuclei. This represents a unique and distinctive advantage that differentiates ultra high from lower field imaging and that can foster a tighter link between fMRI and computational modeling of neural networks. So far, functional brain mapping at submillimeter scale has focused on the processing of sensory information and on well-known systems for which extensive information is available from invasive recordings in animals. It remains an open challenge to extend this methodology to uniquely human functions and, more generally, to systems for which animal models may be problematic. To succeed, the possibility to acquire high-resolution functional data with large spatial coverage, the availability of computational models of neural processing as well as accurate biophysical modeling of neurovascular coupling at mesoscopic scale all appear necessary.

Exploring structure and function of sensory cortex with 7T MRI

In this paper, we present an overview of 7T magnetic resonance imaging (MRI) studies of the detailed function and anatomy of sensory areas of the human brain. We discuss the motivation for the studies, with particular emphasis on increasing the spatial resolution of functional MRI (fMRI) using reduced field-of-view (FOV) data acquisitions. MRI at ultra-high-field (UHF) - defined here as 7T and above - has several advantages over lower field strengths. The intrinsic signal-to-noise ratio (SNR) of images is higher at UHF, and coupled with the increased blood-oxygen-level-dependent (BOLD) signal change, this results in increased BOLD contrast-to-noise ratio (CNR), which can be exploited to improve spatial resolution or detect weaker signals. Additionally, the BOLD signal from the intra-vascular (IV) compartment is relatively diminished compared to lower field strengths. Together, these properties make 7T functional MRI an attractive proposition for high spatial specificity measures. But with the advantages come some challenges. For example, increased vulnerability to susceptibility-induced geometric distortions and signal loss in EPI acquisitions tend to be much larger. Some of these technical issues can be addressed with currently available tools and will be discussed. We highlight the key methodological considerations for high resolution functional and structural imaging at 7 T. We then present recent data using the high spatial resolution available at UHF in studies of the visual and somatosensory cortex to highlight promising developments in this area.

Intracortical microstimulation of human somatosensory cortex

Intracortical microstimulation of the somatosensory cortex offers the potential for creating a sensory neuroprosthesis to restore tactile sensation. Whereas animal studies have suggested that both cutaneous and proprioceptive percepts can be evoked using this approach, the perceptual quality of the stimuli cannot be measured in these experiments. We show that microstimulation within the hand area of the somatosensory cortex of a person with long-term spinal cord injury evokes tactile sensations perceived as originating from locations on the hand and that cortical stimulation sites are organized according to expected somatotopic principles. Many of these percepts exhibit naturalistic characteristics (including feelings of pressure), can be evoked at low stimulation amplitudes, and remain stable for months. Further, modulating the stimulus amplitude grades the perceptual intensity of the stimuli, suggesting that intracortical microstimulation could be used to convey information about the contact location and pressure necessary to perform dexterous hand movements associated with object manipulation.

Recent applications of UHF-MRI in the study of human brain function and structure: a review

The increased availability of ultra-high-field (UHF) MRI has led to its application in a wide range of neuroimaging studies, which are showing promise in transforming fundamental approaches to human neuroscience. This review presents recent work on structural and functional brain imaging, at 7 T and higher field strengths. After a short outline of the effects of high field strength on MR images, the rapidly expanding literature on UHF applications of blood-oxygenation-level-dependent-based functional MRI is reviewed. Structural imaging is then discussed, divided into sections on imaging weighted by relaxation time, including quantitative relaxation time mapping, phase imaging and quantitative susceptibility mapping, angiography, diffusion-weighted imaging, and finally magnetization-transfer imaging. The final section discusses studies using the high spatial resolution available at UHF to identify explicit links between structure and function. Copyright © 2015 John Wiley & Sons, Ltd.

Inter- and Intradigit Somatotopic Map of High-Frequency Vibration Stimulations in Human Primary Somatosensory Cortex

Although more about the somatotopic mapping of fingers continues to be uncovered, there is lack of mapping attempts regarding the integration of within-finger and across-finger somatotopic coordinates in Broadmann area (BA) 3. This study aimed to address the issue by finding an inter-/intradigit somatotopic map with high-frequency (250 Hz) vibrotactile stimulation. Functional magnetic resonance imaging (fMRI) data were acquired while stimulation was applied to 3 phalanxes (distal [p1], intermediate [p2], and proximal [p3] phalanx) of 4 fingers (index, middle, ring, and little finger) for a total of 12 finger-phalanx combinations for a human. Inter-, intra-, and inter-/intradigit distances were calculated from peak activation coordinates in BA 3 for each combination. With regard to interdigit dimensions, the somatotopic coordinates proceeded in the lateral-to-medial direction for the index, middle, ring, and little fingers consecutively. This trend is comparable to that generated from low-frequency stimulation modalities (flutter stimulation). The somatotopic distances between fingers were greatest when p1 was compared across fingers. From an intradigit perspective, stimulation on p1, p2, and p3 yielded BA 3 peak coordinates aligned along the anterior-to-posterior and inferior-to-superior directions for all fingers. An inter-/intradigit map exhibited a radially propagating trend of distances calculated with respect to index p1 as a reference point; this provided an integrated view of inter- and intradigit somatotopies, which are traditionally discussed separately. We expect such an inter-/intradigit somatotopic map approach to contribute in generating a comprehensive somatotopic model of fingers.

Lines of Baillarger in vivo and ex vivo: Myelin contrast across lamina at 7T MRI and histology

The human cerebral cortex is characterized by a number of features that are not uniformly distributed, such as the presence of multiple cytoarchitectonic elements and of myelinated layers running tangentially to the cortex surface. The presence and absence of these features are the basis of the parcellation of the cerebral cortex in several areas. A number of areas show myelin increases localized within the cortex, e.g., the stria of Gennari located in layer IV of the primary visual cortex. Sub-millimeter MRI can resolve myelin variations across the human cortex and may allow in vivo parcellation of these brain areas. Here, we image within-area myelination. We modified a T1-weighted (T1-w) MPRAGE sequence to enhance myelin visualization within the cortex. First, we acquired images from an ex vivo sample, and compared MRI laminar profiles from calcarine (corresponding to primary visual cortex) and extra-calcarine areas with histology sections from the same locations. Laminar profiles between myelin stained sections and the T1-w images were similar both in calcarine as well as extra-calcarine cortex. In calcarine cortex, the profile reveals the stria of Gennari. In extra-calcarine cortex, a similar profile exists which we suggest corresponds to the lines of Baillarger. Next, we adapted the same sequence to image within-area myelination in vivo. Also in in vivo data, we discriminated similar laminar profiles in calcarine and extra-calcarine cortex, extending into parietal and frontal lobes. We argue that this myelin pattern outside the calcarine cortex represents the lines of Baillarger.

Stroking or Buzzing? A Comparison of Somatosensory Touch Stimuli Using 7 Tesla fMRI

Studying body representations in the brain helps us to understand how we humans relate to our own bodies. The in vivo mapping of the somatosensory cortex, where these representations are found, is greatly facilitated by the high spatial resolution and high sensitivity to brain activation available at ultra-high field. In this study, the use of different stimulus types for somatotopic mapping of the digits at ultra-high field, specifically manual stroking and mechanical stimulation, was compared in terms of sensitivity and specificity of the brain responses. Larger positive responses in digit regions of interest were found for manual stroking than for mechanical stimulation, both in terms of average and maximum t-value and in terms of number of voxels with significant responses to the tactile stimulation. Responses to manual stroking were higher throughout the entire post-central sulcus, but the difference was especially large on its posterior wall, i.e. in Brodmann area 2. During mechanical stimulation, cross-digit responses were more negative than during manual stroking, possibly caused by a faster habituation to the stimulus. These differences indicate that manual stroking is a highly suitable stimulus for fast somatotopic mapping procedures, especially if Brodmann area 2 is of interest.

Accurate cortical tissue classification on MRI by modeling cortical folding patterns

Accurate tissue classification is a crucial prerequisite to MRI morphometry. Automated methods based on intensity histograms constructed from the entire volume are challenged by regional intensity variations due to local radiofrequency artifacts as well as disparities in tissue composition, laminar architecture and folding patterns. Current work proposes a novel anatomy-driven method in which parcels conforming cortical folding were regionally extracted from the brain. Each parcel is subsequently classified using nonparametric mean shift clustering. Evaluation was carried out on manually labeled images from two datasets acquired at 3.0 Tesla (n = 15) and 1.5 Tesla (n = 20). In both datasets, we observed high tissue classification accuracy of the proposed method (Dice index >97.6% at 3.0 Tesla, and >89.2% at 1.5 Tesla). Moreover, our method consistently outperformed state-of-the-art classification routines available in SPM8 and FSL-FAST, as well as a recently proposed local classifier that partitions the brain into cubes. Contour-based analyses localized more accurate white matter-gray matter (GM) interface classification of the proposed framework compared to the other algorithms, particularly in central and occipital cortices that generally display bright GM due to their highly degree of myelination. Excellent accuracy was maintained, even in the absence of correction for intensity inhomogeneity. The presented anatomy-driven local classification algorithm may significantly improve cortical boundary definition, with possible benefits for morphometric inference and biomarker discovery.

Parcellation of the primary cerebral cortices based on local connectivity profiles

Connectivity-based parcellation using diffusion MRI has been extensively used to parcellate subcortical areas and the association cortex. Connectivity profiles are vital for connectivity-based parcellation. Two categories of connectivity profiles are generally utilized, including global connectivity profiles, in which the connectivity information is from the seed to the whole brain, and long connectivity profiles, in which the connectivity information is from the seed to other brain regions after excluding the seed. However, whether global or long connectivity profiles should be applied in parcellating the primary cortex utilizing connectivity-based parcellation is unclear. Many sources of evidence have indicated that the primary cerebral cortices are composed of structurally and functionally distinct subregions. Because the primary cerebral cortices are rich in local anatomic hierarchical connections and possess high degree of local functional connectivity profiles, we proposed that local connectivity profiles, that is the connectivity information within a seed region of interest, might be used for parcellating the primary cerebral cortices. In this study, the global, long, and local connectivity profiles were separately used to parcellate the bilateral M1, A1, S1, and V1. We found that results using the three profiles were all quite consistent with reported cytoarchitectonic evidence. More importantly, the results using local connectivity profiles showed less inter-subject variability than the results using the other two, a finding which suggests that local connectivity profiles are superior to global and long connectivity profiles for parcellating the primary cerebral cortices. This also implies that, depending on the characteristics of specific areas of the cerebral cortex, different connectivity profiles may need to be adopted to parcellate different areas.

New whole-body sensory-motor gradients revealed using phase-locked analysis and verified using multivoxel pattern analysis and functional connectivity

Topographic organization is one of the main principles of organization in the human brain. Specifically, whole-brain topographic mapping using spectral analysis is responsible for one of the greatest advances in vision research. Thus, it is intriguing that although topography is a key feature also in the motor system, whole-body somatosensory-motor mapping using spectral analysis has not been conducted in humans outside M1/SMA. Here, using this method, we were able to map a homunculus in the globus pallidus, a key target area for deep brain stimulation, which has not been mapped noninvasively or in healthy subjects. The analysis clarifies contradictory and partial results regarding somatotopy in the caudal-cingulate zone and rostral-cingulate zone in the medial wall and in the putamen. Most of the results were confirmed at the single-subject level and were found to be compatible with results from animal studies. Using multivoxel pattern analysis, we could predict movements of individual body parts in these homunculi, thus confirming that they contain somatotopic information. Using functional connectivity, we demonstrate interhemispheric functional somatotopic connectivity of these homunculi, such that the somatotopy in one hemisphere could have been found given the connectivity pattern of the corresponding regions of interest in the other hemisphere. When inspecting the somatotopic and nonsomatotopic connectivity patterns, a similarity index indicated that the pattern of connected and nonconnected regions of interest across different homunculi is similar for different body parts and hemispheres. The results show that topographical gradients are even more widespread than previously assumed in the somatosensory-motor system. Spectral analysis can thus potentially serve as a gold standard for defining somatosensory-motor system areas for basic research and clinical applications.

2D phase-sensitive inversion recovery imaging to measure in vivo spinal cord gray and white matter areas in clinically feasible acquisition times

PURPOSE

To present and assess a procedure for measurement of spinal cord total cross-sectional areas (TCA) and gray matter (GM) areas based on phase-sensitive inversion recovery imaging (PSIR). In vivo assessment of spinal cord GM and white matter (WM) could become pivotal to study various neurological diseases, but it is challenging because of insufficient GM/WM contrast provided by conventional magnetic resonance imaging (MRI).

MATERIALS AND METHODS

We acquired 2D PSIR images at 3T at each disc level of the spinal axis in 10 healthy subjects and measured TCA, cord diameters, WM and GM areas, and GM area/TCA ratios. Second, we investigated 32 healthy subjects at four selected levels (C2-C3, C3-C4, T8-T9, T9-T10, total acquisition time <8 min) and generated normative reference values of TCA and GM areas. We assessed test-retest, intra- and interoperator reliability of the acquisition strategy, and measurement steps.

RESULTS

The measurement procedure based on 2D PSIR imaging allowed TCA and GM area assessments along the entire spinal cord axis. The tests we performed revealed high test-retest/intraoperator reliability (mean coefficient of variation [COV] at C2-C3: TCA = 0.41%, GM area = 2.75%) and interoperator reliability of the measurements (mean COV on the 4 levels: TCA = 0.44%, GM area = 4.20%; mean intraclass correlation coefficient: TCA = 0.998, GM area = 0.906).

CONCLUSION

2D PSIR allows reliable in vivo assessment of spinal cord TCA, GM, and WM areas in clinically feasible acquisition times. The area measurements presented here are in agreement with previous MRI and postmortem studies.

Individual fMRI maps of all phalanges and digit bases of all fingers in human primary somatosensory cortex

This study determined the individual maps of all fingers in Brodmann area 3b of the human primary somatosensory cortex in a single fMRI session by tactile stimulation at 19 sites across all phalanges and digit bases of the 5 right-hand digits. To quantify basic features of the digit maps within and across subjects, we applied standard descriptive measures, but also implemented a novel quantitative analysis. This so-called Direction/Order (DiOr) method tested whether subjects exhibited an ordering of peak fMRI representations along their individual direction of alignment through the set of analyzed phalanges and whether these individual directions were similar across subjects. Across-digit analysis demonstrated that for each set of homologous phalanges, the D5-to-D1 representations were successively represented along a common direction of alignment. Hence, the well-known mediolateral D5-to-D1 somatotopy was not only confirmed for the distal phalanges (p1), but could also be shown for the medial (p2) and proximal phalanges (p3). In contrast, the peak activation for the digit bases (p4) only partly elicited that digit succession. Complementary, intra-digit analysis revealed a divergent picture of map topography for the different digits. Within D5 (and in a trend: D4), an ordered p1-to-p3 succession was found across subjects, pointing to a consistent intra-digit somatotopy for D5, with p3 generally found medial-posterior to p1. In contrast, for D1, D2, and D3, most subjects did not present with ordered p1-to-p3 maps nor were directions of alignment similarly oriented between subjects. These digits therefore exhibited highly diverse representation patterns across subjects.

Distinct fine-scale fMRI activation patterns of contra- and ipsilateral somatosensory areas 3b and 1 in humans

Inter-areal and ipsilateral cortical responses to tactile stimulation have not been well described in human S1 cortex. By taking advantage of the high signal-to-noise ratio at 7 T, we quantified blood oxygenation level dependent (BOLD) response patterns and time courses to tactile stimuli on individual distal finger pads at a fine spatial scale, and examined whether there are inter-areal (area 3b versus area 1) and interhemispheric response differences to unilateral tactile stimulation in healthy human subjects. We found that 2-Hz tactile stimulation of individual fingertips evoked detectable BOLD signal changes in both contralateral and ipsilateral area 3b and area 1. Contralateral digit activations were organized in an orderly somatotopic manner, and BOLD responses in area 3b were more digit selective than those in area 1. However, the area of cortex that was responsive to stimulation of a single digit (stimulus-response field) was similar across areas. In the ipsilateral hemisphere, response magnitudes in both areas 3b and 1 were significantly weaker than those of the contralateral hemisphere. Digit activations exhibited no clear somatotopic organizational pattern in either area 3b or area 1, yet digit selectivity was retained in area 1 but not in area 3b. The observation of distinct digit-selective responses of contralateral area 3b versus area 1 supports a higher order function of contralateral area 1 in spatial integration. In contrast, ipsilateral cortices may play a less discriminative role in the perception of unilateral tactile sensation in humans.

Event-related fMRI at 7T reveals overlapping cortical representations for adjacent fingertips in S1 of individual subjects

Recent fMRI studies of the human primary somatosensory cortex have been able to differentiate the cortical representations of different fingertips at a single‐subject level. These studies did not, however, investigate the expected overlap in cortical activation due to the stimulation of different fingers. Here, we used an event‐related design in six subjects at 7 Tesla to explore the overlap in cortical responses elicited in S1 by vibrotactile stimulation of the five fingertips. We found that all parts of S1 show some degree of spatial overlap between the cortical representations of adjacent or even nonadjacent fingertips. In S1, the posterior bank of the central sulcus showed less overlap than regions in the post‐central gyrus, which responded to up to five fingertips. The functional properties of these two areas are consistent with the known layout of cytoarchitectonically defined subareas, and we speculate that they correspond to subarea 3b (S1 proper) and subarea 1, respectively. In contrast with previous fMRI studies, however, we did not observe discrete activation clusters that could unequivocally be attributed to different subareas of S1. Venous maps based on T2*‐weighted structural images suggest that the observed overlap is not driven by extra‐vascular contributions from large veins. Hum Brain Mapp 35:2027–2043, 2014. © 2013 The Authors Human Brain Mapping published by Wiley Periodicals, Inc.