Snapshot whole-brain T1 relaxometry using steady-state prepared spiral multislice variable flip angle imaging.

To determine the in vitro accuracy, test-retest repeatability, and interplatform reproducibility of T1 quantification protocols used for dynamic contrast-enhanced MRI at 1.5 and 3 T. A T1 phantom with 14 samples was imaged at eight centers with a common inversion-recovery spin-echo (IR-SE) protocol and a variable flip angle (VFA) protocol using seven flip angles, as well as site-specific protocols (VFA with different flip angles, variable repetition time, proton density, and Look-Locker inversion recovery). Factors influencing the accuracy (deviation from reference NMR T1 measurements) and repeatability were assessed using general linear mixed models. Interplatform reproducibility was assessed using coefficients of variation. For the common IR-SE protocol, accuracy (median error across platforms = 1.4-5.5%) was influenced predominantly by T1 sample (P < 10-6 ), whereas test-retest repeatability (median error = 0.2-8.3%) was influenced by the scanner (P < 10-6 ). For the common VFA protocol, accuracy (median error = 5.7-32.2%) was influenced by field strength (P = 0.006), whereas repeatability (median error = 0.7-25.8%) was influenced by the scanner (P < 0.0001). Interplatform reproducibility with the common VFA was lower at 3 T than 1.5 T (P = 0.004), and lower than that of the common IR-SE protocol (coefficient of variation 1.5T: VFA/IR-SE = 11.13%/8.21%, P = 0.028; 3 T: VFA/IR-SE = 22.87%/5.46%, P = 0.001). Among the site-specific protocols, Look-Locker inversion recovery and VFA (2-3 flip angles) protocols showed the best accuracy and repeatability (errors < 15%). The VFA protocols with 2 to 3 flip angles optimized for different applications achieved acceptable balance of extensive spatial coverage, accuracy, and repeatability in T1 quantification (errors < 15%). Further optimization in terms of flip-angle choice for each tissue application, and the use of B1 correction, are needed to improve the robustness of VFA protocols for T1 mapping. Magn Reson Med 79:2564-2575, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

To compare the accuracy and precision of four different T1 estimation algorithms for modified Look-Locker inversion recovery (MOLLI). Four T1 estimation algorithms, including the original fit, inversion group (IG) fit, instantaneous signal loss simulation (InSiL), and Bloch equation simulation with slice profile correction (BLESSPC) were studied. T1 estimation accuracy, precision, reproducibility, and sensitivity to heart rate (HR), flip angle (FA), and acquisition scheme (AcS) variations were compared in simulation, phantom, and volunteer studies. T1 estimation accuracy of IG (-2.4% ± 3.9%) and original fit (-3.2% ± 1.4%) were worse than BLESSPC (0.2% ± 1.5%) and InSiL (-0.7% ± 2.1%). The original fit had the best precision for T1 from 409-1884 ms for the same FA (0.67% ± 0.16% versus 0.90% ± 0.23% using IG, 0.78% ± 0.11% using InSiL, 0.77% ± 0.12% using BLESSPC). BLESSPC generated the most consistent in vivo T1 values over different FAs and AcS, and the T1 estimation reproducibility was similar (P > 0.3) among the four methods when FA = 35°. When using FA = 50°, the reproducibility was significantly improved only when using BLESSPC (1.6% ± 0.9 versus 2.6% ± 1.9%, P < 0.05). BLESSPC has superior accuracy and is the least sensitive to FA, HR, and AcS variations. T1 estimation using BLESSPC and FA = 50° is superior to conventional MOLLI with FA = 35° in accuracy and precision. Further clinical studies in varying pathological conditions are warranted to confirm our findings. Magn Reson Med 78:1746-1756, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Variable flip angle T(1) mapping and actual flip-angle imaging B(1) mapping are widely used quantitative MRI methods employing radiofrequency spoiled gradient-echo pulse sequences. Incomplete elimination of the transverse magnetization in these sequences has been found to be a critical source of T(1) and B(1) measurement errors. In this study, comprehensive theoretical analysis of spoiling-related errors in variable flip angle and actual flip-angle imaging methods was performed using the combined isochromat summation and diffusion propagator model and validated by phantom experiments. The key theoretical conclusion is that correct interpretation of spoiling phenomena in fast gradient-echo sequences requires accurate consideration of the diffusion effect. A general strategy for improvement of T(1) and B(1) measurement accuracy was proposed based on the strong spoiling regimen, where diffusion-modulated spatial averaging of isochromats becomes a dominant factor determining magnetization evolution. Practical implementation of strongly spoiled variable flip angle and actual flip-angle imaging techniques requires sufficiently large spoiling gradient areas (A(G)) in combination with optimal radiofrequency phase increments (phi(0)). Optimal regimens providing <2% relative T(1) and B(1) measurement errors in a variety of tissues were theoretically derived for prospective in vivo variable flip angle (pulse repetition time = 15-20 ms, A(G) = 280-450 mT.ms/m, phi(0) = 169 degrees) and actual flip-angle imaging (pulse repetition time(1)/pulse repetition time(2) = 20/100 ms, A(G1)/A(G2) = 450/2250 mT.ms/m, phi(0) = 39 degrees) applications based on 25 mT/m maximal available gradient strength.

