Ultra-high resolution blood volume fMRI and BOLD fMRI in humans at 9.4 T: Capabilities and challenges

Functional mapping of cerebral blood volume (CBV) changes has the potential to reveal brain activity with high localization specificity at the level of cortical layers and columns. Non-invasive CBV imaging using Vascular Space Occupancy (VASO) at ultra-high magnetic field strengths promises high spatial specificity but poses unique challenges in human applications. As such, 9.4 T B1+ and B0 inhomogeneities limit efficient blood tagging, while the specific absorption rate (SAR) constraints limit the application of VASO-specific RF pulses. Moreover, short T2* values at 9.4 T require short readout duration, and long T1 values at 9.4 T can cause blood-inflow contaminations. In this study, we investigated the applicability of layer-dependent CBV-fMRI at 9.4 T in humans. We addressed the aforementioned challenges by combining multiple technical advancements: temporally alternating pTx B1+ shimming parameters, advanced adiabatic RF-pulses, 3D-EPI signal readout, optimized GRAPPA acquisition and reconstruction, and stability-optimized RF channel combination. We found that a combination of suitable advanced methodology alleviates the challenges and potential artifacts, and that VASO fMRI provides reliable measures of CBV change across cortical layers in humans at 9.4 T. The localization specificity of CBV-fMRI, combined with the high sensitivity of 9.4 T, makes this method an important tool for future studies investigating cortical micro-circuitry in humans.

Local excitation universal parallel transmit pulses at 9.4T

PURPOSE

To demonstrate that the concept of "universal pTx pulses" is applicable to local excitation applications.

METHODS

A database of B₀ / B 1 + maps from eight different subjects was acquired at 9.4T. Based on these maps, universal pulses that aim at local excitation of the visual cortex area in the human brain (with a flip angle of 90° or 7°) were calculated. The remaining brain regions should not experience any excitation. The pulses were designed with an extension of the "spatial domain method." A 2D and a 3D target excitation pattern were tested, respectively. The pulse performance was examined on non-database subjects by Bloch simulations and in vivo at 9.4T using a GRE anatomical MRI and a presaturated TurboFLASH B 1 + mapping sequence.

RESULTS

The calculated universal pulses show excellent performance in simulations and in vivo on subjects that were not contained in the design database. The visual cortex region is excited, while the desired non-excitation areas produce the only minimal signal. In simulations, the pulses with 3D target pattern show a lack of excitation uniformity in the visual cortex region; however, in vivo, this inhomogeneity can be deemed acceptable. A reduced field of view application of the universal pulse design concept was performed successfully.

CONCLUSIONS

The proposed design approach creates universal local excitation pulses for a flip angle of 7° and 90°, respectively. Providing universal pTx pulses for local excitation applications prospectively abandons the need for time-consuming subject-specific B₀ / B 1 + mapping and pTx-pulse calculation during the scan session.

Assessment of radio-frequency heating of a parallel transmit coil in a phantom using multi-echo proton resonance frequency shift thermometry

We propose a workflow for validating parallel transmission (pTx) radio-frequency (RF) magnetic field heating patterns using Proton-Resonance Frequency shift (PRF)-based MR thermometry. Electromagnetic (EM) and thermal simulations of a 7 T 8-channel dipole coil were done using commercially available software (Sim4Life) to assess RF heating. The fabrication method for a phantom with electrical properties matched to human tissue is also described, along with methods for its electrical and thermal characterisation. Energy was deposited to specific transmit channels, whilst acquiring 3D PRF data using a pair of interleaved RF shim transmit modes. A multi-echo readout and pre-scan stabilisation protocol were used for increased sensitivity and to correct for measurement-to-measurement instabilities. The electrical properties of the phantom were found to be within 10% of the intended values. Adoption of a 14-min stabilisation scan gave sufficient suppression of any evolving background spatial variation in the B₀ field to achieve <0.001 °C/mm thermometry drift over 10 min of subsequent scanning. Using two RF shim transmit modes enabled full phantom coverage and combining multiple echo times enabled a 13-54% improvement in the RMSE sensitivity to temperature changes. Combining multiple echoes reduced the peak RMSE by 45% and visually reduced measurement-to-measurement instabilities. A reference fibre optic probe showed temperature deviations from the PRF-estimated temperature to be smaller than 0.5 °C. Given the importance of RF safety in pTx applications, this workflow enables accurate validation of RF heating simulations with minimal additional hardware requirements.

