Variable flip angle T1 mapping in the human brain with reduced T2 sensitivity using fast radiofrequency-spoiled gradient echo imaging

Comparison of test-retest reproducibility of DESPOT and 3D-QALAS for water T 1 and T 2 mapping

Purpose

Relaxometry, specifically T 1 and T 2 mapping, has become an essential technique for assessing the properties of biological tissues related to various physiological and pathological conditions. Many techniques are being used to estimate T 1 and T 2 relaxation times, ranging from the traditional inversion or saturation recovery and spin-echo sequences to more advanced methods. Choosing the appropriate method for a specific application is critical since the precision and accuracy of T 1 and T 2 measurements are influenced by a variety of factors including the pulse sequence and its parameters, the inherent properties of the tissue being examined, the MRI hardware, and the image reconstruction. The aim of this study is to evaluate and compare the test-retest reproducibility of two advanced MRI relaxometry techniques (Driven Equilibrium Single Pulse Observation of T 1 and T 2, DESPOT, and 3D Quantification using an interleaved Look-Locker acquisition Sequence with a T 2 preparation pulse, QALAS), for T 1 and T 2 mapping in a healthy volunteer cohort.

Methods

10 healthy volunteers underwent brain MRI at 1.3 mm3 isotropic resolution, acquiring DESPOT and QALAS data (~11.8 and ~5 minutes duration, including field maps, respectively), test-retest with subject repositioning, on a 3.0 Tesla Philips Ingenia Elition scanner. To reconstruct the T 1 and T 2 maps, we used an equation-based algorithm for DESPOT and a dictionary-based algorithm that incorporates inversion efficiency and B 1 -field inhomogeneity for QALAS. The test-retest reproducibility was assessed using the coefficient of variation (CoV), intraclass correlation coefficient (ICC) and Bland-Altman plots.

Results

Our results indicate that both the DESPOT and QALAS techniques demonstrate good levels of test-retest reproducibility for T 1 and T 2 mapping across the brain. Higher whole-brain voxel-to-voxel ICCs are observed in QALAS for T 1 (0.84 ± 0.039) and in DESPOT for T 2 (0.897 ± 0.029). The Bland-Altman plots show smaller bias and variability of T 1 estimates for QALAS (mean of -0.02 s, and upper and lower limits of -0.14 and 0.11 s, 95% CI) than for DESPOT (mean of -0.02 s, and limits of -0.31 and 0.27 s). QALAS also showed less variability (mean 1.08 ms, limits -1.88 to 4.04 ms) for T 2 compared to DESPOT (mean of 2.56 ms, and limits -17.29 to 22.41 ms). The within-subject CoVs for QALAS range from 0.6% (T 2 in CSF) to 5.8% (T 2 in GM), while for DESPOT they range from 2.1% (T 2 in CSF) to 6.7% (T 2 in GM). The between-subject CoVs for QALAS range from 2.5% (T 2 in GM) to 12% (T 2 in CSF), and for DESPOT they range from 3.7% (T 2 in WM) to 9.3% (T 2 in CSF).

Conclusion

Overall, QALAS demonstrated better reproducibility for T 1 and T 2 measurements than DESPOT, in addition to reduced acquisition time.

Intra-scan RF power amplifier drift correction

PURPOSE

The drift in radiofrequency (RF) power amplifiers (RFPAs) is assessed and several contributing factors are investigated. Two approaches for prospective correction of drift are proposed and their effectiveness is evaluated.

METHODS

RFPA drift assessment encompasses both intra-pulse and inter-pulse drift analyses. Scan protocols with varying flip angle (FA), RF length, and pulse repetition time (TR) are used to gauge the influence of these parameters on drift. Directional couplers (DICOs) monitor the forward waveforms of the RFPA outputs. DICOs data is stored for evaluation, allowing calculation of correction factors to adjust RFPAs' transmit voltage. Two correction methods, predictive and run-time, are employed: predictive correction necessitates a calibration scan, while run-time correction calculates factors during the ongoing scan.

