Comparison of left ventricular strains and torsion derived from feature tracking and DENSE CMR

Systematic review and meta-analysis for the value of cardiac magnetic resonance strain to predict cardiac outcomes

Cardiac magnetic resonance (CMR) is the gold standard for the diagnostic classification and risk stratification in most patients with cardiac disorders. The aim of the present study was to investigate the ability of Strain-encoded MR (SENC) for the prediction of major adverse cardiovascular events (MACE). A systematic review and meta-analysis was performed according to the PRISMA Guidelines, including patients with or without cardiovascular disease and asymptomatic individuals. Myocardial strain by HARP were used as pulse sequences in 1.5 T scanners. Published literature in MEDLINE (PubMed) and Cochrane's databases were explored before February 2023 for studies assessing the clinical utility of myocardial strain by Harmonic Phase Magnetic Resonance Imaging (HARP), Strain-encoded MR (SENC) or fast-SENC. In total, 8 clinical trials (4 studies conducted in asymptomatic individuals and 4 in patients with suspected or known cardiac disease) were included in this systematic review, while 3 studies were used for our meta-analysis, based on individual patient level data. Kaplan-Meier analysis and Cox proportional hazard models were used, testing the ability of myocardial strain by HARP and SENC/fast-SENC for the prediction of MACE. Strain enabled risk stratification in asymptomatic individuals, predicting MACE and the development of incident heart failure. Of 1332 patients who underwent clinically indicated CMR, including SENC or fast-SENC acquisitions, 19 patients died, 28 experienced non-fatal infarctions, 52 underwent coronary revascularization and 86 were hospitalized due to heart failure during median 22.4 (17.2-28.5) months of follow-up. SENC/fast-SENC, predicted both all-cause mortality and MACE with high accuracy (HR = 3.0, 95% CI = 1.2-7.6, p = 0.02 and HR = 4.1, 95% CI = 3.0-5.5, respectively, p < 0.001). Using hierarchical Cox-proportional hazard regression models, SENC/fast-SENC exhibited incremental value to clinical data and conventional CMR parameters. Reduced myocardial strain predicts of all-cause mortality and cardiac outcomes in symptomatic patients with a wide range of ischemic or non-ischemic cardiac diseases, whereas in asymptomatic individuals, reduced strain was a precursor of incident heart failure.

Measuring Cardiac Dyssynchrony with DENSE (Displacement Encoding with Stimulated Echoes)-A Systematic Review

Background

In this review, we introduce the displacement encoding with stimulated echoes (DENSE) method for measuring myocardial dyssynchrony using cardiovascular magnetic resonance (CMR) imaging. We provide an overview of research findings related to DENSE from the past two decades and discuss other techniques used for dyssynchrony evaluation. Additionally, the review discusses the potential uses of DENSE in clinical practice.

Methods

A search was conducted to identify relevant articles published from January 2000 through January 2023 using the Scopus, Web of Science, PubMed and Cochrane databases. The following search term was used: (DENSE OR 'displacement encoding with stimulated echoes' OR CURE) AND (dyssynchrony* OR asynchron* OR synchron*) AND (MRI OR 'magnetic resonance' OR CMR).

Results

After removing duplicates, researchers screened a total of 174 papers. Papers that were not related to the topic, reviews, general overview articles and case reports were excluded, leaving 35 articles for further analysis. Of these, 14 studies focused on cardiac dyssynchrony estimation with DENSE, while the remaining 21 studies served as background material. The studies used various methods for presenting synchronicity, such as circumferential uniformity ratio estimate (CURE), CURE-singular value decomposition (SVD), radial uniformity ratio estimate (RURE), longitudinal uniformity ratio estimate (LURE), time to onset of shortening (TOS) and dyssynchrony index (DI). Most of the dyssynchrony studies concentrated on human heart failure, but congenital heart diseases and obesity were also evaluated. The researchers found that DENSE demonstrated high reproducibility and was found useful for detecting cardiac resynchronisation therapy (CRT) responders, optimising CRT device settings and assessing right ventricle synchronicity. In addition, studies showed a correlation between cardiac fibrosis and mechanical dyssynchrony in humans, as well as a decrease in the synchrony of contraction in the left ventricle in obese mice.

Conclusions

DENSE shows promise as a tool for quantifying myocardial function and dyssynchrony, with advantages over other cardiac dyssynchrony evaluation methods. However, there remain challenges related to DENSE due to the relatively time-consuming imaging and analysis process. Improvements in imaging and analysing technology, as well as possible artificial intelligence solutions, may help overcome these challenges and lead to more widespread clinical use of DENSE.

The Interfield Strength Agreement of Left Ventricular Strain Measurements at 1.5 T and 3 T Using Cardiac MRI Feature Tracking

BACKGROUND

Left ventricular (LV) strain measurements can be derived using cardiac MRI from routinely acquired balanced steady-state free precession (bSSFP) cine images.

PURPOSE

To compare the interfield strength agreement of global systolic strain, peak strain rates and artificial intelligence (AI) landmark-based global longitudinal shortening at 1.5 T and 3 T.

STUDY TYPE

Prospective.

SUBJECTS

A total of 22 healthy individuals (mean age 36 ± 12 years; 45% male) completed two cardiac MRI scans at 1.5 T and 3 T in a randomized order within 30 minutes. FIELD STRENGTH/SEQUENCE: bSSFP cine images at 1.5 T and 3 T.

ASSESSMENT

Two software packages, Tissue Tracking (cvi42, Circle Cardiovascular Imaging) and QStrain (Medis Suite, Medis Medical Imaging Systems), were used to derive LV global systolic strain in the longitudinal, circumferential and radial directions and peak (systolic, early diastolic, and late diastolic) strain rates. Global longitudinal shortening and mitral annular plane systolic excursion (MAPSE) were measured using an AI deep neural network model.

STATISTICAL TESTS

Comparisons between field strengths were performed using Wilcoxon signed-rank test (P value < 0.05 considered statistically significant). Agreement was determined using intraclass correlation coefficients (ICCs) and Bland-Altman plots.

