The possible influence of third-order shim coils on gradient-magnet interactions: an inter-field and inter-site study

OBJECTIVE

To assess the possible influence of third-order shim coils on the behavior of the gradient field and in gradient-magnet interactions at 7 T and above.

MATERIALS AND METHODS

Gradient impulse response function measurements were performed at 5 sites spanning field strengths from 7 to 11.7 T, all of them sharing the same exact whole-body gradient coil design. Mechanical fixation and boundary conditions of the gradient coil were altered in several ways at one site to study the impact of mechanical coupling with the magnet on the field perturbations. Vibrations, power deposition in the He bath, and field dynamics were characterized at 11.7 T with the third-order shim coils connected and disconnected inside the Faraday cage.

RESULTS

For the same whole-body gradient coil design, all measurements differed greatly based on the third-order shim coil configuration (connected or not). Vibrations and gradient transfer function peaks could be affected by a factor of 2 or more, depending on the resonances. Disconnecting the third-order shim coils at 11.7 T also suppressed almost completely power deposition peaks at some frequencies.

DISCUSSION

Third-order shim coil configurations can have major impact in gradient-magnet interactions with consequences on potential hardware damage, magnet heating, and image quality going beyond EPI acquisitions.

In vivo B 1 + enhancement of calf MRI at 7 T via optimized flexible metasurfaces

PURPOSE

Ultrahigh field (≥7 T) MRI is at the cutting edge of medical imaging, enabling enhanced spatial and spectral resolution as well as enhanced susceptibility contrast. However, transmit (B 1 + $$ {\mathrm{B}}_1^{+} $$ ) field inhomogeneity due to standing wave effects caused by the shortened RF wavelengths at 7 T is still a challenge to overcome. Novel hardware methods such as dielectric pads have been shown to improve theB 1 + $$ {\mathrm{B}}_1^{+} $$ field inhomogeneity but are currently limited in their corrective effect by the range of high-permittivity materials available and have a fixed shelf life. In this work, an optimized metasurface design is presented that demonstrates in vivo enhancement of theB 1 + $$ {\mathrm{B}}_1^{+} $$ field.

METHODS

A prototype metasurface was optimized by an empirical capacitor sweep and by varying the period size. Phantom temperature experiments were performed to evaluate potential metasurface heating effects during scanning. Lastly, in vivo gradient echo images andB 1 + $$ {\mathrm{B}}_1^{+} $$ maps were acquired on five healthy subjects on a 7 T system. Dielectric pads were also used as a comparison throughout the work as a standard comparison.

RESULTS

The metasurfaces presented here enhanced the average relative SNR of the gradient echo images by a factor of 2.26 compared to the dielectric pads factor of 1.61. AverageB 1 + $$ {\mathrm{B}}_1^{+} $$ values reflected a similar enhancement of 27.6% with the metasurfaces present versus 8.9% with the dielectric pads.

CONCLUSION

The results demonstrate that metasurfaces provide superior performance to dielectric padding as shown byB 1 + $$ {\mathrm{B}}_1^{+} $$ maps reflecting their direct effects and resulting enhancements in image SNR at 7 T.

Simultaneous multi-region detection of GABA+ and Glx using 3D spatially resolved SLOW-editing and EPSI-readout at 7T

GABA+ and Glx (glutamate and glutamine) are widely studied metabolites, yet the commonly used magnetic resonance spectroscopy (MRS) techniques have significant limitations, including sensitivity to B0 and B1+-inhomogeneities, limited bandwidth of MEGA-pulses, high SAR which is accentuated at 7T. To address these limitations, we propose SLOW-EPSI method, employing a large 3D MRSI coverage and achieving a high resolution down to 0.26 ml. Simulation results demonstrate the robustness of SLOW-editing for both GABA+ and Glx against B0 and B1+-inhomogeneities within the range of [-0.3, +0.3] ppm and [40 %, 250 %], respectively. Two protocols, both utilizing a 70 mm thick FOV slab, were employed to target distinct brain regions in vivo, differentiated by their orientation: transverse and tilted. Protocol 1 (n = 11) encompassed 5 locations (cortical gray matter, white matter, frontal lobe, parietal lobe, and cingulate gyrus). Protocol 2 (n = 5) involved 9 locations (cortical gray matter, white matter, frontal lobe, occipital lobe, cingulate gyrus, caudate nucleus, hippocampus, putamen, and inferior thalamus). Quantitative analysis of GABA+ and Glx was conducted in a stepwise manner. First, B1+/B1--inhomogeneities were corrected using water reference data. Next, GABA+ and Glx values were calculated employing spectral fitting. Finally, the GABA+ level for each selected region was compared to the global Glx within the same subject, generating the GABA+/Glx_global ratio. Our findings from two protocols indicate that the GABA+/Glx_global level in cortical gray matter was approximately 16 % higher than in white matter. Elevated GABA+/Glx_global levels acquired with protocol 2 were observed in specific regions such as the caudate nucleus (0.118±0.067), putamen (0.108±0.023), thalamus (0.092±0.036), and occipital cortex (0.091±0.010), when compared to the cortical gray matter (0.079±0.012). Overall, our results highlight the effectiveness of SLOW-EPSI as a robust and efficient technique for accurate measurements of GABA+ and Glx at 7T. In contrast to previous SVS and 2D-MRSI based editing sequences with which only one or a limited number of brain regions can be measured simultaneously, the method presented here measures GABA+ and Glx from any brain area and any arbitrarily shaped volume that can be flexibly selected after the examination. Quantification of GABA+ and Glx across multiple brain regions through spectral fitting is achievable with a 9-minute acquisition. Additionally, acquisition times of 18-27 min (GABA+) and 9-18 min (Glx) are required to generate 3D maps, which are constructed using Gaussian fitting and peak integration.

Double cross magnetic wall decoupling for quadrature transceiver RF array coils using common-mode differential-mode resonators

In contrast to linearly polarized RF coil arrays, quadrature transceiver coil arrays are capable of improving signal-to-noise ratio (SNR), spatial resolution, and parallel imaging performance. Owing to a reduced excitation power, a low specific absorption rate can also be obtained using quadrature RF coils. However, due to the complex nature of their structure and their electromagnetic properties, it is challenging to achieve sufficient electromagnetic decoupling while designing multichannel quadrature RF coil arrays, particularly in ultra-high fields. In this work, we proposed a double-cross magnetic wall decoupling for quadrature transceiver RF arrays and implemented the decoupling method on common-mode differential mode quadrature (CMDM) quadrature transceiver arrays at an ultrahigh field of 7 T. The proposed magnetic decoupling wall, comprised of two intrinsically decoupled loops, is used to reduce the mutual coupling between all the multi-mode currents present in the quadrature CMDM array. The decoupling network has no physical connection with the CMDMs' resonators, which provides less design constraint over size-adjustable RF arrays. To validate the feasibility of the proposed cross-magnetic decoupling wall, systematic studies on the decoupling performance based on the impedance of two intrinsic loops are numerically performed. A pair of quadrature transceiver CMDMs is constructed along with the proposed decoupling network, and their scattering matrix is characterized using a network analyzer. The measured results indicate that all the current modes from coupling are simultaneously suppressed using the proposed cross-magnetic wall. Moreover, field distribution and local specific absorption rate (SAR) are numerically obtained for a well-decoupled 8-channel quadrature knee-coil array.

Neuromelanin related ultra-high field signal intensity of the locus coeruleus differs between Parkinson's disease and controls

INTRODUCTION

Neuromelanin related signal changes in catecholaminergic nuclei are considered as a promising MRI biomarker in Parkinson's disease (PD). Until now, most studies have investigated the substantia nigra (SN), while signal changes might be more prominent in the locus coeruleus (LC). Ultra-high field MRI improves the visualisation of these small brainstem regions and might support the development of imaging biomarkers in PD.

OBJECTIVES

To compare signal intensity of the SN and LC on Magnetization Transfer MRI between PD patients and healthy controls (HC) and to explore its association with cognitive performance in PD.

METHODS

This study was conducted using data from the TRACK-PD study, a longitudinal 7T MRI study. A total of 78 early-stage PD patients and 36 HC were included. A mask for the SN and LC was automatically segmented and manually corrected. Neuromelanin related signal intensity of the SN and LC was compared between PD and HC.

RESULTS

PD participants showed a lower contrast-to-noise ratio (CNR) in the right SN (p = 0.029) and left LC (p = 0.027). After adding age as a confounder, the CNR of the right SN did not significantly differ anymore between PD and HC (p = 0.055). Additionally, a significant positive correlation was found between the SN CNR and memory function.

DISCUSSION

This study confirms that neuromelanin related signal intensity of the LC differs between early-stage PD patients and HC. No significant difference was found in the SN. This supports the theory of bottom-up disease progression in PD. Furthermore, loss of SN integrity might influence working memory or learning capabilities in PD patients.

