Robust motion correction for cardiac T1 and ECV mapping using a T1 relaxation model approach

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.

Exercise-induced myocardial T1 increase and right ventricular dysfunction in recreational cyclists: a CMR study

PURPOSE

Although cardiac troponin I (cTnI) increase following strenuous exercise has been observed, the development of exercise-induced myocardial edema remains unclear. Cardiac magnetic resonance (CMR) native T1/T2 mapping is sensitive to the pathological increase of myocardial water content. Therefore, we evaluated exercise-induced acute myocardial changes in recreational cyclists by incorporating biomarkers, echocardiography and CMR.

METHODS

Nineteen male recreational participants (age: 48 ± 5 years) cycled the 'L'étape du tour de France" (EDT) 2021' (175 km, 3600 altimeters). One week before the race, a maximal graded cycling test was conducted to determine individual heart rate (HR) training zones. One day before and 3-6 h post-exercise 3 T CMR and echocardiography were performed to assess myocardial native T1/T2 relaxation times and cardiac function, and blood samples were collected. All participants were asked to cycle 2 h around their anaerobic gas exchange threshold (HR zone 4).

RESULTS

Eighteen participants completed the EDT stage in 537 ± 58 min, including 154 ± 61 min of cycling time in HR zone 4. Post-race right ventricular (RV) dysfunction with reduced strain and increased volumes (p < 0.05) and borderline significant left ventricular global longitudinal strain reduction (p = 0.05) were observed. Post-exercise cTnI (0.75 ± 5.1 ng/l to 69.9 ± 41.6 ng/l; p < 0.001) and T1 relaxation times (1133 ± 48 ms to 1182 ± 46 ms, p < 0.001) increased significantly with no significant change in T2 (p = 0.474). cTnI release correlated with increase in T1 relaxation time (p = 0.002; r = 0.703), post-race RV dysfunction (p < 0.05; r = 0.562) and longer cycling in HR zone 4 (p < 0.05; r = 0.607).

CONCLUSION

Strenuous exercise causes early post-race cTnI increase, increased T1 relaxation time and RV dysfunction in recreational cyclists, which showed interdependent correlation. The long-term clinical significance of these changes needs further investigation.

TRIAL REGISTRATION NUMBERS AND DATE

NCT04940650 06/18/2021. NCT05138003 06/18/2021.

Motion correction for native myocardial T1 mapping using self-supervised deep learning registration with contrast separation

In myocardial T₁ mapping, undesirable motion poses significant challenges because uncorrected motion can affect T₁ estimation accuracy and cause incorrect diagnosis. In this study, we propose and evaluate a motion correction method for myocardial T₁ mapping using self-supervised deep learning based registration with contrast separation (SDRAP). A sparse coding based method was first proposed to separate the contrast component from T₁ -weighted (T1w) images. Then, a self-supervised deep neural network with cross-correlation (SDRAP-CC) or mutual information as the registration similarity measurement was developed to register contrast separated images, after which signal fitting was performed on the motion corrected T1w images to generate motion corrected T₁ maps. The registration network was trained and tested in 80 healthy volunteers with images acquired using the modified Look-Locker inversion recovery (MOLLI) sequence. The proposed SDRAP was compared with the free form deformation (FFD) registration method regarding (1) Dice similarity coefficient (DSC) and mean boundary error (MBE) of myocardium contours, (2) T₁ value and standard deviation (SD) of T₁ fitting, (3) subjective evaluation score for overall image quality and motion artifact level, and (4) computation time. Results showed that SDRAP-CC achieved the highest DSC of 85.0 ± 3.9% and the lowest MBE of 0.92 ± 0.25 mm among the methods compared. Additionally, SDRAP-CC performed the best by resulting in lower SD value (28.1 ± 17.6 ms) and higher subjective image quality scores (3.30 ± 0.79 for overall quality and 3.53 ± 0.68 for motion artifact) evaluated by a cardiologist. The proposed SDRAP took only 0.52 s to register one slice of MOLLI images, achieving about sevenfold acceleration over FFD (3.7 s/slice).

Case Report: Treatment of Hypertrophic Cardiomyopathy With Stereotactic Body Radiotherapy

