Flip angle profile correction for T₁ and T₂ quantification with look-locker inversion recovery 2D steady-state free precession imaging

Quantification correction for free-breathing myocardial T1ρ mapping in mice using a recursively derived description of a T1ρ* relaxation pathway

BACKGROUND

Fast and accurate T1ρ mapping in myocardium is still a major challenge, particularly in small animal models. The complex sequence design owing to electrocardiogram and respiratory gating leads to quantification errors in in vivo experiments, due to variations of the T1ρ relaxation pathway. In this study, we present an improved quantification method for T1ρ using a newly derived formalism of a T1ρ* relaxation pathway.

METHODS

The new signal equation was derived by solving a recursion problem for spin-lock prepared fast gradient echo readouts. Based on Bloch simulations, we compared quantification errors using the common monoexponential model and our corrected model. The method was validated in phantom experiments and tested in vivo for myocardial T1ρ mapping in mice. Here, the impact of the breath dependent spin recovery time Trec on the quantification results was examined in detail.

RESULTS

Simulations indicate that a correction is necessary, since systematically underestimated values are measured under in vivo conditions. In the phantom study, the mean quantification error could be reduced from - 7.4% to - 0.97%. In vivo, a correlation of uncorrected T1ρ with the respiratory cycle was observed. Using the newly derived correction method, this correlation was significantly reduced from r = 0.708 (p < 0.001) to r = 0.204 and the standard deviation of left ventricular T1ρ values in different animals was reduced by at least 39%.

CONCLUSION

The suggested quantification formalism enables fast and precise myocardial T1ρ quantification for small animals during free breathing and can improve the comparability of study results. Our new technique offers a reasonable tool for assessing myocardial diseases, since pathologies that cause a change in heart or breathing rates do not lead to systematic misinterpretations. Besides, the derived signal equation can be used for sequence optimization or for subsequent correction of prior study results.

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 ( B 1 + ) inhomogeneity inherent in these coils. Retrospective B₁ 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 B₁ 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 T₁ contrast were evaluated after correction.

RESULTS

The three described correction methods achieved dramatic improvements in B₁ homogeneity and significantly improved SI quantification and T₁ 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 B₁ 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 T₁ measurements and X-nuclei quantification using surface transceiver RF coils. This approach is applicable to other MR techniques for which no analytical SI exists.

Exploration of New Contrasts, Targets, and MR Imaging and Spectroscopy Techniques for Neuromuscular Disease - A Workshop Report of Working Group 3 of the Biomedicine and Molecular Biosciences COST Action BM1304 MYO-MRI

Neuromuscular diseases are characterized by progressive muscle degeneration and muscle weakness resulting in functional disabilities. While each of these diseases is individually rare, they are common as a group, and a large majority lacks effective treatment with fully market approved drugs. Magnetic resonance imaging and spectroscopy techniques (MRI and MRS) are showing increasing promise as an outcome measure in clinical trials for these diseases. In 2013, the European Union funded the COST (co-operation in science and technology) action BM1304 called MYO-MRI (www.myo-mri.eu), with the overall aim to advance novel MRI and MRS techniques for both diagnosis and quantitative monitoring of neuromuscular diseases through sharing of expertise and data, joint development of protocols, opportunities for young researchers and creation of an online atlas of muscle MRI and MRS. In this report, the topics that were discussed in the framework of working group 3, which had the objective to: Explore new contrasts, new targets and new imaging techniques for NMD are described. The report is written by the scientists who attended the meetings and presented their data. An overview is given on the different contrasts that MRI can generate and their application, clinical needs and desired readouts, and emerging methods.

Investigating and reducing the effects of confounding factors for robust T1 and T2 mapping with cardiac MR fingerprinting

