Advances in spiral fMRI: A high-resolution dataset

We present data collected for the research article "Advances in Spiral fMRI: A High-resolution Study with Single-shot Acquisition" (Kasper et al. 2022). All data was acquired on a 7T ultra-high field MR system (Philips Achieva), equipped with a concurrent magnetic field monitoring setup based on 16 NMR probes. For task-based fMRI, a visual quarterfield stimulation paradigm was employed, alongside physiological monitoring via peripheral recordings. This data collection contains different datasets pertaining to different purposes: (1) Measured magnetic field dynamics (k0, spiral k-space trajectories, 2nd order spherical harmonics, concomitant fields) during ultra-high field fMRI sessions from six subjects, as well as concurrent temperature curves of the gradient coil, to explore MR system and subject-induced variability of field fluctuations and assess the impact of potential correction methods. (2) MR Raw Data, i.e., coil and concurrent encoding magnetic field (trajectory) data, of a single subject, as well as nominal spiral gradient waveforms, precomputed B0 and coil sensitivity maps, to enable testing of alternative image reconstruction approaches for spiral fMRI data. (3) Reconstructed image time series of the same subject alongside behavioral and physiological logfiles, to reproduce the fMRI preprocessing and analysis, as well as figures presented in the research article related to this article, using the published analysis code repository. All data is provided in standardized formats for the respective research area. In particular, ISMRMRD (HDF5) is used to store raw coil data and spiral trajectories, as well as measured field dynamics. Likewise, the NIfTI format is used for all imaging data (anatomical MRI and spiral fMRI, B0 and coil sensitivity maps).

Predicting the retinotopic organization of human visual cortex from anatomy using geometric deep learning

Whether it be in a single neuron or a more complex biological system like the human brain, form and function are often directly related. The functional organization of human visual cortex, for instance, is tightly coupled with the underlying anatomy with cortical shape having been shown to be a useful predictor of the retinotopic organization in early visual cortex. Although the current state-of-the-art in predicting retinotopic maps is able to account for gross individual differences, such models are unable to account for any idiosyncratic differences in the structure-function relationship from anatomical information alone due to their initial assumption of a template. Here we developed a geometric deep learning model capable of exploiting the actual structure of the cortex to learn the complex relationship between brain function and anatomy in human visual cortex such that more realistic and idiosyncratic maps could be predicted. We show that our neural network was not only able to predict the functional organization throughout the visual cortical hierarchy, but that it was also able to predict nuanced variations across individuals. Although we demonstrate its utility for modeling the relationship between structure and function in human visual cortex, our approach is flexible and well-suited for a range of other applications involving data structured in non-Euclidean spaces.

Feasibility of spiral fMRI based on an LTI gradient model

Spiral imaging is very well suited for functional MRI, however its use has been limited by the fact that artifacts caused by gradient imperfections and B₀ inhomogeneity are more difficult to correct compared to EPI. Effective correction requires accurate knowledge of the traversed k-space trajectory. With the goal of making spiral fMRI more accessible, we have evaluated image reconstruction using trajectories predicted by the gradient impulse response function (GIRF), which can be determined in a one-time calibration step. GIRF-predicted reconstruction was tested for high-resolution (0.8 mm) fMRI at 7T. Image quality and functional results of the reconstructions using GIRF-prediction were compared to reconstructions using the nominal trajectory and concurrent field monitoring. The reconstructions using nominal spiral trajectories contain substantial artifacts and the activation maps contain misplaced activation. Image artifacts are substantially reduced when using the GIRF-predicted reconstruction, and the activation maps for the GIRF-predicted and monitored reconstructions largely overlap. The GIRF reconstruction provides a large increase in the spatial specificity of the activation compared to the nominal reconstruction. The GIRF-reconstruction generates image quality and fMRI results similar to using a concurrently monitored trajectory. The presented approach does not prolong or complicate the fMRI acquisition. Using GIRF-predicted trajectories has the potential to enable high-quality spiral fMRI in situations where concurrent trajectory monitoring is not available.

Sensitivity limitations of high-resolution perfusion-based human fMRI at 7 Tesla