To design and evaluate an optimized PETRA (point-wise encoding time reduction with radial acquisition) sequence with long-T2 suppression at 3 Tesla. An adiabatic inversion recovery-based scheme was used to null the long-T2 signal. To minimize scan time, the signal was sampled multiple times after each inversion with variable excitation flip angles designed to yield constant short-T2 signal amplitude. The excitation pulses were phase-modulated, allowing for increased flip angle and higher signal-to-noise ratio (SNR). A fast, noniterative image reconstruction algorithm was designed to minimize image artifacts due to nonuniform excitation profile. Phase-modulated pulse excitation, along with the noniterative reconstruction algorithm, allows the use of larger radiofrequency pulse flip angles, resulting in effective suppression of long-T2 protons and improved image SNR without causing image artifacts. Midtibia images representative of collagen-bound water yielded SNR of 15 at 1-mm isotropic resolution in 6.5 minutes with a standard extremity coil. Further, the technology is shown to be suited for generating multi-angle projection images of bone akin to X-ray images displaying subtle anatomic detail. Optimized long-T2 suppressed PETRA allows imaging of bone matrix water unencumbered by long-T2 soft tissue and pore water protons, opening up new possibilities for anatomic bone imaging at isotropic resolution and quantification in clinically practical scan times. Magn Reson Med 77:989-997, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Motivation: A rapid 3D whole-brain high-resolution T1 quantification using IR is desired but time-consuming. Goal(s): Optimize IR protocols to make the acquisition faster with sufficient accuracy. Approach: Conventional IR was accelerated using a long 3D stack-of-spirals FLASH readout. Its accuracy was evaluated by comparing protocols using varying saturation times with a reference 2D-EPI readout and with different undersampling factors. Results: This study presented a method for 3D whole-brain 1mm-resolution T1 mapping using IR with stack-of-spirals FLASH acquisition under 3.6min, employing a short presaturation time, a long echo train, and a large undersampling factor. Impact: This study demonstrates a 3D whole-brain 1mm-resolution IR-based T1 mapping method under 3.6min. The optimized IR technique was accelerated using a 3D SoS-TFL readout with a short presaturation time, a long echo train, as well as a large undersampling factor.

To study in simulation and in theory the accuracy and precision of dynamic T1 measurements obtained using the previously published single-reference variable flip angle (SR-VFA) technique, with a focus on the effects of dynamic changes in T2 * on the calculation. Monte Carlo simulations were performed over 1000 noisy iterations for the VFA method, the SR-VFA method, and a proposed method, SR-VFA with a T2 * correction (SR-VFA-T2 *). Dynamic T1 estimates were calculated analytically for each method, with signals modeled by the steady-state spoiled gradient echo equation. The mean and standard deviation of these estimates were calculated and compared to truth, while varying repetition time (TR), baseline and dynamic T1 , echo time (TE), baseline and dynamic T2 *, flip angles, and the number of averages on baseline scans. Additionally, the variance of T1 in the SR-VFA and SR-VFA-T2 * methods was derived analytically based on the theory of propagation of errors. This equation was used to produce an inverse-variance weighted linear combination to improve T1 mapping precision in the SR-VFA-T2 * method. Flip angle sensitivity of dynamic T1 precision in the SR-VFA and SR-VFA-T2 * methods was also performed. Substantial bias can be produced by the SR-VFA method when the ratio of the T2 * decay of the dynamic signal versus that of the baseline signals deviates from 1, with a 0.01 deviation leading to approximately a 1% bias in cases of high SNR and TR ≫ T1 . This bias can be corrected by estimating the baseline and dynamic T2 * values in this ratio via multiecho measurements. The bias and precision of the SR-VFA-T2 * method, when normalized to scan time, is found to rival and sometimes improve upon the two flip angle VFA method when an inverse variance weighted linear combination is applied across its multiecho T1 maps. The analytic variance equation presented is found to be accurate within 1% relative to the Monte Carlo simulations over a broad parameter space. Flip angle ranges that maximize SR-VFA and SR-VFA-T2 *T1 precision over a broad parameter space are given, and each is defined relative to TR and T1 . Multiecho SR-VFA-T2 * T1 mapping is found in simulation and theory to be a promising alternative to the VFA method that maintains speed of the SR-VFA method with accuracy and precision similar to the VFA method.

