A comparison of optimization algorithms for localized in vivo B0 shimming

Fourier-based decomposition for simultaneous 2-voxel MRS acquisition with 2SPECIAL

PURPOSE

To simultaneously acquire spectroscopic signals from two MRS voxels using a multi-banded 2 spin-echo, full-intensity acquired localized (2SPECIAL) sequence, and to decompose the signal to their respective regions by a novel voxel-GRAPPA (vGRAPPA) decomposition approach for in vivo brain applications at 7 T.

METHODS

A wideband, uniform rate, smooth truncation (WURST) multi-banded pulse was incorporated into SPECIAL to implement 2SPECIAL for simultaneous multi-voxel spectroscopy (sMVS). To decompose the acquired data, the voxel-GRAPPA decomposition algorithm is introduced, and its performance is compared to the SENSE-based decomposition. Furthermore, the limitations of two-voxel excitation concerning the multi-banded adiabatic inversion pulse, as well as of the combined B₀ shim and B₁ ⁺ adjustments, are evaluated.

RESULTS

It was successfully shown that the 2SPECIAL sequence enables sMVS without a significant loss in SNR while reducing the total scan time by 21.6% compared to two consecutive acquisitions. The proposed voxel-GRAPPA algorithm properly reassigns the signal components to their respective origin region and shows no significant differences to the well-established SENSE-based algorithm in terms of leakage (both <10%) or Cramér-Rao lower bounds (CRLB) for in vivo applications, while not requiring the acquisition of additional sensitivity maps and thus decreasing motion sensitivity.

CONCLUSION

The use of 2SPECIAL in combination with the novel voxel-GRAPPA decomposition technique allows a substantial reduction of measurement time compared to the consecutive acquisition of two single voxels without a significant decrease in spectral quality or metabolite quantification accuracy and thus provides a new option for multiple-voxel applications.

Simultaneous detection of metabolite concentration changes, water BOLD signal and pH changes during visual stimulation in the human brain at 9.4T

This study presents a method to directly link metabolite concentration changes and BOLD response in the human brain during visual stimulation by measuring the water and metabolite signals simultaneously. Therefore, the metabolite-cycling (MC) non-water suppressed semiLASER localization technique was optimized for functional ¹H MRS in the human brain at 9.4 T. Data of 13 volunteers were acquired during a 26:40 min visual stimulation block-design paradigm. Activation-induced BOLD signal was observed in the MC water signal as well as in the NAA-CH₃ and tCr-CH₃ singlets. During stimulation, glutamate concentration increased 2.3 ± 2.0% to a new steady-state, while a continuous increase over the whole stimulation period could be observed in lactate with a mean increase of 35.6 ± 23.1%. These increases of Lac and Glu during brain activation confirm previous findings reported in literature. A positive correlation of the MC water BOLD signal with glutamate and lactate concentration changes was found. In addition, a pH decrease calculated from a change in the ratio of PCr to Cr was observed during brain activation, particularly at the onset of the stimulation.

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.

Shim optimization with region of interest-specific Tikhonov regularization: Application to second-order slice-wise shimming of the brain

PURPOSE

Slice-wise shimming can improve field homogeneity, but suffers from large noise propagation in the shim calculation. Here, we propose a robust shim current optimization for higher-order dynamic shim updating, based on Tikhonov regularization with a variable regularization parameter, λ . THEORY AND METHODS: λ was selected for each slice separately in a fully automatic procedure based on a combination of boundary constraints and an L-curve search algorithm. Shimming performance was evaluated for second order slice-wise shimming of the brain at 7T, by simulation on a database of field maps from 143 subjects, and by direct measurements in 8 subjects.

RESULTS

Simulations yielded on average 36% reduction in the shim current norm for just 0.4 Hz increase in residual field SD as compared to unconstrained unregularized optimization. In vivo results yielded on average 34.0 Hz residual field SD as compared to 34.3 Hz with a constrained unregularized optimization, while simultaneously reducing the shim current norm to 2.8 A from 3.9 A. The proposed regularization also reduced the average step in the shim current between slices.

CONCLUSION