RESULTS

RFPA drift is indeed influenced by the RF duty-cycle, and in the cases examined with a maximum duty-cycle of 66%, the potential drift is approximately 41% or 15%, depending on the specific RFPA revision. Notably, in low transmit voltage scenarios, FA has minimal impact on RFPA drift. The application of predictive and run-time drift correction techniques effectively reduces the average drift from 10.0% to less than 1%, resulting in enhanced MR signal stability.

CONCLUSION

Utilizing DICO recordings and implementing a feedback mechanism enable the prospective correction of RFPA drift. Having a calibration scan, predictive correction can be utilized with fewer complexity; for enhanced performance, a run-time approach can be employed.

Magnetic resonance relaxometry in quantitative imaging of brain gliomas: A literature review

Magnetic resonance (MR) relaxometry is a quantitative imaging method that measures tissue relaxation properties. This review discusses the state of the art of clinical proton MR relaxometry for glial brain tumors. Current MR relaxometry technology also includes MR fingerprinting and synthetic MRI, which solve the inefficiencies and challenges of earlier techniques. Despite mixed results regarding its capability for brain tumor differential diagnosis, there is growing evidence that MR relaxometry can differentiate between gliomas and metastases and between glioma grades. Studies of the peritumoral zones have demonstrated their heterogeneity and possible directions of tumor infiltration. In addition, relaxometry offers T2* mapping that can define areas of tissue hypoxia not discriminated by perfusion assessment. Studies of tumor therapy response have demonstrated an association between survival and progression terms and dynamics of native and contrast-enhanced tumor relaxometric profiles. In conclusion, MR relaxometry is a promising technique for glial tumor diagnosis, particularly in association with neuropathological studies and other imaging techniques.

Quantification of the in vivo brain ultrashort-T2* component in healthy volunteers

PURPOSE

Recent work has shown MRI is able to measure and quantify signals of phospholipid membrane-bound protons associated with myelin in the human brain. This work seeks to develop an improved technique for characterizing this brain ultrashort-T 2 ∗ $$ {\mathrm{T}}_2\ast $$ component in vivo accounting forT 1 $$ {\mathrm{T}}_1 $$ weighting.

METHODS

Data from ultrashort echo time scans from 16 healthy volunteers with variable flip angles (VFA) were collected and fitted into an advanced regression model to quantify signal fraction, relaxation time, and frequency shift of the ultrashort-T 2 ∗ $$ {\mathrm{T}}_2\ast $$ component.

RESULTS

The fitted components show intra-subject differences of different white matter structures and significantly elevated ultrashort-T 2 ∗ $$ {\mathrm{T}}_2\ast $$ signal fraction in the corticospinal tracts measured at 0.09 versus 0.06 in other white matter structures and significantly elevated ultrashort-T 2 ∗ $$ {\mathrm{T}}_2\ast $$ frequency shift in the body of the corpus callosum at- $$ - $$ 1.5 versus- $$ - $$ 2.0 ppm in other white matter structures.

CONCLUSION

The significantly different measured components and measuredT 1 $$ {\mathrm{T}}_1 $$ relaxation time of the ultrashort-T 2 ∗ $$ {\mathrm{T}}_2\ast $$ component suggest that this method is picking up novel signals from phospholipid membrane-bound protons.

Rapid quantitative magnetization transfer imaging: Utilizing the hybrid state and the generalized Bloch model

PURPOSE

To explore efficient encoding schemes for quantitative magnetization transfer (qMT) imaging with few constraints on model parameters.

THEORY AND METHODS

We combine two recently proposed models in a Bloch-McConnell equation: the dynamics of the free spin pool are confined to the hybrid state, and the dynamics of the semi-solid spin pool are described by the generalized Bloch model. We numerically optimize the flip angles and durations of a train of radio frequency pulses to enhance the encoding of three qMT parameters while accounting for all eight parameters of the two-pool model. We sparsely sample each time frame along this spin dynamics with a three-dimensional radial koosh-ball trajectory, reconstruct the data with subspace modeling, and fit the qMT model with a neural network for computational efficiency.

