Myocardial T1-mapping at 3T using saturation-recovery: reference values, precision and comparison with MOLLI

Color Maps: Facilitating the Clinical Impact of Quantitative MRI

Presenting quantitative data using non-standardized color maps potentially results in unrecognized misinterpretation of data. Clinically meaningful color maps should intuitively and inclusively represent data without misleading interpretation. Uniformity of the color gradient for color maps is critically important. Maximal color and lightness contrast, readability for color vision-impaired individuals, and recognizability of the color scheme are highly desirable features. This article describes the use of color maps in five key quantitative MRI techniques: relaxometry, diffusion-weighted imaging (DWI), dynamic contrast-enhanced (DCE)-MRI, MR elastography (MRE), and water-fat MRI. Current display practice of color maps is reviewed and shortcomings against desirable features are highlighted. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 2.

Clinical Applications of Cardiac Magnetic Resonance Parametric Mapping

Cardiovascular magnetic resonance (CMR) imaging is widely regarded as the gold-standard technique for myocardial tissue characterization, allowing for the detection of structural abnormalities such as myocardial fatty replacement, myocardial edema, myocardial necrosis, and/or fibrosis. Historically, the identification of abnormal myocardial regions relied on variations in tissue signal intensity, often necessitating the use of exogenous contrast agents. However, over the past two decades, innovative parametric mapping techniques have emerged, enabling the direct quantitative assessment of tissue magnetic resonance (MR) properties on a voxel-by-voxel basis. These mapping techniques offer significant advantages by providing comprehensive and precise information that can be translated into color-coded maps, facilitating the identification of subtle or diffuse myocardial abnormalities. As unlikely conventional methods, these techniques do not require a substantial amount of structurally altered tissue to be visually identifiable as an area of abnormal signal intensity, eliminating the reliance on contrast agents. Moreover, these parametric mapping techniques, such as T1, T2, and T2* mapping, have transitioned from being primarily research tools to becoming valuable assets in the clinical diagnosis and risk stratification of various cardiac disorders. In this review, we aim to elucidate the underlying physical principles of CMR parametric mapping, explore its current clinical applications, address potential pitfalls, and outline future directions for research and development in this field.

The diagnostic performance of T1 mapping in the assessment of breast lesions: A preliminary study

PURPOSE

To assess T1 mapping performance in distinguishing between benign and malignant breast lesions and to explore its correlation with histopathologic features in breast cancer.

METHODS

This study prospectively enrolled 103 participants with a total of 108 lesions, including 25 benign and 83 malignant lesions. T1 mapping, diffusion-weighted imaging (DWI), and dynamic contrast-enhanced (DCE) were performed. Two radiologists independently outlined the ROIs and analyzed T1 and apparent diffusion coefficient (ADC) values for each lesion, assessing interobserver reliability with the intraclass correlation coefficient (ICC). T1 and ADC values were compared between benign and malignant lesions, across different histopathological characteristics (histological grades, estrogen, progesterone and HER2 receptors expression, Ki67, N status). Receiver operating characteristic (ROC) analysis and Pearson correlation coefficient (ρ) were performed.

RESULTS

T1 values showed statistically significant differences between benign and malignant groups (P < 0.001), with higher values in the malignant (1817.08 ms ± 126.64) compared to the benign group (1429.31 ms ± 167.66). In addition, T1 values significantly increased in the ER (-) group (P = 0.001). No significant differences were found in T1 values among HER2, Ki67, N status, and histological grades groups. Furthermore, T1 values exhibited a significant correlation (ρ) with ER (P < 0.01) and PR (P = 0.03). The AUC for T1 value in distinguishing benign from malignant lesions was 0.69 (95 % CI: 0.55 - 0.82, P = 0.005), and for evaluating ER status, it was 0.75 (95 % CI: 0.62 - 0.87, P = 0.002).

CONCLUSIONS

T1 mapping holds the potential as an imaging biomarker to assist in the discrimination of benign and malignant breast lesions and assessing the ER expression status in breast cancer.

Virtual MOLLI Target: Generative Adversarial Networks Toward Improved Motion Correction in MRI Myocardial T1 Mapping

BACKGROUND

The modified Look-Locker inversion recovery (MOLLI) sequence is commonly used for myocardial T1 mapping. However, it acquires images with different inversion times, which causes difficulty in motion correction for respiratory-induced misregistration to a given target image.

HYPOTHESIS

Using a generative adversarial network (GAN) to produce virtual MOLLI images with consistent heart positions can reduce respiratory-induced misregistration of MOLLI datasets.

STUDY TYPE

Retrospective.

POPULATION

1071 MOLLI datasets from 392 human participants.

FIELD STRENGTH/SEQUENCE

Modified Look-Locker inversion recovery sequence at 3 T.

ASSESSMENT

A GAN model with a single inversion time image as input was trained to generate virtual MOLLI target (VMT) images at different inversion times which were subsequently used in an image registration algorithm. Four VMT models were investigated and the best performing model compared with the standard vendor-provided motion correction (MOCO) technique.

STATISTICAL TESTS

The effectiveness of the motion correction technique was assessed using the fitting quality index (FQI), mutual information (MI), and Dice coefficients of motion-corrected images, plus subjective quality evaluation of T1 maps by three independent readers using Likert score. Wilcoxon signed-rank test with Bonferroni correction for multiple comparison. Significance levels were defined as P < 0.01 for highly significant differences and P < 0.05 for significant differences.

RESULTS

The best performing VMT model with iterative registration demonstrated significantly better performance (FQI 0.88 ± 0.03, MI 1.78 ± 0.20, Dice 0.84 ± 0.23, quality score 2.26 ± 0.95) compared to other approaches, including the vendor-provided MOCO method (FQI 0.86 ± 0.04, MI 1.69 ± 0.25, Dice 0.80 ± 0.27, quality score 2.16 ± 1.01).

DATA CONCLUSION

Our GAN model generating VMT images improved motion correction, which may assist reliable T1 mapping in the presence of respiratory motion. Its robust performance, even with considerable respiratory-induced heart displacements, may be beneficial for patients with difficulties in breath-holding.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 1.

STEAM-SASHA: a novel approach for blood- and fat-suppressed native T1 measurement in the right ventricular myocardium

OBJECTIVE

The excellent blood and fat suppression of stimulated echo acquisition mode (STEAM) can be combined with saturation recovery single-shot acquisition (SASHA) in a novel STEAM-SASHA sequence for right ventricular (RV) native T1 mapping.

MATERIALS AND METHODS

STEAM-SASHA splits magnetization preparation over two cardiac cycles, nulling blood signal and allowing fat signal to decay. Breath-hold T1 mapping was performed in a T1 phantom and twice in 10 volunteers using STEAM-SASHA and a modified Look-Locker sequence at peak systole at 3T. T1 was measured in 3 RV regions, the septum and left ventricle (LV).

RESULTS

In phantoms, MOLLI under-estimated while STEAM-SASHA over-estimated T1, on average by 3.0% and 7.0% respectively, although at typical 3T myocardial T1 (T1 > 1200 ms) STEAM-SASHA was more accurate. In volunteers, T1 was higher using STEAM-SASHA than MOLLI in the LV and septum (p = 0.03, p = 0.006, respectively), but lower in RV regions (p > 0.05). Inter-study, inter-observer and intra-observer coefficients of variation in all regions were < 15%. Blood suppression was excellent with STEAM-SASHA and noise floor effects were minimal.

DISCUSSION

STEAM-SASHA provides accurate and reproducible T1 in the RV with excellent blood and fat suppression. STEAM-SASHA has potential to provide new insights into pathological changes in the RV in future studies.

Prognostic value of native T1 and extracellular volume in patients with immunoglubin light-chain amyloidosis

BACKGROUND

Cardiac involvement in patients with immunoglubin light-chain amyloidosis (AL) is a major determinant of treatment choice and prognosis, and early identification of high-risk patients can initiate intensive treatment strategies to achieve better survival. This study aimed to investigate the prognostic value of native T1 and ECV in patients with AL-cardiac amyloidosis (CA).

METHODS

A total of 38 patients (mean age 59 ± 11 years) with AL diagnosed histopathologically from July 2017 to October 2021 were collected consecutively. All patients were performed 3.0-T cardiac magnetic resonance (CMR) including cine, T1 mapping, and late gadolinium enhancement (LGE). Pre- and post-contrast T1 mapping images were transferred to a dedicated research software package (CVI42 v5.11.3) to create parametric T1 and ECV values. In addition, clinical and laboratory data of all patients were collected, and patients or their family members were regularly followed up by telephone every 3 months. The starting point of follow-up was the time of definitive pathological diagnosis, and the main endpoint was all-cause death. Kaplan-Meier analysis and Cox proportional risk model were used to evaluate the association between native T1 and ECV and death in patients with CA.

RESULTS

After a median follow-up of 27 (16, 37) months, 12 patients with CA died. Kaplan-Meier analysis showed that elevated native T1 and ECV were closely associated with poor prognosis in patients with CA. The survival rate of patients with ECV > 44% and native T1 > 1389ms were significantly lower than that of patients with ECV ≤ 44% and native T1 ≤ 1389ms (Log-rank P < 0.001), and was not associated with the presence of LGE. After adjusting for clinical risk factors and CMR measurements in a stepwise multivariate Cox regression model, ECV [risk ratio (HR):1.37, 95%CI: 1.09-1.73, P = 0.008] and native T1 (HR:1.01, 95%CI: 1.00-1.02, P = 0.037) remained independent predictors of all-cause mortality in patients with CA.

CONCLUSIONS

Both native T1 and ECV were independently prognostic for mortality in patients with CA, and can be used as important indicators for clinical prognosis assessment of AL.

T1 mapping of myocardium in rats using self-gated golden-angle acquisition

PURPOSE

The aim of this study is to design a method of myocardial T1 quantification in small laboratory animals and to investigate the effects of spatiotemporal regularization and the needed acquisition duration.

METHODS

We propose a compressed-sensing approach to T1 quantification based on self-gated inversion-recovery radial two/three-dimensional (2D/3D) golden-angle stack-of-stars acquisition with image reconstruction performed using total-variation spatiotemporal regularization. The method was tested on a phantom and on a healthy rat, as well as on rats in a small myocardium-remodeling study.

RESULTS

The results showed a good match of the T1 estimates with the results obtained using the ground-truth method on a phantom and with the literature values for rats myocardium. The proposed 2D and 3D methods showed significant differences between normal and remodeling myocardium groups for acquisition lengths down to approximately 5 and 15 min, respectively.

CONCLUSIONS

A new 2D and 3D method for quantification of myocardial T1 in rats was proposed. We have shown the capability of both techniques to distinguish between normal and remodeling myocardial tissue. We have shown the effects of image-reconstruction regularization weights and acquisition length on the T1 estimates.