Slice-wise variable Tikhonov regularization yielded reduced current norm and current steps to a negligible cost in field inhomogeneity. The method holds promise to increase the robustness, and thereby the utility, of higher-order shim updating.

Assessment of measurement precision in single-voxel spectroscopy at 7 T: Toward minimal detectable changes of metabolite concentrations in the human brain in vivo

PURPOSE

To introduce a study design and statistical analysis framework to assess the repeatability, reproducibility, and minimal detectable changes (MDCs) of metabolite concentrations determined by in vivo MRS.

METHODS

An unbalanced nested study design was chosen to acquire in vivo MRS data within different repeatability and reproducibility scenarios. A spin-echo, full-intensity acquired localized (SPECIAL) sequence was employed at 7 T utlizing three different inversion pulses: a hyperbolic secant (HS), a gradient offset independent adiabaticity (GOIA), and a wideband, uniform rate, smooth truncation (WURST) pulse. Metabolite concentrations, Cramér-Rao lower bounds (CRLBs) and coefficients of variation (CVs) were calculated. Both Bland-Altman analysis and a restricted maximum-likelihood estimation (REML) analysis were performed to estimate the different variance contributions of the repeatability and reproducibility of the measured concentration. A Bland-Altmann analysis of the spectral shape was performed to assess the variance of the spectral shape, independent of quantification model influences.

RESULTS

For the used setup, minimal detectable changes of brain metabolite concentrations were found to be between 0.40 µmol/g and 2.23 µmol/g. CRLBs account for only 16 % to 74 % of the total variance of the metabolite concentrations. The application of gradient-modulated inversion pulses in SPECIAL led to slightly improved repeatability, but overall reproducibility appeared to be limited by differences in positioning, calibration, and other day-to-day variations throughout different sessions.

CONCLUSION

A framework is introduced to estimate the precision of metabolite concentrations obtained by MRS in vivo, and the minimal detectable changes for 13 metabolite concentrations measured at 7 T using SPECIAL are obtained.

7 Tesla and Beyond: Advanced Methods and Clinical Applications in Magnetic Resonance Imaging

ABSTRACT

Ultrahigh magnetic fields offer significantly higher signal-to-noise ratio, and several magnetic resonance applications additionally benefit from a higher contrast-to-noise ratio, with static magnetic field strengths of B0 ≥ 7 T currently being referred to as ultrahigh fields (UHFs). The advantages of UHF can be used to resolve structures more precisely or to visualize physiological/pathophysiological effects that would be difficult or even impossible to detect at lower field strengths. However, with these advantages also come challenges, such as inhomogeneities applying standard radiofrequency excitation techniques, higher energy deposition in the human body, and enhanced B0 field inhomogeneities. The advantages but also the challenges of UHF as well as promising advanced methodological developments and clinical applications that particularly benefit from UHF are discussed in this review article.

Advanced magnetic resonance spectroscopic neuroimaging: Experts' consensus recommendations

Magnetic resonance spectroscopic imaging (MRSI) offers considerable promise for monitoring metabolic alterations associated with disease or injury; however, to date, these methods have not had a significant impact on clinical care, and their use remains largely confined to the research community and a limited number of clinical sites. The MRSI methods currently implemented on clinical MRI instruments have remained essentially unchanged for two decades, with only incremental improvements in sequence implementation. During this time, a number of technological developments have taken place that have already greatly benefited the quality of MRSI measurements within the research community and which promise to bring advanced MRSI studies to the point where the technique becomes a true imaging modality, while making the traditional review of individual spectra a secondary requirement. Furthermore, the increasing use of biomedical MR spectroscopy studies has indicated clinical areas where advanced MRSI methods can provide valuable information for clinical care. In light of this rapidly changing technological environment and growing understanding of the value of MRSI studies for biomedical studies, this article presents a consensus from a group of experts in the field that reviews the state-of-the-art for clinical proton MRSI studies of the human brain, recommends minimal standards for further development of vendor-provided MRSI implementations, and identifies areas which need further technical development.

Ultrahigh-resolution quantitative spinal cord MRI at 9.4T

PURPOSE

To present the results of the first human spinal cord in vivo MRI scans at 9.4T.

METHODS