RESULTS

Minimal bias was seen in all strain and strain rate measurements between field strengths. Using Tissue Tracking, strain and strain rate values derived from long-axis images showed poor to fair agreement (ICC range 0.39-0.71), whereas global longitudinal shortening and MAPSE showed good agreement (ICC = 0.81 and 0.80, respectively). Measures derived from short-axis images showed good to excellent agreement (ICC range 0.78-0.91). Similar results for the agreement of strain and strain rate measurements were observed with QStrain.

CONCLUSION

The interfield strength agreement of short-axis derived LV strain and strain rate measurements at 1.5 T and 3 T was better than those derived from long-axis images; however, the agreement of global longitudinal shortening and MAPSE was good.

EVIDENCE LEVEL

2 TECHNICAL EFFICACY: Stage 2.

Automated segmentation of long and short axis DENSE cardiovascular magnetic resonance for myocardial strain analysis using spatio-temporal convolutional neural networks

BACKGROUND

Cine Displacement Encoding with Stimulated Echoes (DENSE) facilitates the quantification of myocardial deformation, by encoding tissue displacements in the cardiovascular magnetic resonance (CMR) image phase, from which myocardial strain can be estimated with high accuracy and reproducibility. Current methods for analyzing DENSE images still heavily rely on user input, making this process time-consuming and subject to inter-observer variability. The present study sought to develop a spatio-temporal deep learning model for segmentation of the left-ventricular (LV) myocardium, as spatial networks often fail due to contrast-related properties of DENSE images.

METHODS

2D + time nnU-Net-based models have been trained to segment the LV myocardium from DENSE magnitude data in short- and long-axis images. A dataset of 360 short-axis and 124 long-axis slices was used to train the networks, from a combination of healthy subjects and patients with various conditions (hypertrophic and dilated cardiomyopathy, myocardial infarction, myocarditis). Segmentation performance was evaluated using ground-truth manual labels, and a strain analysis using conventional methods was performed to assess strain agreement with manual segmentation. Additional validation was performed using an externally acquired dataset to compare the inter- and intra-scanner reproducibility with respect to conventional methods.

RESULTS

Spatio-temporal models gave consistent segmentation performance throughout the cine sequence, while 2D architectures often failed to segment end-diastolic frames due to the limited blood-to-myocardium contrast. Our models achieved a DICE score of 0.83 ± 0.05 and a Hausdorff distance of 4.0 ± 1.1 mm for short-axis segmentation, and 0.82 ± 0.03 and 7.9 ± 3.9 mm respectively for long-axis segmentations. Strain measurements obtained from automatically estimated myocardial contours showed good to excellent agreement with manual pipelines, and remained within the limits of inter-user variability estimated in previous studies.

CONCLUSION

Spatio-temporal deep learning shows increased robustness for the segmentation of cine DENSE images. It provides excellent agreement with manual segmentation for strain extraction. Deep learning will facilitate the analysis of DENSE data, bringing it one step closer to clinical routine.

Effects of spatial and temporal resolution on cardiovascular magnetic resonance feature tracking measurements using a simple realistic numerical phantom

OBJECTIVES

To develop a single-slice numerical phantom with known myocardial motion, at several temporal and in-plane spatial resolutions, for testing and comparison of Cardiovascular Magnetic Resonance (CMR) feature tracking (FT) software.

METHODS

The phantom was developed based on CMR acquisitions of one volunteer (acquired cine, tagging cine, T1 map, T2 map, proton density weighted image). The numerical MRI simulator JEMRIS was used, and the phantom was generated at several in-plane spatial resolutions (1.4 × 1.4 mm2 to 3.0 × 3.0 mm2) and temporal resolutions (20 to 40 cardiac phases). Two feature tracking software packages were tested: Medical Image Tracking Toolbox (MITT) and two versions of cvi42 (v5.3.8 and v5.13.7). The effect of resolution on strain results was investigated with reference to ground-truth radial and circumferential strain.

RESULTS

Peak radial strain was consistently undermeasured more for cvi42 v5.13.7 than for v5.3.8. Increased pixel size produced a trend of increased difference from ground-truth peak strain, with the largest changes for cvi42 obtained using v5.13.7 between 1.4 × 1.4 mm2 and 3.0 × 3.0 mm2, at 9.17 percentage points (radial) and 8.42 percentage points (circumferential).

CONCLUSIONS

The results corroborate the presence of intervendor differences in feature tracking results and show the magnitude of strain differences between software versions.

ADVANCES IN KNOWLEDGE

This study shows how temporal and in-plane spatial resolution can affect feature tracking with reference to the ground-truth strain of a numerical phantom. Results reaffirm the need for numerical phantom development for the validation and testing of FT software.

Biventricular shape modes discriminate pulmonary valve replacement in tetralogy of Fallot better than imaging indices

Current indications for pulmonary valve replacement (PVR) in repaired tetralogy of Fallot (rTOF) rely on cardiovascular magnetic resonance (CMR) image-based indices but are inconsistently applied, lead to mixed outcomes, and remain debated. This study aimed to test the hypothesis that specific markers of biventricular shape may discriminate differences between rTOF patients who did and did not require subsequent PVR better than standard imaging indices. In this cross-sectional retrospective study, biventricular shape models were customized to CMR images from 84 rTOF patients. A statistical atlas of end-diastolic shape was constructed using principal component analysis. Multivariate regression was used to quantify shape mode and imaging index associations with subsequent intervention status (PVR, n = 48 vs. No-PVR, n = 36), while accounting for confounders. Clustering analysis was used to test the ability of the most significant shape modes and imaging indices to discriminate PVR status as evaluated by a Matthews correlation coefficient (MCC). Geometric strain analysis was also conducted to assess shape mode associations with systolic function. PVR status correlated significantly with shape modes associated with right ventricular (RV) apical dilation and left ventricular (LV) dilation (p < 0.01), RV basal bulging and LV conicity (p < 0.05), and pulmonary valve dilation (p < 0.01). PVR status also correlated significantly with RV ejection fraction (p < 0.05) and correlated marginally with LV end-systolic volume index (p < 0.07). Shape modes discriminated subsequent PVR better than standard imaging indices (MCC = 0.49 and MCC = 0.28, respectively) and were significantly associated with RV and LV radial systolic strain. Biventricular shape modes discriminated differences between patients who did and did not require subsequent PVR better than standard imaging indices in current use. These regional features of cardiac morphology may provide insight into adaptive vs. maladaptive types of structural remodeling and point toward an improved quantitative, patient-specific assessment tool for clinical use.