Improved reproducibility for myocardial ASL: Impact of physiological and acquisition parameters

PURPOSE

To investigate and mitigate the influence of physiological and acquisition-related parameters on myocardial blood flow (MBF) measurements obtained with myocardial Arterial Spin Labeling (myoASL).

METHODS

A Flow-sensitive Alternating Inversion Recovery (FAIR) myoASL sequence with bSSFP and spoiled GRE (spGRE) readout is investigated for MBF quantification. Bloch-equation simulations and phantom experiments were performed to evaluate how variations in acquisition flip angle (FA), acquisition matrix size (AMS), heart rate (HR) and bloodT 1 $$ {\mathrm{T}}_1 $$ relaxation time (T 1 , B $$ {\mathrm{T}}_{1,B} $$ ) affect quantification of myoASL-MBF. In vivo myoASL-images were acquired in nine healthy subjects. A corrected MBF quantification approach was proposed based on subject-specificT 1 , B $$ {\mathrm{T}}_{1,B} $$ values and, for spGRE imaging, subtracting an additional saturation-prepared baseline from the original baseline signal.

RESULTS

Simulated and phantom experiments showed a strong dependence on AMS and FA (R 2 $$ {R}^2 $$ >0.73), which was eliminated in simulations and alleviated in phantom experiments using the proposed saturation-baseline correction in spGRE. Only a very mild HR dependence (R 2 $$ {R}^2 $$ >0.59) was observed which was reduced when calculating MBF with individualT 1 , B $$ {\mathrm{T}}_{1,B} $$ . For corrected spGRE, in vivo mean global spGRE-MBF ranged from 0.54 to 2.59 mL/g/min and was in agreement with previously reported values. Compared to uncorrected spGRE, the intra-subject variability within a measurement (0.60 mL/g/min), between measurements (0.45 mL/g/min), as well as the inter-subject variability (1.29 mL/g/min) were improved by up to 40% and were comparable with conventional bSSFP.

CONCLUSION

Our results show that physiological and acquisition-related factors can lead to spurious changes in myoASL-MBF if not accounted for. Using individualT 1 , B $$ {\mathrm{T}}_{1,B} $$ and a saturation-baseline can reduce these variations in spGRE and improve reproducibility of FAIR-myoASL against acquisition parameters.

Multi-site comparison of parametric T1 and T2 mapping: healthy travelling volunteers in the Berlin research network for cardiovascular magnetic resonance (BER-CMR)

BACKGROUND

Parametric mapping sequences in cardiovascular magnetic resonance (CMR) allow for non-invasive myocardial tissue characterization. However quantitative myocardial mapping is still limited by the need for local reference values. Confounders, such as field strength, vendors and sequences, make intersite comparisons challenging. This exploratory study aims to assess whether multi-site studies that control confounding factors provide first insights whether parametric mapping values are within pre-defined tolerance ranges across scanners and sites.

METHODS

A cohort of 20 healthy travelling volunteers was prospectively scanned at three sites with a 3 T scanner from the same vendor using the same scanning protocol and acquisition scheme. A Modified Look-Locker inversion recovery sequence (MOLLI) for T1 and a fast low-angle shot sequence (FLASH) for T2 were used. At one site a scan-rescan was performed to assess the intra-scanner reproducibility. All acquired T1- and T2-mappings were analyzed in a core laboratory using the same post-processing approach and software.

RESULTS

After exclusion of one volunteer due to an accidentally diagnosed cardiac disease, T1- and T2-maps of 19 volunteers showed no significant differences between the 3 T sites (mean ± SD [95% confidence interval] for global T1 in ms: site I: 1207 ± 32 [1192-1222]; site II: 1207 ± 40 [1184-1225]; site III: 1219 ± 26 [1207-1232]; p = 0.067; for global T2 in ms: site I: 40 ± 2 [39-41]; site II: 40 ± 1 [39-41]; site III 39 ± 2 [39-41]; p = 0.543).

CONCLUSION

Parametric mapping results displayed initial hints at a sufficient similarity between sites when confounders, such as field strength, vendor diversity, acquisition schemes and post-processing analysis are harmonized. This finding needs to be confirmed in a powered clinical trial. Trial registration ISRCTN14627679 (retrospectively registered).

Improved Simultaneous Multi-slice imaging with Composition of k-space Interpolations (SMS-COOKIE) for myocardial T1 mapping

The aim of this study is to develop and evaluate a regularized Simultaneous Multi-Slice (SMS) reconstruction method for improved Cardiac Magnetic Resonance Imaging (CMR). The proposed reconstruction method, SMS with COmpOsition of k-space IntErpolations (SMS-COOKIE) combines the advantages of Iterative Self-consistent Parallel Imaging Reconstruction (SPIRiT) and split slice-Generalized Autocalibrating Partially Parallel Acquisitions (GRAPPA), while allowing regularization for further noise reduction. The proposed SMS-COOKIE was implemented with and without regularization, and validated using a Saturation Pulse-Prepared Heart rate Independent inversion REcovery (SAPPHIRE) myocardial T1 mapping sequence. The performance of the proposed reconstruction method was compared to ReadOut (RO)-SENSE-GRAPPA and split slice-GRAPPA, on both retrospectively and prospectively three-fold SMS-accelerated data with an additional two-fold in-plane acceleration. All SMS reconstruction methods yielded similar T1 values compared to single band imaging. SMS-COOKIE showed lower spatial variability in myocardial T1 with significant improvement over RO-SENSE-GRAPPA and split slice-GRAPPA (P < 10-4). The proposed method with additional locally low rank (LLR) regularization reduced the spatial variability, again with significant improvement over RO-SENSE-GRAPPA and split slice-GRAPPA (P < 10-4). In conclusion, improved reconstruction quality was achieved with the proposed SMS-COOKIE, which also provided lower spatial variability with significant improvement over split slice-GRAPPA.

Accelerated cardiac T1 mapping with recurrent networks and cyclic, model-based loss

BACKGROUND

Using the spin-lattice relaxation time (T1) as a biomarker, the myocardium can be quantitatively characterized using cardiac T1 mapping. The modified Look-Locker inversion (MOLLI) recovery sequences have become the standard clinical method for cardiac T1 mapping. However, the MOLLI sequences require an 11-heartbeat breath-hold that can be difficult for subjects, particularly during exercise or pharmacologically induced stress. Although shorter cardiac T1 mapping sequences have been proposed, these methods suffer from reduced precision. As such, there is an unmet need for accelerated cardiac T1 mapping.

PURPOSE

To accelerate cardiac T1 mapping MOLLI sequences by using neural networks to estimate T1 maps using a reduced number of T1-weighted images and their corresponding inversion times.

MATERIALS AND METHODS

In this retrospective study, 911 pre-contrast T1 mapping datasets from 202 subjects (128 males, 56 ± 15 years; 74 females, 54 ± 17 years) and 574 T1 mapping post-contrast datasets from 193 subjects (122 males, 57 ± 15 years; 71 females, 54 ± 17 years) were acquired using the MOLLI-5(3)3 sequence and the MOLLI-4(1)3(1)2 sequence, respectively. All acquisition protocols used similar scan parameters:T R = 2.2 ms $TR\; = \;2.2\;{\rm{ms}}$ ,T E = 1.12 ms $TE\; = \;1.12\;{\rm{ms}}$ , andF A = 35 ∘ $FA\; = \;35^\circ $ , gadoteridol (ProHance, Bracco Diagnostics) dose∼ 0.075 mmol / kg $\sim 0.075\;\;{\rm{mmol/kg}}$ . A bidirectional multilayered long short-term memory (LSTM) network with fully connected output and cyclic model-based loss was used to estimate T1 maps from the first three T1-weighted images and their corresponding inversion times for pre- and post-contrast T1 mapping. The performance of the proposed architecture was compared to the three-parameter T1 recovery model using the same reduction of the number of T1-weighted images and inversion times. Reference T1 maps were generated from the scanner using the full MOLLI sequences and the three-parameter T1 recovery model. Correlation and Bland-Altman plots were used to evaluate network performance in which each point represents averaged regions of interest in the myocardium corresponding to the standard American Heart Association 16-segment model. The precision of the network was examined using consecutively repeated scans. Stress and rest pre-contrast MOLLI studies as well as various disease test cases, including amyloidosis, hypertrophic cardiomyopathy, and sarcoidosis were also examined. Paired t-tests were used to determine statistical significance withp < 0.05 $p < 0.05$ .

RESULTS

Our proposed network demonstrated similar T1 estimations to the standard MOLLI sequences (pre-contrast:1260 ± 94 ms $1260 \pm 94\;{\rm{ms}}$ vs.1254 ± 91 ms $1254 \pm 91\;{\rm{ms}}$ withp = 0.13 $p\; = \;0.13$ ; post-contrast:484 ± 92 ms $484 \pm 92\;{\rm{ms}}$ vs.493 ± 91 ms $493 \pm 91\;{\rm{ms}}$ withp = 0.07 $p\; = \;0.07$ ). The precision of standard MOLLI sequences was well preserved with the proposed network architecture (24 ± 28 ms $24 \pm 28\;\;{\rm{ms}}$ vs.18 ± 13 ms $18 \pm 13\;{\rm{ms}}$ ). Network-generated T1 reactivities are similar to stress and rest pre-contrast MOLLI studies (5.1 ± 4.0 % $5.1 \pm 4.0\;\% $ vs.4.9 ± 4.4 % $4.9 \pm 4.4\;\% $ withp = 0.84 $p\; = \;0.84$ ). Amyloidosis T1 maps generated using the proposed network are also similar to the reference T1 maps (pre-contrast:1243 ± 140 ms $1243 \pm 140\;\;{\rm{ms}}$ vs.1231 ± 137 ms $1231 \pm 137\;{\rm{ms}}$ withp = 0.60 $p\; = \;0.60$ ; post-contrast:348 ± 26 ms $348 \pm 26\;{\rm{ms}}$ vs.346 ± 27 ms $346 \pm 27\;{\rm{ms}}$ withp = 0.89 $p\; = \;0.89$ ).

CONCLUSIONS

A bidirectional multilayered LSTM network with fully connected output and cyclic model-based loss was used to generate high-quality pre- and post-contrast T1 maps using the first three T1-weighted images and their corresponding inversion times. This work demonstrates that combining deep learning with cardiac T1 mapping can potentially accelerate standard MOLLI sequences from 11 to 3 heartbeats.

Cardiac phase-resolved late gadolinium enhancement imaging