Influence of adaptive denoising on Diffusion Kurtosis Imaging at 3T and 7T

BACKGROUND AND OBJECTIVE

Choosing the most appropriate denoising method to improve the quality of diagnostic images maximally is key in pre-processing of diffusion MRI images. Recent advancements in acquisition and reconstruction techniques have questioned traditional noise estimation methods favoring adaptive denoising frameworks, circumventing the need to know a priori information that is hardly available in a clinical setting. In this observational study, we compared two innovative adaptive techniques sharing some features, Patch2Self and Nlsam, through application on reference adult data at 3T and 7T. The primary aim was identifying the most effective method in case of Diffusion Kurtosis Imaging (DKI) data - particularly susceptible to noise and signal fluctuations - at 3T and 7T fields. A side goal consisted of investigating the dependence of kurtosis metrics' variability with respect to the magnetic field on the adopted denoising methodology.

METHODS

For comparison purposes, we focused on qualitative and quantitative analysis of DKI data and related microstructural maps before and after applying the two denoising approaches. Specifically, we assessed computational efficiency, preservation of anatomical details via perceptual metrics, consistency of microstructure model fitting, alleviation of degeneracies in model estimation, and joint variability with varying field strength and denoising method.

RESULTS

Accounting for all these factors, Patch2Self framework has turned out to be specifically suitable for DKI data, with improving performance at 7T. Nlsam method is more robust in alleviating degenerate black voxels while introducing some blurring, which in turn is reflected in an overall loss of image sharpness. Regarding the impact of denoising on field-dependent variability, both methods have been shown to make variations from standard to Ultra-High Field more concordant with theoretical evidence, claiming that kurtosis metrics are sensitive to susceptibility-induced background gradients, directly proportional to the magnetic field strength and sensitive to the microscopic distribution of iron and myelin.

CONCLUSIONS

This study serves as a proof-of-concept stressing the need for an accurate choice of a denoising methodology, specifically tailored for the data under analysis and allowing higher spatial resolution acquisition within clinically compatible timings, with all the potential benefits that improving suboptimal quality of diagnostic images entails.

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.

Commissioning of the Iseult CEA 11.7 T whole-body MRI: current status, gradient-magnet interaction tests and first imaging experience

OBJECTIVES

The Iseult MRI is an actively shielded whole-body magnet providing a homogeneous and stable magnetic field of 11.7 T. After nearly 20 years of research and development, the magnet successfully reached its target field strength for the first time in 2019. This article reviews its commissioning status, the gradient-magnet interaction test results and first imaging experience.

MATERIALS AND METHODS

Vibration, acoustics, power deposition in the He bath, and field monitoring measurements were carried out. Magnet safety system was tested against outer magnetic perturbations, and calibrated to define a safe operation of the gradient coil. First measurements using parallel transmission were also performed on an ex-vivo brain to mitigate the RF field inhomogeneity effect.

RESULTS

Acoustics measurements show promising results with sound pressure levels slightly above the enforced limits only at certain frequency intervals. Vibrations of the gradient coil revealed a linear trend with the B₀ field only in the worst case. Field monitoring revealed some resonances at some frequencies that are still under investigation.

DISCUSSION

Gradient-magnet interaction tests at up to 11.7 T are concluded. The scanner is now kept permanently at field and the final calibrations are on-going to pave the road towards the first acquisitions on volunteers.

Improved laminar specificity and sensitivity by combining SE and GE BOLD signals

The most widely used gradient-echo (GE) blood oxygenation level-dependent (BOLD) contrast has high sensitivity, but low specificity due to draining vein contributions, while spin-echo (SE) BOLD approach at ultra-high magnetic fields is highly specific to neural active sites but has lower sensitivity. To obtain high specificity and sensitivity, we propose to utilize a vessel-size-sensitive filter to the GE-BOLD signal, which suppresses macrovascular contributions and to combine selectively retained microvascular GE-BOLD signals with the SE-BOLD signals. To investigate our proposed idea, fMRI with 0.8 mm isotropic resolution was performed on the primary motor and sensory cortices in humans at 7 T by implementing spin- and gradient-echo (SAGE) echo planar imaging (EPI) acquisition. Microvascular-passed sigmoidal filters were designed based upon the vessel-size-sensitive ΔR₂*/ΔR₂ value for retaining GE-BOLD signals originating from venous vessels with ≤ 45 μm and ≤ 65 μm diameter. Unlike GE-BOLD fMRI, the laminar profile of SAGE-BOLD fMRI with the vessel-size-sensitive filter peaked at ∼ 1.0 mm from the surface of the primary motor and sensory cortices, demonstrating an improvement of laminar specificity over GE-BOLD fMRI. Also, the functional sensitivity of SAGE BOLD at middle layers (0.75-1.5 mm) was improved by ∼ 80% to ∼100% when compared with SE BOLD. In summary, we showed that combined GE- and SE-BOLD fMRI with the vessel-size-sensitive filter indeed yielded improved laminar specificity and sensitivity and is therefore an excellent tool for high spatial resolution ultra-high filed (UHF)-fMRI studies for resolving mesoscopic functional units.

Hairpin RF resonators for MR imaging transceiver arrays with high inter-channel isolation and B1 efficiency at ultrahigh field 7 T

Electromagnetic decoupling among a close-fitting or high-density transceiver RF array elements is required to maintain the integrity of the magnetic flux density from individual channel for enhanced performance in detection sensitivity and parallel imaging. High-impedance RF coils have demonstrated to be a prominent design method to circumvent these coupling issues. Yet, inherent characteristics of these coils have ramification on the B₁ field efficiency and SNR. In this work, we propose a hairpin high impedance RF resonator design for highly decoupled multichannel transceiver arrays at ultrahigh magnetic fields. Due to the high impedance property of the hairpin resonators, the proposed transceiver array can provide high decoupling performance without using any dedicated decoupling circuit among the resonant elements. Because of elimination of lumped inductors in the resonator circuit, higher B₁ field efficiency in imaging subjects can be expected. In order to validate the feasibility of the proposed hairpin RF coils, systematical studies on decoupling performance, field distribution, and SNR are performed, and the results are compared with those obtained from existing high-impedance RF coil, e.g., "self-decoupled RF coil". To further investigate its performance, an 8-channel head coil array using the proposed hairpin resonators loaded with a cylindrical phantom is designed, demonstrating a 19 % increase of the B₁⁺ field intensity compared to the self-decoupled coils at 7 T. Furthermore, the characteristics of the hairpin RF coils are evaluated using a more realistic human head voxel model numerically. The proposed hairpin RF coil provides excellent decoupling performance and superior RF magnetic field efficiency compared to the "self-decoupled" high impedance coils. Bench test of a pair of fabricated hairpin coils prove to be in good accordance with numerical results.

Parallel transmit (pTx) with online pulse design for task-based fMRI at 7 T

PURPOSE

Parallel transmission (pTx) is an approach to improve image uniformity for ultra-high field imaging. In this study, we modified an echo planar imaging (EPI) sequence to design subject-specific pTx pulses online. We compared its performance against EPI with conventional circularly polarised (CP) pulses.

METHODS

We compared the pTx-EPI and CP-EPI sequences in a short EPI acquisition protocol and for two different functional paradigms in six healthy volunteers (2 female, aged 23-36 years, mean age 29.2 years). We chose two paradigms that are typically affected by signal dropout at 7 T: a visual objects localiser to determine face/scene selective brain regions and a semantic-processing task.

RESULTS

Across all subjects, pTx-EPI improved whole-brain mean temporal signal-to-noise ratio (tSNR) by 11.0% compared to CP-EPI. We also compared the ability of pTx-EPI and CP-EPI to detect functional activation for three contrasts over the two paradigms: face > object and scene > object for the visual objects localiser and semantic association > pattern matching for the semantic-processing paradigm. Across all three contrasts, pTx-EPI showed higher median z-scores and detected more active voxels in relevant areas, as determined from previous 3 T studies.

CONCLUSION

We have demonstrated a workflow for EPI acquisitions with online per-subject pulse calculations. We saw improved performance in both tSNR and functional acquisitions from pTx-EPI. Thus, we believe that online calculation pTx-EPI is robust enough for future fMRI studies, especially where activation is expected in brain areas liable to significant signal dropout.

Towards an integrated neonatal brain and cardiac examination capability at 7 T: electromagnetic field simulations and early phantom experiments using an 8-channel dipole array

OBJECTIVE

Neonatal brain and cardiac imaging would benefit from the increased signal-to-noise ratio levels at 7 T compared to lower field. Optimal performance might be achieved using purpose designed RF coil arrays. In this study, we introduce an 8-channel dipole array and investigate, using simulations, its RF performances for neonatal applications at 7 T.

METHODS