Quantitative MRI may require correcting for nuisance parameters which can or must be constrained to independently measured or assumed values. The noise and/or bias in these constraints propagate to fitted parameters. For example, the case of refocusing pulse flip angle constraint in multiple spin echo T2 mapping is explored. An analytical expression for the mean-squared error of a parameter of interest was derived as a function of the accuracy and precision of an independent estimate of a nuisance parameter. The expression was validated by simulations and then used to evaluate the effects of flip angle (θ) constraint on the accuracy and precision of T⁁2 for a variety of multi-echo T2 mapping protocols. Constraining θ improved T⁁2 precision when the θ-map signal-to-noise ratio was greater than approximately one-half that of the first spin echo image. For many practical scenarios, constrained fitting was calculated to reduce not just the variance but the full mean-squared error of T⁁2, for bias in θ⁁≲6%. The analytical expression derived in this work can be applied to inform experimental design in quantitative MRI. The example application to T2 mapping provided specific cases, depending on θ⁁ accuracy and precision, in which θ⁁ measurement and constraint would be beneficial to T⁁2 variance or mean-squared error. Magn Reson Med 79:673-682, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

PurposeWe present a new Bloch‐Siegert four Angle Saturation Transfer (BOAST) method for measuring the creatine kinase (CK) first‐order effective rate constant kf in human myocardium at 7 tesla (T). BOAST combines a variant of the four‐angle saturation transfer (FAST) method using amplitude‐modulated radiofrequency pulses, phosphorus Bloch‐Siegert B1+‐mapping to determine the per‐voxel flip angles, and nonlinear fitting to Bloch simulations for postprocessing.MethodsOptimal flip angles and repetition time parameters were determined from Monte Carlo simulations. BOAST was validated in the calf muscle of two volunteers at 3T and 7T. The myocardial CK forward rate constant was then measured in 10 volunteers at 7T in 82 min (after 1H localization).ResultsBOAST kfCK values were 0.281 ± 0.002 s−1 in the calf and 0.35 ± 0.05 s−1 in myocardium. These are consistent with literature values from lower fields. Using a literature values for adenosine triphosphate concentration, we computed CK flux values of 4.55 ± 1.52 mmol kg−1 s−1. The sensitive volume for BOAST depends on the B1 inhomogeneity of the transmit coil.ConclusionBOAST enables measurement of the CK rate constant in the human heart at 7T, with spatial localization in three dimensions to 5.6 mL voxels, using a 10‐cm loop coil. Magn Reson Med 78:20–32, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

To develop an accurate T1 measurement method for short T2 tissues using a combination of a 3-dimensional ultrashort echo time cones actual flip angle imaging technique and a variable repetition time technique (3D UTE-Cones AFI-VTR) on a clinical 3T scanner. First, the longitudinal magnetization mapping function of the excitation pulse was obtained with the 3D UTE-Cones AFI method, which provided information about excitation efficiency and B1 inhomogeneity. Then, the derived mapping function was substituted into the VTR fitting to generate accurate T1 maps. Numerical simulation and phantom studies were carried out to compare the AFI-VTR method with a B1 -uncorrected VTR method, a B1 -uncorrected variable flip angle (VFA) method, and a B1 -corrected VFA method. Finally, the 3D UTE-Cones AFI-VTR method was applied to bovine bone samples (N = 6) and healthy volunteers (N = 3) to quantify the T1 of cortical bone. Numerical simulation and phantom studies showed that the 3D UTE-Cones AFI-VTR technique provides more accurate measurement of the T1 of short T2 tissues than the B1 -uncorrected VTR and VFA methods or the B1 -corrected VFA method. The proposed 3D UTE-Cones AFI-VTR method showed a mean T1 of 240 ± 25 ms for bovine cortical bone and 218 ± 10 ms for the tibial midshaft of human volunteers, respectively, at 3 T. The 3D UTE-Cones AFI-VTR method can provide accurate T1 measurements of short T2 tissues such as cortical bone. Magn Reson Med 80:598-608, 2018. © 2018 International Society for Magnetic Resonance in Medicine.