Late gadolinium enhancement (LGE) with cardiac magnetic resonance (CMR) imaging is the clinical reference for assessment of myocardial scar and focal fibrosis. However, current LGE techniques are confined to imaging of a single cardiac phase, which hampers assessment of scar motility and does not allow cross-comparison between multiple phases. In this work, we investigate a three step approach to obtain cardiac phase-resolved LGE images: (1) Acquisition of cardiac phase-resolved imaging data with varying T₁ weighting. (2) Generation of semi-quantitative T1* maps for each cardiac phase. (3) Synthetization of LGE contrast to obtain functional LGE images. The proposed method is evaluated in phantom imaging, six healthy subjects at 3T and 20 patients at 1.5T. Phantom imaging at 3T demonstrates consistent contrast throughout the cardiac cycle with a coefficient of variation of 2.55 ± 0.42%. In-vivo results show reliable LGE contrast with thorough suppression of the myocardial tissue is healthy subjects. The contrast between blood and myocardium showed moderate variation throughout the cardiac cycle in healthy subjects (coefficient of variation 18.2 ± 3.51%). Images were acquired at 40–60 ms and 80 ms temporal resolution, at 3T and 1.5, respectively. Functional LGE images acquired in patients with myocardial scar visualized scar tissue throughout the cardiac cycle, albeit at noticeably lower imaging resolution and noise resilience than the reference technique. The proposed technique bears the promise of integrating the advantages of phase-resolved CMR with LGE imaging, but further improvements in the acquisition quality are warranted for clinical use.

Clinical evaluation of the Multimapping technique for simultaneous myocardial T1 and T2 mapping

The Multimapping technique was recently proposed for simultaneous myocardial T₁ and T₂ mapping. In this study, we evaluate its correlation with clinical reference mapping techniques in patients with a range of cardiovascular diseases (CVDs) and compare image quality and inter- and intra-observer repeatability. Multimapping consists of an ECG-triggered, 2D single-shot bSSFP readout with inversion recovery and T₂ preparation modules, acquired across 10 cardiac cycles. The sequence was implemented at 1.5T and compared to clinical reference mapping techniques, modified Look-Locker inversion recovery (MOLLI) and T₂ prepared bSSFP with four echo times (T₂bSSFP), and compared in 47 patients with CVD (of which 44 were analyzed). In diseased myocardial segments (defined as the presence of late gadolinium enhancement), there was a high correlation between Multimapping and MOLLI for native myocardium T₁ (r² = 0.73), ECV (r² = 0.91), and blood T₁ (r² = 0.88), and Multimapping and T₂bSSFP for native myocardial T₂ (r² = 0.80). In healthy myocardial segments, a bias for native T₁ (Multimapping = 1,116 ± 21 ms, MOLLI = 1,002 ± 21, P < 0.001), post-contrast T₁ (Multimapping = 479 ± 31 ms, MOLLI = 426 ± 27 ms, 0.001), ECV (Multimapping = 21.5 ± 1.9%, MOLLI = 23.7 ± 2.3%, P = 0.001), and native T₂ (Multimapping = 48.0 ± 3.0 ms, T₂bSSFP = 53.9 ± 3.5 ms, P < 0.001) was observed. The image quality for Multimapping was scored as higher for all mapping techniques (native T₁, post-contrast T₁, ECV, and T₂bSSFP) compared to the clinical reference techniques. The inter- and intra-observer agreements were excellent (intraclass correlation coefficient, ICC > 0.9) for most measurements, except for inter-observer repeatability of Multimapping native T₁ (ICC = 0.87), post-contrast T₁ (ICC = 0.73), and T₂bSSFP native T₂ (ICC = 0.88). Multimapping shows high correlations with clinical reference mapping techniques for T₁, T₂, and ECV in a diverse cohort of patients with different cardiovascular diseases. Multimapping enables simultaneous T₁ and T₂ mapping and can be performed in a short breath-hold, with image quality superior to that of the clinical reference techniques.

Improving accuracy of myocardial T1 estimation in MyoMapNet

PURPOSE

To improve the accuracy and robustness of T₁ estimation by MyoMapNet, a deep learning-based approach using 4 inversion-recovery T₁ -weighted images for cardiac T₁ mapping.

METHODS

MyoMapNet is a fully connected neural network for T₁ estimation of an accelerated cardiac T₁ mapping sequence, which collects 4 T₁ -weighted images by a single Look-Locker inversion-recovery experiment (LL4). MyoMapNet was originally trained using in vivo data from the modified Look-Locker inversion recovery sequence, which resulted in significant bias and sensitivity to various confounders. This study sought to train MyoMapNet using signals generated from numerical simulations and phantom MR data under multiple simulated confounders. The trained model was then evaluated by phantom data scanned using new phantom vials that differed from those used for training. The performance of the new model was compared with modified Look-Locker inversion recovery sequence and saturation-recovery single-shot acquisition for measuring native and postcontrast T₁ in 25 subjects.

RESULTS

In the phantom study, T₁ values measured by LL4 with MyoMapNet were highly correlated with reference values from the spin-echo sequence. Furthermore, the estimated T₁ had excellent robustness to changes in flip angle and off-resonance. Native and postcontrast myocardium T₁ at 3 Tesla measured by saturation-recovery single-shot acquisition, modified Look-Locker inversion recovery sequence, and MyoMapNet were 1483 ± 46.6 ms and 791 ± 45.8 ms, 1169 ± 49.0 ms and 612 ± 36.0 ms, and 1443 ± 57.5 ms and 700 ± 57.5 ms, respectively. The corresponding extracellular volumes were 22.90% ± 3.20%, 28.88% ± 3.48%, and 30.65% ± 3.60%, respectively.

CONCLUSION

Training MyoMapNet with numerical simulations and phantom data will improve the estimation of myocardial T₁ values and increase its robustness to confounders while also reducing the overall T₁ mapping estimation time to only 4 heartbeats.

Saturation-pulse prepared heart-rate independent inversion-recovery (SAPPHIRE) biventricular T1 mapping: inter-field strength, head-to-head comparison of diastolic, systolic and dark-blood measurements

Background

To assess the feasibility of biventricular SAPPHIRE T₁ mapping in vivo across field strengths using diastolic, systolic and dark-blood (DB) approaches.

Methods

10 healthy volunteers underwent same-day non-contrast cardiovascular magnetic resonance at 1.5 Tesla (T) and 3 T. Left and right ventricular (LV, RV) T₁ mapping was performed in the basal, mid and apical short axis using 4-variants of SAPPHIRE: diastolic, systolic, 0th and 2nd order motion-sensitized DB and conventional modified Look-Locker inversion recovery (MOLLI).

Results

LV global myocardial T₁ times (1.5 T then 3 T results) were significantly longer by diastolic SAPPHIRE (1283 ± 11|1600 ± 17 ms) than any of the other SAPPHIRE variants: systolic (1239 ± 9|1595 ± 13 ms), 0th order DB (1241 ± 10|1596 ± 12) and 2nd order DB (1251 ± 11|1560 ± 20 ms, all p < 0.05). In the mid septum MOLLI and diastolic SAPPHIRE exhibited significant T₁ signal contamination (longer T₁) at the blood-myocardial interface not seen with the other 3 SAPPHIRE variants (all p < 0.025). Additionally, systolic, 0th order and 2nd order DB SAPPHIRE showed narrower dispersion of myocardial T₁ times across the mid septum when compared to diastolic SAPPHIRE (interquartile ranges respectively: 25 ms, 71 ms, 73 ms vs 143 ms, all p < 0.05). RV T₁ mapping was achievable using systolic, 0th and 2nd order DB SAPPHIRE but not with MOLLI or diastolic SAPPHIRE. All 4 SAPPHIRE variants showed excellent re-read reproducibility (intraclass correlation coefficients 0.953 to 0.996).

Conclusion

These small-scale preliminary healthy volunteer data suggest that DB SAPPHIRE has the potential to reduce partial volume effects at the blood-myocardial interface, and that systolic SAPPHIRE could be a feasible solution for right ventricular T₁ mapping. Further work is needed to understand the robustness of these sequences and their potential clinical utility.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12880-022-00843-0.

Diffuse Myocardial Fibrosis at Cardiac MRI in Young Adults Born Prematurely: A Cross-sectional Cohort Study

Purpose

To measure native T1 values, a marker of diffuse fibrosis, by using cardiac MRI (CMR) in young adults born prematurely.

Materials and Methods

This secondary analysis of a prospective cohort study included young adults born moderately to extremely preterm and age-matched, term-born participants. CMR was performed with a 3.0-T imager that included cine imaging for the quantification of left ventricular (LV) and right ventricular (RV) volumes and function and native saturation recovery T1 mapping for the assessment of diffuse myocardial fibrosis. Values between preterm and term were compared by using the Student t test. Associations between T1 values and other variables were analyzed by using linear regression and multivariate regression.

Results

Of the 50 young-adult participants, 32 were born preterm (mean age, 25.8 years ± 4.2 [SD]; 23 women) and 18 were born at term (mean age, 26.2 years ± 5.4; 10 women). Native T1 values were significantly higher in participants born preterm than in participants born at term (1477 msec ± 77 vs 1423 msec ± 71, respectively; unadjusted P = .0019). Native T1 values appeared to be positively associated with indexed LV end-diastolic and end-systolic volumes (β = 2.1, standard error = 0.7 and β = 3.8, standard error = 1.2, respectively), the RV end-diastolic volume index (β = 1.3, standard error = 0.6), and the LV mass index (β = 2.5, standard error = 0.9). Higher T1 values may be associated with reduced cardiac systolic strain measures and diastolic strain measures. Five-minute Apgar scores were inversely associated with native T1 values.

Conclusion

Young adults born moderately to extremely preterm exhibited significantly higher native T1 values than age-matched, term-born young adults.

Keywords: MRI, Cardiac, Heart, Left Ventricle, Cardiomyopathies

Clinical trial registration no. NCT03245723

Published under a CC BY 4.0 license

Supplemental material is available for this article.

Case Report: Transient Increase of CMR T1 Mapping Indices in a Patient With COVID-19 mRNA Vaccine Induced Acute Myocarditis

Background

Acute myocarditis is commonly associated with viral infections, including severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Myocarditis following mRNA COVID-19 vaccination has also been reported, however this is rare and usually resolves within days or weeks. We present a case of acute myocarditis reported after vaccination with mRNA-1273 COVID-19 vaccine (Moderna) diagnosed using cardiac magnetic resonance imaging (CMR). This report describes the utility of CMR in the diagnosis and follow-up of such patients using parameters which could suggest the clinical course of myocarditis.

Case Summary