The study of the brain's functional organization at laminar and columnar level of the cortex with blood oxygenation-level dependent (BOLD) functional MRI (fMRI) is affected by the contribution of large veins downstream from the microvascular response to brain activity. Blood volume- and especially perfusion-based techniques may reduce this problem because of their reduced sensitivity to venous effects, but may not allow the same spatial resolution because of smaller signal changes associated with brain activity. Here we investigated the practical resolution limits of perfusion-weighted fMRI in human visual stimulation experiments. For this purpose, we used a highly sensitive, single-shot perfusion labeling (SSPL) technique at 7 T and compared sensitivity to detect visual activation at low (2 mm, n = 10) and high (1 mm, n = 8) nominal isotropic spatial, and 3 s temporal, resolution with BOLD in 5½-minute-long experiments. Despite the smaller absolute signal change with activation, 2 mm resolution SSPL yielded comparable sensitivity to BOLD. This was attributed to a superior suppression of physiological noise with SSPL. However, at 1 mm nominal resolution, SSPL sensitivity fell on average at least 42% below that of BOLD, and detection of visual activation was compromised. This is explained by the fact that at high resolution, with both techniques, typically thermal noise rather than physiological noise dominates sensitivity. The observed sensitivity loss implies that to perform 1-mm resolution, perfusion weighted fMRI with a robustness similar to BOLD, scan times that are almost 3 times longer than the comparable BOLD experiment are required. This is in line with or slightly better than previous comparisons between perfusion-weighted fMRI and BOLD. The lower sensitivity has to be weighed against the spatial fidelity advantages of high-resolution perfusion-weighted fMRI.

Imaging faster neural dynamics with fast fMRI: A need for updated models of the hemodynamic response

Fast fMRI enables the detection of neural dynamics over timescales of hundreds of milliseconds, suggesting it may provide a new avenue for studying subsecond neural processes in the human brain. The magnitudes of these fast fMRI dynamics are far greater than predicted by canonical models of the hemodynamic response. Several studies have established nonlinear properties of the hemodynamic response that have significant implications for fast fMRI. We first review nonlinear properties of the hemodynamic response function that may underlie fast fMRI signals. We then illustrate the breakdown of canonical hemodynamic response models in the context of fast neural dynamics. We will then argue that the canonical hemodynamic response function is not likely to reflect the BOLD response to neuronal activity driven by sparse or naturalistic stimuli or perhaps to spontaneous neuronal fluctuations in the resting state. These properties suggest that fast fMRI is capable of tracking surprisingly fast neuronal dynamics, and we discuss the neuroscientific questions that could be addressed using this approach.

Increased dynamic flexibility in the medial temporal lobe network following an exercise intervention mediates generalization of prior learning

Recent work has conceptualized the brain as a network comprised of groups of sub-networks or modules. "Flexibility" of brain network(s) indexes the dynamic reconfiguration of comprising modules. Using novel techniques from dynamic network neuroscience applied to high-resolution resting-state functional magnetic resonance imaging (fMRI), the present study investigated the effects of an aerobic exercise intervention on the dynamic rearrangement of modular community structure-a measure of neural flexibility-within the medial temporal lobe (MTL) network. The MTL is one of the earliest brain regions impacted by Alzheimer's disease. It is also a major site of neuroplasticity that is sensitive to the effects of exercise. In a two-group non-randomized, repeated measures and matched control design with 34 healthy older adults, we observed an exercise-related increase in flexibility within the MTL network. Furthermore, MTL network flexibility mediated the beneficial effect aerobic exercise had on mnemonic flexibility, as measured by the ability to generalize past learning to novel task demands. Our results suggest that exercise exerts a rehabilitative and protective effect on MTL function, resulting in dynamically evolving networks of regions that interact in complex communication patterns. These reconfigurations may underlie exercise-induced improvements on cognitive measures of generalization, which are sensitive to subtle changes in the MTL.

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

