Characterization of high-resolution Gradient Echo and Spin Echo EPI for fMRI in the human visual cortex at 7 T

Segmented spin-echo echo-planar imaging improves whole-brain BOLD functional MRI in awake pigeon brains

Functional magnetic resonance imaging (fMRI) in awake small animals such as pigeons or songbirds opens a new window into the neural fundaments of cognitive behavior. However, high-field fMRI in the avian brain is challenging due to strong local magnetic field inhomogeneities caused by air cavities in the skull. A spoiled gradient-echo fMRI sequence has already been used to map the auditory network in songbirds, but due to susceptibility artifacts only 50% of the whole brain could be recorded. Since whole-brain fMRI coverage is vital to reveal whole-brain networks, an MRI sequence that is less susceptible to these artifacts was required. This was recently achieved in various bird species by using a rapid acquisition with relaxation enhancement (RARE) sequence. Weak blood oxygen level-dependent (BOLD) sensitivity, low temporal resolution, and heat caused by the long train of RF refocusing pulses are the main limits of RARE fMRI at high magnetic fields. To go beyond some of these limitations, we here describe the implementation of a two-segmented spin-echo echo-planar imaging (SE-EPI). The proposed sequence covers the whole brain of awake pigeons. The sequence was applied to investigate the auditory network in awake pigeons and assessed the relative merits of this method in comparison with the single-shot RARE sequence. At the same imaging resolution but with a volume acquisition of 3 s versus 4 s for RARE, the two-segmented SE-EPI provided twice the strength of BOLD activity compared with the single-shot RARE sequence, while the image signal-to-noise ratio (SNR) and in particular the temporal SNR were very similar for the two sequences. In addition, the activation patterns in two-segmented SE-EPI data are more symmetric and larger than single-shot RARE results. Two-segmented SE-EPI represents a valid alternative to the RARE sequence in avian fMRI research since it yields more than twice the BOLD sensitivity per unit of time with much less energy deposition and better temporal resolution, particularly for event-related experiments.

Improved sensitivity and microvascular weighting of 3T laminar fMRI with GE-BOLD using NORDIC and phase regression

INTRODUCTION

Functional MRI with spatial resolution in the submillimeter domain enables measurements of activation across cortical layers in humans. This is valuable as different types of cortical computations, e.g., feedforward versus feedback related activity, take place in different cortical layers. Laminar fMRI studies have almost exclusively employed 7T scanners to overcome the reduced signal stability associated with small voxels. However, such systems are relatively rare and only a subset of those are clinically approved. In the present study, we examined if the feasibility of laminar fMRI at 3T could be improved by use of NORDIC denoising and phase regression.

METHODS

5 healthy subjects were scanned on a Siemens MAGNETOM Prisma 3T scanner. To assess across-session reliability, each subject was scanned in 3-8 sessions on 3-4 consecutive days. A 3D gradient echo EPI (GE-EPI) sequence was used for BOLD acquisitions (voxel size 0.82 mm isotopic, TR = 2.2 s) using a block design finger tapping paradigm. NORDIC denoising was applied to the magnitude and phase time series to overcome limitations in temporal signal-to-noise ratio (tSNR) and the denoised phase time series were subsequently used to correct for large vein contamination through phase regression.

RESULTS AND CONCLUSION

NORDIC denoising resulted in tSNR values comparable to or higher than commonly observed at 7T. Layer-dependent activation profiles could thus be extracted robustly, within and across sessions, from regions of interest located in the hand knob of the primary motor cortex (M1). Phase regression led to substantially reduced superficial bias in obtained layer profiles, although residual macrovascular contribution remained. We believe the present results support an improved feasibility of laminar fMRI at 3T.

Correcting for Superficial Bias in 7T Gradient Echo fMRI

The arrival of submillimeter ultra high-field fMRI makes it possible to compare activation profiles across cortical layers. However, the blood oxygenation level dependent (BOLD) signal measured by gradient echo (GE) fMRI is biased toward superficial layers of the cortex, which is a serious confound for laminar analysis. Several univariate and multivariate analysis methods have been proposed to correct this bias. We compare these methods using computational simulations of 7T fMRI data from regions of interest (ROI) during a visual attention paradigm. We also tested the methods on a pilot dataset of human 7T fMRI data. The simulations show that two methods–the ratio of ROI means across conditions and a novel application of Deming regression–offer the most robust correction for superficial bias. Deming regression has the additional advantage that it does not require that the conditions differ in their mean activation over voxels within an ROI. When applied to the pilot dataset, we observed strikingly different layer profiles when different attention metrics were used, but were unable to discern any differences in laminar attention across layers when Deming regression or ROI ratio was applied. Our simulations demonstrates that accurate correction of superficial bias is crucial to avoid drawing erroneous conclusions from laminar analyses of GE fMRI data, and this is affirmed by the results from our pilot 7T fMRI data.