Reproducibility of Systolic Strain in Mice Using Cardiac Magnetic Resonance Feature Tracking of Black-Blood Cine Images

PURPOSE

Mouse models are widely utilized to enhance our understanding of cardiac disease. The goal of this study is to investigate the reproducibility of strain parameters that were measured in mice using cardiac magnetic resonance (CMR) feature-tracking (CMR42, Canada).

METHODS

We retrospectively analyzed black-blood CMR datasets from thirteen C57BL/6 B6.SJL-CD45.1 mice (N = 10 female, N = 3 male) that were imaged previously. The circumferential, longitudinal, and radial (Ecc, Ell, and Err, respectively) parameters of strain were measured in the mid-ventricular region of the left ventricle. Intraobserver and interobserver reproducibility were assessed for both the end-systolic (ES) and peak strain.

RESULTS

The ES strain had larger intraclass correlation coefficient (ICC) values when compared to peak strain, for both the intraobserver and interobserver reproducibility studies. Specifically, the intraobserver study showed excellent reproducibility for all three ES strain parameters, namely, Ecc (ICC 0.95, 95% CI 0.83-0.98), Ell (ICC 0.90, 95% CI 0.59-0.97), and Err (ICC 0.92, 95% CI 0.73-0.97). This was also the case for the interobserver study, namely, Ecc (ICC 0.92, 95% CI 0.60-0.98), Ell (ICC 0.76, 95% CI 0.33-0.93), and Err (ICC 0.93, 95% CI 0.68-0.98). Additionally, the coefficient of variation values were all < 10%.

CONCLUSION

The results of this preliminary study showed excellent reproducibility for all ES strain parameters, with good to excellent reproducibility for the peak strain parameters. Moreover, all ES strain parameters had larger ICC values than the peak strain. In general, these results imply that feature-tracking with CMR42 software and black-blood cine images can be reliably used to assess strain patterns in mice.

Right-left ventricular shape variations in tetralogy of Fallot: associations with pulmonary regurgitation

BACKGROUND

Relationships between right ventricular (RV) and left ventricular (LV) shape and function may be useful in determining optimal timing for pulmonary valve replacement in patients with repaired tetralogy of Fallot (rTOF). However, these are multivariate and difficult to quantify. We aimed to quantify variations in biventricular shape associated with pulmonary regurgitant volume (PRV) in rTOF using a biventricular atlas.

METHODS

In this cross-sectional retrospective study, a biventricular shape model was customized to cardiovascular magnetic resonance (CMR) images from 88 rTOF patients (median age 16, inter-quartile range 11.8-24.3 years). Morphometric scores quantifying biventricular shape at end-diastole and end-systole were computed using principal component analysis. Multivariate linear regression was used to quantify biventricular shape associations with PRV, corrected for age, sex, height, and weight. Regional associations were confirmed by univariate correlations with distances and angles computed from the models, as well as global systolic strains computed from changes in arc length from end-diastole to end-systole.

RESULTS

PRV was significantly associated with 5 biventricular morphometric scores, independent of covariates, and accounted for 12.3% of total shape variation (p < 0.05). Increasing PRV was associated with RV dilation and basal bulging, in conjunction with decreased LV septal-lateral dimension (LV flattening) and systolic septal motion towards the RV (all p < 0.05). Increased global RV radial, longitudinal, circumferential and LV radial systolic strains were significantly associated with increased PRV (all p < 0.05).

CONCLUSION

A biventricular atlas of rTOF patients quantified multivariate relationships between left-right ventricular morphometry and wall motion with pulmonary regurgitation. Regional RV dilation, LV reduction, LV septal-lateral flattening and increased RV strain were all associated with increased pulmonary regurgitant volume. Morphometric scores provide simple metrics linking mechanisms for structural and functional alteration with important clinical indices.

Regional end-systolic circumferential strain demonstrates compensatory segmental contractile function in patients with ST-segment elevation myocardial infarction

Myocardial strain has shown tremendous promise as a potential diagnostic tool for characterizing ventricular function. With regards to myocardial infarction, global circumferential strain (CS) can be used to assess overall function, while regional CS can identify local alterations in contractility. Currently, there is a lack of data related to regional strain in patients with ST-segment elevation myocardial infarction (STEMI). Thus, the goal of this study was to quantify regional strain patterns in STEMI and normal control patients, measuring both peak CS and end-systolic (ES) CS in the mid-ventricular region. This was done by conducting cardiac magnetic resonance (CMR) imaging acutely after STEMI patients underwent primary percutaneous coronary intervention. The CMR datasets were then analyzed using feature-tracking of the cine images. The patients were broken into three groups: (1) control patients (N = 18), (2) STEMI patients with ejection fraction (EF) ≥ 50% (N = 20), and (3) STEMI patients with EF < 50% (N = 20). The key result of the analysis was that ES CS detected a significant increase in the magnitude of strain in the non-infarcted tissue of STEMI patients with EF ≥ 50% when compared to STEMI patients with EF < 50%, whereas peak CS did not detect any differences. This implies that the tissue in this region is contracting more strongly compared to non-infarcted tissue in STEMI patients with EF < 50%. Thus, regional ES CS could potentially be utilized as a diagnostic tool for assessing STEMI patients, by detecting regional changes in contractility after PCI, which could assist in treatment planning.

3D-Encoded DENSE MRI with Zonal Excitation for Quantifying Biventricular Myocardial Strain During a Breath-Hold

PURPOSE