RF pulse methods for use with surface coils: Frequency-modulated pulses and parallel transmission

The first use of a surface coil to obtain a 31P NMR spectrum from an intact rat by Ackerman and colleagues initiated a revolution in magnetic resonance imaging (MRI) and spectroscopy (MRS). Today, we take it for granted that one can detect signals in regions external to an RF coil; at the time, however, this concept was most unusual. In the approximately four decade long period since its introduction, this simple idea gave birth to an increasing number of innovations that has led to transformative changes in the way we collect data in an in vivo magnetic resonance experiment, particularly with MRI of humans. These innovations include spatial localization and/or encoding based on the non-uniform B1 field generated by the surface coil, leading to new spectroscopic localization methods, image acceleration, and unique RF pulses that deal with B1 inhomogeneities and even reduce power deposition. Without the surface coil, many of the major technological advances that define the extraordinary success of MRI in clinical diagnosis and in biomedical research, as exemplified by projects like the Human Connectome Project, would not have been possible.

Magnetic resonance imaging at 9.4 T: the Maastricht journey

The 9.4 T scanner in Maastricht is a whole-body magnet with head gradients and parallel RF transmit capability. At the time of the design, it was conceptualized to be one of the best fMRI scanners in the world, but it has also been used for anatomical and diffusion imaging. 9.4 T offers increases in sensitivity and contrast, but the technical ultra-high field (UHF) challenges, such as field inhomogeneities and constraints set by RF power deposition, are exacerbated compared to 7 T. This article reviews some of the 9.4 T work done in Maastricht. Functional imaging experiments included blood oxygenation level-dependent (BOLD) and blood-volume weighted (VASO) fMRI using different readouts. BOLD benefits from shorter T₂* at 9.4 T while VASO from longer T₁. We show examples of both ex vivo and in vivo anatomical imaging. For many applications, pTx and optimized coils are essential to harness the full potential of 9.4 T. Our experience shows that, while considerable effort was required compared to our 7 T scanner, we could obtain high-quality anatomical and functional data, which illustrates the potential of MR acquisitions at even higher field strengths. The practical challenges of working with a relatively unique system are also discussed.

OTUP workflow: target specific optimization of the transmit k-space trajectory for flexible universal parallel transmit RF pulse design

PURPOSE

To optimize transmit k-space trajectories for a wide range of excitation targets and to design "universal pTx RF pulses" based on these trajectories.

METHODS

Transmit k-space trajectories (stack of spirals and SPINS) were optimized to best match different excitation targets using the parameters of the analytical equations of spirals and SPINS. The performances of RF pulses designed based on optimized and non-optimized trajectories were compared. The optimized trajectories were utilized for universal pulse design. The universal pulse performances were compared with subject specific tailored pulse performances. The OTUP workflow (optimization of transmit k-space trajectories and universal pulse calculation) was tested on three test target excitation patterns. For one target (local excitation of a central area in the human brain) the pulses were tested in vivo at 9.4 T.

RESULTS

The workflow produced appropriate transmit k-space trajectories for each test target. Utilization of an optimized trajectory was crucial for the pulse performance. Using unsuited trajectories diminished the performance. It was possible to create target specific universal pulses. However, not every test target is equally well suited for universal pulse design. There was no significant difference in the in vivo performance between subject specific tailored pulses and a universal pulse at 9.4 T.

CONCLUSIONS