Variable flip angle (VFA) and modified Look-Locker inversion recovery (MOLLI) are frequently used for noninvasive evaluation of renal interstitial fibrosis (IF) in chronic kidney disease (CKD). However, controversy remains over which method is preferred. To compare the diagnostic efficacy of VFA and MOLLI for T1 mapping in evaluating renal IF. Prospective. Fifty-one participants with CKD (CKD stage 1-5, 35 males) and 18 healthy volunteers (eight males). 3.0 T, three-dimensional gradient echo sequence for B1+ VFA, and two-dimensional gradient echo sequence for MOLLI. Image quality was assessed on a five-point scale. Cortex and medulla T1 values (cT1 and mT1), corticomedullary T1 value difference (ΔT1, medulla - cortex), and corticomedullary T1 value ratio (ratio T1, cortex:medulla) were compared between VFA and MOLLI as well as between IF grade (0-4) based on biopsy. Intraclass correlation coefficient, Bland-Altman analysis, analysis of variance, Kruskal-Wallis test, correlation analysis, and receiver operating characteristics analysis with the area under the curve (AUC). P-value <0.05 was considered significant. MOLLI provided significantly better image quality compared to VFA. cT1 and mT1 values significantly differed between VFA and MOLLI (cT1-VFA: 1771.4 ± 139.4 msec vs. cT1-MOLLI: 1729.9 ± 132.1 msec; mT1-VFA: 2076.0 [interquartile range (IQR): 2045.9-2129.9] msec vs. mT1-MOLLI: 2039.2 [IQR: 1997.8-2071.6] msec). ΔT1 and ratio T1 values were not different between VFA and MOLLI (ΔT1: 300.8 ± 71.4 vs. 306.0 ± 78.4, respectively, P = 0.33 and ratio T1: 0.85 ± 0.038 vs. 0.85 ± 0.041, respectively, P = 0.064). No difference was observed between T1 variables and T1 mapping methods in diagnosing IF. ΔT1 and ratio T1 were not different between VFA and MOLLI. Both VFA and MOLLI are effective for noninvasive assessment of renal IF. 2 TECHNICAL EFFICACY: Stage 2.

PurposeThis study aimed to investigate the inter- and intraobserver reproducibility of quantitative T1 (qT1) measurements using manual and semiautomatic region of interest (ROI) placements. We hypothesized the usefulness of the semiautomatic method, which utilizes a three-dimensional (3D) anatomical relationship between the myometrium and other tissues, for minimizing ROI placement variation, thereby improving qT1 reproducibility compared to the manual approach. The semiautomatic approach, which considered anatomical relationships, was expected to enhance reproducibility by reducing ROI placement variabilities.Materials and methodsThis study recruited 23 healthy female volunteers. Data with variable flip angle (VFA) and inversion recovery were acquired using 3D-spoiled gradient echo and spin echo sequences, respectively. T1 maps were generated with VFA. Manual and semiautomatic ROI placements were independently conducted. Mean qT1 values were calculated from the T1 maps using the corresponding pixel values of the myometrial ROI. Inter- and intraobserver reproducibility of qT1 values was investigated. The inter- and intraobserver reproducibility of qT1 values was evaluated by calculating the coefficient of variation (CoV). Further, reproducibility was evaluated with inter- and intraobserver errors and intraclass correlation coefficients (ICCs). Bland–Altman analysis was utilized to compare the results, estimate bias, and determine the limits of agreement.ResultsThe mean inter- and intraobserver CoV of the qT1 values for semiautomatic ROI placement was significantly lower than those for manual ROI placement (p < 0.05 and p < 0.01, respectively). ICCs for semiautomatic ROI placement were greater than those for manual ROI placement. Further, the mean inter- and intraobserver errors for semiautomatic ROI placement were significantly lower than those for manual ROI placement (p < 0.05 and p < 0.01, respectively).ConclusionSemiautomatic ROI placement demonstrated high reproducibility of qT1 measurements compared with manual methods. Semiautomatic ROI placement may be useful for evaluating uterine qT1 with high reproducibility.

