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.

Fast Event-Related Mapping of Population Fingertip Tuning Properties in Human Sensorimotor Cortex at 7T.

fMRI studies that investigate somatotopic tactile representations in the human cortex typically use either block or phase-encoded stimulation designs. Event-related (ER) designs allow for more flexible and unpredictable stimulation sequences than the other methods, but they are less efficient. Here, we compared an efficiency-optimized fast ER design (2.8-s average intertrial interval; ITI) to a conventional slow ER design (8-s average ITI) for mapping voxelwise fingertip tactile tuning properties in the sensorimotor cortex of six participants at 7 Tesla. The fast ER design yielded more reliable responses compared with the slow ER design, but with otherwise similar tuning properties. Concatenating the fast and slow ER data, we demonstrate in each individual brain the existence of two separate somatotopically-organized tactile representations of the fingertips, one in the primary somatosensory cortex (S1) on the postcentral gyrus, and the other shared across the motor and premotor cortices on the precentral gyrus. In both S1 and motor representations, fingertip selectivity decreased progressively, from narrowly-tuned Brodmann area (BA) 3b and BA4a, respectively, toward associative parietal and frontal regions that responded equally to all fingertips, suggesting increasing information integration along these two pathways. In addition, fingertip selectivity in S1 decreased from the cortical representation of the thumb to that of the pinky.

Quantitative T1 mapping using multi-slice multi-shot inversion recovery EPI.

An efficient multi-slice inversion-recovery EPI (MS-IR-EPI) sequence for fast, high spatial resolution, quantitative T1 mapping is presented, using a segmented simultaneous multi-slice acquisition, combined with slice order shifting across multiple acquisitions. The segmented acquisition minimises the effective TE and readout duration compared to a single-shot EPI scheme, reducing geometric distortions to provide high quality T1 maps with a narrow point-spread function. The precision and repeatability of MS-IR-EPI T1 measurements are assessed using both T1-calibrated and T2-calibrated ISMRM/NIST phantom spheres at 3 and 7 T and compared with single slice IR and MP2RAGE methods. Magnetization transfer (MT) effects of the spectrally-selective fat-suppression (FS) pulses required for in vivo imaging are shown to shorten the measured in-vivo T1 values. We model the effect of these fat suppression pulses on T1 measurements and show that the model can remove their MT contribution from the measured T1, thus providing accurate T1 quantification. High spatial resolution T1 maps of the human brain generated with MS-IR-EPI at 7 T are compared with those generated with the widely implemented MP2RAGE sequence. Our MS-IR-EPI sequence provides high SNR per unit time and sharper T1 maps than MP2RAGE, demonstrating the potential for ultra-high resolution T1 mapping and the improved discrimination of functionally relevant cortical areas in the human brain.

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.

Imaging human cortical responses to intraneural microstimulation using magnetoencephalography.

The sensation of touch in the glabrous skin of the human hand is conveyed by thousands of fast-conducting mechanoreceptive afferents, which can be categorised into four distinct types. The spiking properties of these afferents in the periphery in response to varied tactile stimuli are well-characterised, but relatively little is known about the spatiotemporal properties of the neural representations of these different receptor types in the human cortex. Here, we use the novel methodological combination of single-unit intraneural microstimulation (INMS) with magnetoencephalography (MEG) to localise cortical representations of individual touch afferents in humans, by measuring the extracranial magnetic fields from neural currents. We found that by assessing the modulation of the beta (13-30 Hz) rhythm during single-unit INMS, significant changes in oscillatory amplitude occur in the contralateral primary somatosensory cortex within and across a group of fast adapting type I mechanoreceptive afferents, which corresponded well to the induced response from matched vibrotactile stimulation. Combining the spatiotemporal specificity of MEG with the selective single-unit stimulation of INMS enables the interrogation of the central representations of different aspects of tactile afferent signalling within the human cortices. The fundamental finding that single-unit INMS ERD responses are robust and consistent with natural somatosensory stimuli will permit us to more dynamically probe the central nervous system responses in humans, to address questions about the processing of touch from the different classes of mechanoreceptive afferents and the effects of varying the stimulus frequency and patterning.

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.

Is Human Auditory Cortex Organization Compatible With the Monkey Model? Contrary Evidence From Ultra-High-Field Functional and Structural MRI.