RESULTS

We extracted qMT parameter maps of the whole brain with an effective resolution of 1.24 mm from a 12.6-min scan. In lesions of multiple sclerosis subjects, we observe a decreased size of the semi-solid spin pool and longer relaxation times, consistent with previous reports.

CONCLUSION

The encoding power of the hybrid state, combined with regularized image reconstruction, and the accuracy of the generalized Bloch model provide an excellent basis for efficient quantitative magnetization transfer imaging with few constraints on model parameters.

Modeling Normal Mouse Uterine Contraction and Placental Perfusion with Non-invasive Longitudinal Dynamic Contrast Enhancement MRI

The placenta is a transient organ critical for fetal development. Disruptions of normal placental functions can impact health throughout an individual's entire life. Although being recognized by the NIH Human Placenta Project as an important organ, the placenta remains understudied, partly because of a lack of non-invasive tools for longitudinally evaluation for key aspects of placental functionalities. Non-invasive imaging that can longitudinally probe murine placental health in vivo are critical to understanding placental development throughout pregnancy. We developed advanced imaging processing schemes to establish functional biomarkers for non-invasive longitudinal evaluation of placental development. We developed a dynamic contrast enhancement magnetic resonance imaging (DCE-MRI) pipeline combined with advanced image process methods to model uterine contraction and placental perfusion dynamics. Our novel imaging pipeline uses subcutaneous administration of gadolinium for steepest-slope based perfusion evaluation. This enables non-invasive longitudinal monitoring. Additionally, we advance the placental perfusion chamber paradigm with a novel physiologically-based threshold model for chamber localization and demonstrate spatially varying placental chambers using multiple functional metrics that assess mouse placental development and continuing remodeling throughout gestation. Lastly, using optic flow to quantify placental motions arisen from uterine contractions in conjunction with time-frequency analysis, we demonstrated that the placenta exhibited asymmetric contractile motion.

Vendor-neutral sequences and fully transparent workflows improve inter-vendor reproducibility of quantitative MRI

PURPOSE

We developed an end-to-end workflow that starts with a vendor-neutral acquisition and tested the hypothesis that vendor-neutral sequences decrease inter-vendor variability of T1, magnetization transfer ratio (MTR), and magnetization transfer saturation-index (MTsat) measurements.

METHODS

We developed and deployed a vendor-neutral 3D spoiled gradient-echo (SPGR) sequence on three clinical scanners by two MRI vendors. We then acquired T1 maps on the ISMRM-NIST system phantom, as well as T1, MTR, and MTsat maps in three healthy participants. We performed hierarchical shift function analysis in vivo to characterize the differences between scanners when the vendor-neutral sequence is used instead of commercial vendor implementations. Inter-vendor deviations were compared for statistical significance to test the hypothesis.

RESULTS

In the phantom, the vendor-neutral sequence reduced inter-vendor differences from 8% to 19.4% to 0.2% to 5% with an overall accuracy improvement, reducing ground truth T1 deviations from 7% to 11% to 0.2% to 4%. In vivo, we found that the variability between vendors is significantly reduced (p = 0.015) for all maps (T1, MTR, and MTsat) using the vendor-neutral sequence.

CONCLUSION

We conclude that vendor-neutral workflows are feasible and compatible with clinical MRI scanners. The significant reduction of inter-vendor variability using vendor-neutral sequences has important implications for qMRI research and for the reliability of multicenter clinical trials.

Repeatability and robustness of MP-GRASP T1 mapping

PURPOSE

To demonstrate the repeatability of fast 3D T1 mapping using Magnetization-Prepared Golden-angle RAdial Sparse Parallel (MP-GRASP) MRI and its robustness to variation of imaging parameters including flip angle and spatial resolution in phantoms and the brain.

THEORY AND METHODS