A human brain coil was used to image the human spinal cord at 9.4T. All anatomical images were acquired with a T₂ *-weighted gradient-echo sequence. A comparison of the influence of four different B₀ shimming routines on the image quality was performed. Intrinsic signal-to-noise-ratio maps were determined using a pseudo-multiple replica approach. Measurements with different echo times were compared and processed to one multiecho data image combination image. Based on the multiecho acquisitions, T₂ *-relaxation time maps were calculated. Algorithmic spinal cord detection and gray matter/white matter segmentation were tested.

RESULTS

An echo time between 9 and 13.8 ms compromised best between gray matter/white matter contrast and image quality. A maximum in-plane resolution of 0.15 × 0.15 mm² was achieved for anatomical images. These images offered excellent image quality and made small structures of the spinal cord visible. The scanner vendor implemented B₀ shimming routine performed best during this work. Intrinsic signal-to-noise-ratio values of between 6600 and 8060 at the upper cervical spinal cord were achieved. Detection and segmentation worked reliably. An average T₂ *-time of 24.88 ms ± 6.68 ms for gray matter and 19.37 ms ± 8.66 ms for white matter was calculated.

CONCLUSION

The proposed human brain coil can be used to image the spinal cord. The maximum in-plane resolution in this work was higher compared with the 7T results from the literature. The 9.4T acquisitions made the small structures of the spinal cord clearly visible.

B0 shimming of the human heart at 7T

PURPOSE

Inhomogeneities of the static magnetic B₀ field are a major limiting factor in cardiac MRI at ultrahigh field (≥ 7T), as they result in signal loss and image distortions. Different magnetic susceptibilities of the myocardium and surrounding tissue in combination with cardiac motion lead to strong spatio-temporal B₀ -field inhomogeneities, and their homogenization (B₀ shimming) is a prerequisite. Limitations of state-of-the-art shimming are described, regional B₀ variations are measured, and a methodology for spherical harmonics shimming of the B₀ field within the human myocardium is proposed.

METHODS

The spatial B₀ -field distribution in the heart was analyzed as well as temporal B₀ -field variations in the myocardium over the cardiac cycle. Different shim region-of-interest selections were compared, and hardware limitations of spherical harmonics B₀ shimming were evaluated by calibration-based B₀ -field modeling. The role of third-order spherical harmonics terms was analyzed as well as potential benefits from cardiac phase-specific shimming.

RESULTS

The strongest B₀ -field inhomogeneities were observed in localized spots within the left-ventricular and right-ventricular myocardium and varied between systolic and diastolic cardiac phases. An anatomy-driven shim region-of-interest selection allowed for improved B₀ -field homogeneity compared with a standard shim region-of-interest cuboid. Third-order spherical harmonics terms were demonstrated to be beneficial for shimming of these myocardial B₀ -field inhomogeneities. Initial results from the in vivo implementation of a potential shim strategy were obtained. Simulated cardiac phase-specific shimming was performed, and a shim term-by-term analysis revealed periodic variations of required currents.

CONCLUSION

Challenges in state-of-the-art B₀ shimming of the human heart at 7 T were described. Cardiac phase-specific shimming strategies were found to be superior to vendor-supplied shimming.

B0 field homogeneity recommendations, specifications, and measurement units for MRI in radiation therapy

PURPOSE

The purpose is: (a) Relate magnetic resonance imaging (MRI) quality recommendations for radiation therapy (RT) to B₀ field homogeneity; (b) Evaluate manufacturer specifications of B₀ homogeneity for 34 commercial whole-body MRI systems based on the MRI quality recommendations and RT application; (c) Measure field homogeneity in five commercial MRI systems and one commercial MRI-Linac used in RT and compare the results with their B₀ homogeneity specifications.

METHODS

