Abstract
Purpose
To assess the effect of applying an automated heart model based measurements of left ventricle (LV) and compare with manual and semi-automated measurements at Cardiovascular MR Imaging.
Materials and Methods
Sixty-two patients who underwent cardiac 1.5T MR imaging were included. Steady state free precession cine images of 20 phases per cardiac cycle were obtained in short axis views and both 2-chamber and 4-chamber views. Epicardial and endocardial contours were drawn in manual, automated, and semi-automated ways. Based on these acquired contour sets, the end-diastolic (ED) and end-systolic (ES) volumes, ejection fraction (EF), systolic volume (SV) and LV mass were calculated and compared.
Results
In EDV and ESV, the differences among three measurement methods were not statistically significant (P = .399 and .145, respectively). However, in EF, SV, and LV mass, the differences were statistically significant (P=.001, <001, <001, respectively) and the measured value from automated method tend to be consistently higher than the values from other two methods.
The accurate measurement of left ventricular (LV) volume and LV mass is important for risk stratification and clinical management, because these parameters are powerful predictors of prognosis related with a variety of cardiovascular disease (1-3).
Cardiac magnetic resonance (CMR) imaging performed using the cine gradient-echo and steady state free precession (SSFP) is the reference standard in the assessment of the LV volume and mass (4-6). With the use of these techniques, an acquisition of the delineations of endocardial border throughout the entire cardiac cycle was facilitated (7). Thus, the measurement by the summation-of-disk section method has been the established reference standard for assessing LV volume and mass. However, short axis (SAX) stacks, considered to be the standard procedure for quantifying LV volume and function, and as such for validation (5), are inevitably leading to problems that underestimated the longitudinal component of cardiac contraction. On the other hand, long axis (LAX) stacks, which have been proposed as an alternatives (8), may be more reliable and reproducible in normal ventricles (9), but more susceptible to the errors in asymmetric diseased ventricles.
Under the circumstances that there is no accepted gold standard for exact LV volume and mass measurements, one of the plausible ways could be the combination of two approaches, i.e., SAX and LAX stack. The automated tool (Argus 4D VF®, Siemens Healthcare, Erlangen, Germany), which is based on a heart-model-based algorithm, tried to combine both methods by employing end-diastole and end-systole long-axis planes into its SAX based LV volume and mass calculation.
Therefore, the aim of this study is to assess the effect of applying an automated heart model based measurements of LV, which contains more components of LAX, and compare the result with conventional manual 2D SAX stack based measurements at Cardiovascular MR Imaging.
This retrospective study was approved by the institutional review board of our hospital, and the requirement for informed consent was waived. 62 patients who underwent cardiac 1.5T MR imaging between March 2011 and July 2011 were included in this study. 51 patients were male and 11 female, with a mean age of 49 years (range 7-77 years, standard deviation 19.9 years). 21 of them had ischemic heart disease, 7 had valvular heart diseases, 14 had cardiomyopathy, 13 had congenital heart diseases, and the remaining 7 had no significant structural heart disease.
All CMR investigations were performed in a 1.5T unit (Sonata, Siemens Healthcare, Erlangen, Germany), using a standardized clinical protocol. Cine images were obtained in 2-chamber, 4-chamber, RV outflow tract and short-axis planes with a temporal resolution that was sufficient to accommodate 20 true phases per cardiac cycle.
Left ventricular parameters were measured from a stack of short-axis and long-axis cine images with commercially available software (Argus 4D VF®, Siemens Healthcare, Erlangen, Germany). They were analyzed to calculate LV volume and mass using a manual contour tracing technique, automated heart-model-based technique, and semi-automated technique.
First, in the manual measurements, the source cine images of SAX stack were manually segmented and contoured in every phase of the cardiac cycle (10), while, in the automated measurements, the basal border of the LV cavity was automatically adapted in every phase of the cardiac cycle, after mitral valve insertion points on long-axis planes were identified in end-diastole and end-systole phase. LV parameters were calculated based on these 2 LAX views and borders generated on SAX stack. Finally, in the semi-automated measurements, the aforementioned sets of contours, which were generated in the automated measurements, were adjusted manually to represent proper blood pool of LV in SAX stack in every phase of a cardiac cycle.
The end diastolic volume (EDV), end systolic volume (ESV), ejection fraction (EF), systolic volume (SV) and end-diastolic LV mass among three measurement methods were compared.
One-way analysis of variance (ANOVA) test was used for the comparison of the measurements data and Kruskal-Wallis test was used for posthoc analysis among three methods. P < .05 was considered to be a significant difference. Statistical analysis was performed using commercial software (SPSS for Windows, version 18.0, SPSS, Chicago, IL).
All loaded images of SAX stack and 2 LAX were sufficient for tracing of the endocardial and epicardial contours with three different methods (example see Fig. 1). Left ventricular parameters were calculated for each contour set.
In EDV and ESV, the differences among three measurement methods were not statistically significant. (P = .399 and .145, respectively). On the other hand, as for the EF, SV, and LV mass measurements, the differences were statically significant (P=.001, <.001, <.001, respectively) and the measured value from semi-automated method tend to be consistently larger than those from the other methods (Table 1).
In EF, the mean value measured from the automated method (62.87±1.41) turned out to be larger value than one from the manual method (55.00±1.61). After the generated contour sets in automated methods were manually adjusted in each phase, the mean EF value was decreased. Therefore, the value from semi-automated measurements (55.94±1.81), was lower than that from the automated method. As a result, the mean measured value from the automated method was significantly higher than those from other two methods (P = .001), while, the values from manual and semi-automated methods were not significantly different in post hoc multiple comparisons. Similarly, in SV, the mean value from the automated one (103.73±3.41) was higher than one from the manual method (85.17±2.76). The value from semi-automated method was lowered (92±3.13) after the adjustment, and the mean measured value from the automated method was significantly higher than those from other two methods (P < .001), while the difference between manual and semi-automated one was not statistically significant. Finally, in LV mass, the mean value from the automated one (151.73±6.18) was higher than one from the manual method (129.06±5.17). The value from semi-automated one was decreased again (118.04±5.20), therefore, the mean measured value from the automated method was significantly larger than those from other two methods (P <.001). Unlike to EF and SV, the difference between manual and semi-automated was statistically significant in post hoc multiple comparison, and the mean value from semi-automated method was lower than one from the manual measurement. Graphical presentations of the result were illustrated in Fig 2.
Because the most validation studies in CMR have been performed based on SAX measurements (11, 12), the correct selection of the basal slice is of importance for accuracy of volume data, which were often became the source of the problem leading to low reproducibility (13). Furthermore, the measurement sorely based on the SAX stacks without the incorporation of LAX components, could be inevitably underestimate the longitudinal part of heart contraction. Argus 4D tried to visit these issues by incorporating end-diastole and end-systole LAX view and by enabling the user to identify the mitral valve insertion point. The primary finding of our study results is two-fold: in comparison with conventional 2-D disc summation method, automated heart-model-based method grossly overestimated EF, SV and LV mass with statistical significance. The other point is that, in semi-automated measurement, which is based on the manual adjusted contour sets of automatically generated ones, the differences with conventional 2-D method were decreased so that the statistical significance went away although the overall tendencies of exaggeration were preserved (Fig. 2).
A first and straightforward explanation lies in the discrepancy to define mitral valve plane and subsequent inclusion of additional ventricular or myocardial volume in the automated method. As mentioned earlier, the choice of the most basal slice in the 2D disc-summation technique is somewhat arbitrary, which is defined as the closest slice to mitral valve, surrounded by at least 50% of myocardium. In the automated method, the most basal contour can be oriented along the anatomic mitral valve plane with long and short axis views, which reduces the risk of missing ventricular volume deteriorating the accuracy of LV volumetric analysis. Considering the most basal part of the heart takes the highest relative to ventricular tissue and cavity, inclusion or exclusion of basal volume constitutes substantial variation of overall volumetric measurements. This might explain why overestimation was so consistent in the automated method (with statistical significance) and in the semi-automated method (without statistical significance), compared with manual 2D method.
The second explanation is the different approaches to papillary muscles in automated contour drawing, compared with conventional 2D or manually adjusted semi-automated one. The SAX disc-summation method includes papillary muscles in the volume of the LV cavity. While, in the manual contouring, the endocardial contours were consistently drawn to include papillary muscle in the volume of LV cavity on every single cardiac phase, the automated contouring has the tendency to draw endocardial border to exclude some portion of papillary muscle in the LV cavity. Since the papillary muscle account for around 8.9% of the LV volume, this difference contouring should result in variation of overall LV volumetric measurements (14). This tendency is more evident in systolic phase which the papillary muscle appears more compact pattern relative to cardiac wall (Fig. 1). Indeed, the volumetric change after the manual correction is much larger in systolic phase (79.22 -> 65.49, 17.3 %) compared with the change in diastolic phase (169 -> 171, 1.1%). With the results that the differences with conventional method became much smaller after manual adjustment of automatically drawn contour sets, our data suggest that a considerable degree of random error persist in automatic drawing and that intelligent manual contour delineation is still superior to automated one.
A final rationale for discrepancies of the results may lie in the longitudinal part of heart contraction. Naturally, SAX based disc summation method would be susceptible to the underestimation of the longitudinal motion and this underestimation could contribute the difference between the results of conventional and semi-automated methods. Persistent overestimation of EF and SV of automated and semi-method, with and without statistical significance respectively, may reflects this underestimation of longitudinal contraction in the measurement based on the 2D SAX disc summation.
There were limitations to our study. First, our study didn't include the reproducibility of the measurements. Although the semi-automated model has less observer dependency owing to its automaticity, this model also has some points that observer's personal tendency, such as the identification of mitral valve insertion points, can make differences. Second, the best measurement tool is probably that which is the most accurate. Because we measured LV parameters in vivo, we were not able to test for accuracy, among three measurement methods. Therefore, we cannot assess whether one of the three measurements was superior to the others. Third, the sample size 62 may be relatively small, considering wide range of the patients age (7-77). ; However, the sample represents a population with wide spectrum of cardiac disorder from congenital heart disease to chronic ischemic heart disease.
In conclusion, automatic heart model-based tool grossly overestimate EF, SV and LV mass compared with conventional 2D or adjusted automated methods. Hence this method should be used with caution, and its benefits should be weighed against its difference with manual or semi-automated measurement. Therefore, even though the method saves a considerable amount of efforts, further manual adjustment should be considered in critical clinical cases.
Figures and Tables
References
1. Levy D, Garrison RJ, Savage DD, Kannel WB, Castelli WP. Prognostic implications of echocardiographically determined left ventricular mass in the Framingham Heart Study. N Engl J Med. 1990; 322:1561–1566.
2. Kondo C, Caputo GR, Semelka R, Foster E, Shimakawa A, Higgins CB. Right and left ventricular stroke volume measurements with velocity-encoded cine MR imaging: in vitro and in vivo validation. AJR Am J Roentgenol. 1991; 157:9–16.
3. Koren MJ, Devereux RB, Casale PN, Savage DD, Laragh JH. Relation of left ventricular mass and geometry to morbidity and mortality in uncomplicated essential hypertension. Ann Intern Med. 1991; 114:345–352.
4. Rerkpattanapipat P, Mazur W, Link KM, Hundley WG. Assessment of cardiac function with MR imaging. Magn Reson Imaging Clin N Am. 2003; 11:67.
5. Francois CJ, Fieno DS, Shors SM, Finn JP. Left ventricular mass: manual and automatic segmentation of true fisp and flash cine MR images in dogs and pigs. Radiology. 2004; 230:389–395.
6. Alfakih K, Reid S, Jones T, Sivananthan M. Assessment of ventricular function and mass by cardiac magnetic resonance imaging. Eur Radiol. 2004; 14:1813–1822.
7. Thiele H, Nagel E, Paetsch I, et al. Functional cardiac MR imaging with steady-state free precession (ssfp) significantly improves endocardial border delineation without contrast agents. J Magn Reson Imaging. 2001; 14:362–367.
8. Bloomgarden DC, Fayad ZA, Ferrari VA, Chin B, Sutton MG, Axel L. Global cardiac function using fast breath-hold MRI: validation of new acquisition and analysis techniques. Magn Reson Med. 1997; 37:683–692.
9. Childs H, Ma L, Ma M, et al. Comparison of long and short axis quantification of left ventricular volume parameters by cardiovascular magnetic resonance, with ex-vivo validation. J Cardiovasc Magn Reson. 2011; 13:40.
10. Hautvast GL, Salton CJ, Chuang ML, Breeuwer M, O'Donnell CJ, Manning WJ. Accurate computer-aided quantification of left ventricular parameters: experience in 1555 cardiac magnetic resonance studies from the framingham heart study. Magn Reson Med. 2012; 67:1478–1486.
11. Katz J, Milliken MC, Stray-Gundersen J, et al. Estimation of human myocardial mass with MR imaging. Radiology. 1988; 169:495–498.
12. Shapiro MD, Nieman K, Nasir K, et al. Utility of cardiovascular magnetic resonance to predict left ventricular recovery after primary percutaneous coronary intervention for patients presenting with acute ST-segment elevation myocardial infarction. Am J Cardiol. 2007; 100:211–216.
13. Dulce MC, Mostbeck GH, Friese KK, Caputo GR, Higgins CB. Quantification of the left ventricular volumes and function with cine MR imaging: comparison of geometric models with three-dimensional data. Radiology. 1993; 188:371–376.
14. Vogel-Claussen J, Finn JP, Gomes AS, et al. Left ventricular papillary muscle mass: relationship to left ventricular mass and volumes by magnetic resonance imaging. J Comput Assist Tomogr. 2006; 30:426–432.