This study aims to improve the accuracy and consistency of T1 and T2 measurements using cardiac MR Fingerprinting (cMRF) by investigating and accounting for the effects of confounding factors including slice profile, inversion and T2 preparation pulse efficiency, and B1+. The goal is to understand how measurements with different pulse sequences are affected by these factors. This can be used to determine which factors must be taken into account for accurate measurements, and which may be mitigated by the selection of an appropriate pulse sequence. Simulations were performed using a numerical cardiac phantom to assess the accuracy of over 600 cMRF sequences with different flip angles, TRs, and preparation pulses. A subset of sequences, including one with the lowest errors in T1 and T2 maps, was used in subsequent analyses. Errors due to non-ideal slice profile, preparation pulse efficiency, and B1+ were quantified in Bloch simulations. Corrections for these effects were included in the dictionary generation and demonstrated in phantom and in vivo cardiac imaging at 3 T. Neglecting to model slice profile and preparation pulse efficiency led to underestimated T1 and overestimated T2 for most cMRF sequences. Sequences with smaller maximum flip angles were less affected by slice profile and B1+. Simulating all corrections in the dictionary improved the accuracy of T1 and T2 phantom measurements, regardless of acquisition pattern. More consistent myocardial T1 and T2 values were measured using different sequences after corrections. Based on these results, a pulse sequence which is minimally affected by confounding factors can be selected, and the appropriate residual corrections included for robust T1 and T2 mapping.

RF slice profile effects in magnetic resonance fingerprinting

The radio frequency (RF) slice profile effects on T1 and T2 estimation in magnetic resonance fingerprinting (MRF) are investigated with respect to time-bandwidth product (TBW), flip angle (FA) level and field inhomogeneities. Signal evolutions are generated incorporating the non-ideal slice selective excitation process using Bloch simulation and matched to the original dictionary with and without the non-ideal slice profile taken into account. For validation, phantom and in vivo experiments are performed at 3T. Both simulations and experiments results show that T1 and T2 error from non-ideal slice profile increases with increasing FA level, off-resonance, and low TBW values. Therefore, RF slice profile effects should be compensated for accurate determination of the MR parameters.

T1 mapping using saturation recovery single-shot acquisition at 3-tesla magnetic resonance imaging in hypertrophic cardiomyopathy: comparison to late gadolinium enhancement

PURPOSE

We evaluated the T1 values of segments and slices and the reproducibility in healthy controls, using saturation recovery single-shot acquisition (SASHA) at 3T magnetic resonance imaging. Moreover, we examined the difference in T1 values between hypertrophic cardiomyopathy (HCM) and healthy controls, and compared those with late gadolinium enhancement (LGE).

MATERIALS AND METHODS

Twenty-one HCM patients and 10 healthy controls underwent T1 mapping before and after contrast administration. T1 values were measured in 12 segments.

RESULTS

Native T1 values were significantly longer in HCM than in healthy controls [1373 ms (1312-1452 ms) vs. 1279 ms (1229-1326 ms); p < 0.0001]. Even in HCM segments without LGE, native T1 values were significantly longer than in healthy control segments [1366 ms (1300-1439 ms) vs. 1279 ms (1229-1326 ms); p < 0.0001]. Using a cutoff value of 1327 ms for septal native T1 values, we differentiated between HCM and healthy controls with 95% sensitivity, 90% specificity, 94% accuracy, and an area under the curve of 0.95.

CONCLUSIONS

Native T1 values using a SASHA at 3T could differentiate HCM from healthy controls. Moreover, native T1 values have the potential to detect abnormal myocardium that cannot be identified adequately by LGE in HCM.

A novel profile/view ordering with a non-convex star shutter for high-resolution 3D volumetric T1 mapping under multiple breath-holds

PURPOSE

To examine a novel non-convex star ordering/shutter for reducing the number of breath-holds in cardiac three-dimensional (3D) T₁ Mapping MRI with multiple breath-holds.

METHODS

