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.

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.

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.

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.

Correspondence of human visual areas identified using functional and anatomical MRI in vivo at 7 T.

PURPOSE: To study the correspondence of anatomically and functionally defined visual areas (primary visual cortex, V1, and motion selective area V5/human MT+) by using structural magnetic resonance imaging (MRI) and functional MRI (fMRI) in vivo at 7 T. MATERIALS AND METHODS: Four subjects participated in this study. High-resolution (≈0.4 mm isotropic) anatomical MRI was used to identify cortical regions based on their distinct cortical lamination. The optimal contrast for identifying heavily myelinated layers within gray matter was quantitatively assessed by comparing T(1)-weighted magnetization-prepared rapid gradient echo (MPRAGE) and T(2)*-weighted, 3D fast-low angle shot (FLASH) imaging. Retinotopic mapping was performed using GE-based fMRI at 1.5 mm isotropic resolution to identify functional areas. RESULTS: T(2)*-weighted FLASH imaging was found to provide a significantly higher contrast-to-noise ratio, allowing visualization of the stria of Gennari in every slice of a volume covering the occipital cortex in each of the four subjects in this study. The independently derived boundary of V1, identified in the same subjects using retinotopic mapping by fMRI, closely matched the border of anatomically defined striate cortex in the human brain. Evidence of banding was also found within the functionally defined V5 area; however, we did not find a good correlation of this area, or the functionally identified subregion (MT), with the banded area. CONCLUSION: High-resolution T(2)*-weighted images acquired at 7 T can be used to identify myelinated bands within cortical gray matter in reasonable measurement times. Regions where a myelinated band was identified show a high degree of overlap with the functionally defined V1 area.