LayNii: A software suite for layer-fMRI

High-resolution fMRI in the sub-millimeter regime allows researchers to resolve brain activity across cortical layers and columns non-invasively. While these high-resolution data make it possible to address novel questions of directional information flow within and across brain circuits, the corresponding data analyses are challenged by MRI artifacts, including image blurring, image distortions, low SNR, and restricted coverage. These challenges often result in insufficient spatial accuracy of conventional analysis pipelines. Here we introduce a new software suite that is specifically designed for layer-specific functional MRI: LayNii. This toolbox is a collection of command-line executable programs written in C/C++ and is distributed opensource and as pre-compiled binaries for Linux, Windows, and macOS. LayNii is designed for layer-fMRI data that suffer from SNR and coverage constraints and thus cannot be straightforwardly analyzed in alternative software packages. Some of the most popular programs of LayNii contain 'layerification' and columnarization in the native voxel space of functional data as well as many other layer-fMRI specific analysis tasks: layer-specific smoothing, model-based vein mitigation of GE-BOLD data, quality assessment of artifact dominated sub-millimeter fMRI, as well as analyses of VASO data.

Effects of phase regression on high-resolution functional MRI of the primary visual cortex

High-resolution functional MRI studies have become a powerful tool to non-invasively probe the sub-millimeter functional organization of the human cortex. Advances in MR hardware, imaging techniques and sophisticated post-processing methods have allowed high resolution fMRI to be used in both the clinical and academic neurosciences. However, consensus within the community regarding the use of gradient echo (GE) or spin echo (SE) based acquisition remains largely divided. On one hand, GE provides a high temporal signal-to-noise ratio (tSNR) technique sensitive to both the macro- and micro-vascular signal while SE based methods are more specific to microvasculature but suffer from lower tSNR and specific absorption rate limitations, especially at high field and with short repetition times. Fortunately, the phase of the GE-EPI signal is sensitive to vessel size and this provides a potential avenue to reduce the macrovascular weighting of the signal (phase regression, Menon 2002). In order to determine the efficacy of this technique at high-resolution, phase regression was applied to GE-EPI timeseries and compared to SE-EPI to determine if GE-EPI's specificity to the microvascular compartment improved. To do this, functional data was collected from seven subjects on a neuro-optimized 7 T system at 800 μm isotropic resolution with both GE-EPI and SE-EPI while observing an 8 Hz contrast reversing checkerboard. Phase data from the GE-EPI was used to create a microvasculature-weighted time series (GE-EPI-PR). Anatomical imaging (MP2RAGE) was also collected to allow for surface segmentation so that the functional results could be projected onto a surface. A multi-echo gradient echo sequence was collected and used to identify venous vasculature. The GE-EPI-PR surface activation maps showed a high qualitative similarity with SE-EPI and also produced laminar activity profiles similar to SE-EPI. When the GE-EPI and GE-EPI-PR distributions were compared to SE-EPI it was shown that GE-EPI-PR had similar distribution characteristics to SE-EPI (p < 0.05) across the top 60% of cortex. Furthermore, it was shown that GE-EPI-PR has a higher contrast-to-noise ratio (0.5 ± 0.2, mean ± std. dev. across layers) than SE-EPI (0.27 ± 0.07) demonstrating the technique has higher sensitivity than SE-EPI. Taken together this evidence suggests phase regression is a useful method in low SNR studies such as high-resolution fMRI.