A 23-year-old male presented in the emergency department with complaints of chest pain radiating to the left arm following vaccination with the second dose of COVID-19 mRNA-1273 vaccine (Moderna). Patient's history revealed an incidence of myocarditis in the past. CMR showed a mid-range left ventricular ejection fraction (38%) and subepicardial late gadolinium enhancement (LGE) in the inferolateral and apical myocardial segments with diffuse elevation of native T1 mapping relaxation times in all myocardial segments. The patient was admitted briefly in the intensive care unit and after a favorable clinical course was discharged from the hospital in stable condition. A follow-up CMR after 3 months revealed normalization of LVEF (57%) and native T1- times in most segments. Scarred myocardium reflecting chronic myocarditis continued to show elevated T1 times.

Conclusions

Our patient presenting with acute myocarditis after recent COVID-19 mRNA vaccination reported a favorable clinical course. CMR revealed increased T1 mapping relaxation times diffusely spread across the myocardium and an impairment of the left ventricular function (LVEF) during the acute phase. However, the LVEF as well as the T1 times normalized at follow-up in all segments except for myocardium affected by chronic myocarditis.

Simultaneous Multi-Slice Cardiac MR Multitasking for Motion-Resolved, Non-ECG, Free-Breathing T1–T2 Mapping

The aim of this study is to simultaneously quantify T1/T2 across three slices of the left-ventricular myocardium without breath-holds or ECG monitoring, all within a 3 min scan. Radial simultaneous multi-slice (SMS) encoding, self-gating, and image reconstruction was incorporated into the cardiovascular magnetic resonance (CMR) Multitasking framework to simultaneously image three short-axis slices. A T2prep-IR FLASH sequence with two flip angles was designed and implemented to allow B1+-robust T1 and T2 mapping. The proposed Multitasking-SMS method was validated in a standardized phantom and 10 healthy volunteers, comparing T1 and T2 measurements and scan-rescan repeatability against corresponding reference methods in one layer of phantom vials and in 16 American Heart Association (AHA) myocardial segments. In phantom, Multitasking-SMS T1/T2 measurements showed substantial correlation (R² > 0.996) and excellent agreement [intraclass correlation coefficients (ICC) ≥ 0.999)] with reference measurements. In healthy volunteers, Multitasking-SMS T1/T2 maps reported similar myocardial T1/T2 values (1,215 ± 91.0/41.5 ± 6.3 ms) to the reference myocardial T1/T2 values (1,239 ± 67.5/42.7 ± 4.1 ms), with P = 0.347 and P = 0.296, respectively. Bland–Altman analyses also demonstrated good in vivo repeatability in both the multitasking and references, with segment-wise coefficients of variation of 4.7% (multitasking T1), 8.9% (multitasking T2), 2.4% [modified look-locker inversion recovery (MOLLI)], and 4.6% (T2-prep FLASH), respectively. In summary, multitasking-SMS is feasible for free-breathing, non-ECG, myocardial T1/T2 quantification in 16 AHA segments over 3 short-axis slices in 3 min. The method shows the great potential for reducing exam time for quantitative CMR without ECG or breath-holds.

Cardiac MR: From Theory to Practice

Cardiovascular disease (CVD) is the leading single cause of morbidity and mortality, causing over 17. 9 million deaths worldwide per year with associated costs of over $800 billion. Improving prevention, diagnosis, and treatment of CVD is therefore a global priority. Cardiovascular magnetic resonance (CMR) has emerged as a clinically important technique for the assessment of cardiovascular anatomy, function, perfusion, and viability. However, diversity and complexity of imaging, reconstruction and analysis methods pose some limitations to the widespread use of CMR. Especially in view of recent developments in the field of machine learning that provide novel solutions to address existing problems, it is necessary to bridge the gap between the clinical and scientific communities. This review covers five essential aspects of CMR to provide a comprehensive overview ranging from CVDs to CMR pulse sequence design, acquisition protocols, motion handling, image reconstruction and quantitative analysis of the obtained data. (1) The basic MR physics of CMR is introduced. Basic pulse sequence building blocks that are commonly used in CMR imaging are presented. Sequences containing these building blocks are formed for parametric mapping and functional imaging techniques. Commonly perceived artifacts and potential countermeasures are discussed for these methods. (2) CMR methods for identifying CVDs are illustrated. Basic anatomy and functional processes are described to understand the cardiac pathologies and how they can be captured by CMR imaging. (3) The planning and conduct of a complete CMR exam which is targeted for the respective pathology is shown. Building blocks are illustrated to create an efficient and patient-centered workflow. Further strategies to cope with challenging patients are discussed. (4) Imaging acceleration and reconstruction techniques are presented that enable acquisition of spatial, temporal, and parametric dynamics of the cardiac cycle. The handling of respiratory and cardiac motion strategies as well as their integration into the reconstruction processes is showcased. (5) Recent advances on deep learning-based reconstructions for this purpose are summarized. Furthermore, an overview of novel deep learning image segmentation and analysis methods is provided with a focus on automatic, fast and reliable extraction of biomarkers and parameters of clinical relevance.

Improved cardiac T1 mapping accuracy and precision with a new hybrid MOLLI and SASHA technique: MOSHA

PURPOSE

To develop and validate a new myocardial T₁ mapping sequence (MOSHA) which is based on a combination of the modified Look-Locker inversion recovery (MOLLI) and the saturation recovery single-shot acquisition (SASHA) sequences.

METHODS

Prior studies have shown that myocardial T₁ mapping by SASHA is more accurate but less precise than MOLLI. A new myocardial T₁ mapping technique (MOSHA) based on single-shot acquisitions is developed by combining the MOLLI and SASHA sequences. Phantom and patient studies on 15 patients (9 males, median age 21 years) were performed to validate and compare MOSHA with the MOLLI and SASHA sequences in terms of accuracy and precision.

RESULTS

In the phantom study, MOSHA was as accurate as SASHA (P-value = 0.88) and as precise as MOLLI (P-value = 0.59). Similar trends were observed in the patient study. Compared to SASHA, MOSHA accuracy was comparable for blood pre-contrast (P-value≥0.10) and post-contrast (P-value≥0.70), and for myocardium pre-contrast (P-value = 0.70) and post-contrast (P-value = 0.09). Compared to MOLLI, MOSHA precision was lower for blood pre-contrast (P-value<0.01) and higher for blood post-contrast (P-value≤0.01), and comparable for myocardium pre-contrast (P-value = 0.24) and post-contrast (P-value = 0.07). Synthetic Extracellular volume fraction (ECV) calculated by MOSHA was more precise than those of SASHA and MOLLI (P-value ≤0.01).

CONCLUSION

In phantom studies and patients, MOSHA has comparable accuracy as SASHA and nearly similar precision as MOLLI for T₁ mapping. Precision of MOSHA was better than MOLLI and SASHA in synthetic ECV measurements. Therefore, it may be a superior choice in clinical practice for a precise and accurate calculation of T₁ and ECV.

Reference Values of Native T1 at 3T Cardiac Magnetic Resonance-Standardization Considerations between Different Vendors

This study aimed at establishing native T1 reference values for a Canon Vantage Galan 3T system and comparing them with previously published values from different vendors. A total of 20 healthy volunteers (55% Women; 33.9 ± 11.1 years) underwent left ventricular T1 mapping at 3T MR. A MOLLI 5(3)3 sequence was used, acquiring three short-axis slices. Native T1 values are shown as means (±standard deviation) and Student’s independent samples t-test was used to test gender differences in T1 values. Pearson’s correlation coefficient analysis was used to compare two processes of T1 analysis. The results show a global native T1 mean value of 1124.9 ± 55.2 ms (exponential analysis), that of women being statistically higher than men (1163 ± 30.5 vs. 1077.9 ± 39.5 ms, respectively; p < 0.001). There were no specific tendencies for T1 times in different ventricular slices. We found a strong correlation (0.977, p < 0.001) with T1 times derived from parametric maps (1136.4 ± 60.2 ms). Native T1 reference values for a Canon 3T scanner were provided, and they are on par with those already reported from other vendors for a similar sequence. We also found a correlation between native T1 and gender, with higher values for women.

Dual flip-angle IR-FLASH with spin history mapping for B1+ corrected T1 mapping: Application to T1 cardiovascular magnetic resonance multitasking

PURPOSE

To develop a single-scan method for B 1 + -corrected T1 mapping and apply it for free-breathing (FB) cardiac MR multitasking without electrocardiogram (ECG) triggering.

METHODS

One dual flip-angle (2FA) inversion recovery (IR)-FLASH scan provides two observations of T 1 ∗ (apparent T1 ) corresponding to two distinct combinations of the nominal FA α and B 1 + . Spatiotemporally coregistered T1 and B 1 + spin history maps are obtained by fitting the 2FA signal model. T1 estimate accuracy and repeatability for single flip-angle (1FA) and 2FA IR-FLASH sequence MR multitasking were evaluated at 3T. A T1 phantom was first imaged on the scanner table, then on two human subjects' thoraxes in both breath-hold (BH) and FB conditions. IR-turbo spin echo (IR-TSE) static phantom T1 measurements served as reference. In 10 healthy subjects, myocardial T1 was evaluated with ECG-free, FB multitasking sequences alongside ECG-triggered BH MOLLI.

RESULTS

For phantom-on-table T1 estimates, 2FA agreed better with IR-TSE (intraclass correlation coefficient [ICC] = 0.996, mean error ± SD = -1.6% ± 1.9%) than did 1FA (ICC = 0.922; mean error ± SD = -4.3% ± 12%). For phantom-on-thorax, 2FA was more repeatable and robust to respiration than 1FA (coefficient of variation [CoV] = 1.2% 2FA, = 11.3% 1FA). In vivo, in intrasession T1 repeatability, 2FA (septal CoV = 2.4%, six-segment CoV = 4.4%) outperformed 1FA (septal CoV = 3.1%, six-segment CoV = 5.5%). In six-segment T1 homogeneity, 2FA (CoV = 7.9%) also outperformed 1FA (CoV = 11.1%).

CONCLUSION

The 2FA IR-FLASH improves T1 estimate accuracy and repeatability over 1FA IR-FLASH, and enables single-scan B 1 + -corrected T1 mapping without BHs or ECG when used with MR multitasking.

Development, validation, qualification, and dissemination of quantitative MR methods: Overview and recommendations by the ISMRM quantitative MR study group

On behalf of the International Society for Magnetic Resonance in Medicine (ISMRM) Quantitative MR Study Group, this article provides an overview of considerations for the development, validation, qualification, and dissemination of quantitative MR (qMR) methods. This process is framed in terms of two central technical performance properties, i.e., bias and precision. Although qMR is confounded by undesired effects, methods with low bias and high precision can be iteratively developed and validated. For illustration, two distinct qMR methods are discussed throughout the manuscript: quantification of liver proton-density fat fraction, and cardiac T₁ . These examples demonstrate the expansion of qMR methods from research centers toward widespread clinical dissemination. The overall goal of this article is to provide trainees, researchers, and clinicians with essential guidelines for the development and validation of qMR methods, as well as an understanding of necessary steps and potential pitfalls for the dissemination of quantitative MR in research and in the clinic.