Patients with hypertrophic cardiomyopathy (HCM), which is characterized by left ventricular hypertrophy, is usually treated with medications such as calcium channel blockers or beta-blockers and invasive treatments such as transcatheter alcohol septal ablation, percutaneous radiofrequency ablation, or heart transplantation. However, non-invasive methods have not been employed for the management of patients with HCM. A 71-year-old male who presented with occasional chest pain for approximately 2 months and had been diagnosed with HCM since he was 39 years old due to occasional fainting was treated with a novel method for HCM using stereotactic body radiotherapy (SBRT). The administration of 25 Gy of radiation as one fraction led to an improvement in his quality of life. No toxicity occurred during or immediately after the treatment. Our observations suggest that SBRT may be a reasonable treatment approach for patients with HCM who are not suitable for surgery.

Myocardial extracellular volume fraction quantification based on T1 mapping at 3 T: quality optimization by contour-based registration and segmental analysis

BACKGROUND

Myocardial extracellular volume fraction (ECV) assessment can be affected by various technical and subject-related factors.

PURPOSE

To evaluate the role of contour-based registration in quantification of ECV and investigate normal segment-based myocardial ECV values at 3T.

MATERIAL AND METHODS

Pre- and post-contrast T1 mapping images of the left ventricular basal, mid-cavity, and apical slices were obtained in 26 healthy volunteers. ECV maps were generated using motion correction with and without contour-based registration. The image quality of all ECV maps was evaluated by a 4-point scale. Slices were dichotomized according to the occurrence of misregistration in the source data. Contour-registered ECVs and standard ECVs were compared within each subgroup using analysis of variance for repeated measurements and generalized linear mixed models.

RESULTS

In all three slices, higher quality of ECV maps were found using contour-registered method than using standard method. Standard ECVs were statistically different from contour-registered ECVs in global (26.8% ± 2.8% vs. 25.8% ± 2.4%; P = 0.001), mid-cavity (25.4% ± 3.1% vs. 24.3% ± 2.5%; P = 0.016), and apical slices (28.7% ± 4.1% vs. 27.2% ± 3.4%; P = 0.010). In the misregistration subgroups, contour-registered ECVs were lower with smaller SDs (basal: 25.2% ± 1.8% vs. 26.7% ± 2.6%; P = 0.038; mid-cavity: 24.4% ± 2.3% vs. 26.8% ± 3.1%; P = 0.012; apical: 27.5% ± 3.6% vs. 29.7% ± 4.5%; P = 0.016). Apical (27.2% ± 3.4%) and basal-septal ECVs (25.6% ± 2.6%) were statistically higher than mid-cavity ECV (24.3% ± 2.5%; both P < 0.001).

CONCLUSION

Contour-based registration can optimize image quality and improve the precision of ECV quantification in cases demonstrating ventricular misregistration among source images.

Motion correction of free-breathing magnetic resonance renography using model-driven registration

INTRODUCTION

Model-driven registration (MDR) is a general approach to remove patient motion in quantitative imaging. In this study, we investigate whether MDR can effectively correct the motion in free-breathing MR renography (MRR).

MATERIALS AND METHODS

MDR was generalised to linear tracer-kinetic models and implemented using 2D or 3D free-form deformations (FFD) with multi-resolution and gradient descent optimization. MDR was evaluated using a kidney-mimicking digital reference object (DRO) and free-breathing patient data acquired at high temporal resolution in multi-slice 2D (5 patients) and 3D acquisitions (8 patients). Registration accuracy was assessed using comparison to ground truth DRO, calculating the Hausdorff distance (HD) between ground truth masks with segmentations and visual evaluation of dynamic images, signal-time courses and parametric maps (all data).

RESULTS

DRO data showed that the bias and precision of parameter maps after MDR are indistinguishable from motion-free data. MDR led to reduction in HD (HDunregistered = 9.98 ± 9.76, HDregistered = 1.63 ± 0.49). Visual inspection showed that MDR effectively removed motion effects in the dynamic data, leading to a clear improvement in anatomical delineation on parametric maps and a reduction in motion-induced oscillations on signal-time courses.

DISCUSSION

MDR provides effective motion correction of MRR in synthetic and patient data. Future work is needed to compare the performance against other more established methods.