Multiple imaging experiments were performed to (1) assess the robustness of MP-GRASP T1 mapping to B1 inhomogeneity using a single tube phantom filled with uniform MnCl2 liquid; (2) compare the repeatability of T1 mapping between MP-GRASP and inversion recovery-based spin-echo (IR-SE; over 12 scans), using a commercial T1MES phantom; (3) evaluate the longitudinal variation of T1 estimation using MP-GRASP with varying imaging parameters, including spatial resolution, flip angle, TR/TE, and acceleration rate, using the T1MES phantom (106 scans performed over a period of 12 months); and (4) evaluate the variation of T1 estimation using MP-GRASP with varying imaging parameters in the brain (24 scans in a single visit). In addition, the accuracy of MP-GRASP T1 mapping was also validated against IR-SE by performing linear correlation and calculating the Lin's concordance correlation coefficient (CCC).

RESULTS

MP-GRASP demonstrates good robustness to B1 inhomogeneity, with intra-slice variability below 1% in the single tube phantom experiment. The longitudinal variability is good both in the phantom (below 2.5%) and in the brain (below 2%) with varying imaging parameters. The T1 values estimated from MP-GRASP are accurate compared to that from the IR-SE imaging (R2 = 0.997, Lin's CCC = 0.996).

CONCLUSION

MP-GRASP shows excellent repeatability of T1 estimation over time, and it is also robust to variation of different imaging parameters evaluated in this study.

Multi-parameter quantitative mapping of R1, R2*, PD, and MTsat is reproducible when accelerated with Compressed SENSE

Multi-parameter mapping (MPM) magnetic resonance imaging (MRI) provides quantitative estimates of the longitudinal and effective transverse relaxation rates R1 and R2*, proton density (PD), and magnetization transfer saturation (MTsat). Thereby, MPM enables better comparability across sites and time than conventional weighted MRI. However, for MPM, several contrasts must be acquired, resulting in prolonged measurement durations and thus preventing MPM's application in clinical routines. State-of-the-art imaging acceleration techniques such as Compressed SENSE (CS), a combination of compressed sensing and sensitivity encoding, can be used to reduce the scan time of MPM. However, the accuracy and precision of the resulting quantitative parameter maps have not been systematically evaluated. In this study, we therefore investigated the effect of CS acceleration on the fidelity and reproducibility of MPM acquisitions. In five healthy volunteers and in a phantom, we compared MPM metrics acquired without imaging acceleration, with the standard acceleration (SENSE factor 2.5), and with Compressed SENSE with acceleration factors 4 and 6 using a 32-channel head coil. We evaluated the reproducibility and repeatability of accelerated MPM using data from three scan sessions in gray and white matter volumes-of-interest (VOIs). Accelerated MPM provided precise and accurate quantitative parameter maps. For most parameters, the results of the CS-accelerated protocols correlated more strongly with the non-accelerated protocol than the standard SENSE-accelerated protocols. Furthermore, for most VOIs and contrasts, coefficients of variation were lower when calculated from data acquired with different imaging accelerations within a single scan session than from data acquired in different scan sessions. These results suggest that MPM with Compressed SENSE acceleration factors up to at least 6 yields reproducible quantitative parameter maps that are highly comparable to those acquired without imaging acceleration. Compressed SENSE can thus be used to considerably reduce the scan duration of R1, R2*, PD, and MTsat mapping, and is highly promising for clinical applications of MPM.

An optimized MP2RAGE sequence for studying both brain and cervical spinal cord in a single acquisition at 3T

Magnetization Prepared 2 Rapid Acquisition Gradient Echo (MP2RAGE) is a T₁ mapping technique that has been used broadly on brain and recently on cervical spinal cord (cSC). The growing interest for combined investigation of brain and SC in numerous pathologies of the central nervous system such as multiple sclerosis (MS), amyotrophic lateral sclerosis (ALS) and traumatic injuries, now brings about the need for optimization with regards to this specific investigation. This implies large spatial coverage with high spatial resolution and short acquisition time, high CNR and low B₁⁺ sensitivity, as well as high reproducibility and robust post-processing tools for T₁ quantification in different regions of brain and SC. In this work, a dedicated protocol (referred to as Pr-BSC) has been optimized for simultaneous brain and cSC T₁ MP2RAGE acquisition at 3T. After computer simulation optimization, the protocol was applied for in vivo validation experiments and compared to previously published state of the art protocols focusing on either the brain (Pr-B) or the cSC (Pr-SC). Reproducibility and in-ROI standard deviations were assessed on healthy volunteers in the perspective of future clinical use. The mean T₁ values, obtained by the Pr-BSC, in brain white, gray and deep gray matters were: (mean ± in-ROI SD) 792 ± 27 ms, 1339 ± 139 ms and 1136 ± 88 ms, respectively. In cSC, T₁ values for white matter corticospinal, posterior sensory, lateral sensory and rubro/reticulospinal tracts were 902 ± 41 ms, 920 ± 35 ms, 903 ± 46 ms, 891 ± 41 ms, respectively, and 954 ± 32 ms for anterior and intermediate gray matter. The Pr-BSC protocol showed excellent agreement with previously proposed Pr-B on brain and Pr-SC on cSC, with very high inter-scan reproducibility (coefficients of variation of 0.52 ± 0.36% and 1.12 ± 0.62% on brain and cSC, respectively). This optimized protocol covering both brain and cSC with a sub-millimetric isotropic spatial resolution in one acquisition of less than 8 min, opens up great perspectives for clinical applications focusing on degenerative tissue such as encountered in MS and ALS.

Imperfect spoiling in variable flip angle T1 mapping at 7T: Quantifying and minimizing impact

Purpose

The variable flip angle (VFA) approach to T₁ mapping assumes perfectly spoiled transverse magnetisation at the end of each repetition time (TR). Despite radiofrequency (RF) and gradient spoiling, this condition is rarely met, leading to erroneous T₁ estimates (T1app). Theoretical corrections can be applied but make assumptions about tissue properties, for example, a global T₂ time. Here, we investigate the effect of imperfect spoiling at 7T and the interaction between the RF and gradient spoiling conditions, additionally accounting for diffusion. We provide guidance on the optimal approach to maximise the accuracy of the T₁ estimate in the context of 3D multi‐echo acquisitions.

Methods

The impact of the spoiling regime was investigated through numerical simulations, phantom and invivo experiments.

Results

The predicted dependence of T1app on tissue properties, system settings, and spoiling conditions was observed in both phantom and in vivo experiments. Diffusion effects modulated the dependence of T1app on both B1+ efficiency and T₂ times.

Conclusion

Error in T1app can be minimized by using an RF spoiling increment and gradient spoiler moment combination that minimizes T₂‐dependence and safeguards image quality. Although the diffusion effect was comparatively small at 7T, correction factors accounting for this effect are recommended.

Click here for author‐reader discussions

The effect of gadolinium-based contrast agent administration on magnetic resonance fingerprinting-based T1 relaxometry in patients with prostate cancer

Magnetic resonance fingerprinting (MRF) is a rapidly developing fast quantitative mapping technique able to produce multiple property maps with reduced sensitivity to motion. MRF has shown promise in improving the diagnosis of clinically significant prostate cancer but requires further validation as part of a prostate multiparametric (mp) MRI protocol. mpMRI protocol mandates the inclusion of dynamic contrast enhanced (DCE) imaging, known for its significant T1 shortening effect. MRF could be used to measure both pre- and post-contrast T1 values, but its utility must be assessed. In this proof-of-concept study, we sought to evaluate the variation in MRF T1 measurements post gadolinium-based contrast agent (GBCA) injection and the utility of such T1 measurements to differentiate peripheral and transition zone tumours from normal prostatic tissue. We found that the T1 variation in all tissues increased considerably post-GBCA following the expected significant T1 shortening effect, compromising the ability of MRF T1 to identify transition zone lesions. We, therefore, recommend performing MRF T1 prior to DCE imaging to maintain its benefit for improving detection of both peripheral and transition zone lesions while reducing additional scanning time. Demonstrating the effect of GBCA on MRF T1 relaxometry in patients also paves the way for future clinical studies investigating the added value of post-GBCA MRF in PCa, including its dynamic analysis as in DCE-MRF.

Controlled saturation magnetization transfer for reproducible multivendor variable flip angle T1 and T2 mapping

Purpose

The widespread clinical application of quantitative MRI has been hindered by a lack of reproducibility across sites and vendors. Previous work has attributed this to incorrect B₁ mapping or insufficient spoiling conditions. We recently proposed the controlled saturation magnetization transfer (CSMT) framework and hypothesized that the lack of reproducibility can also be attributed to magnetization transfer effects. This work seeks to validate this hypothesis and demonstrate that reproducible multivendor single‐pool relaxometry can be achieved with the CSMT approach.

Methods

Three healthy volunteers were scanned on scanners from 3 vendors (GE Healthcare, Philips, Siemens). An extensive set of images necessary for joint T₁ and T₂ estimation were acquired with (1) each vendor default RF pulses and spoiling conditions; (2) harmonized RF spoiling; and (3) harmonized RF spoiling and CSMT pulses. Different subsets of images were used to generate 6 different T₁ and T₂ maps for each subject’s data from each vendor. Cross‐protocol, cross‐vendor, and test/retest variability were estimated.

Results

Harmonized RF spoiling conditions are insufficient to ensure good cross‐vendor reproducibility. Controlled saturation magnetization transfer allows cross‐protocol variability to be reduced from 18.3% to 4.0%. Whole‐brain variability using the same protocol was reduced from a maximum of 19% to 4.5% across sites. Both CSMT and native vendor RF conditions have a reported variability of less than 5% for repeat measures on the same vendor.

Conclusion

Magnetization transfer effects are a major contributor to intersite/intrasite variability of T₁ and T₂ estimation. Controlled saturation magnetization transfer stabilizes these effects, paving the way for the use of single‐pool T₁ and T₂ as a reliable source for clinical diagnosis across sites.

hMRI - A toolbox for quantitative MRI in neuroscience and clinical research

Neuroscience and clinical researchers are increasingly interested in quantitative magnetic resonance imaging (qMRI) due to its sensitivity to micro-structural properties of brain tissue such as axon, myelin, iron and water concentration. We introduce the hMRI-toolbox, an open-source, easy-to-use tool available on GitHub, for qMRI data handling and processing, presented together with a tutorial and example dataset. This toolbox allows the estimation of high-quality multi-parameter qMRI maps (longitudinal and effective transverse relaxation rates R1 and R2⋆, proton density PD and magnetisation transfer MT saturation) that can be used for quantitative parameter analysis and accurate delineation of subcortical brain structures. The qMRI maps generated by the toolbox are key input parameters for biophysical models designed to estimate tissue microstructure properties such as the MR g-ratio and to derive standard and novel MRI biomarkers. Thus, the current version of the toolbox is a first step towards in vivo histology using MRI (hMRI) and is being extended further in this direction. Embedded in the Statistical Parametric Mapping (SPM) framework, it benefits from the extensive range of established SPM tools for high-accuracy spatial registration and statistical inferences and can be readily combined with existing SPM toolboxes for estimating diffusion MRI parameter maps. From a user's perspective, the hMRI-toolbox is an efficient, robust and simple framework for investigating qMRI data in neuroscience and clinical research.

Age-Related Changes in Tissue Value Properties in Children: Simultaneous Quantification of Relaxation Times and Proton Density Using Synthetic Magnetic Resonance Imaging

OBJECTIVES

The properties of brain tissue undergo dynamic changes during maturation. T1 relaxation time (T1), T2 relaxation time (T2), and proton density (PD) are now simultaneously quantifiable within a clinically acceptable time, using a synthetic magnetic resonance imaging (MRI) sequence. This study aimed to provide age-specific reference values for T1, T2, and PD in children, using synthetic MRI.

MATERIALS AND METHODS

We included 89 children (median age, 18 months; range, 34 weeks of gestational age to 17 years) who underwent quantitative MRI, using a multidynamic, multiecho sequence on 3 T MRI, between December 2015 and November 2016, and had no abnormal MRI/neurologic assessment findings. T1, T2, and PD were simultaneously measured in each of the 22 defined white matter and gray matter regions of interest. The measured values were plotted against age, and a curve fitting model that best explained the age dependence of tissue values was identified. Age-specific regional tissue values were calculated using a fit equation.

RESULTS

The tissue values of all brain regions, except cortical PD, decreased with increasing age, and the robust negative association was best explained by modified biexponential model of the form Tissue values = T1 × exp (-C1 × age) + T2 × exp (-C2 × age). The quality of fit to the modified biexponential model was high in white matter and deep gray matter (white matter, R = 97%-99% [T1], 88%-95% [T2], 88%-97% [PD]; deep gray matter, R = 96%-97% [T1], 96% [T2], 49%-88% [PD]; cortex, 70%-83% [T1], 87%-90% [T2], 5%-27% [PD]). The white matter and deep gray matter changed the most dynamically within the first year of life.

CONCLUSIONS

Our study provides age-specific regional reference values, from the neonate to adolescent, of T1, T2, and PD, which could be objective tools for assessment of normal/abnormal brain development using synthetic MRI.

NOVIFAST: A Fast Algorithm for Accurate and Precise VFA MRI Mapping

In quantitative magnetic resonance mapping, the variable flip angle (VFA) steady state spoiled gradient recalled echo (SPGR) imaging technique is popular as it provides a series of high resolution weighted images in a clinically feasible time. Fast, linear methods that estimate maps from these weighted images have been proposed, such as DESPOT1 and iterative re-weighted linear least squares. More accurate, non-linear least squares (NLLS) estimators are in play, but these are generally much slower and require careful initialization. In this paper, we present NOVIFAST, a novel NLLS-based algorithm specifically tailored to VFA SPGR mapping. By exploiting the particular structure of the SPGR model, a computationally efficient, yet accurate and precise map estimator is derived. Simulation and in vivo human brain experiments demonstrate a twenty-fold speed gain of NOVIFAST compared with conventional gradient-based NLLS estimators while maintaining a high precision and accuracy. Moreover, NOVIFAST is eight times faster than the efficient implementations of the variable projection (VARPRO) method. Furthermore, NOVIFAST is shown to be robust against initialization.

Accuracy, repeatability, and interplatform reproducibility of T1 quantification methods used for DCE-MRI: Results from a multicenter phantom study

PURPOSE

To determine the in vitro accuracy, test-retest repeatability, and interplatform reproducibility of T₁ quantification protocols used for dynamic contrast-enhanced MRI at 1.5 and 3 T.

METHODS

A T₁ 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 T₁ measurements) and repeatability were assessed using general linear mixed models. Interplatform reproducibility was assessed using coefficients of variation.

