A critical assessment of data quality and venous effects in sub-millimeter fMRI

Ultra-high-resolution fMRI of Human Ventral Temporal Cortex Reveals Differential Representation of Categories and Domains

Human ventral temporal cortex (VTC) is critical for visual recognition. It is thought that this ability is supported by large-scale patterns of activity across VTC that contain information about visual categories. However, it is unknown how category representations in VTC are organized at the submillimeter scale and across cortical depths. To fill this gap in knowledge, we measured BOLD responses in medial and lateral VTC to images spanning 10 categories from five domains (written characters, bodies, faces, places, and objects) at an ultra-high spatial resolution of 0.8 mm using 7 Tesla fMRI in both male and female participants. Representations in lateral VTC were organized most strongly at the general level of domains (e.g., places), whereas medial VTC was also organized at the level of specific categories (e.g., corridors and houses within the domain of places). In both lateral and medial VTC, domain-level and category-level structure decreased with cortical depth, and downsampling our data to standard resolution (2.4 mm) did not reverse differences in representations between lateral and medial VTC. The functional diversity of representations across VTC partitions may allow downstream regions to read out information in a flexible manner according to task demands. These results bridge an important gap between electrophysiological recordings in single neurons at the micron scale in nonhuman primates and standard-resolution fMRI in humans by elucidating distributed responses at the submillimeter scale with ultra-high-resolution fMRI in humans.SIGNIFICANCE STATEMENT Visual recognition is a fundamental ability supported by human ventral temporal cortex (VTC). However, the nature of fine-scale, submillimeter distributed representations in VTC is unknown. Using ultra-high-resolution fMRI of human VTC, we found differential distributed visual representations across lateral and medial VTC. Domain representations (e.g., faces, bodies, places, characters) were most salient in lateral VTC, whereas category representations (e.g., corridors/houses within the domain of places) were equally salient in medial VTC. These results bridge an important gap between electrophysiological recordings in single neurons at a micron scale and fMRI measurements at a millimeter scale.

Layer-dependent activity in human prefrontal cortex during working memory

Working memory involves storing and/or manipulating previously encoded information over a short-term delay period, which is typically followed by a behavioral response based on the remembered information. Although working memory tasks often engage dorsolateral prefrontal cortex, few studies have investigated whether their subprocesses are localized to different cortical depths in this region, and none have done so in humans. Here we use high-resolution functional MRI to interrogate the layer specificity of neural activity during different periods of a delayed-response task in dorsolateral prefrontal cortex. We detect activity time courses that follow the hypothesized patterns: namely, superficial layers are preferentially active during the delay period, specifically in trials requiring manipulation (rather than mere maintenance) of information held in working memory, and deeper layers are preferentially active during the response. Results demonstrate that layer-specific functional MRI can be used in higher-order brain regions to noninvasively map cognitive processing in humans.

Evaluating the Utility of EPIK in a Finger Tapping fMRI Experiment using BOLD Detection and Effective Connectivity

EPI with Keyhole (EPIK) is a hybrid imaging technique that overcomes many of the performance disadvantages associated with EPI. Previously, EPIK was shown to provide a higher temporal resolution and fewer image distortions than EPI whilst maintaining comparable performance for the detection of BOLD-based signals. This work carefully examines the putative enhanced sensitivity of EPIK in a typical fMRI setting by using a robust fMRI paradigm - visually guided finger tapping - to demonstrate the advantages of EPIK for fMRI at 3 T. The data acquired were directly compared to the community standard fMRI protocol using single-shot EPI to ascertain a clear comparison. Each sequence was optimised to offer its highest possible spatial resolution for a given set of imaging conditions, i.e., EPIK and EPI achieved an in-planar resolution of 2.08 × 2.08 mm² with 32 slices and 3.13 × 3.13 mm² with 36 slices, respectively. EPIK demonstrated a number of clear improvements, such as superior spatial resolution with favourable robustness against susceptibility artefacts. Both imaging sequences revealed robust activation within primary motor, premotor and visual regions, although significantly higher BOLD amplitudes were detected using EPIK within the primary and supplementary motor areas. Dynamic causal modelling, in combination with Bayesian model selection, identified identical winning models for EPIK and EPI data. Coupling parameters reflecting task-related modulations and the connectivity of fixed connections were comparably robust for both sequences. However, fixed connections from the left motor cortex to the right visual cortex were estimated as being significantly more robust for EPIK data.