Improving the robustness of MOLLI T1 maps with a dedicated motion correction algorithm

Myocardial tissue T1 constitutes a reliable indicator of several heart diseases related to extracellular changes (e.g. edema, fibrosis) as well as fat, iron and amyloid content. Magnetic resonance (MR) T1-mapping is typically achieved by pixel-wise exponential fitting of a series of inversion or saturation recovery measurements. Good anatomical alignment between these measurements is essential for accurate T1 estimation. Motion correction is recommended to improve alignment. However, in the case of inversion recovery sequences, this correction is compromised by the intrinsic contrast variation between frames. A model-based, non-rigid motion correction method for MOLLI series was implemented and validated on a large database of cardiac clinical cases (n = 186). The method relies on a dedicated similarity metric that accounts for the intensity changes caused by T1 magnetization relaxation. The results were compared to uncorrected series and to the standard motion correction included in the scanner. To automate the quantitative analysis of results, a custom data alignment metric was defined. Qualitative evaluation was performed on a subset of cases to confirm the validity of the new metric. Motion correction caused noticeable (i.e. > 5%) performance degradation in 12% of cases with the standard method, compared to 0.3% with the new dedicated method. The average alignment quality was 85% ± 9% with the default correction and 90% ± 7% with the new method. The results of the qualitative evaluation were found to correlate with the quantitative metric. In conclusion, a dedicated motion correction method for T1 mapping MOLLI series has been evaluated on a large database of clinical cardiac MR cases, confirming its increased robustness with respect to the standard method implemented in the scanner.

Native T1 mapping and extracellular volume fraction for differentiation of myocardial diseases from normal CMR controls in routine clinical practice

Background

This study aimed to determine native T1 and extracellular volume fraction (ECV) in distinct types of myocardial disease, including amyloidosis, dilated cardiomyopathy (DCM), hypertrophic cardiomyopathy (HCM), myocarditis and coronary artery disease (CAD), compared to controls.

Methods

We retrospectively enrolled patients with distinct types of myocardial disease, CAD patients, and control group (no known heart disease and negative CMR study) who underwent 3.0 Tesla CMR with routine T1 mapping. The region of interest (ROI) was drawn in the myocardium of the mid left ventricular (LV) short axis slice and at the interventricular septum of mid LV slice. ECV was calculated by actual hematocrit (Hct) and synthetic Hct. T1 mapping and ECV was compared between myocardial disease and controls, and between CAD and controls. Diagnostic yield and cut-off values were assessed.

Results

A total of 1188 patients were enrolled. The average T1 values in the control group were 1304 ± 42 ms at septum, and 1294 ± 37 ms at mid LV slice. The average T1 values in patients with myocardial disease and CAD were significantly higher than in controls (1441 ± 72, 1349 ± 59, 1345 ± 59, 1355 ± 56, and 1328 ± 54 ms for septum of amyloidosis, DCM, HCM, myocarditis, and CAD). Native T1 of the mid LV level and ECV at septum and mid LV with actual and synthetic Hct of patients with myocardial disease or CAD were significantly higher than in controls.

Conclusions

Although native T1 and ECV of patients with cardiomyopathy and CAD were significantly higher than controls, the values overlapped. The greatest clinical utilization was found for the amyloidosis group.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12872-021-02086-3.

Endurance exercise and the risk of cardiovascular pathology in men: a comparison between lifelong and late-onset endurance training and a non-athletic lifestyle - rationale and design of the Master@Heart study, a prospective cohort trial

Introduction

Low and moderate endurance exercise is associated with better control of cardiovascular risk factors, a decreased risk of coronary artery disease and atrial fibrillation (AF). There is, however, a growing proportion of individuals regularly performing strenuous and prolonged endurance exercise in which the health benefits have been challenged. Higher doses of endurance exercise have been associated with a greater coronary atherosclerotic plaque burden, risk of AF and myocardial fibrosis (MF).

Methods and analysis

Master@Heart is a multicentre prospective cohort study aiming to assess the incidence of coronary atherosclerosis, AF and MF in lifelong endurance athletes compared to late-onset endurance athletes (initiation of regular endurance exercise after the age of 30 years) and healthy non-athletes.