Advances in hardware, pulse sequences, and reconstruction techniques have made it possible to perform functional magnetic resonance imaging (fMRI) at sub-millimeter resolution while maintaining high spatial coverage and acceptable signal-to-noise ratio. Here, we examine whether sub-millimeter fMRI can be used as a routine method for obtaining accurate measurements of fine-scale local neural activity. We conducted fMRI in human visual cortex during a simple event-related visual experiment (7 T, gradient-echo EPI, 0.8-mm isotropic voxels, 2.2-s sampling rate, 84 slices), and developed analysis and visualization tools to assess the quality of the data. Our results fall along three lines of inquiry. First, we find that the acquired fMRI images, combined with appropriate surface-based processing, provide reliable and accurate measurements of fine-scale blood oxygenation level dependent (BOLD) activity patterns. Second, we show that the highly folded structure of cortex causes substantial biases on spatial resolution and data visualization. Third, we examine the well-recognized issue of venous contributions to fMRI signals. In a systematic assessment of large sections of cortex measured at a fine scale, we show that time-averaged T₂*-weighted EPI intensity is a simple, robust marker of venous effects. These venous effects are unevenly distributed across cortex, are more pronounced in gyri and outer cortical depths, and are, to a certain degree, in consistent locations across subjects relative to cortical folding. Furthermore, we show that these venous effects are strongly correlated with BOLD responses evoked by the experiment. We conclude that sub-millimeter fMRI can provide robust information about fine-scale BOLD activity patterns, but special care must be exercised in visualizing and interpreting these patterns, especially with regards to the confounding influence of the brain's vasculature. To help translate these methodological findings to neuroscience research, we provide practical suggestions for both high-resolution and standard-resolution fMRI studies.

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.

Simulating laminar neuroimaging data for a visual delayed match-to-sample task

Invasive electrophysiological and neuroanatomical studies in nonhuman mammalian experimental preparations have helped elucidate the lamina (layer) dependence of neural computations and interregional connections. Noninvasive functional neuroimaging can, in principle, resolve cortical laminae (layers), and thus provide insight into human neural computations and interregional connections. However human neuroimaging data are noisy and difficult to interpret; biologically realistic simulations can aid experimental interpretation by relating the neuroimaging data to simulated neural activity. We illustrate the potential of laminar neuroimaging by upgrading an existing large-scale, multiregion neural model that simulates a visual delayed match-to-sample (DMS) task. The new laminar-based neural unit incorporates spiny stellate, pyramidal, and inhibitory neural populations which are divided among supragranular, granular, and infragranular laminae (layers). We simulated neural activity which is translated into local field potential-like data used to simulate conventional and laminar fMRI activity. We implemented the laminar connectivity schemes proposed by Felleman and Van Essen (Cerebral Cortex, 1991) for interregional connections. The hemodynamic model that we employ is a modified version of one due to Heinzle et al. (Neuroimage, 2016) that incorporates the effects of draining veins. We show that the laminar version of the model replicates the findings of the existing model. The laminar model shows the finer structure in fMRI activity and functional connectivity. Laminar differences in the magnitude of neural activities are a prominent finding; these are also visible in the simulated fMRI. We illustrate differences between task and control conditions in the fMRI signal, and demonstrate differences in interregional laminar functional connectivity that reflect the underlying connectivity scheme. These results indicate that multi-layer computational models can aid in interpreting layer-specific fMRI, and suggest that increased use of laminar fMRI could provide unique and fundamental insights to human neuroscience.

APOE ε4 status in healthy older African Americans is associated with deficits in pattern separation and hippocampal hyperactivation

African Americans are 1.4 times more likely than European Americans to carry the apolipoprotein E (APOE) ε4 allele, a risk factor for Alzheimer's disease (AD). However, little is known about the neural correlates of cognitive function in older African Americans and how they relate to genetic risk for AD. In particular, no past study on African Americans has examined the effect of APOE ε4 status on pattern separation-mnemonic discrimination performance and its corresponding neural computations in the hippocampus. Previous work using the mnemonic discrimination paradigm has localized increased activation in the DG/CA3 hippocampal subregions as being correlated with discrimination deficits. In a case-control high-resolution functional magnetic resonance imaging study of 30 healthy African Americans, aged 60 years and older, we observed APOE ε4-related impairments in mnemonic discrimination, coincident with dysfunctional hyperactivation in the DG/CA3, and CA1 regions, despite no evidence of structural differences in the hippocampus between carriers and noncarriers. Our results add to the growing body of evidence that deficits in pattern separation may be an early marker for AD-related neuronal dysfunction.