RESULTS

For the common IR-SE protocol, accuracy (median error across platforms = 1.4-5.5%) was influenced predominantly by T₁ 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%).

CONCLUSIONS

The VFA protocols with 2 to 3 flip angles optimized for different applications achieved acceptable balance of extensive spatial coverage, accuracy, and repeatability in T₁ quantification (errors < 15%). Further optimization in terms of flip-angle choice for each tissue application, and the use of B₁ correction, are needed to improve the robustness of VFA protocols for T₁ mapping. Magn Reson Med 79:2564-2575, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Joint system relaxometry (JSR) and Crámer-Rao lower bound optimization of sequence parameters: A framework for enhanced precision of DESPOT T1 and T2 estimation

Purpose

This study aims to increase the precision of single‐compartment DESPOT relaxometry by two means: (i) a joint system relaxometry (JSR) approach that estimates parameters in a single step using all available data; and (ii) optimizing acquisition parameters by deploying a robust design tool based on the Crámer‐Rao lower bound (CRLB).

Methods

Following the development of the analysis and design capabilities, phantom and four in vivo subject experiments were performed to compare directly the precision achieved with DESPOT and JSR estimation using published protocols and protocols designed using a proposed CRLB framework.

Results

Experimental data demonstrate JSR's ability to decrease relaxometry estimation variance. Phantom results show 72 to 77% improvement using the same data as conventional DESPOT. This is further improved to 81 to 87% using optimal parameters. Both experiments show systematic bias depending on the acquisition parameters used, which are shown to be highly reproducible and to vary with different magnetization transfer conditions.