Right ventricular (RV) function is increasingly recognized for its prognostic value in many disease states. As with the left ventricle (LV), strain-based measurements may have better prognostic value than typical chamber volumes or ejection fraction. Complete functional characterization of the RV requires high-resolution, 3D displacement tracking methods, which have been prohibitively challenging to implement. Zonal excitation during Displacement ENcoding with Stimulated Echoes (DENSE) magnetic resonance imaging (MRI) has helped reduce scan time for 2D LV strain quantification. We hypothesized that zonal excitation could alternatively be used to reproducibly acquire higher resolution, 3D-encoded DENSE images for quantification of bi-ventricular strain within a single breath-hold.

METHODS

We modified sequence parameters for a 3D zonal excitation DENSE sequence to achieve in-plane resolution < 2 mm and acquired two sets of images in eight healthy adult male volunteers with median (IQR) age 32.5 (32.0-33.8) years. We assessed the inter-test reproducibility of this technique, and compared computed strains and torsion with previously published data.

RESULTS

Data for one subject was excluded based on image artifacts. Reproducibility for LV (CoV: 6.1-9.0%) and RV normal strains (CoV: 6.3-8.2%) and LV torsion (CoV = 7.1%) were all very good. Reproducibility of RV torsion was lower (CoV = 16.7%), but still within acceptable limits. Computed global strains and torsion were within reasonable agreement with published data, but further studies in larger cohorts are needed to confirm.

CONCLUSION

Reproducible acquisition of 3D-encoded biventricular myocardial strain data in a breath-hold is feasible using DENSE with zonal excitation.

A comprehensive comparison between shortest-path HARP refinement, SinMod, and DENSEanalysis processing tools applied to CSPAMM and DENSE images

We addressed comprehensively the performance of Shortest-Path HARP Refinement (SP-HR), SinMod, and DENSEanalysis using 2D slices of synthetic CSPAMM and DENSE images with realistic contrasts obtained from 3D phantoms. The three motion estimation techniques were interrogated under ideal and no-ideal conditions (with MR induced artifacts, noise, and through-plane motion), considering several resolutions and noise levels. Under noisy conditions, and for isotropic pixel sizes of 1.5 mm and 3.0 mm in CSPAMM and DENSE images respectively, the nRMSE obtained for the circumferential and radial strain components were 10.7 ± 10.8% and 25.5 ± 14.8% using SP-HR, 11.9 ± 2.5% and 29.3 ± 6.5% using SinMod, and 6.4 ± 2.0% and 18.2 ± 4.6% using DENSEanalysis. Overall, the results showed that SP-HR tends to fail for large tissue motions, whereas SinMod and DENSEanalysis gave accurate displacement and strain field estimations, being the last which performed the best.

Regional Brain Tissue Displacement and Strain is Elevated in Subjects with Chiari Malformation Type I Compared to Healthy Controls: A Study Using DENSE MRI

While the degree of cerebellar tonsillar descent is considered the primary radiologic marker of Chiari malformation type I (CMI), biomechanical forces acting on the brain tissue in CMI subjects are less studied and poorly understood. In this study, regional brain tissue displacement and principal strains in 43 CMI subjects and 25 controls were quantified using a magnetic resonance imaging (MRI) methodology known as displacement encoding with stimulated echoes (DENSE). Measurements from MRI were obtained for seven different brain regions-the brainstem, cerebellum, cingulate gyrus, corpus callosum, frontal lobe, occipital lobe, and parietal lobe. Mean displacements in the cerebellum and brainstem were found to be 106 and 64% higher, respectively, for CMI subjects than controls (p < .001). Mean compression and extension strains in the cerebellum were 52 and 50% higher, respectively, in CMI subjects (p < .001). Brainstem mean extension strain was 41% higher in CMI subjects (p < .001), but no significant difference in compression strain was observed. The other brain structures revealed no significant differences between CMI and controls. These findings demonstrate that brain tissue displacement and strain in the cerebellum and brainstem might represent two new biomarkers to distinguish between CMI subjects and controls.

Impact of pulmonary valve replacement on left ventricular rotational mechanics in repaired tetralogy of Fallot

BACKGROUND

In repaired tetralogy of Fallot (rTOF), abnormal left ventricular (LV) rotational mechanics are associated with adverse clinical outcomes. We performed a comprehensive analysis of LV rotational mechanics in rTOF patients using cardiac magnetic resonance (CMR) prior to and following surgical pulmonary valve replacement (PVR).

METHODS

In this single center retrospective study, we identified rTOF patients who (1) had both a CMR ≤ 1 year before PVR and ≤ 5 years after PVR, (2) had no other intervening procedure between CMRs, (3) had a body surface area > 1.0 m2 at CMR, and (4) had images suitable for feature tracking analysis. These subjects were matched to healthy age- and sex-matched control subjects. CMR feature tracking analysis was performed on a ventricular short-axis stack of balanced steady-state free precession images. Measurements included LV basal and apical rotation, twist, torsion, peak systolic rates of rotation and torsion, and timing of events. Associations with LV torsion were assessed.

RESULTS

A total of 60 rTOF patients (23.6 ± 7.9 years, 52% male) and 30 healthy control subjects (20.8 ± 3.1 years, 50% male) were included. Compared with healthy controls, rTOF patients had lower apical and basal rotation, twist, torsion, and systolic rotation rates, and these parameters peaked earlier in systole. The only parameters that were correlated with LV torsion were right ventricular (RV) end-systolic volume (r = - 0.28, p = 0.029) and RV ejection fraction (r = 0.26, p = 0.044). At a median of 1.0 year (IQR 0.5-1.7) following PVR, there was no significant change in LV rotational parameters versus pre-PVR despite reductions in RV volumes, RV mass, pulmonary regurgitation, and RV outflow tract obstruction.

CONCLUSION

In this comprehensive study of CMR-derived LV rotational mechanics in rTOF patients, rotation, twist, and torsion were diminished compared to controls and did not improve at a median of 1 year after PVR despite favorable RV remodeling.

