Only-train-once MR fingerprinting for B0 and B1 inhomogeneity correction in quantitative magnetization-transfer contrast.

PURPOSE: To develop a fast, deep-learning approach for quantitative magnetization-transfer contrast (MTC)-MR fingerprinting (MRF) that simultaneously estimates multiple tissue parameters and corrects the effects of B0 and B1 variations. METHODS: An only-train-once recurrent neural network was designed to perform the fast tissue-parameter quantification for a large range of different MRF acquisition schedules. It enabled a dynamic scan-wise linear calibration of the scan parameters using the measured B0 and B1 maps, which allowed accurate, multiple-tissue parameter mapping. MRF images were acquired from 8 healthy volunteers at 3 T. Estimated parameter maps from the MRF images were used to synthesize the MTC reference signal (Zref ) through Bloch equations at multiple saturation power levels. RESULTS: The B0 and B1 errors in MR fingerprints, if not corrected, would impair the tissue quantification and subsequently corrupt the synthesized MTC reference images. Bloch equation-based numerical phantom studies and synthetic MRI analysis demonstrated that the proposed approach could correctly estimate water and semisolid macromolecule parameters, even with severe B0 and B1 inhomogeneities. CONCLUSION: The only-train-once deep-learning framework can improve the reconstruction accuracy of brain-tissue parameter maps and be further combined with any conventional MRF or CEST-MRF method.

B1 inhomogeneity correction of RARE MRI at low SNR: Quantitative in vivo 19 F MRI of mouse neuroinflammation with a cryogenically-cooled transceive surface radiofrequency probe.

PURPOSE: Low SNR in fluorine-19 (19 F) MRI benefits from cryogenically-cooled transceive surface RF probes (CRPs), but strong B1 inhomogeneities hinder quantification. Rapid acquisition with refocused echoes (RARE) is an SNR-efficient method for MRI of neuroinflammation with perfluorinated compounds but lacks an analytical signal intensity equation to retrospectively correct B1 inhomogeneity. Here, a workflow was proposed and validated to correct and quantify 19 F-MR signals from the inflamed mouse brain using a 19 F-CRP. METHODS: In vivo 19 F-MR images were acquired in a neuroinflammation mouse model with a quadrature 19 F-CRP using an imaging setup including 3D-printed components to acquire co-localized anatomical and 19 F images. Model-based corrections were validated on a uniform 19 F phantom and in the neuroinflammatory model. Corrected 19 F-MR images were benchmarked against reference images and overlaid on in vivo 1 H-MR images. Computed concentration uncertainty maps using Monte Carlo simulations served as a measure of performance of the B1 corrections. RESULTS: Our study reports on the first quantitative in vivo 19 F-MR images of an inflamed mouse brain using a 19 F-CRP, including in vivo T1 calculations for 19 F-nanoparticles during pathology and B1 corrections for 19 F-signal quantification. Model-based corrections markedly improved 19 F-signal quantification from errors > 50% to < 10% in a uniform phantom (p < 0.001). Concentration uncertainty maps ex vivo and in vivo yielded uncertainties that were generally < 25%. Monte Carlo simulations prescribed SNR ≥ 10.1 to reduce uncertainties < 10%, and SNR ≥ 4.25 to achieve uncertainties < 25%. CONCLUSION: Our model-based correction method facilitated 19 F signal quantification in the inflamed mouse brain when using the SNR-boosting 19 F-CRP technology, paving the way for future low-SNR 19 F-MRI applications in vivo.

A model-based framework for correcting B1+ inhomogeneity effects in magnetization transfer saturation and inhomogeneous magnetization transfer saturation maps.