Conclusions

Compared with DESPOT, JSR produces reproducible relaxation maps with improved precision. Further improvement was achieved using CRLB as a protocol design tool. With this combined approach, it is possible to achieve submillimeter maps of ρ,T1,T2, and B0 in an 11‐min examination, making the approach appealing for potential clinical use. Magn Reson Med 79:234–245, 2018. © 2017 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.

Rapid and robust variable flip angle T1 mapping using interleaved two-dimensional multislice spoiled gradient echo imaging

PURPOSE

Conventional T₁ mapping using three-dimensional (3D) radiofrequency (RF) spoiled gradient echo (SPGR) imaging with short repetition times (TR) is adversely affected by incomplete spoiling (i.e. residual T₂ dependency). In this work, an optimized interleaved 2D multislice SPGR sequence scheme and an adapted postprocessing procedure are evaluated for highly T₂ -insensitive T₁ quantification of human brain tissues.

METHODS

An efficient 2D multislice SPGR protocol including a relatively long TR of 200 ms is investigated with careful consideration of cross talk and magnetization transfer effects. Based on the derived scan protocol, T₁ is quantified from the signal ratio of two SPGR datasets acquired at different flip angles. The effect of nonideal RF excitation profiles is incorporated into the SPGR signal model by performing Bloch simulations.