Comparing cardiovascular magnetic resonance strain software packages by their abilities to discriminate outcomes in patients with heart failure with preserved ejection fraction

BACKGROUND

Cardiovascular magnetic resonance (CMR) myocardial strain analysis using feature tracking (FT) is an increasingly popular method to assess cardiac function. However, different software packages produce different strain values from the same images and there is little guidance regarding which software package would be the best to use. We explored a framework under which different software packages could be compared and used based on their abilities to differentiate disease from health and differentiate disease severity based on outcome.

METHOD

To illustrate this concept, we compared 4-chamber left ventricular (LV) peak longitudinal strain (GLS) analyzed from retrospective electrocardiogram gated cine imaging performed on 1.5 T CMR scanners using three CMR post-processing software packages in their abilities to discriminate a group of 45 patients with heart failure with preserved ejection fraction (HFpEF) from 26 controls without cardiovascular disease and to discriminate disease severity based on outcomes. The three different post-processing software used were SuiteHeart, cvi42, and DRA-Trufistrain.

RESULTS

All three software packages were able to distinguish HFpEF patients from controls. 4-chamber peak GLS by SuiteHeart was shown to be a better discriminator of adverse outcomes in HFpEF patients than 4-chamber GLS derived from cvi42 or DRA-Trufistrain.

CONCLUSION

We illustrated a framework to compare feature tracking GLS derived from different post-processing software packages. Publicly available imaging data sets with outcomes would be important to validate the growing number of CMR-FT software packages.

Multivendor comparison of global and regional 2D cardiovascular magnetic resonance feature tracking strains vs tissue tagging at 3T

BACKGROUND

Cardiovascular magnetic resonance (CMR) 2D feature tracking (FT) left ventricular (LV) myocardial strain has seen widespread use to characterize myocardial deformation. Yet, validation of CMR FT measurements remains scarce, particularly for regional strain. Therefore, we aimed to perform intervendor comparison of 3 different FT software against tagging.

METHODS

In 61 subjects (18 healthy subjects, 18 patients with chronic myocardial infarction, 15 with dilated cardiomyopathy, and 10 with LV hypertrophy due to hypertrophic cardiomyopathy or aortic stenosis) were prospectively compared global (G) and regional transmural peak-systolic Lagrangian longitudinal (LS), circumferential (CS) and radial strains (RS) by 3 FT software (cvi42, Segment, and Tomtec) among each other and with tagging at 3T. We also evaluated the ability of regional LS, CS, and RS by different FT software vs tagging to identify late gadolinium enhancement (LGE) in the 18 infarct patients.

RESULTS

GLS and GCS by all 3 software had an excellent agreement among each other (ICC = 0.94-0.98 for GLS and ICC = 0.96-0.98 for GCS respectively) and against tagging (ICC = 0.92-0.94 for GLS and ICC = 0.88-0.91 for GCS respectively), while GRS showed inconsistent agreement between vendors (ICC 0.10-0.81). For regional LS, the agreement was good (ICC = 0.68) between 2 vendors but less vs the 3rd (ICC 0.50-0.59) and moderate to poor (ICC 0.44-0.47) between all three FT software and tagging. Also, for regional CS agreement between 2 software was higher (ICC = 0.80) than against the 3rd (ICC = 0.58-0.60), and both better agreed with tagging (ICC = 0.70-0.72) than the 3rd (ICC = 0.57). Regional RS had more variation in the agreement between methods ranging from good (ICC = 0.75) to poor (ICC = 0.05). Finally, the accuracy of scar detection by regional strains differed among the 3 FT software. While the accuracy of regional LS was similar, CS by one software was less accurate (AUC 0.68) than tagging (AUC 0.80, p < 0.006) and RS less accurate (AUC 0.578) than the other two (AUC 0.76 and 0.73, p < 0.02) to discriminate segments with LGE.

CONCLUSIONS

We confirm good agreement of CMR FT and little intervendor difference for GLS and GCS evaluation, with variable agreement for GRS. For regional strain evaluation, intervendor difference was larger, especially for RS, and the diagnostic performance varied more substantially among different vendors for regional strain analysis.

HARP-I: A Harmonic Phase Interpolation Method for the Estimation of Motion From Tagged MR Images

We proposed a novel method called HARP-I, which enhances the estimation of motion from tagged Magnetic Resonance Imaging (MRI). The harmonic phase of the images is unwrapped and treated as noisy measurements of reference coordinates on a deformed domain, obtaining motion with high accuracy using Radial Basis Functions interpolations. Results were compared against Shortest Path HARP Refinement (SP-HR) and Sine-wave Modeling (SinMod), two harmonic image-based techniques for motion estimation from tagged images. HARP-I showed a favorable similarity with both methods under noise-free conditions, whereas a more robust performance was found in the presence of noise. Cardiac strain was better estimated using HARP-I at almost any motion level, giving strain maps with less artifacts. Additionally, HARP-I showed better temporal consistency as a new method was developed to fix phase jumps between frames. In conclusion, HARP-I showed to be a robust method for the estimation of motion and strain under ideal and non-ideal conditions.

Layer-Specific Strain Is Preload Dependent: Comparison between Speckle-Tracking Echocardiography and Cardiac Magnetic Resonance Feature-Tracking

BACKGROUND

Speckle-tracking echocardiographic (STE) imaging and cardiac magnetic resonance feature-tracking (CMR-FT) are novel imaging techniques enabling layer-specific quantification of myocardial deformation. Conventional echocardiographic parameters are load dependent, but few studies have investigated the effects of loading conditions on STE and CMR-FT layer-specific strain and the interchangeability of the two modalities. The aim of this study was to evaluate the effects of acute preload augmentation by saline infusion on STE and CMR-FT longitudinal and circumferential layer-specific strain parameters and their intermodal agreement.

METHODS