PURPOSE: In this work, we propose that Δ B1+ -induced errors in magnetization transfer (MT) saturation (MTsat ) maps can be corrected with use of an R1 and B1+ map and through numerical simulations of the sequence. THEORY AND METHODS: One healthy subject was scanned at 3.0T using a partial quantitative MT protocol to estimate the relationship between observed R1 (R1,obs ) and apparent bound pool size ( M0,appB ) in the brain. MTsat values were simulated for a range of B1+ , R1,obs , and M0,appB . An equation was fit to the simulated MTsat , then a linear relationship between R1,obs and M0,appB was generated. These results were used to generate correction factor maps for the MTsat acquired from single-point data. The proposed correction was compared to an empirical correction factor with different MT-preparation schemes. RESULTS: M0,appB was highly correlated with R1,obs (r > 0.96), permitting the use of R1,obs to estimate M0,appB for B1+ correction. All B1+ corrected MTsat maps displayed a decreased correlation with B1+ compared to uncorrected MTsat and MTsat corrected with an empirical factor in the corpus callosum. There was good agreement between the proposed approach and the empirical correction with radiofrequency saturation at 2 kHz, with larger deviations seen when using saturation pulses further off-resonance and in inhomogeneous (ih) MTsat maps. CONCLUSION: The proposed correction decreases the dependence of MTsat on B1+ inhomogeneities. Furthermore, this flexible framework permits the use of different saturation protocols, making it useful for correcting B1+ inhomogeneities in ihMT.

Specialized computational methods for denoising, B1 correction, and kinetic modeling in hyperpolarized 13 C MR EPSI studies of liver tumors.

PURPOSE: To develop a novel post-processing pipeline for hyperpolarized (HP) 13 C MRSI that integrates tensor denoising and B1+ correction to measure pyruvate-to-lactate conversion rates (kPL ) in patients with liver tumors. METHODS: Seven HP 13 C MR scans of progressing liver tumors were acquired using a custom 13 C surface transmit/receive coil and the echo-planar spectroscopic imaging (EPSI) data analysis included B0 correction, tensor rank truncation, and zero- and first-order phase corrections to recover metabolite signals that would otherwise be obscured by spectral noise as well as a correction for inhomogeneous transmit ( B1+ ) using a B1+ map aligned to the coil position for each patient scan. Processed HP data and corrected flip angles were analyzed with an inputless two-site exchange model to calculate kPL . RESULTS: Denoising averages SNR increases of pyruvate, lactate, and alanine were 37.4-, 34.0-, and 20.1-fold, respectively, with lactate and alanine dynamics most noticeably recovered and better defined. In agreement with Monte Carlo simulations, over-flipped regions underestimated kPL and under-flipped regions overestimated kPL . B1+ correction addressed this issue. CONCLUSION: The new HP 13 C EPSI post-processing pipeline integrated tensor denoising and B1+ correction to measure kPL in patients with liver tumors. These technical developments not only recovered metabolite signals in voxels that did not receive the prescribed flip angle, but also increased the extent and accuracy of kPL estimations throughout the tumor and adjacent regions including normal-appearing tissue and additional lesions.

B1 inhomogeneity correction of RARE MRI with transceive surface radiofrequency probes.

PURPOSE: The use of surface radiofrequency (RF) coils is common practice to boost sensitivity in (pre)clinical MRI. The number of transceive surface RF coils is rapidly growing due to the surge in cryogenically cooled RF technology and ultrahigh-field MRI. Consequently, there is an increasing need for effective correction of the excitation field ( B1+ ) inhomogeneity inherent in these coils. Retrospective B1 correction permits quantitative MRI, but this usually requires a pulse sequence-specific analytical signal intensity (SI) equation. Such an equation is not available for fast spin-echo (Rapid Acquisition with Relaxation Enhancement, RARE) MRI. Here we present, test, and validate retrospective B1 correction methods for RARE. METHODS: We implemented the commonly used sensitivity correction and developed an empirical model-based method and a hybrid combination of both. Tests and validations were performed with a cryogenically cooled RF probe and a single-loop RF coil. Accuracy of SI quantification and T1 contrast were evaluated after correction. RESULTS: The three described correction methods achieved dramatic improvements in B1 homogeneity and significantly improved SI quantification and T1 contrast, with mean SI errors reduced from >40% to >10% following correction in all cases. Upon correction, images of phantoms and mouse heads demonstrated homogeneity comparable to that of images acquired with a volume resonator. This was quantified by SI profile, SI ratio (error < 10%), and percentage of integral uniformity (PIU > 80% in vivo and ex vivo compared to PIU > 87% with the reference RF coil). CONCLUSION: This work demonstrates the efficacy of three B1 correction methods tailored for transceive surface RF probes and RARE MRI. The corrected images are suitable for quantification and show comparable results between the three methods, opening the way for T1 measurements and X-nuclei quantification using surface transceiver RF coils. This approach is applicable to other MR techniques for which no analytical SI exists.