Measurement of sub-zero temperatures in MRI using T1 temperature sensitive soft silicone materials: Applications for MRI-guided cryosurgery

PURPOSE

One standard method, proton resonance frequency shift, for measuring temperature using magnetic resonance imaging (MRI), in MRI-guided surgeries, fails completely below the freezing point of water. Because of this, we have developed a new methodology for monitoring temperature with MRI below freezing. The purpose of this paper is to show that a strong temperature dependence of the nuclear relaxation time T₁ in soft silicone polymers can lead to temperature-dependent changes of MRI intensity acquired with T₁ weighting. We propose the use of silicone filaments inserted in tissue for measuring temperature during MRI-guided cryoablations.

METHODS

The temperature dependence of T₁ in bio-compatible soft silicone polymers was measured using nuclear magnetic resonance spectroscopy and MRI. Phantoms, made of bulk silicone materials and put in an MRI-compatible thermal container with dry ice, allowed temperature measurements ranging from -60°C to + 20°C. T₁ -weighted gradient echo images of the phantoms were acquired at spatially uniform temperatures and with a gradient in temperature to determine the efficacy of using these materials as temperature indicators in MRI. Ex vivo experiments on silicone rods, 4 mm in diameter, inserted in animal tissue were conducted to assess the practical feasibility of the method.

RESULTS

Measurements of nuclear relaxation times of protons in soft silicone polymers show a monotonic, nearly linear, change with temperature (R²  > 0.98) and have a significant correlation with temperature (Pearson's r > 0.99, p < 0.01). Similarly, the intensity of the MR images in these materials, taken with a gradient echo sequence, are also temperature dependent. There is again a monotonic change in MRI intensity that correlates well with the measured temperature (Pearson's r < -0.98 and p < 0.01). The MRI experiments show that a temperature change of 3°C can be resolved in a distance of about 2.5 mm. Based on MRI images and external sensor calibrations for a sample with a gradient in temperature, temperature maps with 3°C isotherms are created for a bulk phantom. Experiments demonstrate that these changes in MRI intensity with temperature can also be seen in 4 mm silicone rods embedded in ex vivo animal tissue.

CONCLUSIONS

We have developed a new method for measuring temperature in MRI that potentially could be used during MRI-guided cryoablation operations, reducing both procedure time and cost, and making these surgeries safer.

Free-breathing simultaneous T1 , T2 , and T2∗ quantification in the myocardium

PURPOSE

To implement a free-breathing sequence for simultaneous quantification of T 1 , T 2 , and T 2 ∗ for comprehensive tissue characterization of the myocardium in a single scan using a multi-gradient-echo readout with saturation and T 2 preparation pulses.

METHODS

In the proposed Saturation And T 2 -prepared Relaxometry with Navigator-gating (SATURN) technique, a series of multi-gradient-echo (GRE) images with different magnetization preparations was acquired during free breathing. A total of 35 images were acquired in 26.5 ± 14.9 seconds using multiple saturation times and T 2 preparation durations and with imaging at 5 echo times. Bloch simulations and phantom experiments were used to validate a 5-parameter fit model for accurate relaxometry. Free-breathing simultaneous T 1 , T 2 , and T 2 ∗ measurements were performed in 10 healthy volunteers and 2 patients using SATURN at 3T and quantitatively compared to conventional single-parameter methods such as SASHA for T 1 , T 2 -prepared bSSFP, and multi-GRE for T 2 ∗ .

RESULTS

Simulations confirmed accurate fitting with the 5-parameter model. Phantom measurements showed good agreement with the reference methods in the relevant range for in vivo measurements. Compared to single-parameter methods comparable accuracy was achieved. SATURN produced in vivo parameter maps that were visually comparable to single-parameter methods. No significant difference between T 1 , T 2 , and T 2 ∗ times acquired with SATURN and single-parameter methods was shown in quantitative measurements (SATURN T 1 = 1573 ± 86 ms , T 2 = 33.2 ± 3.6 ms , T 2 ∗ = 25.3 ± 6.1 ms ; conventional methods: T 1 = 1544 ± 107 ms , T 2 = 33.2 ± 3.6 ms , T 2 ∗ = 23.8 ± 5.5 ms ; P > . 2 ) CONCLUSION: SATURN enables simultaneous quantification of T 1 , T 2 , and T 2 ∗ in the myocardium for comprehensive tissue characterization with co-registered maps, in a single scan with good agreement to single-parameter methods.

Ferumoxytol-enhanced magnetic resonance T1 reactivity for depiction of myocardial hypoperfusion

Myocardial T1 reactivity, defined as the relative change in T1 between rest and vasodilator-induced stress, has been proposed as a magnetic resonance imaging (MRI) biomarker of tissue perfusion. We hypothesize that the superparamagnetic iron-oxide nanoparticle, ferumoxytol, sensitizes T1 to changes in the intramyocardial vascular compartment and improves the sensitivity and specificity of T1 reactivity as an imaging biomarker of tissue perfusion. We aim to assess the diagnostic performance of ferumoxytol-enhanced (FE) myocardial T1 reactivity in swine models of myocardial hypoperfusion. We induced acute myocardial hypoperfusion in 13 swine via percutaneous, transcatheter deployment of a 3D printed intracoronary stenosis implant into the left anterior descending coronary artery. We performed native and FE adenosine stress testing using 5(3)3(3)3 MOLLI and SASHA T1 mapping sequences with bSSFP readout on a clinical 3.0 T magnet. MOLLI T1 maps were fitted using both the conventional MOLLI and the Instantaneous Signal Loss (InSiL) T1-fitting algorithms. Regardless of the MOLLI or SASHA pulse sequence or T1-fitting algorithm, ferumoxytol contrast increased the dynamic range of T1 reactivity in both the remote and ischemic myocardial regions. Relative to remote myocardium, native and FE T1 reactivity were blunted in ischemic myocardium (p < 0.05) with InSiL-MOLLI, MOLLI and SASHA. An InSiL-MOLLI-derived FE T1 reactivity threshold of -4.65% had 73.3% sensitivity and 96.2% specificity for prediction of regional wall motion abnormalities (AUC 0.915, 95% CI 0.786-0.979), whereas a SASHA-derived FE T1 reactivity threshold of -5.25% had 75.0% sensitivity and 95.2% specificity (AUC 0.905, 95% CI 0.751-0.979). Ferumoxytol significantly increased the dynamic range of T1 reactivity as a measure of myocardial hypoperfusion in vasodilator stress T1 mapping studies. FE T1 reactivity maps can be used to quantitatively distinguish ischemic and remote myocardium with high specificity in swine models of acute myocardial hypoperfusion.

Cardiac stress T1-mapping response and extracellular volume stability of MOLLI-based T1-mapping methods

Stress and rest T1-mapping may assess for myocardial ischemia and extracellular volume (ECV). However, the stress T1 response is method-dependent, and underestimation may lead to misdiagnosis. Further, ECV quantification may be affected by time, as well as the number and dosage of gadolinium (Gd) contrast administered. We compared two commonly available T1-mapping approaches in their stress T1 response and ECV measurement stability. Healthy subjects (n = 10, 50% female, 35 ± 8 years) underwent regadenoson stress CMR (1.5 T) on two separate days. Prototype ShMOLLI 5(1)1(1)1 sequence was used to acquire consecutive mid-ventricular T1-maps at rest, stress and post-Gd contrast to track the T1 time evolution. For comparison, standard MOLLI sequences were used: MOLLI 5(3)3 Low (256 matrix) & High (192 matrix) Heart Rate (HR) to acquire rest and stress T1-maps, and MOLLI 4(1)3(1)2 Low & High HR for post-contrast T1-maps. Stress and rest myocardial blood flow (MBF) maps were acquired after IV Gd contrast (0.05 mmol/kg each). Stress T1 reactivity (delta T1) was defined as the relative percentage increase in native T1 between rest and stress. Myocardial T1 values for delta T1 (dT1) and ECV were calculated. Residuals from the identified time dependencies were used to assess intra-method variability. ShMOLLI achieved a greater stress T1 response compared to MOLLI Low and High HR (peak dT1 = 6.4 ± 1.7% vs. 4.8 ± 1.3% vs. 3.8 ± 1.0%, respectively; both p < 0.0001). ShMOLLI dT1 correlated strongly with stress MBF (r = 0.77, p < 0.001), compared to MOLLI Low HR (r = 0.65, p < 0.01) and MOLLI High HR (r = 0.43, p = 0.07). ShMOLLI ECV was more stable to gadolinium dose with less time drift (0.006-0.04% per minute) than MOLLI variants. Overall, ShMOLLI demonstrated less intra-individual variability than MOLLI variants for stress T1 and ECV quantification. Power calculations indicate up to a fourfold (stress T1) and 7.5-fold (ECV) advantage in sample-size reduction using ShMOLLI. Our results indicate that ShMOLLI correlates strongly with increased MBF during regadenoson stress and achieves a significantly higher stress T1 response, greater effect size, and greater ECV measurement stability compared with the MOLLI variants tested.

Impact of the Choice of Native T1 in Pixelwise Myocardial Blood Flow Quantification

BACKGROUND

Quantification of myocardial blood flow (MBF) from dynamic contrast-enhanced (DCE) MRI can be performed using a signal intensity model that incorporates T1 values of blood and myocardium.

PURPOSE

To assess the impact of T1 values on pixelwise MBF quantification, specifically to evaluate the influence of 1) study population-averaged vs. subject-specific, 2) diastolic vs. systolic, and 3) regional vs. global myocardial T1 values.

STUDY TYPE

Prospective.

SUBJECTS

Fifteen patients with chronic coronary heart disease.

FIELD STRENGTH/SEQUENCE

3T; modified Look-Locker inversion recovery for T1 mapping and saturation recovery gradient echo for DCE imaging, both acquired in a mid-ventricular short-axis slice in systole and diastole.

ASSESSMENT

MBF was estimated using Fermi modeling and signal intensity nonlinearity correction with different T1 values: study population-averaged blood and myocardial, subject-specific systolic and diastolic, and segmental T1 values. Myocardial segments with perfusion deficits were identified visually from DCE series.

STATISTICAL TESTS

The relationships between MBF parameters derived by different methods were analyzed by Bland-Altman analysis; corresponding mean values were compared by t-test.

RESULTS