To apply quantitative whole-brain T1 -rho (T1ρ ) and T2 imaging to the detection and quantification of brain changes resulting from multiple sclerosis (MS). Twenty-three MS patients with clinically isolated syndrome (10) and relapsing remitting MS (13) phenotypes, compared with 24 age-matched healthy controls were imaged at 3 Tesla. An axial T1ρ -weighted three-dimensional turbo spin echo sequence with a variable flip angle and fluid suppression was used. Spin-lock times of 0, 20, 40, 60, 80, and 100 ms were used. Corresponding T2 maps were also acquired. Whole brain white matter (WM) T1ρ maps were elevated compared with controls (P = 0.002). WM lesion T1ρ and T2 values were highly correlated (r = 0.83), but T1ρ demonstrated 25% better contrast to noise ratio (P < 0.001). WM lesion T1ρ correlated with disease duration. Gray matter T1ρ was negatively correlated with the Expanded Disability Status Scale, r = -0.45, P = 0.03. Normal appearing gray matter and cortical gray matter lesions were negatively correlated on T1ρ , but not on T2 (rT1ρ = -0.63, pT1ρ = 0.03; rT2 = -0.17, pT2 = 0.6). T1ρ MRI demonstrates enhanced lesion contrast compared with T2 , and in some cases may provide complementary information. T1ρ may provide a useful measure of demyelinating processes in MS.

Variable flip angles are used in steady-state free precession (SSFP) acquisitions (e.g., time-of-flight) but to a lesser extent than in spin echo acquisitions. In balanced steady-state free precession, imaging is often assumed to occur during the steady state, which has been well described in the literature. However, in many cases, imaging occurs during the transient stage, and the use of variable flip angles can improve signal and thus image quality. Here, we present the calculation of flip angles in transient balanced steady-state free precession to generate a predefined signal profile. The signal profile was iteratively optimized to maximize the integral of the signal versus time curve. The key contribution of this work is the formulation of the flip angle as a deterministic function of the preceding and desired magnetization. Catalyzation schemes, e.g., Kaiser-windowed ramp, can be combined with variable flip angles balanced steady-state free precession to reduce signal oscillations. A uniform signal profile was used as an example to demonstrate the variable flip angle algorithm. Accuracy of the algorithm and Bloch simulations were verified with MRI phantom acquisitions. Renal angiograms were acquired using an inflow-based balanced steady-state free precession MR angiography technique; improved small-vessel depiction was observed in volunteer examinations.

Comparison of DCE-MRI parametric mapping using MP2RAGE and variable flip angle T1 mapping

A new time-efficient and accurate technique for simultaneous mapping of T(1) and B(1) is proposed based on a combination of the actual flip angle (FA) imaging and variable FA methods. Variable FA-actual FA imaging utilizes a single actual FA imaging and one or more spoiled gradient-echo acquisitions with a simultaneous nonlinear fitting procedure to yield accurate T(1)/B(1) maps. The advantage of variable FA-actual FA imaging is high accuracy at either short T(1) times or long repetition times in the actual FA imaging sequence. Simulations show this method is accurate to 0.03% in FA and 0.07% in T(1) for ratios of repetition time to T1 time over the range of 0.01-0.45. We show for the case of brain imaging that it is sufficient to use only one small FA spoiled gradient-echo acquisition, which results in reduced spoiling requirements and a significant scan time reduction compared to the original variable FA method. In vivo validation yielded high-quality 3D T(1) maps and T(1) measurements within 10% of previously published values and within a clinically acceptable scan time. The variable FA-actual FA imaging method will increase the accuracy and clinical feasibility of many quantitative MRI methods requiring T(1)/B(1) mapping such as dynamic contrast enhanced perfusion and quantitative magnetization transfer imaging.