The 8-channel dipole array was designed and evaluated for neonatal brain/cardiac configurations in terms of SAR efficiency (ratio between transmit-field and maximum specific-absorption-rate level) using adjusted dielectric properties for neonate. A birdcage coil operating in circularly polarized mode was simulated for comparison. Validation of the simulation model was performed on phantom for the coil array.

RESULTS

The 8-channel dipole array demonstrated up to 46% higher SAR efficiency levels compared to the birdcage coil in neonatal configurations, as the specific-absorption-rate levels were alleviated. An averaged normalized root-mean-square-error of 6.7% was found between measured and simulated transmit field maps on phantom.

CONCLUSION

The 8-channel dipole array design integrated for neonatal brain and cardiac MR was successfully demonstrated, in simulation with coverage of the baby and increased SAR efficiency levels compared to the birdcage. We conclude that the 8Tx-dipole array promises safe operating procedures for MR imaging of neonatal brain and heart at 7 T.

On the relationship between GABA+ and glutamate across the brain

Equilibrium between excitation and inhibition (E/I balance) is key to healthy brain function. Conversely, disruption of normal E/I balance has been implicated in a range of central neurological pathologies. Magnetic resonance spectroscopy (MRS) provides a non-invasive means of quantifying in vivo concentrations of excitatory and inhibitory neurotransmitters, which could be used as diagnostic biomarkers. Using the ratio of excitatory and inhibitory neurotransmitters as an index of E/I balance is common practice in MRS work, but recent studies have shown inconsistent evidence for the validity of this proxy. This is underscored by the fact that different measures are often used in calculating E/I balance such as glutamate and Glx (glutamate and glutamine). Here we used a large MRS dataset obtained at ultra-high field (7 T) measured from 193 healthy young adults and focused on two brain regions - prefrontal and occipital cortex - to resolve this inconsistency. We find evidence that there is an inter-individual common ratio between GABA+ (γ-aminobutyric acid and macromolecules) and Glx in the occipital, but not prefrontal cortex. We further replicate the prefrontal result in a legacy dataset (n = 78) measured at high-field (3 T) strength. By contrast, with ultra-high field MRS data, we find extreme evidence that there is a common ratio between GABA+ and glutamate in both prefrontal and occipital cortices, which cannot be explained by participant demographics, signal quality, fractional tissue volume, or other metabolite concentrations. These results are consistent with previous electrophysiological and theoretical work supporting E/I balance. Our findings indicate that MRS-detected GABA+ and glutamate (but not Glx), are a reliable measure of E/I balance .

Nonparametric 5D D-R2 distribution imaging with single-shot EPI at 21.1 T: Initial results for in vivo rat brain

In vivo human diffusion MRI is by default performed using single-shot EPI with greater than 50-ms echo times and associated signal loss from transverse relaxation. The individual benefits of the current trends of increasing B₀ to boost SNR and employing more advanced signal preparation schemes to improve the specificity for selected microstructural properties eventually may be cancelled by increased relaxation rates at high B₀ and echo times with advanced encoding. Here, initial attempts to translate state-of-the-art diffusion-relaxation correlation methods from 3 T to 21.1 T are made to identify hurdles that need to be overcome to fulfill the promises of both high SNR and readily interpretable microstructural information.

Physical limits to human brain B0 shimming with spherical harmonics, engineering implications thereof

OBJECTIVE

As the MRI main magnetic field rises for improved signal-to-noise ratio, susceptibility-induced B0-inhomogeneity increases proportionally, aggravating related artifacts. Considering only susceptibility disparities between air and biological tissue, we explore the topological conditions for which perfect shimming could be performed in a Region of Interest (ROI) such as the human brain or part thereof.

MATERIALS AND METHODS

After theoretical considerations for perfect shimming, spherical harmonic (SH) shimming simulations of very high degree are performed, based on a 100-subject database of 1.7-mm-resolved brain fieldmaps acquired at 3T . In addition to the whole brain, shimmed ROIs include slabs targeting the prefrontal cortex, both or single temporal lobes, or spheres in the frontal brain above the nasal sinus.

RESULTS AND DISCUSSION

We show "perfect" SH shimming is possible only if the ROI can be contained in a sphere that does not enclose sources of magnetic field inhomogeneity, which are gathered at the air-tissue interface. We establish a [Formula: see text]Hz inhomogeneity hard shim limit at 7T for whole brain SH shimming, that can only be attained at shimming degree higher than 90. On the other hand, under limited power and SH degree resources, 3D region-specific shimming is shown to greatly improve homogeneity in critical zones such as the prefrontal cortex and around ear canals.

A novel male Japanese quail structural connectivity atlas using ultra-high field diffusion MRI at 11.7 T

The structural connectivity of animal brains can be revealed using post-mortem diffusion-weighted magnetic resonance imaging (MRI). Despite the existence of several structural atlases of avian brains, few of them address the bird’s structural connectivity. In this study, a novel atlas of the structural connectivity is proposed for the male Japanese quail (Coturnix japonica), aiming at investigating two lines divergent on their emotionality trait: the short tonic immobility (STI) and the long tonic immobility (LTI) lines. The STI line presents a low emotionality trait, while the LTI line expresses a high emotionality trait. 21 male Japanese quail brains from both lines were scanned post-mortem for this study, using a preclinical Bruker 11.7 T MRI scanner. Diffusion-weighted MRI was performed using a 3D segmented echo planar imaging (EPI) pulsed gradient spin-echo (PGSE) sequence with a 200 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document}μm isotropic resolution, 75 diffusion-encoding directions and a b-value fixed at 4500 s/mm². Anatomical MRI was likewise performed using a 2D anatomical T₂-weighted spin-echo (SE) sequence with a 150 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upmu$$\end{document}μm isotropic resolution. This very first anatomical connectivity atlas of the male Japanese quail reveals 34 labeled fiber tracts and the existence of structural differences between the connectivity patterns characterizing the two lines. Thus, the link between the male Japanese quail’s connectivity and its underlying anatomical structures has reached a better understanding.

Emerging methods and applications of ultra-high field MR spectroscopic imaging in the human brain

Magnetic Resonance Spectroscopic Imaging (MRSI) of the brain enables insights into the metabolic changes and fluxes in diseases such as tumors, multiple sclerosis, epilepsy, or hepatic encephalopathy, as well as insights into general brain functionality. However, the routine application of MRSI is mostly hampered by very low signal-to-noise ratios (SNR) due to the low concentrations of metabolites, about 10000 times lower than water. Furthermore, MRSI spectra have a dense information content with many overlapping metabolite resonances, especially for proton MRSI. MRI scanners at ultra-high field strengths, like 7 T or above, offer the opportunity to increase SNR, as well as the separation between resonances, thus promising to solve both challenges. Yet, MRSI at ultra-high field strengths is challenged by decreased B0- and B1-homogeneity, shorter T2 relaxation times, stronger chemical shift displacement errors, and aggravated lipid contamination. Therefore, to capitalize on the advantages of ultra-high field strengths, these challenges must be overcome. This review focuses on the challenges MRSI of the human brain faces at ultra-high field strength, as well as the possible applications to this date.

Automatic segmentation of Deep Grey Nuclei using a high-resolution 7T MRI Atlas - quantification of T1 values in healthy volunteers

We present a new consensus atlas of deep grey nuclei obtained by shape-based averaging of manual segmentation of two experienced neuroradiologists and optimized from 7T MP2RAGE images acquired at (0.6mm)³ in 60 healthy subjects. A group-wise normalization method was used to build a high-contrast and high-resolution T₁ -weighted brain template (0.5mm)³ using data from 30 out of the 60 controls. Delineation of 24 deep grey nuclei per hemisphere, including the claustrum and twelve thalamic nuclei, was then performed by two expert neuroradiologists and reviewed by a third neuroradiologist according to tissue contrast and external references based on the Morel atlas. Corresponding deep grey matter structures were also extracted from the Morel and CIT168 atlases. The data-derived, Morel and CIT168 atlases were all applied at the individual level using non-linear registration to fit the subject reference and to extract absolute mean quantitative T₁ values derived from the 3D-MP2RAGE volumes, after correction for residual B₁ ⁺ biases. Three metrics (The Dice and the volumetric similarity coefficients, and a novel Hausdorff distance) were used to estimate the inter-rater agreement of manual MRI segmentation and inter-atlas variability, and these metrics were measured to quantify biases due to image registration and their impact on the measurements of the quantitative T₁ values was highlighted. This represents a fully-automated segmentation process permitting the extraction of unbiased normative T₁ values in a population of young healthy controls as a reference for characterizing subtle structural alterations of deep grey nuclei relevant to a range of neurological diseases.

Concurrent CBF and bold FMRI with dual-echo spiral simultaneous multi-slice acquisitions at 7T