Using subject-specific diastolic T1 values, global diastolic MBF was 0.61 ± 0.13 mL/(min·g). It did not differ from global MBF derived from the study population-averaged T1 (P = 0.88), but the standard deviation of differences was large (0.07 mL/(min·g), 11% of mean MBF). Global diastolic and systolic MBF did not differ (P = 0.12), whereas global diastolic MBF using systolic (0.62 ± 0.13 mL/(min·g)) and diastolic T1 values differed (P < 0.05). If regional instead of global T1 values were used, segmental MBF was lower in segments with perfusion deficits (bias = -0.03 mL/(min·g), -7% of mean MBF, P < 0.05) but higher in segments without perfusion deficits (bias = 0.01 mL/(min·g), 1% of mean MBF, P < 0.05).

DATA CONCLUSION

Whereas cardiac phase-specific T1 values have a minor impact on MBF estimates, subject-specific and myocardial segment-specific T1 values substantially affect MBF quantification.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY STAGE: 3.

Myocardial T1 mapping using SMART1 Map and MOLLI mapping in asymptomatic patients with recent extracardiac sarcoidosis

INTRODUCTION

Sarcoidosis is a systemic granulomatous disease affecting in particular the respiratory tract. Cardiac magnetic resonance (CMR), including a measurement of T₁ relaxation time, could potentially detect early stadia of sarcoidosis of the heart. The study aims to assess T₁ mapping in the detection of early cardiac involvement in asymptomatic patients with sarcoidosis.

METHODS

One hundred and twenty patients with extracardiac sarcoidosis and without any heart disease history were included. One hundred and thirteen of them underwent a CMR examination. The mean time from the diagnosis of sarcoidosis was 0.8 (0.2-3.3) years. Cine images for the assessment of left ventricular (LV) functional parameters and pre- and post-contrast saturation method using adaptive recovery times for cardiac T₁ mapping (SMART₁ Map) and modified Look-Locker inversion recovery (MOLLI) images were acquired for the assessment of native T₁ relaxation time and extracellular volume (ECV). The measured parameters were compared between sarcoidosis patients and 22 controls.

RESULTS

The sarcoidosis patients had normal global and regional systolic LV function-LV ejection fraction 65 ± 5% versus 66 ± 7% (p NS). The mean native T₁ relaxation times were not prolonged-1465 ± 93 ms versus 1480 ± 88 ms (p NS) measured by SMART₁ Map and 1317 ± 60 ms versus 1313 ± 83 ms (p NS) measured using a MOLLI sequence. Similarly, the mean ECV values did not increase-16.9 ± 3.9% versus 17.9 ± 3.7% (p NS) measured by SMART₁ Map and 30.9 ± 2.9% versus 31.6 ± 8.3% (p NS) measured using a MOLLI sequence.

CONCLUSION

Myocardial native T₁ relaxation times were not prolonged and ECV was not increased in asymptomatic patients with extracardiac sarcoidosis.

Standardization of T1-mapping in cardiovascular magnetic resonance using clustered structuring for benchmarking normal ranges

BACKGROUND

Cardiovascular magnetic resonance T1-mapping is increasingly used for tissue characterization, commonly based on Modified Look-Locker Inversion recovery (MOLLI). However, there are numerous MOLLI variants with differing normal ranges. This lack of standardization presents confusion and difficulty in inter-center comparisons, hindering widespread adoption of T1-mapping.

METHODS

To address this, we performed a structured literature search for native left ventricular myocardial T1-mapping in healthy humans measured using MOLLI variants at 1.5 and 3 Tesla, across scanner vendors. We then used k-means clustering to structure normal MOLLI-T1 values according to magnetic field strength, and investigated correlations between common imaging parameters: repetition time (TR), echo time (TE), flip angle (FA).

RESULTS

We analyzed data from 2207 healthy controls in 76 independent reports. Normal MOLLI-T1 standard deviations varied by 11-fold, and dependencies on TE, TR, and FA differed between 1.5 T and 3 T, thwarting meaningful T1 standardization even within a single field strength, including the use of Z-score. However, divergent MOLLI-T1 norms may be structured using data clustering. For 1.5 T, two clusters emerged: Cluster1_1.5T: T1 = 958 ± 16 ms (n = 1280); Cluster2_1.5T: T1 = 1027 ± 19 ms (n = 386). For 3 T, three clusters emerged: Cluster1_3T: T1 = 1160 ± 21 ms (n = 330); Cluster2_3T: T1 = 1067 ± 18 ms (n = 178); Cluster3_3T: T1 = 1227 ± 19 ms (n = 41). We then propose the concept of an online calculator for assigning local norms to a known MOLLI-T1 cluster, allowing benchmarking against published norms.

CONCLUSION

Clustered structuring allows T1 standardization of widely-divergent MOLLI variants, benchmarking local norms (usually based on smaller samples) against published norms (larger samples). This may increase confidence and quality control in method implementation, facilitating wider clinical adoption of T1-mapping.

Magnetic resonance fingerprinting of the pancreas at 1.5 T and 3.0 T

Magnetic resonance imaging of the pancreas is increasingly used as an important diagnostic modality for characterisation of pancreatic lesions. Pancreatic MRI protocols are mostly qualitative due to time constraints and motion sensitivity. MR Fingerprinting is an innovative acquisition technique that provides qualitative data and quantitative parameter maps from a single free-breathing acquisition with the potential to reduce exam times. This work investigates the feasibility of MRF parameter mapping for pancreatic imaging in the presence of free-breathing exam. Sixteen healthy participants were prospectively imaged using MRF framework. Regions-of-interest were drawn in multiple solid organs including the pancreas and T1 and T2 values determined. MRF T1 and T2 mapping was performed successfully in all participants (acquisition time:2.4-3.6 min). Mean pancreatic T1 values were 37-43% lower than those of the muscle, spleen, and kidney at both 1.5 and 3.0 T. For these organs, the mean pancreatic T2 values were nearly 40% at 1.5 T and < 12% at 3.0 T. The feasibility of MRF at 1.5 T and 3 T was demonstrated in the pancreas. By enabling fast and free-breathing quantitation, MRF has the potential to add value during the clinical characterisation and grading of pathological conditions, such as pancreatitis or cancer.

Cardiac cine magnetic resonance fingerprinting for combined ejection fraction, T1 and T2 quantification

This study introduces a technique called cine magnetic resonance fingerprinting (cine-MRF) for simultaneous T₁ , T₂ and ejection fraction (EF) quantification. Data acquired with a free-running MRF sequence are retrospectively sorted into different cardiac phases using an external electrocardiogram (ECG) signal. A low-rank reconstruction with a finite difference sparsity constraint along the cardiac motion dimension yields images resolved by cardiac phase. To improve SNR and precision in the parameter maps, these images are nonrigidly registered to the same phase and matched to a dictionary to generate T₁ and T₂ maps. Cine images for computing left ventricular volumes and EF are also derived from the same data. Cine-MRF was tested in simulations using a numerical relaxation phantom. Phantom and in vivo scans of 19 subjects were performed at 3 T during a 10.9 seconds breath-hold with an in-plane resolution of 1.6 x 1.6 mm² and 24 cardiac phases. Left ventricular EF values obtained with cine-MRF agreed with the conventional cine images (mean bias -1.0%). Average myocardial T₁ times in diastole/systole were 1398/1391 ms with cine-MRF, 1394/1378 ms with ECG-triggered cardiac MRF (cMRF) and 1234/1212 ms with MOLLI; and T₂ values were 30.7/30.3 ms with cine-MRF, 32.6/32.9 ms with ECG-triggered cMRF and 37.6/41.0 ms with T₂ -prepared FLASH. Cine-MRF and ECG-triggered cMRF relaxation times were in good agreement. Cine-MRF T₁ values were significantly longer than MOLLI, and cine-MRF T₂ values were significantly shorter than T₂ -prepared FLASH. In summary, cine-MRF can potentially streamline cardiac MRI exams by combining left ventricle functional assessment and T₁ -T₂ mapping into one time-efficient acquisition.

Free-breathing whole-heart multi-slice myocardial T1 mapping in 2 minutes

PURPOSE

To develop and validate a saturation-delay-inversion recovery preparation, slice tracking and multi-slice based sequence for measuring whole-heart native T₁ .

METHOD

The proposed free-breathing sequence performs T₁ mapping of multiple left-ventricular slices by slice-interleaved acquisition to collect 10 electrocardiogram-triggered single-shot slice-selective images for each slice. A saturation-delay-inversion recovery pulse is used for T₁ preparation. Prospective slice tracking by the diaphragm navigator and retrospective registration are used to reduce through-plane and in-plane motion, respectively. The proposed sequence was validated in both phantom and human subjects (12 healthy subjects and 15 patients who were referred for a clinical cardiac MR exam) and compared with saturation recovery single-shot acquisition (SASHA) and modified Look-Locker inversion recovery (MOLLI).

RESULTS

Phantom T₁ measured by the proposed sequence had excellent agreement (R²  = 0.99) with the ground-truth T₁ and was insensitive to heart rate. In both healthy subjects and patients, the proposed sequence yielded nine left-ventricular T₁ maps per volume in less than 2 minutes (healthy volunteers: 1.8 ± 0.4 minutes; patients: 1.9 ± 0.2 minutes). The average T₁ of whole left ventricle for all healthy subjects and patients were 1560 ± 61 and 1535 ± 49 ms by SASHA, 1208 ± 42 and 1233 ± 56 ms by MOLLI5(3)3, and 1397 ± 34 and 1433 ± 56 ms by the proposed sequence, respectively. The corresponding coefficient of variation of T₁ were 6.2 ± 1.4% and 5.8 ± 1.6%, 5.3 ± 1.1% and 5.1 ± 0.8%, and 4.9 ± 0.8% and 4.5 ± 0.8%, respectively.

CONCLUSION

The proposed sequence enables quantification of whole heart T₁ with good accuracy and precision in less than 2 minutes during free breathing.

Towards measuring the effect of flow in blood T 1 assessed in a flow phantom and in vivo

