Search

Introduction: Functional magnetic resonance imaging (fMRI) can be used to assess language and memory function as part of pre-surgical decision making in refractory epilepsy. Although language paradigms are well established, memory paradigms are not widely used in clinical practice due to a lack of evidence for robust and reliable methods. Here, we aim to investigate the clinical utility of the Home Town Walk (HTW) paradigm for personalized treatment decisions in medial temporal lobe epilepsy., Methods: A cohort of 123 consecutive patients having HTW-fMRI as part of routine MRI scans over a 7.5 year period were included in this retrospective study. Of these, 111 patients underwent repeated HTW-fMRI in two scanning sessions one to three days apart. fMRI analysis was performed at the time of the scans using clinically approved software and retrospectively validated using FSL. We assessed the test-retest within subject reliability of activations within the posterior parahippocampal gyri (pPHG) at the individual subject level., Results and Discussion: Activations within the pPHG region were observed for 101 patients (91%) in at least one of the fMRI sessions and for 88 patients (79%) in both fMRI sessions, with 82 patients showing overlapping unilateral or bilateral activations and 8 further patients showing overlapping activations in one of the hemispheres but not the other. Reproducibility was evaluated using metrics based on the concordance ratios for size (R size ) and location (R overlap ) within the pPHG region, as well as the lateralization index (LI) metric to reflect the asymmetry of hemispheric activations, which is of crucial relevance to inform surgery. Test-retest reliability of visuospatial memory LIs, assessed by an intra-class correlation coefficient (ICC) yielded a value of 0.76, indicating excellent between session stability of memory lateralization., Conclusion: The HTW-fMRI paradigm shows reproducible activations in the medial temporal lobes of individual epilepsy patients sufficient to consistently lateralize visuospatial memory function, demonstrating the clinical utility of HTW memory fMRI and its potential for application in the pre-surgical assessment of people with temporal lobe epilepsy., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Sanchez Panchuelo, Flintham, Wesolowski, Jalali, Herbert, Samarasekera, Bagshaw, Chelvarajah, Davies and Sawlani.)

An efficient multi-slice inversion-recovery EPI (MS-IR-EPI) sequence for fast, high spatial resolution, quantitative T 1 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 T 1 maps with a narrow point-spread function. The precision and repeatability of MS-IR-EPI T 1 measurements are assessed using both T 1 -calibrated and T 2 -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 T 1 values. We model the effect of these fat suppression pulses on T 1 measurements and show that the model can remove their MT contribution from the measured T 1 , thus providing accurate T 1 quantification. High spatial resolution T 1 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 T 1 maps than MP2RAGE, demonstrating the potential for ultra-high resolution T 1 mapping and the improved discrimination of functionally relevant cortical areas in the human brain., Competing Interests: Declaration of Competing Interest None., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

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.

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., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© 2018 The Authors. Human Brain Mapping published by Wiley Periodicals, Inc.)

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.

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., (Copyright © 2017 Elsevier Inc. All rights reserved.)

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.

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.

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., (Copyright © 2014. Published by Elsevier Inc.)

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., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Previous studies have demonstrated that the perceived direction of motion of a visual stimulus can be decoded from the pattern of functional magnetic resonance imaging (fMRI) responses in occipital cortex using multivariate analysis methods (Kamitani and Tong, 2006). One possible mechanism for this is a difference in the sampling of direction selective cortical columns between voxels, implying that information at a level smaller than the voxel size might be accessible with fMRI. Alternatively, multivariate analysis methods might be driven by the organization of neurons into clusters or even orderly maps at a much larger scale. To assess the possible sources of the direction selectivity observed in fMRI data, we tested how classification accuracy varied across different visual areas and subsets of voxels for classification of motion-direction. To enable high spatial resolution functional MRI measurements (1.5mm isotropic voxels), data were collected at 7T. To test whether information about the direction of motion is represented at the scale of retinotopic maps, we looked at classification performance after combining data across different voxels within visual areas (V1-3 and MT+/V5) before training the multivariate classifier. A recent study has shown that orientation biases in V1 are both necessary and sufficient to explain classification of stimulus orientation (Freeman et al., 2011). Here, we combined voxels with similar visual field preference as determined in separate retinotopy measurements and observed that classification accuracy was preserved when averaging in this 'retinotopically restricted' way, compared to random averaging of voxels. This insensitivity to averaging of voxels (with similar visual angle preference) across substantial distances in cortical space suggests that there are large-scale biases at the level of retinotopic maps underlying our ability to classify direction of motion., (Copyright © 2012 Elsevier Inc. All rights reserved.)

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.