A total of 80 subjects, including 41 control subjects (mean age, 40 ± 12 years; 49% men) and 39 patients with cardiac disease (mean age, 47 ± 15 years; 92% men) were examined using STE and CMR-FT layer-specific strain analysis before and after saline infusion (median, 2.0 L) with quantification of transmural global longitudinal strain (GLS), epicardial GLS, endocardial GLS, transmural global circumferential strain (GCS), epicardial GCS, and endocardial GCS in addition to epicardial-endocardial gradients. Bland-Altman plots and Pearson correlation coefficients were used to evaluate agreement between the two modalities across all strain parameters.

RESULTS

Acute saline infusion increased all STE and CMR-FT layer-specific strain parameters in both groups. STE and CMR-FT GLS increased by 1.4 ± 1.5% and 1.5 ± 2.0% (P < .001) in control subjects and by 0.9 ± 1.8% and 0.9 ± 1.9% (P < .001) in patients with cardiac disease. STE and CMR-FT GCS increased by 2.0 ± 2.2% and 1.8 ± 2.3% (P < .001) in control subjects and by 1.8 ± 2.3% and 1.7 ± 3.6% in patients with cardiac disease (P < .001 and P = .03). STE longitudinal strain correlated strongly with corresponding CMR-FT longitudinal strain (GLS, epicardial GLS, and endocardial GLS: r = 0.81, r = 0.82, and r = 0.81, respectively) despite poor intermodal agreement (bias ± limits of agreement, -2.84 ± 4.06%, 0.16 ± 3.68%, and 2.33 ± 3.52%, respectively) whereas GCS, epicardial GCS, and endocardial GCS correlated weakly between the two modalities (r = 0.28, r = 0.19, and r = 0.34, respectively) and displayed poor intermodal agreement (bias ± limits of agreement, -1.33 ± 6.86%, 4.43 ± 6.49%, and -9.92 ± 8.55%, respectively).

CONCLUSIONS

STE and CMR-FT longitudinal and circumferential layer-specific strain parameters are preload dependent in both control subjects and patients with cardiac disease. STE and CMR-FT longitudinal layer-specific strain parameters are strongly correlated, whereas circumferential layer-specific strain parameters are weakly correlated. STE and CMR-FT longitudinal and circumferential strain should not be used interchangeably, because of poor intermodal agreement.

Accuracy of cardiac-induced brain motion measurement using displacement-encoding with stimulated echoes (DENSE) magnetic resonance imaging (MRI): A phantom study

PURPOSE

The goal of this study was to determine the accuracy of displacement-encoding with stimulated echoes (DENSE) MRI in a tissue motion phantom with displacements representative of those observed in human brain tissue.

METHODS

The phantom was comprised of a plastic shaft rotated at a constant speed. The rotational motion was converted to a vertical displacement through a camshaft. The phantom generated repeatable cyclical displacement waveforms with a peak displacement ranging from 92 µm to 1.04 mm at 1-Hz frequency. The surface displacement of the tissue was obtained using a laser Doppler vibrometer (LDV) before and after the DENSE MRI scans to check for repeatability. The accuracy of DENSE MRI displacement was assessed by comparing the laser Doppler vibrometer and DENSE MRI waveforms.

RESULTS

Laser Doppler vibrometer measurements of the tissue motion demonstrated excellent cycle-to-cycle repeatability with a maximum root mean square error of 9 µm between the ensemble-averaged displacement waveform and the individual waveforms over 180 cycles. The maximum difference between DENSE MRI and the laser Doppler vibrometer waveforms ranged from 15 to 50 µm. Additionally, the peak-to-peak difference between the 2 waveforms ranged from 1 to 18 µm.

CONCLUSION

Using a tissue phantom undergoing cyclical motion, we demonstrated the percent accuracy of DENSE MRI to measure displacement similar to that observed for in vivo cardiac-induced brain tissue.

Regional myocardial strain by cardiac magnetic resonance feature tracking for detection of scar in ischemic heart disease

BACKGROUND

Although cardiac magnetic resonance (CMR) can accurately quantify global left ventricular strain using feature tracking (FT), it has been suggested that FT cannot reliably quantify regional strain. We aimed to determine whether abnormalities in regional strain measured using FT can be detected within areas of myocardial scar and to determine the extent to which the regional strain measurement is impacted by LV ejection fraction (EF).

METHODS

We retrospectively studied 96 patients (46 with LVEF ≤ 40%, 50 with LVEF > 40%) with coronary artery disease and a late gadolinium enhancement (LGE) pattern consistent with myocardial infarction, who underwent CMR imaging (1.5T). Regional peak systolic longitudinal and circumferential strains (RLS, RCS) were measured within LGE and non-LGE areas. Linear regression analysis was performed for strain in both areas against LVEF to determine whether the relationship between strain and LGE holds across the LV function spectrum. Receiver-operating curve (ROC) analysis was performed in 33 patients (derivation cohort) to optimize strain cutoff, which was tested in the remaining 63 patients (validation cohort) for its ability to differentiate LGE from non-LGE areas.

RESULTS

Both RLS and RCS magnitudes were reduced in LGE areas: RLS = -10.4 ± 6.2% versus -21.0 ± 8.5% (p < 0.001); RCS = -10.4 ± 6.0% versus -18.9 ± 8.6%, respectively (p < 0.001), but there was considerable overlap between LGE and non-LGE areas. Linear regression revealed that it was partially driven by the natural dependence between strain and EF, suggesting that EF-corrected strain cutoff is needed to detect LGE. ROC analysis showed the ability of both RLS and RCS to differentiate LGE from non-LGE areas: area under curve 0.95 and 0.89, respectively. In the validation cohort, optimal cutoffs of RLS/EF = 0.36 and RCS/EF = 0.37 yielded sensitivity, specificity and accuracy 0.74-0.78.

CONCLUSION

Abnormalities in RLS and RCS within areas of myocardial scar can be detected using CMR-FT; however, LVEF must be accounted for.

Altered regional myocardial velocities by tissue phase mapping and feature tracking in pediatric patients with hypertrophic cardiomyopathy

BACKGROUND