Measurement of the blood T ₁ time using conventional myocardial T ₁ mapping methods has gained clinical significance in the context of extracellular volume (ECV) mapping and synthetic hematocrit (Hct). However, its accuracy is potentially compromised by in-flow of non-inverted/non-saturated spins and in-flow of spins which are not partially saturated from previous imaging pulses. Bloch simulations were used to analyze various flow effects separately. T ₁ measurements of gadolinium doped water were performed using a flow phantom with adjustable flow velocities at 3 T. Additionally, in vivo blood T ₁ measurements were performed in 6 healthy subjects (26 ± 5 years, 2 female). To study the T ₁ time as a function of the instantaneous flow velocity, T ₁ times were evaluated in an axial imaging slice of the descending aorta. Velocity encoded cine measurements were performed to quantify the flow velocity throughout the cardiac cycle. Simulation results show more than 30% loss in accuracy for 10% non-prepared in-flowing spins. However, in- and out-flow to the imaging plane only demonstrated minor impact on the T ₁ time. Phantom T ₁ times were decreased by up to 200 ms in the flow phantom, due to in-flow of non-prepared spins. High flow velocities cause in-flow of spins that lack partial saturation from the imaging pulses but only lead to negligible T ₁ time deviation (less than 30 ms). In vivo measurements confirm a substantial variation of the T ₁ time depending on the flow velocity. The highest aortic T ₁ times are observed at the time point of minimal flow with increased flow velocity leading to reduction of the measured T ₁ time by up to [Formula: see text] at peak velocity. In this work we attempt to dissect the effects of flow on T ₁ times, by using simulations, well-controlled, simplified phantom setup and the linear flow pattern in the descending aorta in vivo.

Single-shot EPI for ASL-CMR

PURPOSE

To evaluate single-shot echo planar imaging (SS-EPI), as an alternative to snapshot balanced steady state free precession (bSSFP) imaging, for arterial-spin-labeled cardiac MR (ASL-CMR). This study presents a practical implementation SS-EPI tailored to the needs of ASL-CMR at 3T and demonstrates sequential multi-slice ASL with no increase in scan time.

METHODS

Reduced field of view SS-EPI was performed using a 2DRF pulse. A spin-echo was used with crushers optimized to maximize blood suppression and minimize myocardial signal loss, based on experiments in 4 healthy volunteers. SS-EPI was evaluated against the widely used bSSFP reference method in single-slice ASL-CMR in 4 healthy volunteers, during both systole and diastole. Sequential multi-slice ASL-CMR with SS-EPI was demonstrated during diastole (3 slices: basal, mid, and apical short-axis) and during systole (2 slices: mid and apical short-axis), in 3 volunteers.

RESULTS

Global myocardial perfusion for diastolic SS-EPI (1.66 ± 0.73 mL/g/min) and systolic SS-EPI (1.50 ± 0.36 mL/g/min) were found to be statistically equivalent (2 one-sided test with a difference of 0.4 mL/g/min) to diastolic bSSFP (duration of 1 cardiac cycle, 1.60 ± 0.80 mL/g/min) with P-values of 0.022 and 0.031, respectively. Global myocardial perfusion for sequential multi-slice experiments was 1.64 ± 0.47, 1.34 ± 0.29, and 1.88 ± 0.58 for basal, mid, and apical SAX slices during diastole and was 1.61 ± 0.35, and 1.66 ± 0.49 for mid and apical slice during systole. These values are comparable to published ASL-CMR and positron emission tomography studies.

CONCLUSION

SS-EPI is a promising alternative to bSSFP imaging for ASL-CMR and can potentially improve the spatial coverage of ASL-CMR by 3-fold during diastole and 2-fold during systole, without increasing scan time.

Feasibility of coronary endothelial function assessment using arterial spin labeled CMR

Coronary endothelial dysfunction (CED) is an independent predictor of cardiovascular disease, but its assessment has been limited to invasive coronary angiography. Myocardial perfusion imaging using arterial spin labeled (ASL) cardiac magnetic resonance (CMR) may be an effective non-invasive alternative for detection of CED. Thirty-four patients were recruited: 10 healthy volunteers, 13 at high-risk for coronary artery disease (CAD), and 11 with established CAD. ASL-CMR was performed continuously in a single mid-short axis slice during rest, stress, and recovery. Stress was induced with sustained isometric handgrip exercise, an endothelial dependent stressor. Myocardial perfusion (MP) during rest, peak stress, and recovery were calculated and compared. After excluding subjects unable to complete the protocol or who exhibited poor data quality, 6 healthy, 10 high-risk, and 7 CAD patients were included in the analysis. Average MP (ml/g/min) was 1.31 ± 1.23, 1.61 ± 1.12, and 1.40 ± 0.97 at rest, and 1.64 ± 1.49, 2.31 ± 1.61, and 2.84 ± 1.77 during stress, for the CAD, high-risk and healthy group, respectively. The average MP response (MP_stress - MP_rest , ml/g/min) was 0.32 ± 1.93, 0.69 ± 1.34, and 1.44 ± 1.46 for CAD, high-risk and healthy group, respectively. MP during handgrip stress was significantly lower for both the CAD (p = 0.0005) and high-risk groups (p = 0.05) compared to the healthy volunteers. In only the healthy subjects, MP was significantly higher in stress compared to rest (p = 0.0002). Participants with CAD had significantly lower MP response compared to healthy volunteers, as detected by ASL-CMR. These findings support the feasibility of ASL-CMR for non-invasive assessment of CED.

Low-Rank Tensor Models for Improved Multi-Dimensional MRI: Application to Dynamic Cardiac T 1 Mapping

Multi-dimensional, multi-contrast magnetic resonance imaging (MRI) has become increasingly available for comprehensive and time-efficient evaluation of various pathologies, providing large amounts of data and offering new opportunities for improved image reconstructions. Recently, a cardiac phase-resolved myocardial T 1 mapping method has been introduced to provide dynamic information on tissue viability. Improved spatio-temporal resolution in clinically acceptable scan times is highly desirable but requires high acceleration factors. Tensors are well-suited to describe inter-dimensional hidden structures in such multi-dimensional datasets. In this study, we sought to utilize and compare different tensor decomposition methods, without the use of auxiliary navigator data. We explored multiple processing approaches in order to enable high-resolution cardiac phase-resolved myocardial T 1 mapping. Eight different low-rank tensor approximation and processing approaches were evaluated using quantitative analysis of accuracy and precision in T 1 maps acquired in six healthy volunteers. All methods provided comparable T 1 values. However, the precision was significantly improved using local processing, as well as a direct tensor rank approximation. Low-rank tensor approximation approaches are well-suited to enable dynamic T 1 mapping at high spatio-temporal resolutions.

Differentiation between Fabry disease and hypertrophic cardiomyopathy with cardiac T1 mapping

PURPOSE

To evaluate the potential of non-contrast myocardial T1 mapping on cardiovascular magnetic resonance examination (CMR) in differentiating patients with Fabry disease (FD) from those with hypertrophic cardiomyopathy (HCM) and healthy control subjects.

MATERIALS AND METHODS

Seventeen patients with FD (8 men, 9 women; mean age, 48 ±18 [SD] years; [range: 19-73 years]; 53% with left ventricular hypertrophy [LVH]) were matched with 36 patients with hypertrophic cardiomyopathy (HCM) (22 men, 14 women; mean age, 57±16 [SD] years; [range: 22-85 years]) and 70 healthy control subjects (34 men, 36 women; mean age, 38 ±15 [SD] years; [range: 18-65 years]). Cardiac T1 mapping was performed using the modified Look-Locker inversion (MOLLI®) sequence on a 1.5-T magnet. T1 values were calculated, on midventricular section, for septal left ventricular segments (S8-S9) and all mid-ventricular ones (global T1 values; S7-S12). Statistical analysis included unpaired Mann-Whitney test, receiver operating characteristic curve and likelihood ratios.

RESULTS

Septal native T1 values were significantly decreased in patients with FD (889±61 [SD] ms; range: 784-980ms) compared to those with HCM (995±48 [SD] ms; range: 935-1125ms) (P<0.001) and versus healthy controls (965±29 [SD] ms; range: 910-1028ms) (P<0.001). Global native T1 values were also significantly decreased in patients with FD (891±49 [SD] ms; range 794-970ms) compared to those with HCM (995±34 [SD] ms; range: 952-1086ms) (P<0.001) and versus healthy controls (966±27 [SD] ms; range: 920-1042ms) (P<0.001). A septal left ventricular native T1 cutoff value of 940ms could distinguish FD from HCM with 88% sensitivity (95% CI: 73-100%) and 92% specificity (95% CI: 83-100%). Positive likelihood ratio was 11, negative likelihood ratio was 0.12. Compared to controls, the same threshold could distinguish FD with 88% sensitivity (95% CI: 73-100%) and 86% specificity (95% CI: 78-94%). Positive likelihood ratio was 6.3, negative likelihood ratio was 0.14. T1 value was abnormal in 4 of 8 (50%) of FD patients who did not have LVH.

CONCLUSION

Native T1 values are significantly lower in patients with FD by comparison with those with HCM and healthy volunteers.

Early detection of heart function abnormality by native T1: a comparison of two T1 quantification methods

Objective

To compare the robustness of native T1 mapping using mean and median pixel-wise quantification methods.

Methods

Fifty-seven consecutive patients without overt signs of heart failure were examined in clinical routine for suspicion of cardiomyopathy. MRI included the acquisition of native T1 maps by a motion-corrected modified Look-Locker inversion recovery sequence at 1.5 T. Heart function status according to four established volumetric left ventricular (LV) cardio MRI parameter thresholds was used for retrospective separation into subgroups of normal (n = 26) or abnormal heart function (n = 31). Statistical normality of pixel-wise T1 was tested on each myocardial segment and mean and median segmental T1 values were assessed.

Results

Segments with normally distributed pixel-wise T1 (57/58%) showed no difference between mean and median quantification in either patient group, while differences were highly significant (p < 0.001) for the respective 43/42% non-normally distributed segments. Heart function differentiation between two patient groups was significant in 14 myocardial segments (p < 0.001–0.040) by median quantification compared with six (p < 0.001–0.042) by using the mean. The differences by median quantification were observed between the native T1 values of the three coronary artery territories of normal heart function patients (p = 0.023) and insignificantly in the abnormal patients (p = 0.053).

Conclusion

Median quantification increases the robustness of myocardial native T1 definition, regardless of statistical normality of the data. Compared with the currently prevailing method of mean quantification, differentiation between LV segments and coronary artery territories is better and allows for earlier detection of heart function impairment.

Key Points

• Median pixel-wise quantification of native T1 maps is robust and can be applied regardless of the statistical distribution of data points.

• Median quantification is more sensitive to early heart function abnormality compared with mean quantification.

• The new method yields significant native T1 value differentiation between the three coronary artery territories.

Electronic supplementary material

The online version of this article (10.1007/s00330-019-06364-9) contains supplementary material, which is available to authorized users.

T1- and ECV-mapping in clinical routine at 3 T: differences between MOLLI, ShMOLLI and SASHA

BACKGROUND

T1 mapping sequences such as MOLLI, ShMOLLI and SASHA make use of different technical approaches, bearing strengths and weaknesses. It is well known that obtained T1 relaxation times differ between the sequence techniques as well as between different hardware. Yet, T1 quantification is a promising tool for myocardial tissue characterization, disregarding the absence of established reference values. The purpose of this study was to evaluate the feasibility of native and post-contrast T1 mapping methods as well as ECV maps and its diagnostic benefits in a clinical environment when scanning patients with various cardiac diseases at 3 T.