CEST MRI quantification procedures for breast cancer treatment-related lymphedema therapy evaluation.

PURPOSE: To quantify chemical exchange saturation transfer contrast in upper extremities of participants with lymphedema before and after standardized lymphatic mobilization therapy using correction procedures for B0 and B1 heterogeneity, and T1 relaxation. METHODS: Females with (n = 12) and without (n = 17) breast cancer treatment-related lymphedema (BCRL) matched for age and body mass index were scanned at 3.0T MRI. B1 efficiency and T1 were calculated in series with chemical exchange saturation transfer in bilateral axilla (B1 amplitude = 2µT, Δω = ±5.5 ppm, slices = 9, spatial resolution = 1.8 × 1.47 × 5.5 mm3 ). B1 dispersion measurements (B1 = 1-3 µT; increment = 0.5 µT) were performed in controls (n = 6 arms in 3 subjects). BCRL participants were scanned pre- and post-manual lymphatic drainage (MLD) therapy. Chemical exchange saturation transfer amide proton transfer (APT) and nuclear Overhauser effect (NOE) metrics corrected for B1 efficiency were calculated, including proton transfer ratio (PTR'), magnetization transfer ratio asymmetry (MTRasymmetry') , and apparent exchange-dependent relaxation (AREX'). Nonparametric tests were used to evaluate relationships between metrics in BCRL participants pre- versus post-MLD (two-sided P < 0.05 required for significance). RESULTS: B1 dispersion experiments showed nonlinear dependence of Z-values on B1 efficiency in the upper extremities; PTR' showed < 1% mean fractional difference between subject-specific and group-level correction procedures. PTR'APT significantly correlated with T1 (Spearman's rho = 0.57, P < 0.001) and body mass index (Spearman's rho = -0.37, P = 0.029) in controls and with lymphedema stage (Spearman's rho = 0.48, P = 0.017) in BCRL participants. Following MLD therapy, PTR'APT significantly increased in the affected arm of BCRL participants (pre- vs. post-MLD: 0.41 ± 0.05 vs. 0.43 ± 0.03, P = 0.02), consistent with treatment effects from mobilized lymphatic fluid. CONCLUSION: Chemical exchange saturation transfer metrics, following appropriate correction procedures, respond to lymphatic mobilization therapies and may have potential for evaluating treatments in participants with secondary lymphedema.

Slice profile and B1 corrections in 2D magnetic resonance fingerprinting.

PURPOSE: The goal of this study is to characterize and improve the accuracy of 2D magnetic resonance fingerprinting (MRF) scans in the presence of slice profile (SP) and B1 imperfections, which are two main factors that affect quantitative results in MRF. METHODS: The SP and B1 imperfections are characterized and corrected separately. The SP effect is corrected by simulating the radiofrequency pulse in the dictionary, and the B1 is corrected by acquiring a B1 map using the Bloch-Siegert method before each scan. The accuracy, precision, and repeatability of the proposed method are evaluated in phantom studies. The effects of both SP and B1 imperfections are also illustrated and corrected in the in vivo studies. RESULTS: The SP and B1 corrections improve the accuracy of the T1 and T2 values, independent of the shape of the radiofrequency pulse. The T1 and T2 values obtained from different excitation patterns become more consistent after corrections, which leads to an improvement of the robustness of the MRF design. CONCLUSION: This study demonstrates that MRF is sensitive to both SP and B1 effects, and that corrections can be made to improve the accuracy of MRF with only a 2-s increase in acquisition time. Magn Reson Med 78:1781-1789, 2017. © 2017 International Society for Magnetic Resonance in Medicine.

Perturbation of longitudinal relaxation rate in rotating frame (PLRF) analysis for quantification of chemical exchange saturation transfer signal in a transient state.