Hypertrophic cardiomyopathy (HCM) is associated with heart failure, atrial fibrillation and sudden death. Reduced myocardial function has been reported in HCM despite normal left ventricular (LV) ejection fraction. Additionally, LV fibrosis is associated with elevated T1 and might be an outcome predictor.

OBJECTIVE

To systematically compare tissue phase mapping and feature tracking for assessing regional LV function in children and young adults with HCM and pediatric controls, and to evaluate structure-function relationships among myocardial velocities, LV wall thickness and myocardial T1.

MATERIALS AND METHODS

Seventeen pediatric patients with HCM and 21 age-matched controls underwent cardiac MRI including standard cine imaging, tissue phase mapping (two-dimensional cine phase contrast with three-directional velocity encoding), and modified Look-Locker inversion recovery to calculate native global LV T1. Maximum LV wall thickness was measured on cine images. LV radial, circumferential and long-axis myocardial velocity time courses, as well as global and segmental systolic and diastolic peak velocities, were quantified from tissue phase mapping and feature tracking.

RESULTS

Both tissue phase mapping and feature tracking detected significantly decreased global and segmental diastolic radial and long-axis peak velocities (by 12-51%, P<0.001-0.05) in pediatric patients with HCM vs. controls. Feature tracking peak velocities were lower than directly measured tissue phase mapping velocities (mean bias = 0.3-2.9 cm/s). Diastolic global peak velocities correlated moderately with global T1 (r = -0.57 to -0.72, P<0.01) and maximum wall thickness (r = -0.37 to -0.61, P<0.05).

CONCLUSION

Both tissue phase mapping and feature tracking detected myocardial velocity changes in children and young adults with HCM vs. controls. Associations between impaired diastolic LV velocities and elevated T1 indicate structure-function relationships in HCM.

Journal of Cardiovascular Magnetic Resonance: 2017/2018 in review