Intracortical smoothing of small-voxel fMRI data can provide increased detection power without spatial resolution losses compared to conventional large-voxel fMRI data

Continued improvement in MRI acquisition technology has made functional MRI (fMRI) with small isotropic voxel sizes down to 1 mm and below more commonly available. Although many conventional fMRI studies seek to investigate regional patterns of cortical activation for which conventional voxel sizes of 3 mm and larger provide sufficient spatial resolution, smaller voxels can help avoid contamination from adjacent white matter (WM) and cerebrospinal fluid (CSF), and thereby increase the specificity of fMRI to signal changes within the gray matter. Unfortunately, temporal signal-to-noise ratio (tSNR), a metric of fMRI sensitivity, is reduced in high-resolution acquisitions, which offsets the benefits of small voxels. Here we introduce a framework that combines small, isotropic fMRI voxels acquired at 7 T field strength with a novel anatomically-informed, surface mesh-navigated spatial smoothing that can provide both higher detection power and higher resolution than conventional voxel sizes. Our smoothing approach uses a family of intracortical surface meshes and allows for kernels of various shapes and sizes, including curved 3D kernels that adapt to and track the cortical folding pattern. Our goal is to restrict smoothing to the cortical gray matter ribbon and avoid noise contamination from CSF and signal dilution from WM via partial volume effects. We found that the intracortical kernel that maximizes tSNR does not maximize percent signal change (ΔS/S), and therefore the kernel configuration that optimizes detection power cannot be determined from tSNR considerations alone. However, several kernel configurations provided a favorable balance between boosting tSNR and ΔS/S, and allowed a 1.1-mm isotropic fMRI acquisition to have higher performance after smoothing (in terms of both detection power and spatial resolution) compared to an unsmoothed 3.0-mm isotropic fMRI acquisition. Overall, the results of this study support the strategy of acquiring voxels smaller than the cortical thickness, even for studies not requiring high spatial resolution, and smoothing them down within the cortical ribbon with a kernel of an appropriate shape to achieve the best performance-thus decoupling the choice of fMRI voxel size from the spatial resolution requirements of the particular study. The improvement of this new intracortical smoothing approach over conventional surface-based smoothing is expected to be modest for conventional resolutions, however the improvement is expected to increase with higher resolutions. This framework can also be applied to anatomically-informed intracortical smoothing of higher-resolution data (e.g. along columns and layers) in studies with prior information about the spatial structure of activation.

The morphology of the human cerebrovascular system

While several methodologies exist for quantifying gray and white matter properties in humans, relatively little is known regarding the spatial organization and the intersubject variability of cerebral vessels. To resolve this, we developed a fast, open-source processing algorithm using advanced vessel segmentation schemes and iterative nonlinear registration to isolate, extract, and quantify cerebral vessels in susceptibility weighting imaging (SWI) and time-of-flight angiography (TOF-MRA) datasets acquired in a large cohort (n = 42) of healthy individuals. From this, whole-brain venous and arterial probabilistic maps were generated along with the computation of regional densities and diameters within regions based on popular anatomical and functional atlases. The results show that cerebral vasculature is highly heterogeneous, displaying disproportionally large vessel densities in brain areas such as the anterior and posterior cingulate, cuneus, precuneus, parahippocampus, insula, and temporal gyri. On average, venous densities were slightly higher and less variable across subjects than arterial. Moreover, regional variations in both venous and arterial density were significantly correlated to cortical thickness (R = 0.42). This publicly available new atlas of the human cerebrovascular system provides a first step toward quantifying morphological changes in the diseased brain and serving as a potential regression tool in fMRI analysis.