Forging a path to mesoscopic imaging success with ultra-high field functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) studies with ultra-high field (UHF, 7+ Tesla) technology enable the acquisition of high-resolution images. In this work, we discuss recent achievements in UHF fMRI at the mesoscopic scale, on the order of cortical columns and layers, and examine approaches to addressing common challenges. As researchers push to smaller and smaller voxel sizes, acquisition and analysis decisions have greater potential to degrade spatial accuracy, and UHF fMRI data must be carefully interpreted. We consider the impact of acquisition decisions on the spatial specificity of the MR signal with a representative dataset with 0.8 mm isotropic resolution. We illustrate the trade-offs in contrast with noise ratio and spatial specificity of different acquisition techniques and show that acquisition blurring can increase the effective voxel size by as much as 50% in some dimensions. We further describe how different sources of degradations to spatial resolution in functional data may be characterized. Finally, we emphasize that progress in UHF fMRI depends not only on scientific discovery and technical advancement, but also on informal discussions and documentation of challenges researchers face and overcome in pursuit of their goals. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

New acquisition techniques and their prospects for the achievable resolution of fMRI

This work reviews recent advances in technologies for functional magnetic resonance imaging (fMRI) of the human brain and highlights the push for higher functional specificity based on increased spatial resolution and specific MR contrasts to reveal previously undetectable functional properties of small-scale cortical structures. We discuss how the combination of MR hardware, advanced acquisition techniques and various MR contrast mechanisms have enabled recent progress in functional neuroimaging. However, these advanced fMRI practices have only been applied to a handful of neuroscience questions to date, with the majority of the neuroscience community still using conventional imaging techniques. We thus discuss upcoming challenges and possibilities for fMRI technology development in human neuroscience. We hope that readers interested in functional brain imaging acquire an understanding of current and novel developments and potential future applications, even if they don't have a background in MR physics or engineering. We summarize the capabilities of standard fMRI acquisition schemes with pointers to relevant literature and comprehensive reviews and introduce more recent developments.

Interhemispheric functional connectivity in the zebra finch brain, absent the corpus callosum in normal ontogeny

Bilaterally symmetric intrinsic brain activity (homotopic functional connectivity; FC) is a fundamental feature of the mammalian brain's functional architecture. In mammals, homotopic FC is primarily mediated by the corpus callosum (CC), a large interhemispheric white matter tract thought to balance the bilateral coordination and hemispheric specialization critical for many complex brain functions, including human language. The CC first emerged with the Eutherian (placental) mammals ∼160 MYA and is not found among other vertebrates. Despite this, other vertebrates also exhibit complex brain functions requiring hemispheric specialization and coordination. For example, the zebra finch (Taeniopygia guttata) songbird learns to sing from tutors much as humans acquire speech and must balance hemispheric specialization and coordination to successfully learn and produce song. We therefore tested whether the zebra finch also exhibits homotopic FC, despite lacking the CC. Resting-state fMRI analyses demonstrated widespread homotopic FC throughout the zebra finch brain across development, including within a network required for learned song that lacks direct interhemispheric structural connectivity. The presence of homotopic FC in a non-Eutherian suggests that ancestral pathways, potentially including indirect connectivity via the anterior commissure, are sufficient for maintaining a homotopic functional architecture, an insight with broad implications for understanding interhemispheric coordination across phylogeny.

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.

Improving fMRI in signal drop-out regions at 7 T by using tailored radio-frequency pulses: application to the ventral occipito-temporal cortex

OBJECTIVE

Signal drop-off occurs in echo-planar imaging in inferior brain areas due to field gradients from susceptibility differences between air and tissue. Tailored-RF pulses based on a hyperbolic secant (HS) have been shown to partially recover signal at 3 T, but have not been tested at higher fields.

MATERIALS AND METHODS

The aim of this study was to compare the performance of an optimized tailored-RF gradient-echo echo-planar imaging (TRF GRE-EPI) sequence with standard GRE-EPI at 7 T, in a passive viewing of faces or objects fMRI paradigm in healthy subjects.

RESULTS

Increased temporal-SNR (tSNR) was observed in the middle and inferior temporal lobes and orbitofrontal cortex of all subjects scanned, but elsewhere tSNR decreased relative to the standard acquisition. In the TRF GRE-EPI, increased functional signal was observed in the fusiform, lateral occipital cortex, and occipital pole, regions known to be part of the visual pathway involved in face-object perception.

CONCLUSION

This work highlights the potential of TRF approaches at 7 T. Paired with a reversed-gradient distortion correction to compensate for in-plane susceptibility gradients, it provides an improved acquisition strategy for future neurocognitive studies at ultra-high field imaging in areas suffering from static magnetic field inhomogeneities.