The increased blood oxygenation level-dependent contrast-to-noise ratio at ultrahigh field (7 T) has been exploited in a comparison of the spatial location and strength of activation in high-resolution (1.5  mm isotropic) gradient echo (GE) and spin echo (SE), echo planar imaging data acquired during the execution of a simple motor task in five subjects. SE data were acquired at six echo times from 30 to 55  ms. Excellent fat suppression was achieved in the SE echo planar images using slice-selective gradient reversal. Threshold-free cluster enhancement was used to define regions of interest (ROIs) containing voxels showing significant stimulus-locked signal changes from the GE and average SE data. These were used to compare the signal changes and spatial locations of activated regions in SE and GE data. T(2) and T(2)* values were measured, with means of 48.3 ± 1.1  ms and 36.5 ± 3.4  ms in the SE ROI. In addition, we identified a dark band in SE images of the motor cortex corresponding to a region in which T(2) and T(2)* were significantly lower than in the surrounding grey matter. The fractional SE signal change in the ROI was found to vary linearly as a function of TE, with a slope that was dependent on the particular ROI assessed: the mean ΔR(2) value was found to be 0.85 ± 0.11  s(-1) for the SE ROI and -0.37 ± 0.05  s(-1) for the GE ROI. The fractional signal change relative to the shortest TE revealed that the largest signal change occurred at a TE of 45  ms outside of the dark band. At this TE, the ratio of the fractional signal change in GE and SE data was found to be 0.48 ± 0.05. Phase maps produced from high-resolution GE images spanning the right motor cortex were used to identify veins. The GE ROI was found to contain 18% more voxels overlying the venous mask than the SE ROI., (Copyright © 2011 John Wiley & Sons, Ltd.)

Functional magnetic resonance imaging (fMRI) is now routinely used to map the topographic organization of human visual cortex. Mapping the detailed topography of somatosensory cortex, however, has proven to be more difficult. Here we used the increased blood-oxygen-level-dependent contrast-to-noise ratio at ultra-high field (7 Tesla) to measure the topographic representation of the digits in human somatosensory cortex at 1 mm isotropic resolution in individual subjects. A "traveling wave" paradigm was used to locate regions of cortex responding to periodic tactile stimulation of each distal phalangeal digit. Tactile stimulation was applied sequentially to each digit of the left hand from thumb to little finger (and in the reverse order). In all subjects, we found an orderly map of the digits on the posterior bank of the central sulcus (postcentral gyrus). Additionally, we measured event-related responses to brief stimuli for comparison with the topographic mapping data and related the fMRI responses to anatomical images obtained with an inversion-recovery sequence. Our results have important implications for the study of human somatosensory cortex and underscore the practical utility of ultra-high field functional imaging with 1 mm isotropic resolution for neuroscience experiments. First, topographic mapping of somatosensory cortex can be achieved in 20 min, allowing time for further experiments in the same session. Second, the maps are of sufficiently high resolution to resolve the representations of all five digits and third, the measurements are robust and can be made in an individual subject. These combined advantages will allow somatotopic fMRI to be used to measure the representation of digits in patients undergoing rehabilitation or plastic changes after peripheral nerve damage as well as tracking changes in normal subjects undergoing perceptual learning.