Measurement of cerebral blood flow (CBF) using the Arterial Spin Labeling (ASL) technique is a desirable fMRI approach due to the higher specificity of CBF to the site of neural activation. However, ASL has inherent limitations, such as a low signal-to-noise ratio (SNR) and low coverage/resolution due to the limited readout window following the labeling. Recently, ASL has been implemented at ultra-high field (UHF) strengths in an attempt to mitigate the SNR challenges. Even though ASL intrinsically allows concurrent acquisition of CBF and BOLD contrasts, a compromise in the echo time (TE) for either of the contrasts is inevitable with single-echo acquisitions. Long durations of the Cartesian EPI readout do not allow for multi-echo acquisitions for resolutions ≤2 mm where both contrasts can be acquired at their optimal TE at UHF. With its higher acquisition efficiency, single-shot spiral imaging provides a promising alternative to EPI, and with a dual-echo, out-in trajectory allows both CBF and BOLD contrasts to be acquired at their respective optimal TE. In this work, we implemented a dual-echo spiral out-in ASL sequence with simultaneous multi-slice (SMS) readout for increased coverage, and validated its application to fMRI with a visuomotor paradigm. Conventional Cartesian EPI acquisitions with matched parameters served as a reference. The dual-echo spiral ASL acquisitions resulted in robust CBF and BOLD activations maps. The absolute and relative CBF changes measured with the dual-echo spiral readout were in agreement with previous reports in the literature as well as the reference Cartesian acquisitions. The BOLD response amplitude was higher compared to the Cartesian acquisitions, attributable to a more optimal TE of the second echo. In conclusion, dual-echo spiral out-in SMS acquisition shows promise for concurrent acquisitions of BOLD and non-BOLD contrasts that require a short TE, with no loss in temporal resolution.

Simultaneous pure T2 and varying T2'-weighted BOLD fMRI using Echo Planar Time-resolved Imaging for mapping cortical-depth dependent responses

Spin-echo (SE) BOLD fMRI has high microvascular specificity, and thus provides a more reliable means to localize neural activity compared to conventional gradient-echo BOLD fMRI. However, the most common SE BOLD acquisition method, SE-EPI, is known to suffer from T₂′ contrast contamination with undesirable draining vein bias. To address this, in this study, we extended a recently developed distortion/blurring-free multi-shot EPI technique, Echo-Planar Time-resolved Imaging (EPTI), to cortical-depth dependent SE-fMRI at 7T to test whether it could provide purer SE BOLD contrast with minimal T₂′ contamination for improved neuronal specificity. From the same acquisition, the time-resolved feature of EPTI also provides a series of asymmetric SE (ASE) images with varying T₂′ weightings, and enables extraction of data equivalent to conventional SE EPI with different echo train lengths (ETLs). This allows us to systematically examine how T₂′-contribution affects different SE acquisition strategies using a single dataset. A low-rank spatiotemporal subspace reconstruction was implemented for the SE-EPTI acquisition, which incorporates corrections for both shot-to-shot phase variations and dynamic B₀ drifts. SE-EPTI was used in a visual task fMRI experiment to demonstrate that i) the pure SE image provided by EPTI results in the highest microvascular specificity; ii) the ASE EPTI series, with a graded introduction of T₂′ weightings at time points farther away from the pure SE, show a gradual sensitivity increase along with increasing draining vein bias; iii) the longer ETL seen in conventional SE EPI acquisitions will induce more draining vein bias. Consistent results were observed across multiple subjects, demonstrating the robustness of the proposed technique for SE-BOLD fMRI with high specificity.

Epilepsy under the scope of ultra-high field MRI

Ultra-high field magnetic resonance imaging (UHF-MRI) is capable of unraveling anatomical structures in a submillimeter range. In addition, its high resonance regime allows the quantification of constitutive molecules in a spatially sensitive manner, a crucial capability for determining the extent and localization of a probable epileptogenic region or the severity of the epilepsy. The main technical challenges for data acquisition under UHF are to produce a strong, homogeneous transverse field, while keeping the tissue power deposition within the safe regulatory guidelines. The nonuniformities caused by destructive and constructive interferences at UHFs required new technologies to accelerate and increase yield regarding time spent and quality achieved. Image quality is the paramount contribution of UHF high-resolution imaging, which is capable to disclose fine details of the hippocampal formation and its surroundings and their changes in the course of epilepsy. Other sequences like diffusion tensor imaging (DTI) and multiecho susceptibility imaging at 7 T in vivo can assist the creation of normative atlases of the hippocampal subfields or the reconstruction of the highly arborized cerebral blood vessels. In our review, we specify the impact of these advanced relevant techniques onto the study of epilepsy. In this context, we focused onto high field high-resolution scanners and clinically-enriched decision-making. Studies on focal dysplasias correlating ex vivo high-resolution imaging with specific histological and ultrastructural patterns showed that white matter hyperintensities were related to a demyelination process and other alterations. Preliminary results correlating thick serial sections through bioptic epileptogenic tissue could extend the strategy to localize degenerated tissue sectors, correlate nature and extent of tissue loss with preoperative diagnosis and postoperative outcome. Finally, this protocol will provide the neurosurgeon with a detailed depiction of the removed pathologic tissue and possible adverse effects by the pathologic tissue left in situ. This article is part of the special issue "NEWroscience 2018".

Improvement of sensitivity and specificity for laminar BOLD fMRI with double spin-echo EPI in humans at 7 T

Mapping mesoscopic cortical functional units such as columns or laminae is increasingly pursued by ultra-high field (UHF) functional magnetic resonance imaging (fMRI). The most popular approach for high-resolution fMRI is currently gradient-echo (GE) blood oxygenation level-dependent (BOLD) fMRI. However, its spatial accuracy is reduced due to its sensitivity to draining vessels, including pial veins, whereas spin-echo (SE) BOLD signal is expected to have higher spatial accuracy, albeit with lower sensitivity than the GE-BOLD signal. Here, we introduce a new double spin-echo (dSE) echo-planar imaging (EPI) method to improve the sensitivity of SE-BOLD contrast by averaging two spin-echoes using three radiofrequency pulses. Human fMRI experiments were performed with slices perpendicular to the central sulcus between motor and sensory cortices at 7 T during fist-clenching with touching. First, we evaluated the feasibility of single-shot dSE-EPI for BOLD fMRI with 1.5 mm isotropic resolution and found that dSE-BOLD fMRI has higher signal-to-noise ratio (SNR), temporal SNR (tSNR), and higher functional sensitivity than conventional SE-BOLD fMRI. Second, to investigate the laminar specificity of dSE-BOLD fMRI, we implemented a multi-shot approach to achieve 0.8-mm isotropic resolution with sliding-window reconstruction. Unlike GE-BOLD fMRI, the cortical profile of dSE-BOLD fMRI peaked at ~ 1.0 mm from the surface of the primary motor and sensory cortices, demonstrating an improvement of laminar specificity in humans over GE-BOLD fMRI. The proposed multi-shot dSE-EPI method is viable for high spatial resolution UHF-fMRI studies in the pursuit of resolving mesoscopic functional units.

Ultra-high field imaging in Major Depressive Disorder: a review of structural and functional studies

BACKGROUND

Major depressive disorder (MDD) is a severe and pervasive psychiatric condition with a lifetime prevalence of 15-25%. Numerous Magnetic Resonance Imaging (MRI) studies employing scans at field strengths of 1.5T or 3T have been carried out in the last decades, providing an unprecedented insight into the neural correlates of MDD. However, in recent years, MRI technology has largely progressed and the use of scans at ultra-high field (≥ 7T) has improved the sensitivity and the resolution of MR images. In this context, with this review we aim to summarize evidence of structural and functional brain mechanisms underlying MDD obtained with ultra-high field MRI.

METHODS

We conducted a search on PubMed, Scopus and Web of Science of neuroimaging studies on MDD patients, which employed ultra-high field MRI. We detected six structural MRI studies, two Diffusion Tensor Imaging (DTI) studies and five functional MRI (fMRI) studies.

RESULTS

Overall, the MRI and DTI studies showed volumetric and structural connectivity alterations in the hippocampus and, to a lesser extent, in the amygdala. In contrast, more heterogeneous results were reported by fMRI studies, which, though, described functional abnormalities in the cingulate cortex, thalamus and several other brain areas.

LIMITATIONS

The small sample size and the heterogeneity in patients' samples, processing and study design limit the conclusion of the present review.

CONCLUSIONS

Studies employing scans at ultra-high magnetic field may provide a useful contribution to the mixed body of literature on MDD. This preliminary but promising evidence confirms the importance of performing ultra-high field MRI investigations in order to detect and better characterize subtle brain abnormalities in MDD.

Individualized cognitive neuroscience needs 7T: Comparing numerosity maps at 3T and 7T MRI