PURPOSE: To develop a novel analytical method for quantification of chemical exchange saturation transfer (CEST) in the transient state. The proposed method aims to reduce the effects of non-chemical-exchange (non-CE) parameters on the CEST signal, emphasizing the effect of chemical exchange. METHODS: The difference in the longitudinal relaxation rate in the rotating frame ( ΔR1ρ) was calculated based on perturbation of the Z-value by R1ρ, and a saturation-pulse-amplitude-compensated exchange-dependent relaxation rate (SPACER) was determined with a high-exchange-rate approximation. In both phantom and human subject experiments, MTRasym (representative of the traditional CEST index), ΔR1ρ, and SPACER were measured, evaluated, and compared by altering the non-CE parameters in a transient-state continuous-wave CEST sequence. RESULTS: In line with the theoretical expectation, our experimental data demonstrate that the effects of the non-CE parameters can be more effectively reduced using the proposed indices (  ΔR1ρ and SPACER) than using the traditional CEST index ( MTRasym). CONCLUSION: The proposed method allows for the chemical exchange weight to be better emphasized in the transient-state CEST signal, which is beneficial, in practice, for quantifying the CEST signal. Magn Reson Med 78:1711-1723, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

On the transmit field inhomogeneity correction of relaxation-compensated amide and NOE CEST effects at 7 T.

High field MRI is beneficial for chemical exchange saturation transfer (CEST) in terms of high SNR, CNR, and chemical shift dispersion. These advantages may, however, be counter-balanced by the increased transmit field inhomogeneity normally associated with high field MRI. The relatively high sensitivity of the CEST contrast to B1 inhomogeneity necessitates the development of correction methods, which is essential for the clinical translation of CEST. In this work, two B1 correction algorithms for the most studied CEST effects, amide-CEST and nuclear Overhauser enhancement (NOE), were analyzed. Both methods rely on fitting the multi-pool Bloch-McConnell equations to the densely sampled CEST spectra. In the first method, the correction is achieved by using a linear B1 correction of the calculated amide and NOE CEST effects. The second method uses the Bloch-McConnell fit parameters and the desired B1 amplitude to recalculate the CEST spectra, followed by the calculation of B1 -corrected amide and NOE CEST effects. Both algorithms were systematically studied in Bloch-McConnell equations and in human data, and compared with the earlier proposed ideal interpolation-based B1 correction method. In the low B1 regime of 0.15-0.50 µT (average power), a simple linear model was sufficient to mitigate B1 inhomogeneity effects on a par with the interpolation B1 correction, as demonstrated by a reduced correlation of the CEST contrast with B1 in both the simulations and the experiments.

Correction of B1-inhomogeneities for relaxation-compensated CEST imaging at 7 T.

Chemical exchange saturation transfer (CEST) imaging of endogenous agents in vivo is influenced by direct water proton saturation (spillover) and semi-solid macromolecular magnetization transfer (MT). Lorentzian fit isolation and application of the inverse metric yields the pure CEST contrast AREX, which is less affected by these processes, but still depends on the measurement technique, in particular on the irradiation amplitude B1 of the saturation pulses. This study focuses on two well-known CEST effects in the slow exchange regime originating from amide and aliphatic protons resonating at 3.5 ppm or -3.5 ppm from water protons, respectively. A B1-correction of CEST contrasts is crucial for the evaluation of data obtained in clinical studies at high field strengths with strong B1-inhomogeneities. Herein two approaches for B1-inhomogeneity correction, based on either CEST contrasts or Z-spectra, are investigated. Both rely on multiple acquisitions with different B1-values. One volunteer was examined with eight different B1-values to optimize the saturation field strength and the correction algorithm. Histogram evaluation allowed quantification of the quality of the B1-correction. Finally, the correction was applied to CEST images of a patient with oligodendroglioma WHO grade 2, and showed improvement of the image quality compared with the non-corrected CEST images, especially in the tumor region.

Improved quantitative (19) F MR molecular imaging with flip angle calibration and B1 -mapping compensation.