It is commonly assumed that the human auditory cortex is organized similarly to that of macaque monkeys, where the primary region, or "core," is elongated parallel to the tonotopic axis (main direction of tonotopic gradients), and subdivided across this axis into up to 3 distinct areas (A1, R, and RT), with separate, mirror-symmetric tonotopic gradients. This assumption, however, has not been tested until now. Here, we used high-resolution ultra-high-field (7 T) magnetic resonance imaging (MRI) to delineate the human core and map tonotopy in 24 individual hemispheres. In each hemisphere, we assessed tonotopic gradients using principled, quantitative analysis methods, and delineated the core using 2 independent (functional and structural) MRI criteria. Our results indicate that, contrary to macaques, the human core is elongated perpendicular rather than parallel to the main tonotopic axis, and that this axis contains no more than 2 mirror-reversed gradients within the core region. Previously suggested homologies between these gradients and areas A1 and R in macaques were not supported. Our findings suggest fundamental differences in auditory cortex organization between humans and macaques.

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.

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.

An intra-neural microstimulation system for ultra-high field magnetic resonance imaging and magnetoencephalography.

BACKGROUND: Intra-neural microstimulation (INMS) is a technique that allows the precise delivery of low-current electrical pulses into human peripheral nerves. Single unit INMS can be used to stimulate individual afferent nerve fibres during microneurography. Combining this with neuroimaging allows the unique monitoring of central nervous system activation in response to unitary, controlled tactile input, with functional magnetic resonance imaging (fMRI) providing exquisite spatial localisation of brain activity and magnetoencephalography (MEG) high temporal resolution. NEW METHOD: INMS systems suitable for use within electrophysiology laboratories have been available for many years. We describe an INMS system specifically designed to provide compatibility with both ultra-high field (7T) fMRI and MEG. Numerous technical and safety issues are addressed. The system is fully analogue, allowing for arbitrary frequency and amplitude INMS stimulation. RESULTS: Unitary recordings obtained within both the MRI and MEG screened-room environments are comparable with those obtained in 'clean' electrophysiology recording environments. Single unit INMS (current <7µA, 200µs pulses) of individual mechanoreceptive afferents produces appropriate and robust responses during fMRI and MEG. COMPARISON WITH EXISTING METHOD(S): This custom-built MRI- and MEG-compatible stimulator overcomes issues with existing INMS approaches; it allows well-controlled switching between recording and stimulus mode, prevents electrical shocks because of long cable lengths, permits unlimited patterns of stimulation, and provides a system with improved work-flow and participant comfort. CONCLUSIONS: We demonstrate that the requirements for an INMS-integrated system, which can be used with both fMRI and MEG imaging systems, have been fully met.

Mapping quantal touch using 7 Tesla functional magnetic resonance imaging and single-unit intraneural microstimulation.

Using ultra-high field 7 Tesla (7T) functional magnetic resonance imaging (fMRI), we map the cortical and perceptual responses elicited by intraneural microstimulation (INMS) of single mechanoreceptive afferent units in the median nerve, in humans. Activations are compared to those produced by applying vibrotactile stimulation to the unit's receptive field, and unit-type perceptual reports are analyzed. We show that INMS and vibrotactile stimulation engage overlapping areas within the topographically appropriate digit representation in the primary somatosensory cortex. Additional brain regions in bilateral secondary somatosensory cortex, premotor cortex, primary motor cortex, insula and posterior parietal cortex, as well as in contralateral prefrontal cortex are also shown to be activated in response to INMS. The combination of INMS and 7T fMRI opens up an unprecedented opportunity to bridge the gap between first-order mechanoreceptive afferent input codes and their spatial, dynamic and perceptual representations in human cortex.

7 tesla FMRI reveals systematic functional organization for binocular disparity in dorsal visual cortex.

The binocular disparity between the views of the world registered by the left and right eyes provides a powerful signal about the depth structure of the environment. Despite increasing knowledge of the cortical areas that process disparity from animal models, comparatively little is known about the local architecture of stereoscopic processing in the human brain. Here, we take advantage of the high spatial specificity and image contrast offered by 7 tesla fMRI to test for systematic organization of disparity representations in the human brain. Participants viewed random dot stereogram stimuli depicting different depth positions while we recorded fMRI responses from dorsomedial visual cortex. We repeated measurements across three separate imaging sessions. Using a series of computational modeling approaches, we report three main advances in understanding disparity organization in the human brain. First, we show that disparity preferences are clustered and that this organization persists across imaging sessions, particularly in area V3A. Second, we observe differences between the local distribution of voxel responses in early and dorsomedial visual areas, suggesting different cortical organization. Third, using modeling of voxel responses, we show that higher dorsal areas (V3A, V3B/KO) have properties that are characteristic of human depth judgments: a simple model that uses tuning parameters estimated from fMRI data captures known variations in human psychophysical performance. Together, these findings indicate that human dorsal visual cortex contains selective cortical structures for disparity that may support the neural computations that underlie depth perception.