The field of cognitive neuroscience is weighing evidence about whether to move from the current standard field strength of 3 Tesla (3T) to ultra-high field (UHF) of 7T and above. The present study contributes to the evidence by comparing a computational cognitive neuroscience paradigm at 3T and 7T. The goal was to evaluate the practical effects, i.e. model predictive power, of field strength on a numerosity task using accessible pre-processing and analysis tools. Previously, using 7T functional magnetic resonance imaging and biologically-inspired analyses, i.e. population receptive field modelling, we discovered topographical organization of numerosity-selective neural populations in human parietal cortex. Here we show that these topographic maps are also detectable at 3T. However, averaging of many more functional runs was required at 3T to reliably reconstruct numerosity maps. On average, one 7T run had about four times the model predictive power of one 3T run. We believe that this amount of scanning would have made the initial discovery of the numerosity maps on 3T highly infeasible in practice. Therefore, we suggest that the higher signal-to-noise ratio and signal sensitivity of UHF MRI is necessary to build mechanistic models of the organization and function of our cognitive abilities in individual participants.

Ultra-high field (7T) functional magnetic resonance imaging in amyotrophic lateral sclerosis: a pilot study

•

Participants with ALS exhibited impaired function between the cortex and cerebellum.

•

The cerebellum is associated with complex motor and cognitive processing tasks.

•

These findings add to the growing number of ALS reports implicating the cerebellum.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease of the central nervous system that results in a progressive loss of motor function and ultimately death. It is critical, yet also challenging, to develop non-invasive biomarkers to identify, localize, measure and/or track biological mechanisms implicated in ALS. Such biomarkers may also provide clues to identify potential molecular targets for future therapeutic trials. Herein we report on a pilot study involving twelve participants with ALS and nine age-matched healthy controls who underwent high-resolution resting state functional magnetic resonance imaging at an ultra-high field of 7 Tesla. A group-level whole-brain analysis revealed a disruption in long-range functional connectivity between the superior sensorimotor cortex (in the precentral gyrus) and bilateral cerebellar lobule VI. Post hoc analyses using atlas-derived left and right cerebellar lobule VI revealed decreased functional connectivity in ALS participants that predominantly mapped to bilateral postcentral and precentral gyri. Cerebellar lobule VI is a transition zone between anterior motor networks and posterior non-motor networks in the cerebellum, and is associated with a wide range of key functions including complex motor and cognitive processing tasks. Our observation of the involvement of cerebellar lobule VI adds to the growing number of studies implicating the cerebellum in ALS. Future avenues of scientific investigation should consider how high-resolution imaging at 7T may be leveraged to visualize differences in functional connectivity disturbances in various genotypes and phenotypes of ALS along the ALS-frontotemporal dementia spectrum.

Ultra-high field magnetic resonance imaging in human epilepsy: A systematic review

RATIONALE

Resective epilepsy surgery is an evidence-based curative treatment option for patients with drug-resistant focal epilepsy. The major preoperative predictor of a good surgical outcome is detection of an epileptogenic lesion by magnetic resonance imaging (MRI). Application of ultra-high field (UHF) MRI, i.e. field strengths ≥ 7 Tesla (T), may increase the sensitivity to detect such a lesion.

METHODS

A keyword search strategy was submitted to Pubmed, EMBASE, Cochrane Database and clinicaltrials.gov to select studies on UHF MRI in patients with epilepsy. Follow-up study selection and data extraction were performed following PRISMA guidelines. We focused on I) diagnostic gain of UHF- over conventional MRI, II) concordance of MRI-detected lesion, seizure onset zone and surgical decision-making, and III) postoperative histopathological diagnosis and seizure outcome.

RESULTS

Sixteen observational cohort studies, all using 7T MRI were included. Diagnostic gain of 7T over conventional MRI ranged from 8% to 67%, with a pooled gain of 31%. Novel techniques to visualize pathological processes in epilepsy and lesion detection are discussed. Seizure freedom was achieved in 73% of operated patients; no seizure outcome comparison was made between 7T MRI positive, 7T negative and 3T positive patients. 7T could influence surgical decision-making, with high concordance of lesion and seizure onset zone. Focal cortical dysplasia (54%), hippocampal sclerosis (12%) and gliosis (8.1%) were the most frequently diagnosed histopathological entities.

SIGNIFICANCE

UHF MRI increases, yet variably, the sensitivity to detect an epileptogenic lesion, showing potential for use in clinical practice. It remains to be established whether this results in improved seizure outcome after surgical treatment. Prospective studies with larger cohorts of epilepsy patients, uniform scan and sequence protocols, and innovative post-processing technology are equally important as further increasing field strengths. Besides technical ameliorations, improved correlation of imaging features with clinical semiology, histopathology and clinical outcome has to be established.

MRI phase offset correction method impacts quantitative susceptibility mapping

Individual channel ultra-high field (7T) phase images have to be phase offset corrected prior to the mapping of magnetic susceptibility of tissue. Whilst numerous methods have been proposed for gradient recalled echo MRI phase offset correction, it remains unclear how they affect quantitative magnetic susceptibility values derived from phase images. Methods already proposed either employ a single or multiple echo time MRI data. In terms of the latter, offsets can be derived using an ultra-short echo time acquisition, or by estimating the offset based on two echo points with the assumption of linear phase evolution with echo time. Our evaluation involved 32 channel multi-echo time 7T GRE (Gradient Recalled Echo) and ultra-short echo time PETRA (Pointwise Encoding Time Reduction with Radial Acquisition) MRI data collected for a susceptibility phantom and three human brains. The combined phase images generated using four established offset correction methods (two single and two multiple echo time) were analysed, followed by an assessment of quantitative susceptibility values obtained for a phantom and human brains. The effectiveness of each method in removing the offsets was shown to reduce with increased echo time, decreased signal intensity and reduced overlap in coil sensitivity profiles. Quantitative susceptibility values and how they change with echo time were found to be method specific. Phase offset correction methods based on single echo time data have a tendency to produce more accurate and less noisy quantitative susceptibility maps in comparison with methods employing multiple echo time data.

Reduced striatal GABA in unmedicated children with ADHD at 7T

Attention-deficit hyperactive disorder (ADHD) is characterized by inattention and increased impulsive and hypermotoric behaviors.Despite the high prevalence and impact of ADHD, little is known about the underlying neurophysiology of ADHD. The main inhibitory and excitatory neurotransmitters γ-aminobutyric acid (GABA) and glutamate are receiving increased attention in ADHD and can be measured using Magnetic Resonance Spectroscopy (MRS). However, MRS studies in ADHD are limited. We measured GABA and glutamate in young unmedicated participants, utilizing high magnetic field strength. Fifty unmedicated children (26 with ADHD, 24 controls) aged 5-9 years completed MRS at 7T and behavioral testing. GABA and glutamate were measured in dorsolateral prefrontal cortex (DLPFC), anterior cingulate cortex (ACC), premotor cortex (PMC), and striatum, and estimated using LCModel. Children with ADHD showed poorer inhibitory control and significantly reduced GABA/Cr in the striatum, but not in ACC, DLPFC, or PMC regions. There were no significant group differences for Glu/Cr levels, or correlations with behavioral manifestations of ADHD. The primary finding of this study is a reduction of striatal GABA levels in unmedicated children with ADHD at 7T. These findings provide guidance for future studies or interventions. Reduced striatal GABA may be a marker for specific GABA-related treatment for ADHD.

A probabilistic atlas of finger dominance in the primary somatosensory cortex

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.

A Data-Driven Multi-scale Technique for fMRI Mapping of the Human Somatosensory Cortex

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.

High resolution time-of-flight MR-angiography at 7 T exploiting VERSE saturation, compressed sensing and segmentation

BACKGROUND

3D Time-of-Flight (TOF) MR-angiography (MRA) substantially benefits from ultra-high magnetic field strengths (≥7 T) due to increased Signal-to-Noise ratio and improved contrast. However, high-resolution TOF-MRA usually requires long acquisition times. In addition, specific absorption rate constraints limit the choice of optimal pulse sequence parameters, especially if venous saturation is employed.

PURPOSE

To implement and evaluate an arterial TOF-MRA for accelerated high-resolution angiography at ultra-high magnetic field strength.

FIELD STRENGTHS/SEQUENCE

7 T modified gradient-echo TOF sequence including venous saturation using Variable-Rate Selective Excitation (VERSE), Compressed Sensing (CS) and sparse application of saturation pulses, called segmentation, were included for acceleration.

ASSESSMENT

To analyze the acceleration techniques all volunteers were examined with the same protocols. CS with different sampling patterns and regularization factors as well as segmentation were applied for acceleration. For comparison, conventional acceleration techniques were applied (GRAPPA PAT 3 and Partial Fourier (6/8 in slice/phase encoding)). Images were co-registered and 40 mm thick transversal maximum intensity projections were created to calculate the relative number of vessels. To analyze the visibility of small vessels, the lenticulostriate arteries (LSA) were examined. This was done via multiscale vessel enhancement filtering in a VOI and quantification via Fiji ImageJ as well as qualitatively evaluation by two radiologists. Additionally, the venous/arterial vessel-to-background ratios (vVBR/aVBR) were calculated for chosen protocols.