METHODS

Native and post-contrast T1 mapping data acquired on a 3 T full-body scanner using the three pulse sequences 5(3)3 MOLLI, ShMOLLI and SASHA in 19 patients with clinical indication for contrast enhanced MRI were compared. We analyzed global and segmental T1 relaxation times as well as respective extracellular volumes and compared the emerged differences between the used pulse sequences.

RESULTS

T1 times acquired with MOLLI and ShMOLLI exhibited systematic T1 deviation compared to SASHA. Myocardial MOLLI T1 times were 19% lower and ShMOLLI T1 times 25% lower compared to SASHA. Native blood T1 times from MOLLI were 13% lower than SASHA, while post-contrast MOLLI T1-times were only 5% lower. ECV values exhibited comparably biased estimation with MOLLI and ShMOLLI compared to SASHA in good agreement with results reported in literature. Pathology-suspect segments were clearly differentiated from remote myocardium with all three sequences.

CONCLUSION

Myocardial T1 mapping yields systematically biased pre- and post-contrast T1 times depending on the applied pulse sequence. Additionally calculating ECV attenuates this bias, making MOLLI, ShMOLLI and SASHA better comparable. Therefore, myocardial T1 mapping is a powerful clinical tool for classification of soft tissue abnormalities in spite of the absence of established reference values.

FASt single-breathhold 2D multislice myocardial T1 mapping (FAST1) at 1.5T for full left ventricular coverage in three breathholds

Background

Conventional myocardial T₁ mapping techniques such as modified Look–Locker inversion recovery (MOLLI) generate one T₁ map per breathhold. T₁ mapping with full left ventricular coverage may be desirable when spatial T₁ variations are expected. This would require multiple breathholds, increasing patient discomfort and prolonging scan time.

Purpose

To develop and characterize a novel FASt single‐breathhold 2D multislice myocardial T₁ mapping (FAST1) technique for full left ventricular coverage.

Study Type

Prospective.

Population/Phantom

Numerical simulation, agarose/NiCl₂ phantom, 9 healthy volunteers, and 17 patients.

Field Strength/Sequence

1.5T/FAST1.

Assessment

Two FAST1 approaches, FAST1‐BS and FAST1‐IR, were characterized and compared with standard 5‐(3)‐3 MOLLI in terms of accuracy, precision/spatial variability, and repeatability.

Statistical Tests

Kruskal‐Wallis, Wilcoxon signed rank tests, intraclass correlation coefficient analysis, analysis of variance, Student's t‐tests, Pearson correlation analysis, and Bland–Altman analysis.

Results

In simulation/phantom, FAST1‐BS, FAST1‐IR, and MOLLI had an accuracy (expressed as T₁ error) of 0.2%/4%, 6%/9%, and 4%/7%, respectively, while FAST1‐BS and FAST1‐IR had a precision penalty of 1.7/1.5 and 1.5/1.4 in comparison with MOLLI, respectively. In healthy volunteers, FAST1‐BS/FAST1‐IR/MOLLI led to different native myocardial T₁ times (1016 ± 27 msec/952 ±22 msec/987 ± 23 msec, P < 0.0001) and spatial variability (66 ± 10 msec/57 ± 8 msec/46 ± 7 msec, P < 0.001). There were no statistically significant differences between all techniques for T₁ repeatability (P = 0.18). In vivo native and postcontrast myocardial T₁ times in both healthy volunteers and patients using FAST1‐BS/FAST1‐IR were highly correlated with MOLLI (Pearson correlation coefficient ≥0.93).

Data Conclusion

FAST1 enables myocardial T₁ mapping with full left ventricular coverage in three separated breathholds. In comparison with MOLLI, FAST1 yield a 5‐fold increase of spatial coverage, limited penalty of T₁ precision/spatial variability, no significant difference of T₁ repeatability, and highly correlated T₁ times. FAST1‐IR provides improved T₁ precision/spatial variability but reduced accuracy when compared with FAST1‐BS.

Level of Evidence: 1

Technical Efficacy: Stage 3

J. Magn. Reson. Imaging 2020;51:492–504.

Functional LGE Imaging: Cardiac Phase-Resolved Assessment of Focal Fibrosis

Cardiac Magnetic Resonance Imaging (CMR) is a central tool for diagnosis of various ischemic and non-ischemic cardiomyopathies. CMR protocols commonly comprise assessment of functional properties using cardiac phase-resolved CINE MRI and characterization of myocardial viability using late gadolinium enhancement (LGE) imaging. Conventional LGE imaging requires inversion recovery preparation with a specific inversion time to null the healthy myocardium, which restricts the acquisition to a single cardiac phase. In turn, this necessitates separate scans for cardiac function and viability. In this work, we develop a new method for functional LGE imaging in a single breath-hold using a three-step approach: 1) ECG-triggered multi-contrast data is acquired for each cardiac phase, 2) semi-quantitative relaxation maps are generated, 3) LGE imaging contrast is synthesized based on the semi-quantitative maps. The proposed functional LGE method is evaluated in four healthy subject and 20 patients at 1.5T and 3T. Thorough suppression of the healthy myocardium, as well as 40-80ms temporal resolution are achieved, with no visually apparent temporal blurring at tissue interfaces. Functional LGE in patients with focal scar demonstrates robust hyperenhancement in the scar area throughout all cardiac phases, allowing for visual assessment of scar motility. The proposed technique bears the potential to simplify and speedup common cardiac imaging protocols, while enabling improved data fusion of functional and viability information for improved evaluation of CMR.

Liver MR relaxometry at 3T - segmental normal T1 and T2* values in patients without focal or diffuse liver disease and in patients with increased liver fat and elevated liver stiffness

Magnetic resonance (MR) T₁ and T₂* mapping allows quantification of liver relaxation times for non-invasive characterization of diffuse liver disease. We hypothesized that liver relaxation times are not only influenced by liver fibrosis, inflammation and fat, but also by air in liver segments adjacent to the lung – especially in MR imaging at 3T. A total of 161 study participants were recruited, while 6 patients had to be excluded due to claustrophobia or technically uninterpretable MR elastography. Resulting study population consisted of 12 healthy volunteers and 143 patients who prospectively underwent multiparametric MR imaging at 3T. Of those 143 patients, 79 had normal liver stiffness in MR elastography (shear modulus <2.8 kPa, indicating absence of fibrosis) and normal proton density fat fraction (PDFF < 10%, indicating absence of steatosis), defined as reference population. T₁ relaxation times in these patients were significantly shorter in liver segments adjacent to the lung than in those not adjacent to the lung (p < 0.001, mean of differences 33 ms). In liver segments not adjacent to the lung, T₁ allowed to differentiate significantly between the reference population and patients with steatosis and/or fibrosis (p ≤ 0.011), while there was no significant difference of T₁ between the reference population and healthy volunteers. In conclusion, we propose to measure T₁ relaxation times in liver segments not adjacent to the lung. Otherwise, we recommend taking into account slightly shorter T₁ values in liver segments adjacent to the lung.

Kidney transplantation is associated with reduced myocardial fibrosis. A cardiovascular magnetic resonance study with native T1 mapping

BACKGROUND

The measurement of native T1 through cardiovascular magnetic resonance (CMR) is a noninvasive method of assessing myocardial fibrosis without gadolinium contrast. No studies so far have evaluated native T1 after renal transplantation. The primary aim of the current study is to assess changes in the myocardium native T1 6 months after renal transplantation.

METHODS

We prospectively evaluated 44 renal transplant patients with 3 T CMR exams: baseline at the beginning of transplantation and at 6 months after transplantation.

RESULTS

The native T1 time was measured in the midventricular septum and decreased significantly from 1331 ± 52 ms at the baseline to 1298 ± 42 ms 6 months after transplantation (p = 0.001). The patients were split into two groups through a two-step cluster algorithm: In cluster-1 (n = 30) the left ventricular (LV) mass index and the prevalence of diabetes were lower. In cluster-2 (n = 14) the LV mass index and diabetes prevalence were higher. Decrease in native T1 values was significant only in the patients in cluster-1 (p = 0.001).

CONCLUSIONS

The native myocardial T1 time decreased significantly 6 months after renal transplant, which may be associated with the regression of the reactive fibrosis. The patients with greater baseline LV mass index and the diabetic group did not reach a significant decrease in T1.

Assessment of sympathetic reinnervation after cardiac transplantation using hybrid cardiac PET/MRI: A pilot study

Background

Sympathetic reinnervation after heart transplantation (HTX) is a known phenomenon, which has an impact on patient heart rate variability and exercise capacity. The impact of reinnervation on myocardial structure has not been evaluated yet.

Propose

To evaluate the feasibility of simultaneous imaging of cardiac reinnervation and cardiac structure using a hybrid PET/MRI system.

Study type

Prospective / pilot study.

Subjects

Ten patients, 4–21 years after cardiac transplantation.

Field Strength/Sequence

3 T hybrid PET/MRI system. Cine SSFP, T₁ mapping (modified Look–Locker inversion recovery sequence) pre/postcontrast as well as dynamic [¹¹C]meta‐hydroxyephedrine ([¹¹C]mHED) PET.

Assessment

All MRI and PET parameters were evaluated by experienced readers using dedicated postprocessing software packages for cardiac MRI and PET. For all parameters a 16‐segment model for the left ventricle was applied.

Statistical Tests

Mann–Whitney U‐test; Spearman correlations.

Results

Thirty‐six of 160 myocardial segments showed evidence of reinnervation by PET. On a segment‐based analysis, mean native T₁ relaxation times were nonsignificantly altered in segments with evidence of reinnervation (1305 ± 151 msec vs. 1270 ± 112 msec; P = 0.1), whereas mean extracellular volume (ECV) was significantly higher in segments with evidence of reinnervation (35.8 ± 11% vs. 30.9 ± 7%; P = 0.019). There were no significant differences in wall motion (WM) and wall thickening (WT) between segments with or without reinnervation (mean WM: 7.6 ± 4 mm vs. group B: 9.3 ± 7 mm [P = 0.13]; WT: 79 ± 63% vs. 94 ± 74% [P = 0.27]) under resting conditions.

Data Conclusion

The assessment of cardiac reinnervation using a hybrid PET/MRI system is feasible. Segments with evidence of reinnervation by PET showed nonsignificantly higher T₁ relaxation times and a significantly higher ECV, suggesting a higher percentage of diffuse fibrosis in these segments, without impairment of rest WM and WT.

Level of Evidence: 3

Technical Efficacy: Stage 3

J. Magn. Reson. Imaging 2019;50:1326–1335.