PURPOSE: To improve (19) F flip angle calibration and compensate for B1 inhomogeneities in quantitative (19) F MRI of sparse molecular epitopes with perfluorocarbon (PFC) nanoparticle (NP) emulsion contrast agents. MATERIALS AND METHODS: Flip angle sweep experiments on PFC-NP point source phantoms with three custom-designed (19) F/(1) H dual-tuned coils revealed a difference in required power settings for (19) F and (1) H nuclei, which was used to calculate a calibration ratio specific for each coil. An image-based correction technique was developed using B1 -field mapping on (1) H to correct for (19) F and (1) H images in two phantom experiments. RESULTS: Optimized (19) F peak power differed significantly from that of (1) H power for each coil (P < 0.05). A ratio of (19) F/(1) H power settings yielded a coil-specific and spatially independent calibration value (surface: 1.48 ± 0.06; semicylindrical: 1.71 ± 0.02, single-turn-solenoid: 1.92 ± 0.03). (1) H-image-based B1 correction equalized the signal intensity of (19) F images for two identical (19) F PFC-NP samples placed in different parts of the field, which were offset significantly by ~66% (P < 0.001), before correction. CONCLUSION: (19) F flip angle calibration and B1 -mapping compensations to the (19) F images employing the more abundant (1) H signal as a basis for correction resulted in a significant change in the quantification of sparse (19) F MR signals from targeted PFC NP emulsions.

Combined estimation of B1and T1for dynamic contrast-enhanced MRI by accounting for incomplete spoiling of transverse magnetization.

Objective.The variable flip angle (VFA) method for longitudinal relaxation time (T1) measurement is inherently sensitive to inaccuracies in the radiofRequency transmit field (B1) and incomplete spoiling of transverse magnetization. The objective of this study is to devise a computational method that addresses the problems of incomplete spoiling andB1inhomogeneity in the estimation ofT1using VFA method.Approach. Using an analytical expression of the gradient echo signal with account of incomplete spoiling, we first showed that ill-posedness in the simultaneous estimation ofB1andT1can be lifted with the use of flip angles larger than the Ernst angle. We then devised a nonlinear optimization method based on this signal model of incomplete spoiling for simultaneous estimation ofB1andT1.Main results. We evaluated the proposed method on a graded-concentration phantom to show that the derivedT1estimates offers an improvement over the regular VFA method and compares well with reference values measured by inversion recovery. Reduction of the number of flip angles from 17 to 5 yielded consistent results indicating that the proposed method is numerically stable.T1estimates derived from in-vivo brain imaging were consistent with literature values for gray and white matter tissues.Significance. Contrary to the common notion thatB1correction in the VFA method forT1mapping should be performed separately, we show that combined estimation ofB1andT1is feasible by the proposed method simply with the acquisition of 5 flip angles, as demonstrated on both phantom and in-vivo imaging data.

Fast WASABI post-processing: Access to rapid B0 and B1 correction in clinical routine for CEST MRI.

CEST MRI methods, such as APT and NOE imaging reveal biomarkers with significant diagnostic potential due to their ability to access molecular tissue information. Regardless of the technique used, CEST MRI data are affected by static magnetic B0 and radiofrequency B1 field inhomogeneities that degrade their contrast. For this reason, the correction of B0 field-induced artefacts is essential, whereas accounting for B1 field inhomogeneities have shown significant improvements in image readability. In a previous work, an MRI protocol called WASABI was presented, which can map simultaneously B0 and B1 field inhomogeneities, while maintaining the same sequence and readout types as used for CEST MRI. Despite the highly satisfactory quality of B0 and B1 maps computed from the WASABI data, the post-processing method is based on an exhaustive search of a four-parameter space and an additional four-parameter non-linear model fitting step. This leads to long post-processing times that are prohibitive in clinical practice. This work provides a new method for fast post-processing of WASABI data with outstanding acceleration of the parameter estimation procedure and without compromising its stability. The resulting computational acceleration makes the WASABI technique suitable for clinical use. The stability of the method is demonstrated on phantom data and clinical 3 Tesla in vivo data.

Corrigendum: The Influence of Radio-Frequency Transmit Field Inhomogeneities on the Accuracy of G-ratio Weighted Imaging.

[This corrects the article DOI: 10.3389/fnins.2021.674719.].

The Influence of Radio-Frequency Transmit Field Inhomogeneities on the Accuracy of G-ratio Weighted Imaging.