RESULTS

For the acceleration of a high resolution TOF-MRA (0.31 mm isotropic), under-sampling of 9.6 showed aliasing artifacts, whereas 7.2 showed no aliasing. The regularization factor R had a strong impact on the image quality according to smoothing (R = 0.01 to R = 0.005) and noise (R = 0.0005 to R = 0.00005). With the alternating sampling patterns it was shown that the k-space center should not be under-sampled too much. Additionally segmentation could be verified to be feasible for stronger acceleration with sufficient venous suppression.

CONCLUSION

The combination of several independent techniques (VERSE, CS with acceleration factor 7.2, R = 0.001, Poisson disc radius of 80%, 3 segments) enables the application of high-resolution (0.31 mm isotropic) TOF-MRA with venous saturation at 7 T in clinical time settings (TA ≈ 5 min) and within the SAR limits.

Multi-slice passband bSSFP for human and rodent fMRI at ultra-high field

Balanced steady-state free precession (bSSFP) can be used as an alternative to gradient-echo (GE) EPI for BOLD functional MRI when image distortions and signal drop-outs are severe such as at ultra-high field. However, 3D-bSSFP acquisitions have distinct drawbacks on either human or animal MR systems. On clinical scanners, 3D imaging is suboptimal for localized fMRI applications. It can also display distortions when acceleration methods such as spiral read-outs are used, and, compared to multi-slice acquisitions, suffers from increased sensitivity to motion or physiological noise which further results in blurring. On pre-clinical systems, 3D acquisitions have low temporal resolution due to limited acceleration options, while single slice often results in insufficient coverage. The aim of the present study was to implement a multi-slice bSSFP acquisition with Cartesian read-out to obtain non-distorted BOLD fMRI activation maps in the human and rat brain at ultra-high field. We show that, when using a new pseudo-steady-state, the bSSFP signal characteristics are preserved. In the human brain at 7 T, we demonstrate that both task- and resting-state fMRI can be performed with multi-slice bSSFP, with a temporal SNR that matches that of 3D-bSSFP, resulting in - at least - equal performance. In the rat brain at 14 T, we show that the multi-slice bSSFP protocol has similar sensitivity to gradient-echo EPI for task fMRI, while benefitting from much reduced distortions and drop-outs. The advantages of passband bSSFP at 14 T in comparison with GE-EPI are expected to be even more marked for mouse brain.

Investigating the variability of cardiac pulse artifacts across heartbeats in simultaneous EEG-fMRI recordings: A 7T study

Electroencephalography (EEG) recordings performed in magnetic resonance imaging (MRI) scanners are affected by complex artifacts caused by heart function, often termed pulse artifacts (PAs). PAs can strongly compromise EEG data quality, and remain an open problem for EEG-fMRI. This study investigated the properties and mechanisms of PA variability across heartbeats, which has remained largely unaddressed to date, and evaluated its impact on PA correction approaches. Simultaneous EEG-fMRI was performed at 7T on healthy participants at rest or under visual stimulation, with concurrent recordings of breathing and cardiac activity. PA variability was found to contribute to EEG variance with more than 500 μV² at 7T, which extrapolates to 92 μV² at 3T. Clustering analyses revealed that PA variability not only is linked to variations in head position/orientation, as previously hypothesized, but also, and more importantly, to the respiratory cycle and to heart rate fluctuations. The latter mechanisms are associated to short-timescale variability (even across consecutive heartbeats), and their importance varied across EEG channels. In light of this PA variability, three PA correction techniques were compared: average artifact subtraction (AAS), optimal basis sets (OBS), and an approach based on K-means clustering. All methods allowed the recovery of visual evoked potentials from the EEG data; nonetheless, OBS and K-means tended to outperform AAS, likely due to the inability of the latter in modeling short-timescale variability. Altogether, these results offer novel insights into the dynamics and underlying mechanisms of the pulse artifact, with important consequences for its correction, relevant to most EEG-fMRI applications.

Pros and cons of ultra-high-field MRI/MRS for human application

Magnetic resonance imaging and spectroscopic techniques are widely used in humans both for clinical diagnostic applications and in basic research areas such as cognitive neuroimaging. In recent years, new human MR systems have become available operating at static magnetic fields of 7 T or higher (≥300 MHz proton frequency). Imaging human-sized objects at such high frequencies presents several challenges including non-uniform radiofrequency fields, enhanced susceptibility artifacts, and higher radiofrequency energy deposition in the tissue. On the other side of the scale are gains in signal-to-noise or contrast-to-noise ratio that allow finer structures to be visualized and smaller physiological effects to be detected. This review presents an overview of some of the latest methodological developments in human ultra-high field MRI/MRS as well as associated clinical and scientific applications. Emphasis is given to techniques that particularly benefit from the changing physical characteristics at high magnetic fields, including susceptibility-weighted imaging and phase-contrast techniques, imaging with X-nuclei, MR spectroscopy, CEST imaging, as well as functional MRI. In addition, more general methodological developments such as parallel transmission and motion correction will be discussed that are required to leverage the full potential of higher magnetic fields, and an overview of relevant physiological considerations of human high magnetic field exposure is provided.

Depth-dependent intracortical myelin organization in the living human brain determined by in vivo ultra-high field magnetic resonance imaging

BACKGROUND

Intracortical myelin is a key determinant of neuronal synchrony and plasticity that underpin optimal brain function. Magnetic resonance imaging (MRI) facilitates the examination of intracortical myelin but presents with methodological challenges. Here we describe a whole-brain approach for the in vivo investigation of intracortical myelin in the human brain using ultra-high field MRI.

METHODS

Twenty-five healthy adults were imaged in a 7 Tesla MRI scanner using diffusion-weighted imaging and a T₁-weighted sequence optimized for intracortical myelin contrast. Using an automated pipeline, T₁ values were extracted at 20 depth-levels from each of 148 cortical regions. In each cortical region, T₁ values were used to infer myelin concentration and to construct a non-linearity index as a measure the spatial distribution of myelin across the cortical ribbon. The relationship of myelin concentration and the non-linearity index with other neuroanatomical properties were investigated. Five patients with multiple sclerosis were also assessed using the same protocol as positive controls.

RESULTS

Intracortical T₁ values decreased between the outer brain surface and the gray-white matter boundary following a slope that showed a slight leveling between 50% and 75% of cortical depth. Higher-order regions in the prefrontal, cingulate and insular cortices, displayed higher non-linearity indices than sensorimotor regions. Across all regions, there was a positive association between T₁ values and non-linearity indices (P < 10^-5). Both T₁ values (P < 10^-5) and non-linearity indices (P < 10^-15) were associated with cortical thickness. Higher myelin concentration but only in the deepest cortical levels was associated with increased subcortical fractional anisotropy (P = 0.05).

CONCLUSIONS

We demonstrate the usefulness of an automatic, whole-brain method to perform depth-dependent examination of intracortical myelin organization. The extracted metrics, T₁ values and the non-linearity index, have characteristic patterns across cortical regions, and are associated with thickness and underlying white matter microstructure.

Visualizing the Human Subcortex Using Ultra-high Field Magnetic Resonance Imaging

With the recent increased availability of ultra-high field (UHF) magnetic resonance imaging (MRI), substantial progress has been made in visualizing the human brain, which can now be done in extraordinary detail. This review provides an extensive overview of the use of UHF MRI in visualizing the human subcortex for both healthy and patient populations. The high inter-subject variability in size and location of subcortical structures limits the usability of atlases in the midbrain. Fortunately, the combined results of this review indicate that a large number of subcortical areas can be visualized in individual space using UHF MRI. Current limitations and potential solutions of UHF MRI for visualizing the subcortex are also discussed.

Signal compartments in ultra-high field multi-echo gradient echo MRI reflect underlying tissue microstructure in the brain

Gradient recalled echo magnetic resonance imaging (GRE-MRI) at ultra-high field holds great promise for new contrast mechanisms and delineation of putative tissue compartments that contribute to the multi-echo GRE-MRI signal may aid structural characterization. Several studies have adopted the three water-pool compartment model to study white matter brain regions, associating individual compartments with myelin, axonal and extracellular water. However, the number and identifiability of GRE-MRI signal compartments has not been fully explored. We undertook this task for human brain imaging data. Multiple echo time GRE-MRI data were acquired in five healthy participants, specific anatomical structures were segmented in each dataset (substantia nigra, caudate, insula, putamen, thalamus, fornix, internal capsule, corpus callosum and cerebrospinal fluid), and the signal fitted with models comprising one to six signal compartments using a complex-valued plane wave formulation. Information criteria and cluster analysis methods were used to ascertain the number of distinct compartments within the signal from each structure and to determine their respective frequency shifts. We identified five principal signal compartments with different relative contributions to each structure's signal. Voxel-based maps of the volume fraction of each of these compartments were generated and demonstrated spatial correlation with brain anatomy.