RESULTS

Simulations showed that the parameters of the SPGR protocol (such as TR and the spoiler gradient moments) guarantee virtually complete spoiling. This result was confirmed by T₁ measurements both in vitro using a 2% agar probe doped with 0.1 mM Gd (Gadovist) and in vivo in the human brain.

CONCLUSION

The derived 2D multislice SPGR protocol offers efficient, highly reproducible, and in particular T₂ -insensitive T₁ quantification of human brain tissues. Magn Reson Med 77:1606-1611, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Data-driven optimized flip angle selection for T1 estimation from spoiled gradient echo acquisitions

PURPOSE

Define criteria for selection of optimal flip angle sets for T1 estimation and evaluate effects on T1 mapping.

THEORY AND METHODS

Flip angle sets for spoiled gradient echo-based T1 mapping were selected by minimizing T1 estimate variance weighted by the joint density of M0 and T1 in an initial acquisition. The effect of optimized flip angle selection on T1 estimate error was measured using simulations and experimental data in the human and rat brain.

RESULTS

For two-point acquisitions, optimized angle sets were similar to those proposed by other groups and, therefore, performed similarly. For multipoint acquisitions, optimal angle sets for T1 mapping in the brain consisted of a repetition of two angles. Implementation of optimal angles reduced T1 estimate variance by 30-40% compared with a multipoint acquisition using a range of angles. Performance of the optimal angle set was equivalent to that of a repetition of the two-angle set selected using criteria proposed by other researchers.

CONCLUSION

Repetition of two carefully selected flip angles notably improves the precision of resulting T1 estimates compared with acquisitions using a range of flip angles. This work provides a flexible and widely applicable optimization method of particular use for those who repeatedly perform T1 estimation. Magn Reson Med 76:792-802, 2016. © 2015 Wiley Periodicals, Inc.