Magnetic resonance imaging quality recommendations for spatial integrity, image blurring, fat saturation, and null banding in RT were developed based on the literature. Guaranteed (maximum) and typical B₀ field homogeneity specifications for various diameter spherical volumes (DSVs) were provided by GE, Philips, Siemens, and Canon. For each system, the DSV that conforms to each MRI quality recommendation and anatomical RT application was estimated based on the manufacturer specifications. B₀ field homogeneity was measured on six MRI systems including Philips (1.5 T), Siemens (1.5 and 3 T), and ViewRay MRI (0.35 T) systems using 24 and 35 cm DSV spherical phantoms. Two measurement techniques were used: (a) MRI using phase contrast field mapping to measure peak-to-peak (pk-pk), volume root mean square (VRMS), and standard deviation (SD); and (b) Magnetic resonance (MR) spectroscopy by acquiring a volumetric free induction decay (FID) to measure full width at half maximum (FWHM). The measurements were used to assess: (a) conformance with the manufacturer specifications; and (b) the relationship between the various field homogeneity measurement units. Measurements were made with and without gradient shimming (gradshim) or second-order active shimming. Multiple comparisons, analysis of variance (ANOVA), and Pearson correlations were performed to assess the dependence of pk-pk, VRMS, SD, and FWHM measurements of field homogeneity on shim volume, level of shim, and MRI system.

RESULTS

For a 40 cm DSV, the B₀ homogeneity specifications ranged from 0.35 to 5 ppm (median = 0.75 ppm) VRMS for 1.5 T systems and 0.2 to 1.4 ppm (median = 0.5 ppm) VRMS for 3 T systems. The usable DSVs ranged from 16 to 49 cm (median = 35 cm) based on the image quality recommendations and the manufacturer specifications. There was general compliance between the six measured field homogeneities and manufacturer specifications although signal dephasing was observed in two systems at < 35 cm DSV. The relationships between pk-pk, VRMS, SD, and FWHM varied based on MRI system, shim volume, and quality of shim. However, VRMS and SD measurements were highly correlated.

CONCLUSIONS

The delineation of the diseased lesion from organs at risk is the main priority for RT. Therefore, field homogeneity performance for RT must minimize image blurring and image artifacts (null bands and signal dephasing) while optimizing spatial integrity and fat saturation. Based on the specifications and recommendations for field homogeneity, some MRI systems are not well suited to meet the strict demands of RT particularly for the large imaging volumes used in body MRI. VRMS and SD measurements of B₀ field homogeneity tend to be more stable and sensitive to field inhomogeneities in RT applications than pk-pk and FWHM.

Overdiscrete echo-planar spectroscopic imaging with correlated higher-order phase correction

PURPOSE

To introduce a robust methodology for fast ¹ H MRSI of the brain at 3T with improved SNR and reduced phase-related artifacts.

METHOD

An accelerated acquisition scheme using echo-planar spectroscopic imaging (EPSI) was combined with the overdiscrete reconstruction framework. This approach enables the interleaved acquisition of a water reference scan at each phase encoding step, maximizing its correlation with the water-suppressed measurement. Moreover, a generalized high-order phase correction was incorporated into the reconstruction pipeline. The spatial-temporal phase correction term was estimated from the reference scan and interpolated to high resolution using a polynomial basis. The method was implemented at 3T and validated with phantom and in vivo experiments.

RESULTS

The methodology showed the elimination of spectral artifacts generated by phase disturbances and achieved mean SNR gains in vivo of 3.18 and 1.19 compared to standard reconstructions with corrections performed at nominal and high resolution, respectively. EPSI scans with interleaved water acquisition showed to be robust to system instabilities and potentially to patient motion. Moreover, phase distortions were effectively corrected in a single step, avoiding additional reference measurements and post-processing steps.

CONCLUSION

The overdiscrete reconstruction framework with high-order phase correction allowed to effectively correct for distortions, related to B₀ inhomogeneities, B₀ drift, eddy currents, and system vibrations. Furthermore, the presented reconstruction method, combined with EPSI acquisitions, demonstrated improved measurement stability, substantial SNR enhancement, better spectral linewidth, and effective artifact removal.

Methodological consensus on clinical proton MRS of the brain: Review and recommendations