G-ratio weighted imaging is a non-invasive, in-vivo MRI-based technique that aims at estimating an aggregated measure of relative myelination of axons across the entire brain white matter. The MR g-ratio and its constituents (axonal and myelin volume fraction) are more specific to the tissue microstructure than conventional MRI metrics targeting either the myelin or axonal compartment. To calculate the MR g-ratio, an MRI-based myelin-mapping technique is combined with an axon-sensitive MR technique (such as diffusion MRI). Correction for radio-frequency transmit (B1+) field inhomogeneities is crucial for myelin mapping techniques such as magnetization transfer saturation. Here we assessed the effect of B1+ correction on g-ratio weighted imaging. To this end, the B1+ field was measured and the B1+ corrected MR g-ratio was used as the reference in a Bland-Altman analysis. We found a substantial bias (≈-89%) and error (≈37%) relative to the dynamic range of g-ratio values in the white matter if the B1+ correction was not applied. Moreover, we tested the efficiency of a data-driven B1+ correction approach that was applied retrospectively without additional reference measurements. We found that it reduced the bias and error in the MR g-ratio by a factor of three. The data-driven correction is readily available in the open-source hMRI toolbox (www.hmri.info) which is embedded in the statistical parameter mapping (SPM) framework.

Fast B1+ mapping using three consecutive RF pulses and balanced gradients for improved bSSFP imaging.

PURPOSE: To develop a B1+ mapping during the transient phase of balanced steady state free precession (bSSFP) imaging which can be used for subsequent B1+ inhomogeneity compensation. METHODS: Two images with different flip angles (FA) are acquired using single-shot spiral technique during the transient phase of bSSFP with three consecutive RF pulses and balanced gradients. Under the assumptions that the transmit (B1+) field varies slowly in spatial domain and T1 and T2 relaxation effects are negligible during 2·TR, B1+ was estimated using the two magnitude images and bSSFP data was sequentially acquired. B1+ estimation error due to the assumptions and other factors such as FA and off-resonance were assessed using Bloch simulation. Phantom and in vivo experiments were performed with α-2α-3α scheme. RESULTS: The simulation results indicated that the proposed method was less sensitive to T1 relaxation and B1+ mapping FA (α) of approximately 60° produced minimum estimation error. The B1+-induced intensity variation was reduced with the proposed method in the phantom experiment. For both the phantom and in vivo experiments, the estimated B1+ map showed comparable to the conventional B1+ map using spin-echo DAM. CONCLUSION: B1+ map was estimated during the transient phase of bSSFP and subsequently compensated bSSFP images. There was no scan time increment and hence the technique can be used in a prescan manner for B1+ mapping or shimming.

Minimum acquisition methods for simultaneously imaging T(1), T(2), and proton density with B(1) correction and no spin-echoes.

The spin lattice (T(1)) and spin-spin (T(2)) relaxation times, along with the proton density (PD) contain almost all of the information that (1)H MRI routinely uses in clinical diagnosis and research, but are seldom imaged directly. Here, three methods for directly imaging T(1), T(2), and PD with the least possible number of acquisitions - three, are presented. All methods utilize long 0° self-refocusing adiabatic pre-pulses instead of spin-echoes to encode the T(2) information prior to a conventional gradient-echo MRI sequence. T(1) information is encoded by varying the flip-angle (FA) in the 'Dual-τ Dual-FA' and 'Four-FA' methods, or the sequence repetition period, TR, in the 'Dual-τ Dual-TR' method. Inhomogeneity in the FA distribution and slice-selection profile are recognized as the main error sources for T(1) measurements. The former is remedied by integrating an extra FA-dependent acquisition into the 'Four-FA' method to provide self-corrected T(1), T(2), PD, and FA in just four acquisitions - again, the minimum possible. Slice profile errors - which manifest as differences between 2D and 3D T(1) measurements, can be addressed by Bloch equation analysis and experimental calibration. All three methods are validated in phantom studies, and the 'Dual-τ Dual-FA' and 'Four-FA' methods are validated in human brain studies using standard partial saturation and spin-echo methods for reference. The new methods offer a minimum-acquisition option for imaging single-component T(1), T(2), and PD. 'Four-FA' performs best overall in accuracy, with high efficiency per unit accuracy vs. existing methods when B(1)-inhomogeneity is appropriately addressed.