Pulse sequences and parallel imaging for high spatiotemporal resolution MRI at ultra-high field

The SNR and CNR benefits of ultra-high field (UHF) have helped push the envelope of achievable spatial resolution in MRI. For applications based on susceptibility contrast where there is a large CNR gain, high quality sub-millimeter resolution imaging is now being routinely performed, particularly in fMRI and phase imaging/QSM. This has enabled the study of structure and function of very fine-scale structures in the brain. UHF has also helped push the spatial resolution of many other MRI applications as will be outlined in this review. However, this push in resolution comes at a cost of a large encoding burden leading to very lengthy scans. Developments in parallel imaging with controlled aliasing and the move away from 2D slice-by-slice imaging to much more SNR-efficient simultaneous multi-slice (SMS) and 3D acquisitions have helped address this issue. In particular, these developments have revolutionized the efficiency of UHF MRI to enable high spatiotemporal resolution imaging at an order of magnitude faster acquisition. In addition to describing the main approaches to these techniques, this review will also outline important key practical considerations in using these methods in practice. Furthermore, new RF pulse design to tackle the B1+ and SAR issues of UHF and the increased SAR and power requirement of SMS RF pulses will also be touched upon. Finally, an outlook into new developments of smart encoding in more dimensions, particularly through using better temporal/across-contrast encoding and reconstruction will be described. Just as controlled aliasing fully exploits spatial encoding in parallel imaging to provide large multiplicative gains in accelerations, the complimentary use of these new approaches in temporal and across-contrast encoding are expected to provide exciting opportunities for further large gains in efficiency to further push the spatiotemporal resolution of MRI.

Surface-based characteristics of the cerebellar cortex visualized with ultra-high field MRI

Although having a relatively homogeneous cytoarchitectonic organization, the cerebellar cortex is a heterogeneous region characterized by different amounts of myelin, iron and protein expression profiles. In this study, we used quantitative T₁ and T₂* mapping at ultra-high field (7T) MRI to investigate the tissue characteristics of the cerebellar gray matter surface and its layers. Detailed subject-specific surfaces were generated at three different cortical depths and averaged across subjects to create averaged T₁- and T₂*-maps on the cerebellar surface. T₁ surfaces showed an alternation of lower and higher T₁ values when going from the median to the lateral part of the cerebellar hemispheres. In addition, longer T₁ values were observed in the more superficial gray matter layers. T₂*-maps showed a similar longitudinal pattern, but no change related to the cortical depths. These patterns are possibly due to variations in the level of myelination, iron and zebrin protein expression.

Quantification of perivascular spaces at 7T: A potential MRI biomarker for epilepsy

PURPOSE

7T (7T) magnetic resonance imaging (MRI) facilitates the visualization of the brain with resolution and contrast beyond what is available at conventional clinical field strengths, enabling improved detection and quantification of small structural features such as perivascular spaces (PVSs). The distribution of PVSs, detected in vivo at 7T, may act as a biomarker for the effects of epilepsy. In this work, we systematically quantify the PVSs in the brains of epilepsy patients and compare them to healthy controls.

METHODS

T₂-weighted turbo spin echo images were obtained at 7T on 21 epilepsy patients and 17 healthy controls. For all subjects, PVSs were manually marked on Osirix image analysis software. Marked PVSs with diameter≥0.5mm were then mapped by hemisphere and lobe. The asymmetry index (AI) was calculated for each region and the maximum asymmetry index (|AI_max|) was reported for each subject. The asymmetry in epilepsy subjects was compared to that of controls, and the region with highest asymmetry was compared to the suspected seizure onset zone.

RESULTS

There was a significant difference between the |AI_max| in epilepsy subjects and in controls (p=0.016). In 72% of patients, the region or lobe of the brain showing maximum PVS asymmetry was the same as the region containing the suspected seizure onset zone.

CONCLUSION

These findings suggest that epilepsy may be associated with significantly asymmetric distribution of PVSs in the brain. Furthermore, the region of maximal asymmetry of the PVSs may help provide localization or confirmation of the seizure onset zone.

Laminar fMRI: What can the time domain tell us?

The rapid developments in functional MRI (fMRI) acquisition methods and hardware technologies in recent years, particularly at high field (≥7 T), have enabled unparalleled visualization of functional detail at a laminar or columnar level, bringing fMRI close to the intrinsic resolution of brain function. These advances highlight the potential of high resolution fMRI to be a valuable tool to study the fundamental processing performed in cortical micro-circuits, and their interactions such as feedforward and feedback processes. Notably, because fMRI measures neuronal activity via hemodynamics, the ultimate resolution it affords depends on the spatial specificity of hemodynamics to neuronal activity at a detailed spatial scale, and by the evolution of this specificity over time. Several laminar (≤1 mm spatial resolution) fMRI studies have examined spatial characteristics of the measured hemodynamic signals across cortical depth, in light of understanding or improving the spatial specificity of laminar fMRI. Few studies have examined temporal features of the hemodynamic response across cortical depth. Temporal features of the hemodynamic response offer an additional means to improve the specificity of fMRI, and could help target neuronal processes and neurovascular coupling relationships across laminae, for example by differences in the onset times of the response across cortical depth. In this review, we discuss factors that affect the timing of neuronal and hemodynamic responses across laminae, touching on the neuronal laminar organization, and focusing on the laminar vascular organization. We provide an overview of hemodynamics across the cortical vascular tree based on optical imaging studies, and review temporal aspects of hemodynamics that have been examined across cortical depth in high spatiotemporal resolution fMRI studies. Last, we discuss the limits and potential of high spatiotemporal resolution fMRI to study laminar neurovascular coupling and neuronal processes.

A purpose-built neck coil for black-blood DANTE-prepared carotid artery imaging at 7T

Atherosclerotic plaques in the bifurcation of the carotid arteries can pose a significant health risk due to possible plaque rupture and subsequent stroke. The assessment of plaques, and evaluation of the risk they pose, can be performed with Black-Blood (BB) vessel wall magnetic resonance imaging. However, resolution at standard clinical field strengths (up to 3T) is limited, hampering reliable assessment and diagnosis. The aim of this study was to investigate the benefits of 7T MRI using a BB application that has been successful at clinical field strengths. Therefore, for BB imaging, each sequence was preceded with 'Delay Alternating with Nutation for Tailored Excitation' (DANTE) preparation pulses for blood signal suppression. A coil comprising a 4-channel Tx array was designed and built to provide the required excitation coverage for the DANTE train; and a 4-channel Rx array was constructed to target the carotid bifurcation. Human and phantom results showed satisfactory blood suppression and comparable SNR and CNR to 3T, therefore demonstrating the feasibility of the application at 7T. However, the imposed SAR restrictions led to long scan times and subsequent motion artifacts. Thus, more accurate local SAR supervision schemes are required which could lead to a further improvement of BB DANTE vessel wall imaging at 7T.

Diffusion MRI of the human brain at ultra-high field (UHF): A review

The continued drive towards MRI scanners operating at increasingly higher main magnetic fields is primarily motivated by the maxim that more teslas mean more signal and lead to better images. This promise of increased signal, which cannot easily be achieved in other ways, encourages efforts to overcome the inextricable technical challenges which accompany this endeavor. Unlike for many applications, however, diffusion imaging is not currently able to directly reap these potential signal gains - at the time of writing it seems fair to say that, for matched gradient and RF hardware, the majority of diffusion images acquired at 7T, while comparable in quality to those achievable at 3T, do not demonstrate a clear advantage over what can be obtained at lower field. This does not mean that diffusion imaging at UHF is not a worthwhile pursuit - but more a reflection of the fact that the associated challenges are manifold - and converting the potential of higher field strengths into 'better' diffusion imaging is by no means a straightforward task. This article attempts to summarize the specific reasons that make diffusion imaging at UHF more complicated than one might expect, and to highlight the range of developments that have already been made which have enabled diffusion images of excellent quality to be acquired at 7T.

The clinical relevance of distortion correction in presurgical fMRI at 7T