Analysis strategies for high-resolution UHF-fMRI data

Functional MRI (fMRI) benefits from both increased sensitivity and specificity with increasing magnetic field strength, making it a key application for Ultra-High Field (UHF) MRI scanners. Most UHF-fMRI studies utilize the dramatic increases in sensitivity and specificity to acquire high-resolution data reaching sub-millimeter scales, which enable new classes of experiments to probe the functional organization of the human brain. This review article surveys advanced data analysis strategies developed for high-resolution fMRI at UHF. These include strategies designed to mitigate distortion and artifacts associated with higher fields in ways that attempt to preserve spatial resolution of the fMRI data, as well as recently introduced analysis techniques that are enabled by these extremely high-resolution data. Particular focus is placed on anatomically-informed analyses, including cortical surface-based analysis, which are powerful techniques that can guide each step of the analysis from preprocessing to statistical analysis to interpretation and visualization. New intracortical analysis techniques for laminar and columnar fMRI are also reviewed and discussed. Prospects for single-subject individualized analyses are also presented and discussed. Altogether, there are both specific challenges and opportunities presented by UHF-fMRI, and the use of proper analysis strategies can help these valuable data reach their full potential.

Reduction of across-run variability of temporal SNR in accelerated EPI time-series data through FLEET-based robust autocalibration

Temporal signal-to-noise ratio (tSNR) is a key metric for assessing the ability to detect brain activation in fMRI data. A recent study has shown substantial variation of tSNR between multiple runs of accelerated EPI acquisitions reconstructed with the GRAPPA method using protocols commonly used for fMRI experiments. Across-run changes in the location of high-tSNR regions could lead to misinterpretation of the observed brain activation patterns, reduced sensitivity of the fMRI studies, and biased results. We compared conventional EPI autocalibration (ACS) methods with the recently-introduced FLEET ACS method, measuring their tSNR variability, as well as spatial overlap and displacement of high-tSNR clusters across runs in datasets acquired from human subjects at 7T and 3T. FLEET ACS reconstructed data had higher tSNR levels, as previously reported, as well as better temporal consistency and larger overlap of the high-tSNR clusters across runs compared with reconstructions using conventional multi-shot (ms) EPI ACS data. tSNR variability across two different runs of the same protocol using ms-EPI ACS data was about two times larger than for the protocol using FLEET ACS for acceleration factors (R) 2 and 3, and one and half times larger for R=4. The level of across-run tSNR consistency for data reconstructed with FLEET ACS was similar to within-run tSNR consistency. The displacement of high-tSNR clusters across two runs (inter-cluster distance) decreased from ∼8mm in the time-series reconstructed using conventional ms-EPI ACS data to ∼4mm for images reconstructed using FLEET ACS. However, the performance gap between conventional ms-EPI ACS and FLEET ACS narrowed with increasing parallel imaging acceleration factor. Overall, the FLEET ACS method provides a simple solution to the problem of varying tSNR across runs, and therefore helps ensure that an assumption of fMRI analysis-that tSNR is largely consistent across runs-is met for accelerated acquisitions.

Amygdala subnuclei response and connectivity during emotional processing

The involvement of the human amygdala in emotion-related processing has been studied using functional magnetic resonance imaging (fMRI) for many years. However, despite the amygdala being comprised of several subnuclei, most studies investigated the role of the entire amygdala in processing of emotions. Here we combined a novel anatomical tracing protocol with event-related high-resolution fMRI acquisition to study the responsiveness of the amygdala subnuclei to negative emotional stimuli and to examine intra-amygdala functional connectivity. The greatest sensitivity to the negative emotional stimuli was observed in the centromedial amygdala, where the hemodynamic response amplitude elicited by the negative emotional stimuli was greater and peaked later than for neutral stimuli. Connectivity patterns converge with extant findings in animals, such that the centromedial amygdala was more connected with the nuclei of the basal amygdala than with the lateral amygdala. Current findings provide evidence of functional specialization within the human amygdala.