There were 89 articles published in the Journal of Cardiovascular Magnetic Resonance (JCMR) in 2017, including 76 original research papers, 4 reviews, 5 technical notes, 1 guideline, and 3 corrections. The volume was down slightly from 2017 with a corresponding 15% decrease in manuscript submissions from 405 to 346 and thus reflects a slight increase in the acceptance rate from 25 to 26%. The decrease in submissions for the year followed the initiation of the increased author processing charge (APC) for Society for Cardiovascular Magnetic Resonance (SCMR) members for manuscripts submitted after June 30, 2018. The quality of the submissions continues to be high. The 2018 JCMR Impact Factor (which is published in June 2019) was slightly lower at 5.1 (vs. 5.46 for 2017; as published in June 2018. The 2018 impact factor means that on average, each JCMR published in 2016 and 2017 was cited 5.1 times in 2018. Our 5 year impact factor was 5.82.In accordance with Open-Access publishing guidelines of BMC, the JCMR articles are published on-line in a continuus fashion in the chronologic order of acceptance, with no collating of the articles into sections or special thematic issues. For this reason, over the years, the Editors have felt that it is useful for the JCMR audience to annually summarize the publications into broad areas of interest or themes, so that readers can view areas of interest in a single article in relation to each other and contemporaneous JCMR publications. In this publication, the manuscripts are presented in broad themes and set in context with related literature and previously published JCMR papers to guide continuity of thought within the journal. In addition, as in the past two years, I have used this publication to also convey information regarding the editorial process and as a "State of our JCMR."This is the 12th year of JCMR as an open-access publication with BMC (formerly known as Biomed Central). The timing of the JCMR transition to the open access platform was "ahead of the curve" and a tribute to the vision of Dr. Matthias Friedrich, the SCMR Publications Committee Chair and Dr. Dudley Pennell, the JCMR editor-in-chief at the time. The open-access system has dramatically increased the reading and citation of JCMR publications and I hope that you, our authors, will continue to send your very best, high quality manuscripts to JCMR for consideration. It takes a village to run a journal and I thank our very dedicated Associate Editors, Guest Editors, Reviewers for their efforts to ensure that the review process occurs in a timely and responsible manner. These efforts have allowed the JCMR to continue as the premier journal of our field. This entire process would also not be possible without the dedication and efforts of our managing editor, Diana Gethers. Finally, I thank you for entrusting me with the editorship of the JCMR as I begin my 4th year as your editor-in-chief. It has been a tremendous experience for me and the opportunity to review manuscripts that reflect the best in our field remains a great joy and highlight of my week!

CMR DENSE and the Seattle Heart Failure Model Inform Survival and Arrhythmia Risk After CRT

OBJECTIVES

This study sought to determine if combining the Seattle Heart Failure Model (SHFM-D) and cardiac magnetic resonance (CMR) provides complementary prognostic data for patients with cardiac resynchronization therapy (CRT) defibrillators.

BACKGROUND

The SHFM-D is among the most widely used risk stratification models for overall survival in patients with heart failure and implantable cardioverter-defibrillators (ICDs), and CMR provides highly detailed information regarding cardiac structure and function.

METHODS

CMR Displacement Encoding with Stimulated Echoes (DENSE) strain imaging was used to generate the circumferential uniformity ratio estimate with singular value decomposition (CURE-SVD) circumferential strain dyssynchrony parameter, and the SHFM-D was determined from clinical parameters. Multivariable Cox proportional hazards regression was used to determine adjusted hazard ratios and time-dependent areas under the curve for the primary endpoint of death, heart transplantation, left ventricular assist device, or appropriate ICD therapies.

RESULTS

The cohort consisted of 100 patients (65.5 [interquartile range 57.7 to 72.7] years; 29% female), of whom 47% had the primary clinical endpoint and 18% had appropriate ICD therapies during a median follow-up of 5.3 years. CURE-SVD and the SHFM-D were independently associated with the primary endpoint (SHFM-D: hazard ratio: 1.47/SD; 95% confidence interval: 1.06 to 2.03; p = 0.02) (CURE-SVD: hazard ratio: 1.54/SD; 95% confidence interval: 1.12 to 2.11; p = 0.009). Furthermore, a favorable prognostic group (Group A, with CURE-SVD <0.60 and SHFM-D <0.70) comprising approximately one-third of the patients had a very low rate of appropriate ICD therapies (1.5% per year) and a greater (90%) 4-year survival compared with Group B (CURE-SVD ≥0.60 or SHFM-D ≥0.70) patients (p = 0.02). CURE-SVD with DENSE had a stronger correlation with CRT response (r = -0.57; p < 0.0001) than CURE-SVD with feature tracking (r = -0.28; p = 0.004).

CONCLUSIONS

A combined approach to risk stratification using CMR DENSE strain imaging and a widely used clinical risk model, the SHFM-D, proved to be effective in this cohort of patients referred for CRT defibrillators. The combined use of CMR and clinical risk models represents a promising and novel paradigm to inform prognosis and device selection in the future.

Echocardiography and cardiovascular magnetic resonance based evaluation of myocardial strain and relationship with late gadolinium enhancement

OBJECTIVES

We sought to: (1) determine the agreement in cardiovascular magnetic resonance (CMR) and speckle tracking echocardiography (STE) derived strain measurements, (2) compare their reproducibility, (3) determine which approach is best related to CMR late gadolinium enhancement (LGE).

BACKGROUND

While STE-derived strain is routinely used to assess left ventricular (LV) function, CMR strain measurements are not yet standardized. Strain can be measured using dedicated pulse sequences (strain-encoding, SENC), or post-processing of cine images (feature tracking, FT). It is unclear whether these measurements are interchangeable, and whether strain can be used as an alternative to LGE.

METHODS

Fifty patients underwent 2D echocardiography and 1.5 T CMR. Global longitudinal strain (GLS) was measured by STE (Epsilon), FT (NeoSoft) and SENC (Myocardial Solutions) and circumferential strain (GCS) by FT and SENC.

RESULTS

GLS showed good inter-modality agreement (r-values: 0.71-0.75), small biases (< 1%) but considerable limits of agreement (- 7 to 8%). The agreement between the CMR techniques was better for GLS than GCS (r = 0.81 vs 0.67; smaller bias). Repeated measurements showed low intra- and inter-observer variability for both GLS and GCS (intraclass correlations 0.86-0.99; coefficients of variation 3-13%). LGE was present in 22 (44%) of patients. Both SENC- and FT-derived GLS and GCS were associated with LGE, while STE-GLS was not. Irrespective of CMR technique, this association was stronger for GCS (AUC 0.77-0.78) than GLS (AUC 0.67-0.72) and STE-GLS (AUC = 0.58).

CONCLUSION

There is good inter-technique agreement in strain measurements, which were highly reproducible, irrespective of modality or analysis technique. GCS may better reflect the presence of underlying LGE than GLS.

Cardiac magnetic resonance feature tracking for quantifying right ventricular deformation in type 2 diabetes mellitus patients

To determine the feasibility of deformation analysis in the right ventricle (RV) using cardiovascular magnetic resonance myocardial feature tracking (CMR-FT) in type 2 diabetes mellitus (T2DM) patients. We enrolled 104 T2DM patients, including 14 with impaired right ventricular ejection fraction (RVEF) and 90 with preserved RVEF, and 26 healthy controls in this prospective study. CMR was used to determine RV feature-tracking parameters. RV strain parameters were compared among the controls, patients with preserved and reduced RVEF. Binary logistic regression was used to predict RV dysfunction. Receiver operating characteristic analysis was used to assess the diagnostic accuracy. The agreement was tested by Bland-Altman analysis. Compared with controls, longitudinal and circumferential global peak strain (PS) and PS at mid-ventricular, apical slices were significantly decreased in T2DM patients with or without reduced RVEF (p < 0.05). Within the T2DM patients, the global longitudinal PS (GLPS) and the longitudinal PS at mid-ventricular segments were significantly reduced in the reduced RVEF group than in preserved RVEF groups (p < 0.05). GLPS was an independent predictor of RV dysfunction (odds ratio: 1.246, 95% CI: 1.037-1.496; p = 0.019). The GLPS demonstrated greater diagnostic accuracy (area under curve: 0.716) to predict RV dysfunction. On Bland-Altman analysis, global circumferential PS and GLPS had the best intra- and inter-observer agreement, respectively. In T2DM patients, CMR-FT could quantify RV deformation and identify subclinical RV dysfunction in those with normal RVEF. Further, RV strain parameters are potential predictors for RV dysfunction in T2DM patients.

Reproducibility and Agreement of Tissue Tracking versus Feature Tracking for Strain Measurement on Cardiac MR Images in Patients with Repaired Tetralogy of Fallot

PURPOSE

To evaluate differences in strain measurements by using tissue-tracking (TT) and feature-tracking (FT) software on cardiovascular MR images in patients with repaired tetralogy of Fallot (TOF).

MATERIALS AND METHODS

In this retrospective cross-sectional study of 25 patients with repaired TOF (median age, 33.1 years; interquartile range, 25.3-38.3 years) from 2008 through 2014, left ventricular (LV) and right ventricular (RV) global circumferential and longitudinal strain were measured from cardiac MR images by using TT and FT software. Time to process was measured from opening the study to acceptance of contours. Intra- and interobserver reproducibility were evaluated with Bland-Altman analysis, coefficient of variation, and intraclass correlation coefficient.

RESULTS

Time to process was slightly longer for TT (10.2 minutes ± 3.1 [standard deviation] vs 9.0 minutes ± 1.7, P = .04). Fewer patients required contour revision with TT than with FT. Both TT and FT measurements had similar moderate-to-strong correlations with LV and RV ejection fractions; correlation of RV longitudinal strain with RV ejection fraction did not reach significance by using either method. With the exception of LV circumferential strain, strain measurements were lower with FT relative to TT. Intra- and interobserver reproducibility were lower with FT for longitudinal strain measurements.

CONCLUSION

TT and FT have systematic differences in strain values and reproducibility, particularly for longitudinal strain. Software-specific normative data are necessary, as are studies to evaluate correlation with clinical outcomes for each modality.© RSNA, 2019.