Presurgical planning with fMRI benefits from increased reliability and the possibility to reduce measurement time introduced by using ultra-high field. Echo-planar imaging suffers, however, from geometric distortions which scale with field strength and potentially give rise to clinically significant displacement of functional activation. We evaluate the effectiveness of a dynamic distortion correction (DDC) method based on unmodified single-echo EPI in the context of simulated presurgical planning fMRI at 7T and compare it with static distortion correction (SDC). The extent of distortion in EPI and activation shifts are investigated in a group of eleven patients with a range of neuropathologies who performed a motor task. The consequences of neglecting to correct images for susceptibility-induced distortions are assessed in a clinical context. It was possible to generate time series of EPI-based field maps which were free of artifacts in the eloquent brain areas relevant to presurgical fMRI, despite the presence of signal dropouts caused by pathologies and post-operative sites. Distortions of up to 5.1mm were observed in the primary motor cortex in raw EPI. These were accurately corrected with DDC and slightly less accurately with SDC. The dynamic nature of distortions in UHF clinical fMRI was demonstrated via investigation of temporal variation in voxel shift maps, confirming the potential inadequacy of SDC based on a single reference field map, particularly in the vicinity of pathologies or in the presence of motion. In two patients, the distortion correction was potentially clinically significant in that it might have affected the localization or interpretation of activation and could thereby have influenced the treatment plan. Distortion correction is shown to be effective and clinically relevant in presurgical planning at 7T.

Fast triple-spin-echo Dixon (FTSED) sequence for water and fat imaging

A number of 'Dixon' techniques based on fast spin echo (FSE) sequence have been proposed and successfully used in many branches of medicine. Some require only one scan, but most of them need multiple scans and long scan times. This article describes a new fast triple-spin-echo Dixon (FTSED) technique suitable for ultra-high field MRI, in which three specific time shifts are introduced in the echo train; thus, three images with defined water-fat phase-differences (0, π, 2π) are encoded in the phase of the acquired images without extreme restrictions upon the echo duration. The water and fat images are then calculated by iterative least-squares estimation method. The sequence was successfully implemented at a 9.4T ultra-high field MRI system and tested on a phantom and a rat.

A method for the dynamic correction of B0-related distortions in single-echo EPI at 7T

We propose a method to calculate field maps from the phase of each EPI in an fMRI time series. These field maps can be used to correct the corresponding magnitude images for distortion caused by inhomogeneity in the static magnetic field. In contrast to conventional static distortion correction, in which one ‘snapshot’ field map is applied to all subsequent fMRI time points, our method also captures dynamic changes to B₀ which arise due to motion and respiration. The approach is based on the assumption that the non-B₀-related contribution to the phase measured by each radio-frequency coil, which is dominated by the coil sensitivity, is stable over time and can therefore be removed to yield a field map from EPI.

Our solution addresses imaging with multi-channel coils at ultra-high field (7 T), where phase offsets vary rapidly in space, phase processing is non-trivial and distortions are comparatively large. We propose using dual-echo gradient echo reference scan for the phase offset calculation, which yields estimates with high signal-to-noise ratio. An extrapolation method is proposed which yields reliable estimates for phase offsets even where motion is large and a tailored phase unwrapping procedure for EPI is suggested which gives robust results in regions with disconnected tissue or strong signal decay.

Phase offsets are shown to be stable during long measurements (40 min) and for large head motions. The dynamic distortion correction proposed here is found to work accurately in the presence of large motion (up to 8.1°), whereas a conventional method based on single field map fails to correct or even introduces distortions (up to 11.2 mm). Finally, we show that dynamic unwarping increases the temporal stability of EPI in the presence of motion.

Our approach can be applied to any EPI measurements without the need for sequence modification.

Pulsed arterial spin labelling at ultra-high field with a B 1 (+) -optimised adiabatic labelling pulse

OBJECTIVE

Arterial spin labelling (ASL) techniques benefit from the increased signal-to-noise ratio and the longer T 1 relaxation times available at ultra-high field. Previous pulsed ASL studies at 7 T concentrated on the superior regions of the brain because of the larger transmit radiofrequency inhomogeneity experienced at ultra-high field that hinders an adequate inversion of the blood bolus when labelling in the neck. Recently, researchers have proposed to overcome this problem with either the use of dielectric pads, through dedicated transmit labelling coils, or special adiabatic inversion pulses.

MATERIALS AND METHODS

We investigate the performance of an optimised time-resampled frequency-offset corrected inversion (TR-FOCI) pulse designed to cause inversion at much lower peak B 1 (+) . In combination with a PICORE labelling, the perfusion signal obtained with this pulse is compared against that obtained with a FOCI pulse, with and without dielectric pads.

RESULTS

Mean grey matter perfusion with the TR-FOCI was 52.5 ± 10.3 mL/100 g/min, being significantly higher than the 34.6 ± 2.6 mL/100 g/min obtained with the FOCI pulse. No significant effect of the dielectric pads was observed.

CONCLUSION

The usage of the B 1 (+) -optimised TR-FOCI pulse results in a significantly higher perfusion signal. PICORE-ASL is feasible at ultra-high field with no changes to operating conditions.

Improving the clinical potential of ultra-high field fMRI using a model-free analysis method based on response consistency

Objective

To develop an analysis method that is sensitive to non-model-conform responses often encountered in ultra-high field presurgical planning fMRI. Using the consistency of time courses over a number of experiment repetitions, it should exclude low quality runs and generate activation maps that reflect the reliability of responses.

Materials and methods

7 T fMRI data were acquired from six healthy volunteers: three performing purely motor tasks and three a visuomotor task. These were analysed with the proposed approach (UNBIASED) and the GLM.

Results

UNBIASED results were generally less affected by false positive results than the GLM. Runs that were identified as being of low quality were confirmed to contain little or no activation. In two cases, regions were identified as activated in UNBIASED but not GLM results. Signal changes in these areas were time-locked to the task, but were delayed or transient.

Conclusion

UNBIASED is shown to be a reliable means of identifying consistent task-related signal changes regardless of response timing. In presurgical planning, UNBIASED could be used to rapidly generate reliable maps of the consistency with which eloquent brain regions are activated without recourse to task timing and despite modified hemodynamics.

Electronic supplementary material

The online version of this article (doi:10.1007/s10334-016-0533-8) contains supplementary material, which is available to authorized users.

Presurgical brain mapping in epilepsy using simultaneous EEG and functional MRI at ultra-high field: feasibility and first results

OBJECTIVES

The aim of this study was to demonstrate that eloquent cortex and epileptic-related hemodynamic changes can be safely and reliably detected using simultaneous electroencephalography (EEG)-functional magnetic resonance imaging (fMRI) recordings at ultra-high field (UHF) for clinical evaluation of patients with epilepsy.

MATERIALS AND METHODS

Simultaneous EEG-fMRI was acquired at 7 T using an optimized setup in nine patients with lesional epilepsy. According to the localization of the lesion, mapping of eloquent cortex (language and motor) was also performed in two patients.

RESULTS

Despite strong artifacts, efficient correction of intra-MRI EEG could be achieved with optimized artifact removal algorithms, allowing robust identification of interictal epileptiform discharges. Noise-sensitive topography-related analyses and electrical source localization were also performed successfully. Localization of epilepsy-related hemodynamic changes compatible with the lesion were detected in three patients and concordant with findings obtained at 3 T. Local loss of signal in specific regions, essentially due to B 1 inhomogeneities were found to depend on the geometric arrangement of EEG leads over the cap.

CONCLUSION

These results demonstrate that presurgical mapping of epileptic networks and eloquent cortex is both safe and feasible at UHF, with the benefits of greater spatial resolution and higher blood-oxygenation-level-dependent sensitivity compared with the more traditional field strength of 3 T.

(31)P CSI of the human brain in healthy subjects and tumor patients at 9.4 T with a three-layered multi-nuclear coil: initial results

OBJECTIVE

Investigation of the feasibility and performance of phosphorus ((31)P) magnetic resonance spectroscopic imaging (MRSI) at 9.4 T with a three-layered phosphorus/proton coil in human normal brain tissue and tumor.

MATERIALS AND METHODS

A multi-channel (31)P coil was designed to enable MRSI of the entire human brain. The performance of the coil was evaluated by means of electromagnetic field simulations and actual measurements. A 3D chemical shift imaging approach with a variable repetition time and flip angle was used to increase the achievable signal-to-noise ratio of the acquired (31)P spectra. The impact of the resulting k-space modulation was investigated by simulations. Three tumor patients and three healthy volunteers were scanned and differences between spectra from healthy and cancerous tissue were evaluated qualitatively.

RESULTS

The high sensitivity provided by the 27-channel (31)P coil allowed acquiring CSI data in 22 min with a nominal voxel size of 15 × 15 × 15 mm(3). Shimming and anatomical localization could be performed with the integrated four-channel proton dipole array. The amplitudes of the phosphodiesters and phosphoethanolamine appeared reduced in tumorous tissue for all three patients. A neutral or slightly alkaline pH was measured within the brain lesions.

CONCLUSION

These initial results demonstrate that (31)P 3D CSI is feasible at 9.4 T and could be performed successfully in healthy subjects and tumor patients in under 30 min.