Quantification of late gadolinium enhanced CMR in viability assessment in chronic ischemic heart disease: a comparison to functional outcome
Journal of Cardiovascular Magnetic Resonance volume 11, Article number: 6 (2009)
Quantification of late gadolinium enhanced cardiovascular magnetic resonance (LGE CMR) by objective window setting increases reproducibility and facilitates multicenter comparison and cooperation. So far, quantification methods or models have only been validated to postmortem animal studies. This study was undertaken to evaluate quantification of LGE in relation to the clinical standard of viability, i.e. functional outcome after revascularization.
Thirty-eight patients with chronic ischemic myocardial dysfunction underwent cine and LGE 1 month before and cine CMR 6 months after coronary revascularization. Enhancement was quantified by thresholding window setting at: 2-8SD above mean signal intensity of a remote normal region, and according to the full width at half maximum method (FWHM). Dysfunctional segments were divided in 5 groups according to segmental extent of enhancement (SEE): SEE 1 – no enhancement to SEE 5 – 76–100% with each quantification method.
Quantification methods had a strong influence on SEE and total infarct size. Multilevel analysis showed that thresholding contrast images at 6SD best predicted segmental functional outcome after revascularization, but the difference with other methods was small and non-significant.
Simple thresholding techniques strongly influence global and segmental extent of LGE, but have relatively little influence on the accuracy to predict segmental functional improvement after revascularization.
Revascularization of dysfunctional but viable myocardium may lead to reversed remodelling, improved regional and global function and better prognosis in patients with chronic ischemic heart disease . The diagnostic accuracy of imaging modalities to predict functional outcome is influenced by the definition of disease (what is viable). Although visual or qualitative analysis may provide satisfactory results, standardization and quantification of these definitions increases reproducibility and reliability in follow-up studies, and facilitates comparison between different centers.
Late gadolinium enhanced cardiovascular magnetic resonance (LGE CMR) accurately visualizes the transmural extent of ischemia-related scar and has been shown to predict the likelihood of functional improvement after revascularization [2, 3]. Several methods have been proposed to differentiate enhanced, non-viable from non-enhancing, viable myocardium, all using the in-slice signal intensity of infarcted or remote myocardium, and ranging from simple thresholding to more complex computer algorithms. We have previously shown that the use of common thresholds based on the suppressed signal of remote myocardium may lead to considerable overestimation of the infarct size . However, in this study, we used visual estimation as the reference standard. Although a number of experimental studies have used ex-vivo imaging or 2,3,5-triphenyltetrazolium choride (TTC) staining to determine the optimal threshold of enhancement, so far, no study has used the clinically useful standard of viability i.e. functional outcome after revascularization [5–7].
Therefore, the aim of this study was to evaluate the relation between quantification of LGE and functional outcome after revascularization in patients with chronic ischemic myocardial dysfunction. To quantify LGE, we chose simple thresholding techniques that are easily applicable in any clinical or research situation.
All patients with known coronary artery disease and regional wall motion abnormalities on echocardiography or left ventricular (LV) angiography, without CMR contraindications, who were scheduled to undergo surgical or percutaneous revascularisation, were study candidates. The Committee on Research Involving Human Subjects of the VU University Medical Centre, Amsterdam, approved the study protocol. All patients gave written informed consent.
Forty-seven patients were initially included in this study protocol. After revascularization, 7 patients were excluded because of left ventricular aneurysmectomy (1), electrocardiographic and/or biochemical evidence of peri-procedural myocardial infarction (defined as post procedural peak CK-MB > 3 upper limit of normal) (4), pacemaker implantation (1), and incomplete data (1). During analysis, 2 more patients were excluded because of absence of wall motion abnormalities at baseline and non-diagnostic image quality, leaving 38 patients as the final study group. All patients were in stable clinical condition at the time of both CMR examinations without clinical evidence of ischemic events during the study period.
CMR scans were acquired at 4 ± 4 weeks before and 30 ± 4 weeks after revascularisation. All scans were performed on a 1.5T scanner (Sonata, Siemens, Erlangen, Germany) with the patient in a supine position using a four-element phased array cardiac receiver coil. ECG-gated cine images were acquired using a breath-hold segmented steady-state free precession sequence (true FISP; echo time/repetition time of 1.2/3.2 ms; resolution of 1.3 × 1.8 × 5 mm). Per patient eight to ten short-axis views were obtained every 10 mm starting from the mitral valve insertion and covering the entire left ventricle. Ten to 15 minutes after injection of a gadolinium-based contrast agent (Magnevist, Schering AG, Berlin, Germany; 0,2 mmol/kg) contrast-enhanced images were acquired in the same orientation as the cine images using a 2D-segmented inversion recovery gradient-echo pulse sequence triggered to end-diastole (repetition time/echo time = 9.6/4.4 ms, flip angle 25°, inversion time set to suppress signal from remote myocardium, matrix 208 × 256 and a typical voxel size of 1.6 × 1.3 × 5.0).
All data were analysed on a separate workstation using a dedicated software package (MASS v15, 2008, Medis, Leiden, The Netherlands).
Segmental wall thickness was measured at end-systole and end-diastole after manual tracing of endocardial end epicardial borders in stop-frame images, carefully excluding trabeculations and papillary muscles. The observer (OB) was blinded to other patient or imaging data such as extent of coronary artery disease, use of medication and results of the LGE analysis. Baseline and follow-up slices were analysed separately after registration using the scanner slice position and various anatomical landmarks such as right ventricle septal insertion sites, papillary muscle location, and trabecularization patterns in the right and left ventricles. Although the observer was unaware of the timing of the study, blinding to pre-/postoperative status was impossible because of the artifacts related to sternal wires. For segmental analysis, the 2 most basal and apical slices were excluded because of the left ventricular outflow tract and partial volume effects, respectively. The remaining slices were divided into 6 segments each, starting at the inferior insertion of the right ventricle to the septum.
Segmental wall thickening (SWT) in millimetres was calculated as: end systolic wall thickness minus end diastolic wall thickness. Segments with SWT < 3 mm (mean – 2SD) were considered dysfunctional . Functional improvement was defined as an increase in SWT of ≥ 1.5 mm compared to baseline. Intraobserver and interobserver variability of SWT were 0.0 ± 0.4 mm (mean difference between 2 measurements (OB); intraclass correlation coefficient = 0.97, 95% confidence interval 0.86 – 0.99), and 0.1 ± 0.7 mm (mean difference between values of observer 1 (OB) and 2 (AMB); intraclass correlation coefficient = 0.89, 95% confidence interval 0.55–0.98), as reported previously .
Left ventricular end-diastolic volume (EDV) and end-systolic volumes (ESV) were determined by planimetry of all short-axis images in each patient and indexed to body surface area. Left ventricular ejection fraction (EF, %) was calculated as (LVEDV - LVESV)/LVEDV * 100%.
Endocardial and epicardial contours were manually traced, again avoiding papillary muscles and trabecularizations. Enhanced regions were then determined in the following ways:
I. After thresholding signal intensity at 2–8 standard deviations (SD) above the mean signal intensity of remote normal myocardium in the same slice. A region of interest of 0.5–1 cm2 was manually drawn and placed in remote normal myocardium (defined as normal function without enhancement on visual assessment). If a slice contained no remote normal myocardium, mean and SD of the nearest slice with normal remote myocardium was used; if the 2 neighbouring slices contained normal myocardium, mean and SD were averaged.
II. After thresholding signal intensity using the full-width at half-maximum (FWHM) method that defines the enhanced area by using 50% of the maximum signal found within the enhanced area. The maximum signal was found by computer-assisted window thresholding of the enhanced area. If no enhancement was found in a slice, the maximum signal of the nearest slice with enhancement was used; if the 2 neighbouring slices showed enhancement, maximum signals were averaged.
All areas of enhancement were quantified by computer-assisted planimetry on each of the short axis images and segmental extent of enhancement (SEE) was expressed as percentage of segmental area. Obvious artifacts such as caused by motion were excluded by highlighting them using a tool from the software package. Other small isolated regions of enhancement that were clearly not of ischemic origin (like small subepicardial spots) were also excluded from analysis. Total infarct size was calculated by summation of all slice volumes of enhancement, i.e. including the slices that were excluded for segmental analysis.
The paired sample t test (with Bonferroni correction) and the independent samples t test were used to compare means within the study group or between subgroups. To account for the non-independence of the data, we used multilevel logistic regression (MlwiN, version 1.02.0007, Centre for Multilevel Modelling, London, United Kingdom) to analyse the relation between SEE and the likelihood of improvement. Details of the multilevel analysis have been published previously . All dysfunctional segments were assigned to one of the following groups according to SEE: 1 – 0%, 2 – 1 to 25%, 3 – 26 to 50%, 4 – 51 to 75%, and 5 – 76 to 100% enhancement. The regression coefficients were used to calculate odds ratios that express the likelihood of improvement of SEE-groups 2–5 relative to SEE-group 1. To compare the different quantification methods, multilevel analysis was repeated assuming a linear correlation between SEE and likelihood of improvement, which generates one regression coefficient per method.
Receiver operating characteristics (ROC) analysis was used to find the SEE cut-off with highest diagnostic accuracy to predict segmental functional improvement.
Improvement in ejection fraction was defined as an increase from baseline to follow-up of ≥ 5%. Since the number of viable segments per patient is not normally distributed, the Mann-Whitney test was used to compare these between patients with and without improved ejection fraction.
All values are expressed as mean ± SD. P-values < 0.05 were considered statistically significant.
Patient characteristics are provided in table 1.
Extent of enhancement
The extent of enhancement according to quantification method is shown in table 2. Both total infarct size and segmental extent of enhancement were strongly correlated between all methods (p < 0.001), except total infarct size between FWHM and 2SD and 3SD (non-significant after Bonferroni correction). Total infarct size decreased with increasing number of SD's (all steps p < 0.001). Total infarct size according to FWHM differed significantly from 2–4 SD (p < 0.001). Mean SEE decreased with increasing number of SD's and FWHM (all p < 0.001).
A total of 1122 segments were analysed (29.6 ± 5.5/patient). Of these, 628 were dysfunctional (mean SWT 1.1 ± 1.1 mm). Segmental function improved in 174 segments (mean SWT 0.9 ± 1.2 to 3.5 ± 1.5 mm). The likelihood of improvement was inversely related to the SEE according to all quantification methods. The results of the multilevel analysis are displayed in table 3. The odds ratios were inverted so that they reflect the number of times that a certain SEE group is less likely to improve compared to SEE group 1 in that category. Quantifying enhancement using 6SD best predicted segmental functional improvement at follow-up: segments in SEE-groups 2, 3, 4 and 5 were 2.7, 4.7, 5.2 and 14.8 less likely to improve than segments in SEE-group 1 (table 3).
Multilevel analysis was repeated assuming a linear relationship between SEE-groups and likelihood of improvement. This analysis generates one odds ratio that represents the likelihood of improvement from one SEE-group to the next. 6SD had the highest predictive power: SEE-group x was 2.04 (95% confidence intervals 1.62–2.58) less likely to improve than the SEE-group x-1. 2SD had the lowest power with (inverted) odds ratio 1.66 (1.31–2.08). Figure 1 shows 95%-confidence intervals for each method, demonstrating that there is considerable overlap.
Left ventricular EDV and ESV were 123 ± 36 ml/m2 and 79 ± 35 ml/m2 at baseline, and decreased significantly to 108 ± 36 ml/m2 (p < 0.001) and 68 ± 35 ml/m2 (p = 0.002), respectively. Mean EF improved from 38.2 ± 11.4% at baseline to 40.4 ± 12% at follow-up (p = 0.067).
ROC analysis showed that the 6SD method optimally predicted segmental improvement at < 10% SEE (AUC 0.700 (0.654–0.746), sensitivity 70%, specificity 65%) (table 4). Twenty-two patients showed significant (≥ 5%) improvement in EF. Although the number of dysfunctional segments with < 10% SEE-6SD was higher in patients with improved ejection fraction compared to patients without improvement (median 7.5, range 1–24, vs 4.0, range 0–14), there was considerable overlap, and no significant difference could be found (Mann-Whitney-test, p = 0.251).
Our results show that quantification methods have a strong influence on total and regional infarct extent. Although analysis using 6SD above signal intensity of a remote normal region resulted in the best prediction of segmental functional improvement after revascularization, differences with the other methods were small and non-significant.
Both SD- and FWHM-methods are based on experience rather than on firm pathophysiological evidence. SD uses the (purposely suppressed) signal of a remote normal region, whereas FWHM uses the peak signal within the enhanced region. FWHM has been suggested on theoretical grounds to be less sensitive to partial volume effects. Amado et al compared total infarct size according to the SD- and FWHM-thresholding methods and postmortem data in an animal occlusion-reperfusion study . They found that FWHM correlated best with TTC-staining with a small overestimation (significance not reported), whereas no significant correlation and considerable overestimation was found with 6SD. Heiberg et al recently proposed a quantification method to compensate for partial volume effects in thicker (8–10) mm) slices by weighting pixels according to their signal intensity . The standard LGE protocol in our institution uses 5 mm slice thickness, which maintains excellent image quality while minimizing partial volume effects.
Hsu et al developed a computer algorithm based on feature analysis and combined thresholding involving sequential use of the 2SD- and FWHM-method . The algorithm was validated in a canine study, and compared to visual assessment and 2SD- and FWHM-thresholding using ex-vivo and in-vivo LGE. The authors found that infarct size measured by their model showed no difference with TTC-staining, and was more accurate than the other methods, that all overestimated the histopathological infarct size. The algorithm was subsequently tested in 20 patients with acute or chronic infarction, and resulted in significantly lower values for both regional and global infarct size when compared to visual assessment, although this study lacked a reference standard . Hsu et al used phase-sensitive reconstruction in both studies and also used 8 mm slice thickness, which makes comparison of their findings to our study difficult and again stresses the importance of standardization.
In a previous study we quantified enhancement in 15 patients with chronic ischemic heart disease. In comparison to visual assessment, we found significant overestimation with 2–4SD of both segmental extent and total infarct size . The use of 6SD led to a significant underestimation, whereas 5SD showed no difference.
In the current study we used functional outcome after revascularisation as the clinical reference standard of viability. We found that, although the extent of enhancement was profoundly influenced by the quantification method used, there was no significant difference in their accuracy to predict segmental functional improvement at follow-up. 6SD showed the highest accuracy, followed closely by FWHM, 4SD and 5SD. Thresholding at 2SD and 3SD above remote signal intensity was less accurate, but there was still considerable overlap of confidence intervals with 6SD. Our data suggest that simple thresholding methods using the signal intensity of remote or of the enhanced area have relatively little effect on the assessment of regional viability and the prediction of likelihood of improvement. The lack of a clear difference may seem surprising considering the large differences that the use of the various methods caused in total infarct size and segmental extent of enhancement. However, with each method the inverse relation between SEE and likelihood of improvement remains intact, which makes it hard to detect small differences in predictive power. A larger study group and longer follow-up may be required to strengthen the statistical position of 6SD .
All current LGE quantification methods are limited by the fact that essential steps like the delineation of myocardial contours, the drawing a region of interest in remote myocardium and the exclusion of artifacts are all (still) done manually.
The applicability of our results is limited to patients with chronic ischemic heart disease, although they are probably valid in acute myocardial infarction, since both animal and clinical studies have shown that enhancement reflects necrosis at all stages after infarction [10–12]. However, this remains to be established. Also, our results cannot be extrapolated to other techniques (e.g. 3D-acquisition) or scanner types (especially with different magnetic field strengths).
In conclusion, in this study we evaluated the relation between simple thresholding techniques to quantify LGE images and functional outcome after revascularization in patients with chronic ischemic myocardial dysfunction. Although quantitative analysis with window setting thresholded at mean + 6SD of the signal of a remote normal region best predicted segmental functional improvement, there was no significant difference with the other SD or FWHM methods. Further study is needed to evaluate SD and FWHM in larger groups with a longer follow-up and compare them to the more complex algorithms that have not yet been tested in human viability studies.
Camici PG, Prasad SK, Rimoldi OE: Stunning, Hibernation, and Assessment of Myocardial Viability. Circulation. 2008, 117: 103-114. 10.1161/CIRCULATIONAHA.107.702993.
Kim RJ, Wu E, Rafael A, Chen EL, Parker MA, Simonetti O, Klocke FJ, Bonow RO, Judd RM: The Use of Contrast-Enhanced Magnetic Resonance Imaging to Identify Reversible Myocardial Dysfunction. N Engl J Med. 2000, 343: 1445-1453. 10.1056/NEJM200011163432003.
Bondarenko O, Beek AM, Twisk JW, Visser CA, van Rossum AC: Time course of functional recovery after revascularization of hibernating myocardium: a contrast-enhanced cardiovascular magnetic resonance study. Eur Heart J. 2008, 29: 2000-5. 10.1093/eurheartj/ehn266.
Bondarenko O, Beek AM, Hofman MB, Kühl HP, Twisk JW, van Dockum WG, Visser CA, van Rossum AC: Standardizing the Definition of Hyperenhancement in the Quantitative Assessment of Infarct Size and Myocardial Viability Using Delayed Contrast-Enhanced CMR. J Cardiovasc Magn Reson. 2005, 7: 481-485. 10.1081/JCMR-200053623.
Amado LC, Gerber BL, Gupta SN, Rettmann DW, Szarf G, Schock R, Nasir K, Kraitchman DL, Lima JA: Accurate and objective infarct sizing by contrast-enhanced magnetic resonance imaging in a canine myocardial infarction model. J Am Coll Cardiol. 2004, 44: 2383-2389. 10.1016/j.jacc.2004.09.020.
Hsu LY, Natanzon A, Kellman P, Hirsch GA, Aletras AH, Arai AE: Quantitative myocardial infarction on delayed enhancement MRI. Part I: Animal validation of an automated feature analysis and combined thresholding infarct sizing algorithm. J Magn Reson Imaging. 2006, 23: 298-308. 10.1002/jmri.20496.
Heiberg E, Ugander M, Engblom H, Götberg M, Olivecrona GK, Erlinge D, Arheden H: Automated Quantification of Myocardial Infarction from MR Images by Accounting for Partial Volume Effects: Animal, Phantom, and Human Study. Radiology. 2008, 246: 581-588. 10.1148/radiol.2461062164.
Simonetti OP, Kim RJ, Fieno DS, Hillenbrand HB, Wu E, Bundy JM, Finn JP, Judd RM: An improved MR imaging technique for the visualization of myocardial infarction. Radiology. 2001, 218: 215-223.
Hsu LY, Ingkanisorn WP, Kellman P, Aletras AH, Arai AE: Quantitative myocardial infarction on delayed enhancement MRI. Part II: Clinical application of an automated feature analysis and combined thresholding infarct sizing algorithm. J Magn Reson Imaging. 2006, 23: 309-314. 10.1002/jmri.20495.
Beek AM, Kühl HP, Bondarenko O, Twisk JW, Hofman MB, van Dockum WG, Visser CA, van Rossum AC: Delayed contrast-enhanced magnetic resonance imaging for the prediction of regional functional improvement after acute myocardial infarction. J Am Coll Cardiol. 2003, 42: 895-901. 10.1016/S0735-1097(03)00835-0.
Fieno DS, Kim RJ, Chen EL, Lomasney JW, Klocke FJ, Judd RM: Contrast-enhanced magnetic resonance imaging of myocardium at risk: distinction between reversible and irreversible injury throughout infarct healing. J Am Coll Cardiol. 2000, 36: 1985-1991. 10.1016/S0735-1097(00)00958-X.
Selvanayagam JB, Kardos A, Francis JM, Wiesmann F, Petersen SE, Taggart DP, Neubauer S: Value of Delayed-Enhancement Cardiovascular Magnetic Resonance Imaging in Predicting Myocardial Viability After Surgical Revascularization. Circulation. 2004, 110: 1535-1541. 10.1161/01.CIR.0000142045.22628.74.
This study was supported by Netherlands Heart Foundation grant 2001.158.
Prof. Jos W.R. Twisk, PhD, Department of Clinical Epidemiology and Biostatistics is kindly acknowledged for his contribution to the statistical analysis.
Mary Belderok is kindly acknowledged for her contribution to the organization of the study and the data collection.
The authors declare that they have no competing interests.
AMB participated in the design of the study, data acquisition, image analysis, statistical analysis and drafting of the manuscript. OB and FA participated in the data acquisition and the image analysis. ACvR participated in the design of the study and drafting of the manuscript. All authors have read and approved the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Beek, A.M., Bondarenko, O., Afsharzada, F. et al. Quantification of late gadolinium enhanced CMR in viability assessment in chronic ischemic heart disease: a comparison to functional outcome. J Cardiovasc Magn Reson 11, 6 (2009). https://doi.org/10.1186/1532-429X-11-6
- Cardiovascular Magnetic Resonance
- Infarct Size
- Standard Deviation
- Quantification Method
- Remote Myocardium