Assessing the spatial precision of SE and GE-BOLD contrast at 7 Tesla.

Spin echo (SE) EPI offers an alternative to standard gradient echo (GE) EPI for functional MRI. SE-EPI offers improved spatial specificity, since signal changes originate from the microvasculature, but its lower functional sensitivity has limited the usage of this sequence in fMRI experiments. Differential fMRI paradigms, in which two closely matched stimulus conditions are used, can suppress the contribution from veins, thus also offering improved spatial specificity compared to conventional block or event-related designs with long "rest" periods. In this study, we employed a differential fMRI paradigm to stimulate bands of primary visual cortex with pre-defined widths by using visual stimuli comprised of complementary rings of contrast-reversing checkerboard patterns (8 Hz). This paradigm was used to investigate the spatial specificity of GE and SE-BOLD contrast at 7T. Results show that the contrast-to-noise ratio (CNR) is larger for GE-EPI data than for the SE-EPI data for band widths in the range 1.7-6.6 mm, however as the width of the band decreases the CNR for GE and SE sequences converges. These results suggest that when using a differential mapping paradigm, GE-BOLD contrast is better for studying functional features that are larger than ~1.5 mm in size.

Neuroimaging paradigms for tonotopic mapping (II): the influence of acquisition protocol.

Numerous studies on the tonotopic organisation of auditory cortex in humans have employed a wide range of neuroimaging protocols to assess cortical frequency tuning. In the present functional magnetic resonance imaging (fMRI) study, we made a systematic comparison between acquisition protocols with variable levels of interference from acoustic scanner noise. Using sweep stimuli to evoke travelling waves of activation, we measured sound-evoked response signals using sparse, clustered, and continuous imaging protocols that were characterised by inter-scan intervals of 8.8, 2.2, or 0.0 s, respectively. With regard to sensitivity to sound-evoked activation, the sparse and clustered protocols performed similarly, and both detected more activation than the continuous method. Qualitatively, tonotopic maps in activated areas proved highly similar, in the sense that the overall pattern of tonotopic gradients was reproducible across all three protocols. However, quantitatively, we observed substantial reductions in response amplitudes to moderately low stimulus frequencies that coincided with regions of strong energy in the scanner noise spectrum for the clustered and continuous protocols compared to the sparse protocol. At the same time, extreme frequencies became over-represented for these two protocols, and high best frequencies became relatively more abundant. Our results indicate that although all three scanning protocols are suitable to determine the layout of tonotopic fields, an exact quantitative assessment of the representation of various sound frequencies is substantially confounded by the presence of scanner noise. In addition, we noticed anomalous signal dynamics in response to our travelling wave paradigm that suggest that the assessment of frequency-dependent tuning is non-trivially influenced by time-dependent (hemo)dynamics when using sweep stimuli.

Functional quantitative susceptibility mapping (fQSM).

Blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) is a powerful technique, typically based on the statistical analysis of the magnitude component of the complex time-series. Here, we additionally interrogated the phase data of the fMRI time-series and used quantitative susceptibility mapping (QSM) in order to investigate the potential of functional QSM (fQSM) relative to standard magnitude BOLD fMRI. High spatial resolution data (1mm isotropic) were acquired every 3 seconds using zoomed multi-slice gradient-echo EPI collected at 7 T in single orientation (SO) and multiple orientation (MO) experiments, the latter involving 4 repetitions with the subject's head rotated relative to B0. Statistical parametric maps (SPM) were reconstructed for magnitude, phase and QSM time-series and each was subjected to detailed analysis. Several fQSM pipelines were evaluated and compared based on the relative number of voxels that were coincidentally found to be significant in QSM and magnitude SPMs (common voxels). We found that sensitivity and spatial reliability of fQSM relative to the magnitude data depended strongly on the arbitrary significance threshold defining "activated" voxels in SPMs, and on the efficiency of spatio-temporal filtering of the phase time-series. Sensitivity and spatial reliability depended slightly on whether MO or SO fQSM was performed and on the QSM calculation approach used for SO data. Our results present the potential of fQSM as a quantitative method of mapping BOLD changes. We also critically discuss the technical challenges and issues linked to this intriguing new technique.

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.

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