The primary endpoint is the incidence of mixed coronary plaques. Secondary endpoints include coronary calcium scores, coronary stenosis >50%, the prevalence of calcified and soft plaques and AF and MF presence. Tertiary endpoints include ventricular arrhythmias, left and right ventricular function at rest and during exercise, arterial stiffness and carotid artery intima media thickness.

Two hundred male lifelong athletes, 200 late-onset athletes and 200 healthy non-athletes aged 45–70 will undergo comprehensive cardiovascular phenotyping using CT, coronary angiography, echocardiography, cardiac MRI, 12-lead ECG, exercise ECG and 24-hour Holter monitoring at baseline. Follow-up will include online tracking of sports activities, telephone calls to assess clinical events and a 7-day ECG recording after 1 year.

Ethics and dissemination

Local ethics committees approved the Master@Heart study. The trial was launched on 18 October 2018, recruitment is complete and inclusions are ongoing.

Trial registration number

NCT03711539.

Free breathing lung T1 mapping using image registration in patients with idiopathic pulmonary fibrosis

PURPOSE

To assess the use of image registration for correcting respiratory motion in free breathing lung T₁ mapping acquisition in patients with idiopathic pulmonary fibrosis (IPF).

THEORY AND METHODS

The method presented used image registration to synthetic images during postprocessing to remove respiratory motion. Synthetic images were generated from a model of the inversion recovery signal of the acquired images that incorporated a periodic lung motion model. Ten healthy volunteers and 19 patients with IPF underwent 2D Look-Locker T₁ mapping acquisition at 1.5T during inspiratory breath-hold and free breathing. Eight healthy volunteers and seven patients with IPF underwent T₁ mapping acquisition during expiratory breath-hold. Fourteen patients had follow-up scanning at 6 months. Dice similarity coefficient (DSC) was used to evaluate registration efficacy.

RESULTS

Image registration increased image DSC (P < .001) in the free breathing inversion recovery images. Lung T₁ measured during a free breathing acquisition was lower in patients with IPF when compared with healthy controls (inspiration: P = .238; expiration: P = .261; free breathing: P = .021). Measured lung T₁ was higher in expiration breath-hold than inspiration breath-hold in healthy volunteers (P < .001) but not in patients with IPF (P = .645). There were no other significant differences between lung T₁ values within subject groups.

CONCLUSIONS

The registration technique significantly reduced motion in the Look-Locker images acquired during free breathing and may improve the robustness of lung T₁ mapping in patients who struggle to hold their breath. Lung T₁ measured during a free breathing acquisition was significantly lower in patients with IPF when compared with healthy controls.

Improved fetal blood oxygenation and placental estimated measurements of diffusion-weighted MRI using data-driven Bayesian modeling

PURPOSE

Motion correction in placental DW-MRI is challenging due to maternal breathing motion, maternal movements, and rapid intensity changes. Parameter estimates are usually obtained using least-squares methods for voxel-wise fitting; however, they typically give noisy estimates due to low signal-to-noise ratio. We introduce a model-driven registration (MDR) technique which incorporates a placenta-specific signal model into the registration process, and we present a Bayesian approach for Diffusion-rElaxation Combined Imaging for Detailed placental Evaluation model to obtain individual and population trends in estimated parameters.

METHODS

MDR exploits the fact that a placenta signal model is available and thus we incorporate it into the registration to generate a series of target images. The proposed registration method is compared to a pre-existing method used for DCE-MRI data making use of principal components analysis. The Bayesian shrinkage prior (BSP) method has no user-defined parameters and therefore measures of parameter variation in a region of interest are determined by the data alone. The MDR method and the Bayesian approach were evaluated on 10 control 4D DW-MRI singleton placental data.

RESULTS

MDR method improves the alignment of placenta data compared to the pre-existing method. It also shows a further reduction of the residual error between the data and the fit. BSP approach showed higher precision leading to more clearly apparent spatial features in the parameter maps. Placental fetal oxygen saturation (FO₂ ) showed a negative linear correlation with gestational age.

CONCLUSIONS

The proposed pipeline provides a robust framework for registering DW-MRI data and analyzing longitudinal changes of placental function.