The proposed workflow further exploited and improved the universal pulse concept by combining it with gradient trajectory optimization for stack of spirals and SPINS. It emphasized the importance of a well suited trajectory for pTx RF pulse design. Universal and tailored pulses performed with a sufficient degree of similarity in simulations and a high degree of similarity in vivo. The implemented OTUP workflow and the B₀ /B₁⁺ map data from 18 subjects measured at 9.4 T are available as open source (https://github.com/ole1965/workflow_OTUP.git).

Dynamic parallel transmit diffusion MRI at 7T

Diffusion MRI (dMRI) is inherently limited by SNR. Scanning at 7 T increases intrinsic SNR but 7 T MRI scans suffer from regions of signal dropout, especially in the temporal lobes and cerebellum. We applied dynamic parallel transmit (pTx) to allow whole-brain 7 T dMRI and compared with circularly polarized (CP) pulses in 6 subjects. Subject-specific 2-spoke dynamic pTx pulses were designed offline for 8 slabs covering the brain. We used vendor-provided B0 and B1+ mapping. Spokes positions were set using the Fourier difference approach, and RF coefficients optimized with a Jacobi-matrix high-flip-angle optimizer. Diffusion data were analyzed with FSL. Comparing whole-brain averages for pTx against CP scans: mean flip angle error improved by 15% for excitation (2-spoke-VERSE 15.7° vs CP 18.4°, P = 0.012) and improved by 14% for refocusing (2-spoke-VERSE 39.7° vs CP 46.2°, P = 0.008). Computed spin-echo signal standard deviation improved by 14% (2-spoke-VERSE 0.185 vs 0.214 CP, P = 0.025). Temporal SNR increased by 5.4% (2-spoke-VERSE 8.47 vs CP 8.04, P = 0.004) especially in the inferior temporal lobes. Diffusion fitting uncertainty decreased by 6.2% for first fibers (2-spoke VERSE 0.0655 vs CP 0.0703, P < 0.001) and 1.3% for second fibers (2-spoke VERSE 0.139 vs CP 0.141, P = 0.01). In conclusion, dynamic parallel transmit improves the uniformity of 7 T diffusion-weighted imaging. In future, less restrictive SAR limits for parallel transmit scans are expected to allow further improvements.

Hybrid B 1 + -shimming and gradient adaptions for improved pseudo-continuous arterial spin labeling at 7 Tesla

PURPOSE

To improve pseudo-continuous arterial spin labeling (pcASL) at 7T by exploiting a hybrid homogeneity- and efficiency-optimized B 1 + -shim with adapted gradient strength as well as background suppression.

METHODS

The following three experiments were performed at 7T, each employing five volunteers: (1) A hybrid (ie, homogeneity-efficiency optimized) B 1 + -shim was introduced and evaluated for variable-rate selective excitation pcASL labeling. Therefore, B 1 + -maps in the V3 segment and time-of-flight images were acquired to identify the feeding arteries. For validation, a gradient-echo sequence was applied in circular polarized (CP) mode and with the hybrid B 1 + -shim. Additionally, the gray matter (temporal) signal-to-noise ratio (tSNR) in pcASL perfusion images was evaluated. (2) Bloch simulations for the pcASL labeling were conducted and validated experimentally, with a focus on the slice-selective gradients. (3) Background suppression was added to the B 1 + -shimmed, gradient-adapted 7T sequence and this was then compared to a matched sequence at 3T.

RESULTS

The B 1 + -shim improved the signal within the labeling plane (23.6%) and the SNR/tSNR increased (+11%/+11%) compared to its value in CP mode; however, the increase was not significant. In accordance with the simulations, the adapted gradients increased the tSNR (35%) and SNR (45%) significantly. Background suppression further improved the perfusion images at 7T, and this protocol performed as well as a resolution-matched protocol at 3T.

CONCLUSION

The combination of the proposed hybrid B 1 + -phase-shim with the adapted slice-selective gradients and background suppression shows great potential for improved pcASL labeling under suboptimal B 1 + conditions at 7T.