Proton MRS (1 H MRS) provides noninvasive, quantitative metabolite profiles of tissue and has been shown to aid the clinical management of several brain diseases. Although most modern clinical MR scanners support MRS capabilities, routine use is largely restricted to specialized centers with good access to MR research support. Widespread adoption has been slow for several reasons, and technical challenges toward obtaining reliable good-quality results have been identified as a contributing factor. Considerable progress has been made by the research community to address many of these challenges, and in this paper a consensus is presented on deficiencies in widely available MRS methodology and validated improvements that are currently in routine use at several clinical research institutions. In particular, the localization error for the PRESS localization sequence was found to be unacceptably high at 3 T, and use of the semi-adiabatic localization by adiabatic selective refocusing sequence is a recommended solution. Incorporation of simulated metabolite basis sets into analysis routines is recommended for reliably capturing the full spectral detail available from short TE acquisitions. In addition, the importance of achieving a highly homogenous static magnetic field (B0 ) in the acquisition region is emphasized, and the limitations of current methods and hardware are discussed. Most recommendations require only software improvements, greatly enhancing the capabilities of clinical MRS on existing hardware. Implementation of these recommendations should strengthen current clinical applications and advance progress toward developing and validating new MRS biomarkers for clinical use.

General technical remarks on 1HMRS translational research in 7T

Purpose

The aim of the work was to share the practical experience of preclinical and clinical proton magnetic resonance spectroscopy (¹HMRS) studies conducted using a 7-Tesla magnetic field strength scanner, taking into account the specificity of both settings in the context of translational research.

Material and methods

¹HMRS volunteer studies conducted using a Discovery 950 GE 7T scanner, were carried out with PRESS sequence, and a VOI measuring 2.0 × 2.0 × 2.0 cm³ placed in the white matter at the parietal occipital lobe. Rodent spectra obtained using a 7T Bruker were measured with PRESS, with a VOI 2.0 × 2.0 × 5.5 mm³ placed over the hippocampus.

Results

¹HMRS data from humans and rats show that the brain spectra obtained in the same field are characterised by a similar neurochemical structure and spectral resolution. Spectra obtained from rats demonstrate the following metabolites: NAA, Glu, Gln, Ins, Cho, Cr, PCr, Tau, GABA, Lac, NAAG, and Asp. In turn, spectra from humans allowed estimation of the following metabolites: Ala, NAA, Glu, Gln, Ins, Cho, Cr, PCr, Tau, GABA, Lac, NAAG, and Asp. Signals from Gln, Glu with chemical shift around 2.4 ppm, from Cr, PCr, and GABA at 3 ppm, and signals from Cho and Tau at approximately 3.2 ppm, can be properly separated and estimated both in humans and in rats.

Conclusions

These results are promising in terms of broadening the knowledge of many neurological diseases by inducing them on animal models and then transferring this knowledge to clinical practice. In spite of this, important distinctions in the technical aspects and methodological differences of high-field ¹HMRS in both preclinical and clinical conditions should be taken into account.

Modeling real shim fields for very high degree (and order) B0 shimming of the human brain at 9.4 T

PURPOSE

To describe the process of calibrating a B₀ shim system using high-degree (or high order) spherical harmonic models of the measured shim fields, to provide a method that considers amplitude dependency of these models, and to show the advantage of very high-degree B₀ shimming for whole-brain and single-slice applications at 9.4 Tesla (T).

METHODS

An insert shim with up to fourth and partial fifth/sixth degree (order) spherical harmonics was used with a Siemens 9.4T scanner. Each shim field was measured and modeled as input for the shimming algorithm. Optimal shim currents can therefore be calculated in a single iteration. A range of shim currents was used in the modeling to account for possible amplitude nonlinearities. The modeled shim fields were used to compare different degrees of whole-brain B₀ shimming on healthy subjects.

RESULTS

The ideal shim fields did not correctly shim the subject brains. However, using the modeled shim fields improved the B₀ homogeneity from 55.1 (second degree) to 44.68 Hz (partial fifth/sixth degree) on the whole brains of 9 healthy volunteers, with a total applied current of 0.77 and 6.8 A, respectively.

CONCLUSIONS

The necessity of calibrating the shim system was shown. Better B₀ homogeneity drastically reduces signal dropout and distortions for echo-planar imaging, and significantly improves the linewidths of MR spectroscopy imaging. Magn Reson Med 79:529-540, 2018. © 2017 International Society for Magnetic Resonance in Medicine.