A novel ordering, Non-Convex Star (NCS) was designed to acquire 3D volumes in a modified look-locker inversion recovery (MOLLI) T₁ mapping sequence to provide more spatial resolution and coverage in fewer breath-holds. The proposed 3D-MOLLI approach using NCS was first validated in two phantoms using artifact power (AP) measurement against the fully sampled phantom. This was followed by an in vivo study in seven swine, in which the T₁ values of the left ventricular (LV) myocardium divided into the American Heart Association (AHA) 16-segment model was compared against the reference multislice two-dimensional (2D) clinical reference and 3D volume without NCS breath-hold reduction.

RESULTS

NCS breath-hold reduction yielded less AP compared with the matched SENSE accelerated phantom volume (P < 0.0005), and was shown to be optimal at 25% fewer breath-holds. Calculated T₁ values from 3D in vivo volumes with/without NCS were comparable in all AHA segments (P = NS), whereas 3D-NCS yielded significantly higher T₁ values than 2D at midslice of the LV myocardium in each AHA segment (P < 0.05).

CONCLUSION

We successfully demonstrate the feasibility of the NCS approach for a 3D T₁ mapping acquisition requiring fewer breath-holds. Magn Reson Med 77:2215-2224, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Myocardial T1 mapping at 3.0 tesla using an inversion recovery spoiled gradient echo readout and bloch equation simulation with slice profile correction (BLESSPC) T1 estimation algorithm

BACKGROUND

To develop an accurate and precise myocardial T1 mapping technique using an inversion recovery spoiled gradient echo readout at 3.0 Tesla (T).

THEORY AND METHODS

The modified Look-Locker inversion-recovery (MOLLI) sequence was modified to use fast low angle shot (FLASH) readout, incorporating a BLESSPC (Bloch Equation Simulation with Slice Profile Correction) T1 estimation algorithm, for accurate myocardial T1 mapping. The FLASH-MOLLI with BLESSPC fitting was compared with different approaches and sequences with regards to T1 estimation accuracy, precision and image artifact based on simulation, phantom studies, and in vivo studies of 10 healthy volunteers and three patients at 3.0 Tesla.

RESULTS

The FLASH-MOLLI with BLESSPC fitting yields accurate T1 estimation (average error = -5.4 ± 15.1 ms, percentage error = -0.5% ± 1.2%) for T1 from 236-1852 ms and heart rate from 40-100 bpm in phantom studies. The FLASH-MOLLI sequence prevented off-resonance artifacts in all 10 healthy volunteers at 3.0T. In vivo, there was no significant difference between FLASH-MOLLI-derived myocardial T1 values and "ShMOLLI+IE" derived values (1458.9 ± 20.9 ms versus 1464.1 ± 6.8 ms, P = 0.50); However, the average precision by FLASH-MOLLI was significantly better than that generated by "ShMOLLI+IE" (1.84 ± 0.36% variance versus 3.57 ± 0.94%, P < 0.001).

CONCLUSION

The FLASH-MOLLI with BLESSPC fitting yields accurate and precise T1 estimation, and eliminates banding artifacts associated with bSSFP at 3.0T.

Flip-angle profile of slice-selective excitation and the measurement of the MR longitudinal relaxation time with steady-state magnetization

In MRI, the flip angle (FA) of slice-selective excitation is not uniform across the slice-thickness dimension. This work investigates the effect of the non-uniform FA profile on the accuracy of a commonly-used method for the measurement, in which the T1 value, i.e., the longitudinal relaxation time, is determined from the steady-state signals of an equally-spaced RF pulse train. By using the numerical solutions of the Bloch equation, it is shown that, because of the non-uniform FA profile, the outcome of the T1 measurement depends significantly on T1 of the specimen and on the FA and the inter-pulse spacing τ of the pulse train. A new method to restore the accuracy of the T1 measurement is described. Different from the existing approaches, the new method also removes the FA profile effect for the measurement of the FA, which is normally a part of the T1 measurement. In addition, the new method does not involve theoretical modeling, approximation, or modification to the underlying principle of the T1 measurement. An imaging experiment is performed, which shows that the new method can remove the FA-, the τ-, and the T1-dependence and produce T1 measurements in excellent agreement with the ones obtained from a gold standard method (the inversion-recovery method).