Ultra-high-field (UHF) MRI is ideally suited for structural and functional imaging of the brain. High-resolution structural MRI can be used to map the anatomical boundaries between functional domains of the brain by identifying changes related to the pattern of myelination within cortical gray matter, opening up the possibility to study the relationship between functional domains and underlying structure in vivo. In a recent study, we demonstrated the correspondence between functional (based on retinotopic mapping) and structural (based on changes in T2(⁎)-weighted images linked to myelination) parcellations of the primary visual cortex (V1) in vivo at 7T (Sanchez-Panchuelo et al., 2012b). Here, we take advantage of the improved BOLD CNR and high spatial resolution achievable at 7T to study regional structural variations across the functionally defined areas within the primary somatosensory cortex (S1) in individual subjects. Using a traveling wave fMRI paradigm to map the internal somatotopic representation of the index, middle, and ring fingers in S1, we were able to identify multiple map reversals at the tip and base, corresponding to the boundaries between Brodmann areas 3a, 3b, 1 and 2. Based on high resolution structural MRI data acquired in the same subjects, we inspected these functionally-parcellated Brodmann areas for differences in cortical thickness and MR contrast measures (magnetization transfer ratio (MTR) and signal intensity in phase sensitive inversion recovery (PSIR) images) that are sensitive to myelination. Consistent area-related differences in cortical thickness and MTR/PSIR measurements were found across subjects. However these measures did not have sufficient sensitivity to allow definition of areal boundaries.

Single-subject fMRI mapping at 7 T of the representation of fingertips in S1: a comparison of event-related and phase-encoding designs.

A desirable goal of functional MRI (fMRI), both clinically and for basic research, is to produce detailed maps of cortical function in individual subjects. Single-subject mapping of the somatotopic hand representation in the human primary somatosensory cortex (S1) has been performed using both phase-encoding and block/event-related designs. Here, we review the theoretical strengths and limits of each method and empirically compare high-resolution (1.5 mm isotropic) somatotopic maps obtained using fMRI at ultrahigh magnetic field (7 T) with phase-encoding and event-related designs in six subjects in response to vibrotactile stimulation of the five fingertips. Results show that the phase-encoding design is more efficient than the event-related design for mapping fingertip-specific responses and in particular allows us to describe a new additional somatotopic representation of fingertips on the precentral gyrus. However, with sufficient data, both designs yield very similar fingertip-specific maps in S1, which confirms that the assumption of local representational continuity underlying phase-encoding designs is largely valid at the level of the fingertips in S1. In addition, it is shown that the event-related design allows the mapping of overlapping cortical representations that are difficult to estimate using the phase-encoding design. The event-related data show a complex pattern of overlapping cortical representations for different fingertips within S1 and demonstrate that regions of S1 responding to several adjacent fingertips can incorrectly be identified as responding preferentially to one fingertip in the phase-encoding data.

Within-digit functional parcellation of Brodmann areas of the human primary somatosensory cortex using functional magnetic resonance imaging at 7 tesla.

The primary somatosensory cortex (S1) can be subdivided cytoarchitectonically into four distinct Brodmann areas (3a, 3b, 1, and 2), but these areas have never been successfully delineated in vivo in single human subjects. Here, we demonstrate the functional parcellation of four areas of S1 in individual human subjects based on high-resolution functional MRI measurements made at 7 T using vibrotactile stimulation. By stimulating four sites along the length of the index finger, we were able to identify and locate map reversals of the base to tip representation of the index finger in S1. We suggest that these reversals correspond to the areal borders between the mirrored representations in the four Brodmann areas, as predicted from electrophysiology measurements in nonhuman primates. In all subjects, maps were highly reproducible across scanning sessions and stable over weeks. In four of the six subjects scanned, four, mirrored, within-finger somatotopic maps defining the extent of the Brodmann areas could be directly observed on the cortical surface. In addition, by using multivariate classification analysis, the location of stimulation on the index finger (four distinct sites) could be decoded with a mean accuracy of 65% across subjects. Our measurements thus show that within-finger topography is present at the millimeter scale in the cortex and is highly reproducible. The ability to identify functional areas of S1 in vivo in individual subjects will provide a framework for investigating more complex aspects of tactile representation in S1.