How accurate is MOLLI T1 mapping in vivo? Validation by spin echo methods

T₁ mapping is a promising quantitative tool for assessing diffuse cardiomyopathies. The purpose of this study is to quantify in vivo accuracy of the Modified Look-Locker Inversion Recovery (MOLLI) cardiac T₁ mapping sequence against the spin echo gold standard, which has not been done previously. T₁ accuracy of MOLLI was determined by comparing with the gold standard inversion recovery spin echo sequence in the calf muscle, and with a rapid inversion recovery fast spin echo sequence in the heart. T₁ values were obtained with both conventional MOLLI fitting and MOLLI fitting with inversion efficiency correction. In the calf (n = 6), conventional MOLLI fitting produced inconsistent T₁ values with error ranging from 8.0% at 90° to 17.3% at 30°. Modified MOLLI fitting with inversion efficiency correction improved error to under 7.4% at all flip angles. In the heart (n = 5), modified MOLLI fitting with inversion correction reduced T₁ error to 5.5% from 14.0% by conventional MOLLI fitting. This study shows that conventional MOLLI fitting can lead to significant in vivo T₁ errors when not accounting for the lower adiabatic inversion efficiency often experienced in vivo.

Accelerated cardiac MR stress perfusion with radial sampling after physical exercise with an MR-compatible supine bicycle ergometer

PURPOSE

To evaluate the feasibility of accelerated cardiac MR (CMR) perfusion with radial sampling using nonlinear image reconstruction after exercise on an MR-compatible supine bike ergometer.

METHODS

Eight healthy subjects were scanned on two separate days using radial and Cartesian CMR perfusion sequences in rest and exercise stress perfusion. Four different methods (standard gridding, conjugate gradient SENSE [CG-SENSE], nonlinear inversion with joint estimation of coil-sensitivity profiles [NLINV] and compressed sensing with a total variation constraint [TV]) were compared for the reconstruction of radial data. Cartesian data were reconstructed using SENSE. All images were assessed by two blinded readers in terms of image quality and diagnostic value.

RESULTS

CG-SENSE and NLINV were scored more favorably than TV (in both rest and stress perfusion cases, P < 0.05) and gridding (for rest perfusion cases, P < 0.05). TV images showed patchy artifacts, which negatively influenced image quality especially in the stress perfusion images acquired with a low number of radial spokes. Although CG-SENSE and NLINV received better scores than Cartesian sampling in both rest and exercise stress perfusion cases, these differences were not statistically significant (P > 0.05).

CONCLUSION

We have demonstrated the feasibility of accelerated CMR perfusion using radial sampling after physical exercise using a supine bicycle ergometer in healthy subjects. For reconstruction of undersampled radial perfusion, CG-SENSE and NLINV resulted in better image quality than standard gridding or TV reconstruction. Further technical improvements and clinical assessment are needed before using this approach in patients with suspected coronary artery disease.

Consideration of slice profiles in inversion recovery Look-Locker relaxation parameter mapping

PURPOSE

To include the flip angle distribution caused by the slice profile into the model used for describing the relaxation curves observed in inversion recovery Look-Locker FLASH T1 mapping for a more accurate determination of the relaxation parameters.

MATERIALS AND METHODS

For each inversion time, the flip angle dependent signal of the mono-exponential relaxation model is integrated across the slice profile. The resulting Consideration of Slice Profiles (CSP) relaxation curves are compared to the mono-exponential signal model in numerical simulations as well as in phantom and in-vivo experiments.

RESULTS

All measured relaxation curves showed systematic deviations from a mono-exponential curve increasing with flip angle and T1 but decreasing with repetition time. Additionally, the accuracy of T1 was found to be largely dependent on the temporal coverage of the relaxation curve. All these systematic errors were largely reduced by the CSP model.

CONCLUSION

The proposed CSP model represents a useful extension of the conventionally used mono-exponential relaxation model. Despite inherent model inaccuracies, the mono-exponential model was found to be sufficient for many T1 mapping situations. However, if only a poor temporal coverage of the relaxation process is achievable or a very precise modeling of the relaxation course is needed as in model-based techniques, the mono-exponential model leads to systematic errors and the CSP model should be used instead.

Fast in vivo quantification of T1 and T2 MRI relaxation times in the myocardium based on inversion recovery SSFP with in vitro validation post Gd-based contrast administration

PURPOSE

A fast clinical imaging technique for quantifying myocardial T1 and T2 relaxation times after Gadolinium (Gd)-based contrast administration within a single breathhold is presented with in vitro validation.

MATERIALS AND METHODS

From signal intensity curves in ECG-gated segmented inversion recovery balanced steady state free precession (IR-bSSFP) images, T1 and T2 values were determined for 24 agarose samples made from solutions of Omniscan (0.25-2 mg/mL) and copper-sulfate (0.52-22.17 mg/mL). T1 and T2 were also measured using turbo spin-echo (TSE) acquisitions and compared with IR-bSSFP results. In vivo T1 and T2 values from post-contrast IR-bSSFP images of five healthy volunteers were determined for (I) the left ventricular wall, (II) the interventricular septum (IVS) and (III) the lateral wall of the left ventricle (LV). Spin system simulations were performed for selected T1 and T2 values.

RESULTS

Good agreement between TSE and IR-bSSFP for T1 for realistic in vivo post-contrast values (below 1,250 ms, R=0.88) and for T2 (entire range, R=0.97) was found. Spin system simulations were in good agreement with measurements. In vivo average T1 was 546±32 ms and average T2 was 59±9 ms.

CONCLUSIONS

A fast imaging protocol for absolute quantification of myocardial T1 and T2 post-contrast is presented, validated in vitro and consecutively applied in vivo in humans.

Accuracy, precision, and reproducibility of four T1 mapping sequences: a head-to-head comparison of MOLLI, ShMOLLI, SASHA, and SAPPHIRE

PURPOSE

To compare accuracy, precision, and reproducibility of four commonly used myocardial T1 mapping sequences: modified Look-Locker inversion recovery (MOLLI), shortened MOLLI (ShMOLLI), saturation recovery single-shot acquisition (SASHA), and saturation pulse prepared heart rate independent inversion recovery (SAPPHIRE).

MATERIALS AND METHODS

This HIPAA-compliant study was approved by the institutional review board. All subjects provided written informed consent. Accuracy, precision, and reproducibility of the four T1 mapping sequences were first compared in phantom experiments. In vivo analysis was performed in seven healthy subjects (mean age ± standard deviation, 38 years ± 19; four men, three women) who were imaged twice on two separate days. In vivo reproducibility of native T1 mapping and extracellular volume (ECV) were measured. Differences between the sequences were assessed by using Kruskal-Wallis and Wilcoxon rank sum tests (phantom data) and mixed-effect models (in vivo data).

RESULTS

T1 mapping accuracy in phantoms was lower with ShMOLLI (62 msec) and MOLLI (44 msec) than with SASHA (13 msec; P < .05) and SAPPHIRE (12 msec; P < .05). MOLLI had similar precision to ShMOLLI (4.0 msec vs 5.6 msec; P = .07) but higher precision than SAPPHIRE (6.8 msec; P = .002) and SASHA (8.7 msec; P < .001). All sequences had similar reproducibility in phantoms (P = .1). The four sequences had similar in vivo reproducibility for native T1 mapping (∼25-50 msec; P > .05) and ECV quantification (∼0.01-0.02; P > .05).

CONCLUSION

SASHA and SAPPHIRE yield higher accuracy, lower precision, and similar reproducibility compared with MOLLI and ShMOLLI for T1 measurement. Different sequences yield different ECV values; however, all sequences have similar reproducibility for ECV quantification.