- Open Access
Evaluation of diastolic function by three-dimensional volume tracking of the mitral annulus with cardiovascular magnetic resonance: comparison with tissue Doppler imaging
Journal of Cardiovascular Magnetic Resonance volume 16, Article number: 71 (2014)
Measurement of mitral annulus (MA) dynamics is an important component of the evaluation of left ventricular (LV) diastolic function; MA velocities are commonly measured using tissue Doppler imaging (TDI). This study aimed to examine the clinical potential of a semi-automated cardiovascular magnetic resonance (CMR) technique for quantifying global LV diastolic function, using 3D volume tracking of the MA with conventional cine-CMR images.
124 consecutive patients with normal ejection fraction underwent both clinically indicated transthoracic echocardiography (TTE) and CMR within 2 months. Interpolated 3D reconstruction of the MA over time was performed with semi-automated atrioventricular junction (AVJ) tracking in long-axis cine-CMR images, producing an MA sweep volume over the cardiac cycle. CMR-based diastolic function was evaluated, using the following parameters: peak volume sweep rates in early diastole (PSRE) and atrial systole (PSRA), PSRE/PSRA ratio, deceleration time of sweep volume (DTSV), and 50% diastolic sweep volume recovery time (DSVRT50); these were compared with TTE diastolic measurements.
Patients with TTE-based diastolic dysfunction (n = 62) showed significantly different normalized MA sweep volume profiles compared to those with TTE-based normal diastolic function (n = 62), including a lower PSRE (5.25 ± 1.38 s−1 vs. 7.72 ± 1.7 s−1), a higher PSRA (6.56 ± 1.99 s−1 vs. 4.67 ± 1.38 s−1), a lower PSRE/PSRA ratio (0.9 ± 0.44 vs. 1.82 ± 0.69), a longer DTSV (144 ± 55 ms vs. 96 ± 37 ms), and a longer DSVRT50 (25.0 ± 11.0% vs. 15.6 ± 4.0%) (all p < 0.05). CMR diastolic parameters were independent predictors of TTE-based diastolic dysfunction after adjusting for left ventricular hypertrophy, hypertension, and coronary artery disease. Good correlations were observed between CMR PSRE/PSRA and early-to-late diastolic annular velocity ratios (e′/a′) measured by TDI (r = 0.756 to 0.828, p < 0.001).
3D MA sweep volumes generated by semi-automated AVJ tracking in routinely acquired CMR images yielded diastolic parameters that were effective in identifying patients with diastolic dysfunction when correlated with TTE-based variables.
Left ventricular (LV) diastolic dysfunction refers to mechanical abnormalities that impede effective volume filling during diastole. This is the main mechanism of heart failure with preserved ejection fraction (HFPEF), which accounts for up to 50% of overall cases of heart failure . Despite the absence of systolic dysfunction, patients with HFPEF still experience high rates of morbidity and mortality . Factors associated with diastolic dysfunction include hypertrophic, infiltrative, and dilated cardiomyopathies, advanced age, coronary artery disease (CAD), and systemic hypertension (HTN) . Moreover, diastolic dysfunction often precedes systolic dysfunction, and is a sensitive indicator of early coronary and hypertensive heart disease . As a result, the characterization of diastolic function, even in asymptomatic subjects, may have important implications for prognosis and treatment strategies.
Transthoracic echocardiography (TTE) is currently the standard method used in the noninvasive evaluation of diastolic function. Commonly, pulsed-wave Doppler of transmitral flow (TMF) is used, with pulmonary venous flow (PVF) and tissue Doppler imaging (TDI), to assess myocardial relaxation and filling pressures. However, TTE measurements have limitations, due to variability in sampling locations and ultrasound beam alignment . In particular, TDI-derived velocities are position dependent, and may not accurately reflect global diastolic function in the presence of regional dysfunction . One study demonstrated that the average of four annular site measurements is more accurate than a single site measurement . While x-ray computed tomography has also been introduced as a potential way to assess diastolic function , this is still investigational.
Cardiovascular magnetic resonance (CMR) has emerged as a promising alternative modality for quantifying diastolic dysfunction . CMR's good spatial resolution, field of view, and range of tissue contrasts allow correlation of cardiac function with morphologic and tissue characteristics. Measurements analogous to TTE measurements, such as TMF, PVF, and myocardial tissue velocities, can be performed with velocity phase-contrast CMR, which has demonstrated excellent agreement with Doppler ,. Similar agreement with TTE has also been reported for CMR derived LV volume-time curves ,. In addition, myocardial tagging  methods have had promising results. However, despite recent advances, CMR-based methods for diastology remain in their infancy. The existing methods require prolonged imaging acquisition and post-processing times, preventing routine use of CMR in assessing diastolic function .
Mitral annular (MA) dynamics has been previously studied with various methods -; however, the relationship between CMR MA motion and diastolic function remains largely unexplored. In a prior study, a 2D manual tracking technique in assessing CMR MA motion had been reported . In this paper, we aimed to expand on this prior technique and identify patients with diastolic dysfunction by using 3D MA sweep volumes calculated from routinely acquired long-axis cine-CMR images, comparing the results to TTE.
Patient population and study design
The research protocol was approved by the local Institutional Review Board. Study subjects were identified from retrospectively reviewed clinical data of all consecutive patients who underwent both clinically indicated TTE and CMR at New York University Langone Medical Center between June 2011 and December 2013. Inclusion criteria were: (1) cine-CMR performed within 2 months of TTE, (2) TTE included assessment of diastolic function, and (3) normal LV systolic function on both CMR and TTE (EF > 50%). Patients were excluded for: (1) incomplete CMR images for MA reconstruction, (2) poor CMR image quality due to arrhythmias (atrial fibrillation, ventricular ectopy), or (3) sinus tachycardia that rendered interpretation of the cardiac phases difficult . Clinical data at the time of imaging were collected, including symptoms, CAD history, presence of structural heart disease on CMR, and cardiovascular risk factors such as HTN, hyperlipidemia, diabetes, and tobacco use. A positive history of CAD was defined as: reports of myocardial infarction, evidence of disease on diagnostic tests, or previous revascularization.
Study subjects were divided into normal diastolic function and diastolic dysfunction groups, using TDI velocities as reference. According to echocardiography guidelines ,, diastolic dysfunction was defined as septal e′ MA peak velocity < 8 cm/s, lateral e′ MA peak velocity < 10 cm/s, or maximal LA volume index ≥ 34 mL/m2. Abnormal LA dilatation was defined as LA volume index ≥ 29 mL/m2.
TTE was performed using conventional equipment: Philips iE33 (Phillips Medical Systems, MA, USA), General Electric Vivid 7 (General Electric Medical Systems, WI, USA), or Siemens SC2000 (Siemens Medical Solutions, CA, USA). Standard apical and parasternal views were obtained. Function and dimensional measurements were performed according to accepted guidelines . TMF velocities were recorded with pulsed-wave Doppler sampling at the mitral valve leaflet tips in the apical 4-chamber view. TDI was performed with pulsed-wave Doppler sampling at the MA junctions of the septal and lateral walls in the apical 4-chamber view.
Conventional cine-CMR (2D steady state free precession pulse sequence) was performed using 1.5 or 3 T MRI systems (Avanto, Tim Trio, Siemens Medical Solutions, Erlangen, Germany) with a 6-element phased array anterior receiver coil, as part of routine clinical protocol. Cine images were acquired in multiplanar short-and long-axis views with retrospective electrocardiographic gating. Cine short-axis views were used for routine evaluation of global cardiac function. Typical imaging parameters were as follows: TR/TE = 2.4/1.4 ms, temporal resolution = 37–63 ms, in-plane spatial resolution = 1.4 mm × 1.4 mm - 1.8 mm × 1.8 mm, flip angle = 51°, slice thickness = 6 mm, receiver bandwidth = 930 Hz/pixel, and 25 reconstructed phases per cardiac cycle.
CMR image analysis
AVJ tracking and interpolated 3D reconstruction of mitral annulus
The atrioventricular junction (AVJ), the septal and lateral junction between the left atrium and ventricle, was tracked in two-, three-, and four-chamber long-axis CMR views over the cardiac cycle . This was performed semi-automatically in lab-written MATLAB (MathWorks Inc., MA, USA) using normalized cross-correlation (NCC), a well-known feature-tracking algorithm , that has been used in MA location tracking in CMR . Briefly, the user initially selects an AVJ point in the image of cardiac phase 1 (Figure 1A). This point becomes the center of a square template (red square) that undergoes NCC with a larger ROI (yellow square) in the image of phase 2 (Figure 1B), producing a correlation coefficient map (empirically, the use of a 20-pixel square in the initial phase and a 40-pixel square in the subsequent phase produced good tracking results that were not sensitive to the precise number of pixels used in this study).
In this calculated correlation map, the location of the maximum-intensity pixel indicates the best corresponding location of the initial template feature in the new ROI, which is thus taken as the AVJ location in phase 2 (Figure 1C). This new point becomes the center of a new template that then undergoes similar calculation of the NCC with the corresponding ROI in the next phase; the process is similarly automatically repeated for the subsequent phases, in order to track the AVJ locations throughout the cardiac cycle. The software also allows for optional interactive user correction of the AVJ location, should there be tracking errors (e.g., due to blurring artifacts). After user correction for a specific phase, if needed, automated AVJ tracking with the NCC method is then performed again for the subsequent phases in the image series, using the corrected location as the starting point.
For each cardiac phase, two AVJ locations were tracked in each of the two-, three-, and four-chamber long-axis CMR views (Figure 2A), creating six independent spatial 3D coordinates within the AVJ that were tracked over the cardiac cycle. Note that 2D image coordinates were transformed into the corresponding 3D space coordinates, using spatial information about the image acquisition locations from DICOM headers. A 3D spline curve was then used to interpolate these 6 distinct 3D spatial coordinates sampled within the MA at each cardiac cycle phase, in order to create a reconstruction of the MA in 3D space (Figure 2B), using solid-modeling software Rhinoceros (McNeel, WA, USA). This is analogous to methods that have been used to reconstruct the 3D shape of the MA from ultrasound images . 3D spline curves have been similarly used to reconstruct the 3D MA structure in prior CMR studies ,. A 3D MA incremental sweep volume (Vn) was then generated for each cardiac phase tn, using the MA areas at tn-1 and tn, and the 3D distance (positive or negative) through which the MA traversed (Figure 2C, D); the net sweep volume at a given cardiac phase was derived from the sum of the incremental volumes starting from end-diastole.
The resultant MA sweep volume curve was manually divided into its corresponding cardiac cycle intervals according to transitions in its slope: systole, early diastole (ED), mid-diastole (MD), and atrial systole (AS) (Figure 3). In addition, its first derivative was calculated to characterize the MA sweep volume rate in each interval. For comparisons between different heart sizes, which may confound the absolute amounts of myocardial relaxation, the sweep volumes (and rates) were normalized to the corresponding end systolic sweep volume (ESSV). To quantify measures related to LV myocardial relaxation, the following CMR diastolic parameters were derived from the curves: (1) peak sweep rates (PSR) in ED and AS; (2) average sweep rates in ED, MD, AS; (3) percentage sweep volume recovery in ED and AS, defined by the proportion of ESSV that the MA had recovered during the given interval; and (4) the ratio of peak sweep rate in ED to peak sweep rate in AS (PSRE/PSRA). Note that PSRE/PSRA is analogous to the e′/a′ ratio used in TDI.
Three time-interval parameters were selected to characterize diastolic function from the sweep volumes: (1) acceleration time (ATSV), measured from ED onset to the time of PSRE; (2) deceleration time (DTSV), measured by linear extrapolation of PSRE to baseline; and (3) 50% diastolic sweep volume recovery (DSVRT50), defined as the time required in diastole for the MA to recover 50% of its end systolic sweep volume (adjusted for RR interval); this index had been introduced in a prior study to describe CMR LV volume filling .
All statistical analyses were performed using SPSS version 20.0 (SPSS Inc., Chicago, Illinois). Continuous variables were presented as mean ± SD, whereas categorical data were presented as frequencies with percentages. Differences between groups were evaluated using either the Student t-test for normally distributed data, or the Mann–Whitney U-test for asymmetrically distributed data; Shapiro-Wilk test was used to assess normality of distributions. Fisher's exact test was used to compare differences in categorical data. In addition, the association between CMR- and TTE-based parameters was assessed using the Spearman's rank correlation test. For classification performance, receiver operating characteristic (ROC) analyses were used to evaluate the ability of CMR parameters to identify patients with TTE-based diastolic dysfunction. Optimal cutoff values for these parameters were obtained, along with their associated sensitivity and specificity. Multivariate logistic regression was used to determine the predictive ability of CMR parameters with respect to TTE-based outcomes, after controlling for age, HTN, left ventricular hypertrophy (LVH), and CAD status. A two-tailed value of p < 0.05 was considered statistically significant.
Intraobserver and interobserver variability were evaluated for 10 randomly selected subjects using intraclass correlation coefficients (ICC), Bland-Altman analyses, and Pearson's r correlation coefficients. For interobserver variability assessment, analyses were performed by two independent CMR readers who were blinded to TTE and each other's results. For intraobserver variability assessment, analyses were performed by one user who re-analyzed the same 10 subjects after a 3 week interval.
A total of 124 consecutive subjects met the selection criteria and were included in the study. The time difference between the dates of CMR and TTE was 16 ± 16 days. Based on TTE findings, subjects were classified as either having normal diastolic function (n = 62) or diastolic dysfunction (n = 62).
Table 1 details the baseline characteristics of the study population. Compared to the normal diastolic function group, those with TTE-based diastolic dysfunction were older in age (61 ± 14 vs. 40 ± 14 years, p < 0.001) and were more likely to be male (p = 0.046). Clinically, the diastolic dysfunction group had increased frequencies of HTN (61% vs. 26%, p < 0.001), hyperlipidemia (56% vs. 27%, p = 0.002), and CAD (35% vs. 11%, p = 0.003). In terms of symptoms, patients with diastolic dysfunction were more likely to complain of dyspnea (31% vs. 15%, p = 0.052). A higher proportion of subjects in the diastolic dysfunction group, compared to those in the normal diastolic function group, had CMR evidence of structural heart disease (76% vs. 32%, p < 0.001); in particular, the prevalence of LVH was significantly higher in the diastolic dysfunction group (39% vs. 10%, p < 0.001). These clinical findings were consistent with prior reports on the prevalence of diastolic dysfunction . Morphologically, TTE-derived LV mass index and LA volume index were both significantly higher in the diastolic dysfunction group.
AVJ tracking and MA reconstruction
Our CMR method successfully generated MA sweep volume curves for all subjects in the study. The approximate processing time for each case was less than 3 minutes. The processing included the initial user delineation of AVJ points, review and possible correction of semi-automated tracking results, MA reconstruction in 3D, and the identification of cardiac intervals within the sweep volume curve. Manual corrections of AVJ locations were primarily necessary in images with blurring artifacts affecting the regions of interest. Approximately 30% of cine image series required at least one manual correction. Cases that did not require manual correction were analyzed in under 1 minute.
Figure 4A illustrates representative patients, comparing the normalized MA sweep volume curves and their derivatives between subjects with normal diastolic function and diastolic dysfunction. These curves show that the MA of a diastolic dysfunction subject requires more time to recover from its end systolic volume state, as quantified by the parameter DSVRT50. Furthermore, the derivative curve provides a visual representation of the normalized MA sweep rate at each cardiac interval. In normal diastolic function, there was typically a higher MA sweep rate during early diastole that was accompanied by a relatively lower sweep rate during atrial systole; this was reversed in diastolic dysfunction. In these cases, peak MA sweep rates (PSRE and PSRA) on CMR were in qualitative relative agreement with the echocardiographically-derived early (e′) and late (a′) TDI annular velocity profiles (Figure 4B). In addition, the utility of DTSV was demonstrated in the sweep volume derivative curve, where the diastolic dysfunction subjects showed a relative flattening of the deceleration component in early diastole.
LV diastolic function
Table 2 describes the CMR and TTE diastolic parameters for the two groups. Compared to patients with TTE-based normal diastolic function, those with TTE-based diastolic dysfunction had slower normalized peak and average MA sweep rates in early diastole (5.3 ± 1.4 s−1 vs. 7.7 ± 1.7 s−1 and 2.5 ± 0.6 s−1 vs. 4.2 ± 1.0 s−1, respectively, both p < 0.001). The diastolic dysfunction group also had a lower PSRE/PSRA ratio (0.9 ± 0.4 vs. 1.8 ± 0.7, p < 0.001), due to their relatively higher peak sweep rate in atrial systole. During early diastole, the MA of hearts with diastolic dysfunction recovered a smaller proportion of their end systolic sweep volume (54 ± 12% vs. 69 ± 10%, p < 0.001), and their MA sweep volume constituted a smaller percentage of their stroke volume (7.6 ± 2.2% vs. 11.2 ± 2.9%, p < 0.001). Diastolic dysfunction was also differentiated by its longer DTSV (144 ± 55 ms vs. 96 ± 37 ms, p < 0.001) and longer DSVRT50 (25 ± 11% vs. 16 ± 4.0%, p < 0.001).
Figure 5 illustrates the correlations of MA velocities between CMR and TDI measured at different locations. Although there was a strong positive correlation between CMR PSRE/PSRA and echo septal e′/a′ (r = 0.756, p < 0.001), the correlation between PSRE/PSRA and lateral e′/a′ was even stronger (r = 0.803, p < 0.001). The strongest correlation, however, was seen between PSRE/PSRA and the average e′/a′ of septal and lateral positions (r = 0.828, p < 0.001).
Classification performance of CMR parameters
Table 3 summarizes the ROC analysis of CMR diastolic parameters in the identification of TTE-based diastolic dysfunction. Based on areas under the curve (AUC), PSRE and PSRE/PSRA ratio demonstrated excellent classification power (both AUC = 0.88), as did DTSV (AUC = 0.83) and DSVRT50 (AUC = 0.79). The optimal cutoff value for PSRE was 6.7 s−1, with a sensitivity of 89% and specificity of 74%; the optimal cutoff for PSRE/PSRA was 1.17, with a sensitivity of 84% and specificity of 81%. In comparison, DTSV and DSVRT50 exhibited slightly lower accuracies, with optimal cutoffs at 101 ms and 19.45%, respectively. A recovery threshold of 50% of end-systolic sweep volume was chosen for the parameter DSVRT50 because this generated the best classification performance compared to other recovery thresholds (thresholds of 40%, 60%, 70%, and 80% produced AUC's ranging from 0.60 to 0.75).
CMR prediction of diastolic dysfunction
Multivariate logistic regression analysis demonstrated that PSRE, PSRE/PSRA, and DSVRT50 were significant independent predictors of TTE-based diastolic dysfunction after controlling for age, HTN, LVH, and CAD status (Table 4). DTSV was near significant as an independent predictor with an odds ratio of 1.12 per 10 ms increment (p = 0.077). In addition, PSRE/PSRA had an odds ratio of 0.86 per 0.1 increment (p = 0.016), whereas PSRE and DSVRT50 had odds ratios of 0.54 (p = 0.007) and 1.14 (p = 0.013), respectively.
CMR diastolic parameters were used as covariates to construct predictive models with logistic regression. A model including both PSRE and PSRE/PSRA showed an improved AUC of 0.895 in predicting TTE-based diastolic dysfunction. When adding in DTSV and DSVRT50, the AUC improved slightly to 0.896. The latter model with all four parameters was used to demonstrate the prevalence of CMR-based diastolic dysfunction in a number of important clinical conditions in our sample population (Figure 6). As shown, the prevalence of CMR-based diastolic dysfunction was almost twice as great in patients with LVH compared to patients without LVH (p = 0.002). Similar results were observed in patients with HTN and CAD. In addition, the same predictive model showed that LA dilatation was seen in 24 of 65 subjects (37%) with CMR-based diastolic dysfunction and 6 of 59 (10%) subjects with CMR-based normal diastolic function (p < 0.001). Based on ROC analyses, LA volume index demonstrated only moderate classification powers in predicting both CMR-based and TTE-based diastolic dysfunction (AUC = 0.64, AUC = 0.67, respectively).
Table 5 shows the intraobserver and interoberserver variabilities using ICC, the Bland-Altman test, and Pearson's r correlation. CMR diastolic parameters demonstrated good consistency in terms of correlations (ICC ranged from 0.88 to 0.97) with minimal degree of biases for both intraobserver and interobserver measurements.
This study demonstrates that CMR diastolic parameters derived from 3D MA sweep volumes were reproducible and could accurately differentiate between patients with normal diastolic function and diastolic dysfunction, as established by TDI. CMR-based measurements of peak sweep rates were also strongly correlated with analogous TDI velocity indices. Diastolic assessments using CMR 3D MA sweep volumes were validated in a diverse population of 124 subjects with normal systolic function measures, which revealed an increased prevalence of CMR based diastolic dysfunction in patients with HTN, CAD, and LVH. This suggests that MA sweep volume may contribute to the evaluation of LV diastolic function, potentially providing additional prognostic information and guidance that could be useful in management before frank heart failure occurs .
Our method addresses several limitations encountered by current CMR methods of assessing diastolic function. Notably, routine use of phase-contrast imaging and tagging have been hindered by the need for additional imaging pulse-sequences and associated acquisition time. In addition, manual calculation of LV volumetric curves requires extensive post-processing effort and has associated subjectivity, while automatic measurements of these curves require proprietary and complex image-processing algorithms . Our novel CMR method may provide a more practical option, since it assesses diastolic function by using conventional cine-CMR long-axis images that are already routinely acquired in clinical CMR examinations, without the need for additional imaging. This also enables users to analyze pre-existing cases.
To facilitate the measurement of MA sweep volumes, AVJ tracking was performed semi-automatically, using a simple NCC feature-tracking algorithm. NCC is an image-processing technique that has been used for many applications in motion tracking. NCC provides a means to assess the degree of similarity ("correlation") between two images, as a function of pixel position . Assuming that the immediate surroundings of a given point in an initial image provide the features of a template, the algorithm "slides" and centers this template at each pixel within a neighborhood of this point in a subsequent image, and a correlation coefficient is calculated between the template and the subsequent image for each such pixel position. This process produces a map of the correlation between the point (and its surroundings) in the initial image and the points within that corresponding neighborhood in the subsequent image; the location of the maximum correlation represents the likely location of the initial point (and its surrounding template) within the subsequent image. In this study, NCC enabled semi-automated AVJ tracking that reduced the long post-processing time previously required with manual AVJ tracking . In addition, the use of reproducibility analyses here showed that the use of NCC limited variability in calculated MA sweep volume measurements between users.
Our results showed that patients with TTE-based diastolic dysfunction had lower normalized peak and average MA sweep rates in early diastole, resulting in the reversal of PSRE/PSRA ratios. Likewise, prior echocardiographic studies have reported the analogous phenomenon of TDI e′/a′ MA velocity reversal in diastolic dysfunction ,. The present study also showed a relative increase in peak MA sweep rate in atrial systole in these patients. This can potentially be explained by the effect of LA dilatation associated with diastolic dysfunction, which may have led to compensatory enhancement of LA contraction, due to activation of the Frank-Starling mechanism . Based on our data, LA dilatation was significantly associated with both CMR-based and TTE-based diastolic dysfunction compared to their respective controls. In addition, ROC analysis showed that the combined use of PSRE and PSRE/PSRA improved performance in identifying patients with TTE-based diastolic dysfunction. The ability to characterize atrial systole with our methods is an advantage compared to techniques such as MR tagging, which cannot reliably produce late diastolic strain rates because of fading signal intensities due to T1 relaxation .
MA sweep mechanics were reported here in terms of "percentage sweep volume recovery". This index showed that patients with normal diastolic function recovered nearly 70% of the end-systolic sweep volume by the end of early diastole, or 10.6 mL in absolute volume, but patients with diastolic dysfunction only recovered 54%, or 6.8 mL. These findings were in rough overall agreement with a prior study that reported an average MA excursion volume of 6 mL in nine healthy subjects . However, this prior study involved subjects with a considerably lower average stroke volume (52 mL vs. 93 mL), which may explain the moderate difference in results. Overall, the CMR parameters discussed here were all consistent in demonstrating blunted early-diastolic MA kinetics in patients with abnormal TDI velocities.
We also investigated the time intervals of the MA sweep volume. The results showed that patients with TTE-based diastolic dysfunction had significantly longer DTSV (144 ± 55 ms) compared to the normal diastolic function group (96 ± 37 ms). ATSV was slightly longer in the diastolic dysfunction group, but the difference was not significant. Both findings were consistent with a published TDI study that examined acceleration and deceleration times pertaining to MA velocities (normal DT = 84 ms, diastolic dysfunction DT = 156 to 168 ms) .
DSVRT50 was an additional CMR diastolic function parameter, which measured the time needed for the MA to recover 50% of its end systolic sweep volume, adjusted for R-R interval. On average, longer DSVRT50 times were observed in patients with TTE-based diastolic dysfunction. This parameter was analogous to an index used in a prior study named "diastolic volume recovery" , which accounted for both heart rate and volume status when assessing diastolic dysfunction in CMR LV volumetric filling. Lastly, it is important to note that PSRE, PSRE/PSRA, and DSVRT50 each independently predicted TTE-based diastolic dysfunction after controlling for age, LVH, HTN, and CAD.
Clinically, the study demonstrated that patients with HTN, LVH, and CAD had increased prevalence of CMR-based diastolic dysfunction. Among the LVH subgroup, a large percentage (20/23 patients) was also diagnosed with hypertrophic cardiomyopathy. Similar results regarding these clinical subgroups have been reported in the past . In general, HTN, LVH, and CAD are all well-known causes of diastolic dysfunction, and often co-exist in patients with cardiovascular disease. These conditions have been implicated in the disruption of myocyte calcium reuptake that governs the active phase of relaxation, as well as mechanisms that increase collagen synthesis and the ventricular stiffness that governs the passive phase of relaxation ,. Because previous echocardiography studies have shown that diastolic dysfunction often develops before any clinically apparent changes of cardiac function in predisposed individuals , such patients undergoing routine CMR may benefit from additional assessment of their diastolic function using the proposed method. In particular, the CMR-derived diastolic function variables may contribute additional value to the conventional CMR characterization of tissue morphology and contrast agent delayed enhancement.
The MA shape was reconstructed in 3D with six imaged AVJ points, derived from three standard LV long-axis planes situated approximately 60 degrees apart in space. By accounting for MA 3D motion in its entirety, in both transverse and longitudinal directions, the proposed method was hypothesized to be a better representation of global diastolic function than TDI, which only measures longitudinal motion. Typically in TDI, longitudinal velocities are measured regionally at the septal and lateral sides of the MA. However, prior studies have reported significant regional variability in LV myocardial velocities . Specifically, early TDI velocities in the lateral, inferior, and posterior locations were often greater than the septal and anterior segments. As a result, varying regional wall motion abnormalities were shown to undermine the correlation between TDI velocities and invasively acquired LV end-diastolic pressure and tau . In current practice, the average of TDI velocities is often used to better characterize global LV diastolic function . In this study, CMR PSRE/PSRA correlated better with the average TDI e′/a′, compared to single-site measurements, thereby supporting the hypothesis that the proposed method inherently better characterizes global LV function. It is interesting to note that CMR PSRE/PSRA correlated better with lateral TDI e′/a′ than it did with septal e′/a′. This finding is consistent with studies showing that, among single site measurements, lateral TDI velocities have the best correlations with global LV filling pressures and invasive indices of LV compliance .
The CMR-derived MA sweep volume was found here to effectively characterize global diastolic function. Prior echocardiography studies have demonstrated that MA excursion volume plays an important role in LV function ,. In addition, reports have long recognized the relationship between LV long-axis function and MA displacement . For example, the longitudinal MA excursion in diastole envelops and effectively transfers blood from the atrium to the ventricle, separately from the flow of blood across the location of the MA. This can occur even while the blood remains relatively stationary in relation to the apex. Similarly, in atrial systole, the MA is pulled away from the apex by the atrial pectinate muscles, "over" the blood, to facilitate further ventricular filling (and also increasing the ventricular pre-stretch, thus augmenting contractility). As a result, MA excursion is a central component of diastole, along with the transmitral pressure gradient.
There were some limitations in our study design. First, this was a retrospective investigation and most patients did not have their imaging studies performed on the same day. Although the subjects were in stable clinical condition between studies, changes in cardiac function between studies were possible. This may explain some discrepancy in diastolic function as identified by CMR and TTE. Other potential reasons for discrepancy include differences in temporal resolution between the methods, as well as sampling location, where regional dysfunction identified by TDI may not correspond with global dysfunction as assessed by MA sweep volume , or vice versa. Second, the study used TTE as the standard reference for diastolic function instead of invasive hemodynamic data. However, routine use of invasive studies is not practical. Furthermore, despite various limitations of TTE , it remains the generally accepted non-invasive reference method for diastolic function assessment. Third, the study lacked follow-up data that might indicate the prognostic significance of the results of the proposed method.
There are also several limitations of the proposed technique. First, the three long-axis images used to reconstruct the MA were acquired at similar, but different, points in time, e.g., with potential differences in respiratory excursion. Second, the accuracy of NCC feature-tracking can be affected by image scaling, rotation, and distortion ; as with other feature-tracking approaches, it is dependent on there being some suitable distinctive image structure around the point being tracked. In conventional long-axis cine-CMR images, significant scaling or rotation were rarely seen in AVJ motion. However, distortion of the AVJ was relatively common during rapid motion, due to resulting blurring artifacts, which necessitated manual correction of the AVJ tracking results in a minority of cases. Lastly, when using routinely acquired long-axis images, only six points could be used to reconstruct the 3D MA. However, a prior animal study was able to effectively characterize MA dynamics by implanting only eight equally spaced radiopaque markers on the MA .
Future studies involving a larger cohort with assessment of the prognostic utility of the information may offer further evaluation of the potential clinical utility of CMR-derived 3D MA sweep volume as a technique for evaluation of diastolic dysfunction. A larger population could also allow stratification of the different diastolic dysfunction stages. In addition, there may be particular value in applying the technique in evaluating specific patient groups with potential ventricular dysfunction, such as athletes, patients with hypertrophic, dilated, or restrictive cardiomyopathies, or patients undergoing anthracycline therapy.
3D MA sweep volumes generated by semi-automated AVJ tracking in routinely acquired CMR images demonstrated diastolic parameters that were effective in identifying patients with TTE-based diastolic dysfunction. These parameters also correlated well with TDI MA velocities, and may serve as a useful tool for characterizing patients with diastolic dysfunction.
VW: conception and design, development and implementation of algorithm and software, acquisition, analysis and interpretation of data, statistical analysis, drafting of manuscript; JC: acquisition, analysis and interpretation of data, statistical analysis, critical revision of manuscript; SC: acquisition, analysis and interpretation of data, development of original software, critical revision of manuscript; SB: development of original software, critical revision of manuscript; LA: conception and design, acquisition, analysis and interpretation of data, critical revision of manuscript; All authors read and approved the final manuscript.
Cardiovascular magnetic resonance
Tissue Doppler imaging
Peak sweep rate
Deceleration time of sweep volume
Acceleration time of sweep volume
50% diastolic sweep volume recovery time
End systolic sweep volume
Intraclass correlation coefficients
Hogg K, Swedberg K, McMurray J: Heart failure with preserved left ventricular systolic function; epidemiology, clinical characteristics, and prognosis. J Am Coll Cardiol. 2004, 43: 317-327. 10.1016/j.jacc.2003.07.046.
The survival of patients with heart failure with preserved or reduced left ventricular ejection fraction: an individual patient data meta-analysis. Eur Heart J. 2012, 33: 1750-1757. 10.1093/eurheartj/ehr254.
Paulus WJ, Tschöpe C, Sanderson JE, Rusconi C, Flachskampf FA, Rademakers FE, Marino P, Smiseth OA, De Keulenaer G, Leite-Moreira AF: How to diagnose diastolic heart failure: a consensus statement on the diagnosis of heart failure with normal left ventricular ejection fraction by the Heart Failure and Echocardiography Associations of the European Society of Cardiology. Eur Heart J. 2007, 28: 2539-2550. 10.1093/eurheartj/ehm037.
Dougherty AH, Naccarelli GV, Gray EL, Hicks CH, Goldstein RA: Congestive heart failure with normal systolic function. Am J Cardiol. 1984, 54: 778-782. 10.1016/S0002-9149(84)80207-6.
Fyrenius A, Wigström L, Bolger AF, Ebbers T, Öhman KP, Karlsson M, Wranne B, Engvall J: Pitfalls in Doppler evaluation of diastolic function: insights from 3-dimensional magnetic resonance imaging. J Am Soc Echocardiogr. 1999, 12: 817-826. 10.1016/S0894-7317(99)70186-0.
Rivas-Gotz C, Manolios M, Thohan V, Nagueh SF: Impact of left ventricular ejection fraction on estimation of left ventricular filling pressures using tissue Doppler and flow propagation velocity. Am J Cardiol. 2003, 91: 780-784. 10.1016/S0002-9149(02)03433-1.
Rivas-Gotz C, Khoury DS, Manolios M, Rao L, Kopelen HA, Nagueh SF: Time interval between onset of mitral inflow and onset of early diastolic velocity by tissue Doppler: a novel index of left ventricular relaxation: experimental studies and clinical application. J Am Coll Cardiol. 2003, 42: 1463-1470. 10.1016/S0735-1097(03)01034-9.
Boogers MJ, Van Werkhoven JM, Schuijf JD, Delgado V, El-Naggar HM, Boersma E, Nucifora G, van der Geest RJ, Paelinck BP, Kroft LJ: Feasibility of Diastolic Function Assessment With Cardiac CT: Feasibility Study in Comparison With Tissue Doppler Imaging. JACC Cardiovasc Imaging. 2011, 4: 246-256. 10.1016/j.jcmg.2010.11.017.
Caudron J, Fares J, Bauer F, Dacher J-N: Evaluation of left ventricular diastolic function with cardiac MR imaging. Radiographics. 2011, 31: 239-259. 10.1148/rg.311105049.
Rathi VK, Doyle M, Yamrozik J, Williams RB, Caruppannan K, Truman C, Vido D, Biederman RW: Routine evaluation of left ventricular diastolic function by cardiovascular magnetic resonance: A practical approach. J Cardiovasc Magn Reson. 2008, 10: A212-10.1186/1532-429X-10-S1-A212.
Paelinck BP, De Roos A, Bax JJ, Bosmans JM, van Der Geest RJ, Dhondt D, Parizel PM, Vrints CJ, Lamb HJ: Feasibility of tissue magnetic resonance imaging: a pilot study in comparison with tissue Doppler imaging and invasive measurement. J Am Coll Cardiol. 2005, 45: 1109-1116. 10.1016/j.jacc.2004.12.051.
Kawaji K, Codella NC, Prince MR, Chu CW, Shakoor A, LaBounty TM, Min JK, Swaminathan RV, Devereux RB, Wang Y, Weinsaft JW: Automated segmentation of routine clinical cardiac magnetic resonance imaging for assessment of left ventricular diastolic dysfunction. Circ Cardiovasc Imaging. 2009, 2: 476-484. 10.1161/CIRCIMAGING.109.879304.
Mendoza DD, Codella NC, Wang Y, Prince MR, Sethi S, Manoushagian SJ, Kawaji K, Min JK, LaBounty TM, Devereux RB: Impact of diastolic dysfunction severity on global left ventricular volumetric filling-assessment by automated segmentation of routine cine cardiovascular magnetic resonance. J Cardiovasc Magn Reson. 2010, 12: 46-10.1186/1532-429X-12-46.
Garot J: The study of diastole by tagged MRI: are we nearly there yet?. Eur Heart J. 2004, 25: 1376-1377. 10.1016/j.ehj.2004.06.019.
Fernández-Pérez G, Duarte R, Corral De la Calle M, Calatayud J, Sánchez González J: Analysis of left ventricular diastolic function using magnetic resonance imaging. Radiología. 2012, 54: 295-305. 10.1016/j.rx.2011.09.018.
Nammas W, El-Okda E: Atrioventricular plane displacement: does it predict in-hospital outcome after acute myocardial infarction?. Eur Rev Med Pharmacol Sci. 2012, 16: 16-21.
Carlhäl C, Wigström L, Heiberg E, Karlsson M, Bolger A, Nylander E: Contribution of mitral annular excursion and shape dynamics to total left ventricular volume change. Am J Physiol Heart Circ Physiol. 2004, 287: H1836-H1841. 10.1152/ajpheart.00103.2004.
Maffessanti F, Gripari P, Pontone G, Andreini D, Bertella E, Mushtaq S, Tamborini G, Fusini L, Pepi M, Caiani EG: Three-dimensional dynamic assessment of tricuspid and mitral annuli using cardiovascular magnetic resonance. Eur Heart J Cardiovasc Imaging. 2013, 14: 986-995. 10.1093/ehjci/jet004.
Saba SG, Chung S, Bhagavatula S, Donnino R, Srichai MB, Saric M, Katz SD, Axel L: A novel and practical cardiovascular magnetic resonance method to quantify mitral annular excursion and recoil applied to hypertrophic cardiomyopathy. J Cardiovasc Magn Reson. 2014, 16: 1-9. 10.1186/1532-429X-16-35.
Nagueh SF, Appleton CP, Gillebert TC, Marino PN, Oh JK, Smiseth OA, Waggoner AD, Flachskampf FA, Pellikka PA, Evangelisa A: Recommendations for the evaluation of left ventricular diastolic function by echocardiography. Eur J Echocardiogr. 2009, 10: 165-193. 10.1093/ejechocard/jep007.
Morris DA, Gailani M, Vaz Pérez A, Blaschke F, Dietz R, Haverkamp W, Özcelik C: Left atrial systolic and diastolic dysfunction in heart failure with normal left ventricular ejection fraction. J Am Soc Echocardiogr. 2011, 24: 651-662. 10.1016/j.echo.2011.02.004.
Lang RM, Bierig M, Devereux RB, Flachskampf FA, Foster E, Pellikka PA, Picard MH, Roman MJ, Seward J, Shanewise J: Recommendations for chamber quantification. Eur J Echocardiogr. 2006, 7: 79-108. 10.1016/j.euje.2005.12.014.
Lewis J: Fast normalized cross-correlation. Vision Interface. 1995, 10: 120-123.
Hii A, Hann CE, Chase JG, Van Houten EE: Fast normalized cross correlation for motion tracking using basis functions. Comput Methods Programs Biomed. 2006, 82: 144-156. 10.1016/j.cmpb.2006.02.007.
Levine RA, Handschumacher MD, Sanfilippo A, Hagege A, Harrigan P, Marshall J, Weyman A: Three-dimensional echocardiographic reconstruction of the mitral valve, with implications for the diagnosis of mitral valve prolapse. Circulation. 1989, 80: 589-598. 10.1161/01.CIR.80.3.589.
Ennis DB, Rudd-Barnard GR, Li B, Fonseca CG, Young AA, Cowan BR, Stewart RA: Changes in mitral annular geometry and dynamics with β-blockade in patients with degenerative mitral valve disease. Circ Cardiovasc Imaging. 2010, 3: 687-693. 10.1161/CIRCIMAGING.110.959171.
Fischer M, Baessler A, Hense H, Hengstenberg C, Muscholl M, Holmer S, Döring A, Broeckel U, Riegger G, Schunkert H: Prevalence of left ventricular diastolic dysfunction in the community. Results from a Doppler echocardiographic-based survey of a population sample. Eur Heart J. 2003, 24: 320-328. 10.1016/S0195-668X(02)00428-1.
Zile MR, Brutsaert DL: New concepts in diastolic dysfunction and diastolic heart failure: Part I diagnosis, prognosis, and measurements of diastolic function. Circulation. 2002, 105: 1387-1393. 10.1161/hc1102.105289.
Sohn DW, Chai IH, Lee DJ, Kim HC, Kim HS, Oh BH, Lee MM, Park YB, Choi YS, Seo JD: Assessment of mitral annulus velocity by Doppler tissue imaging in the evaluation of left ventricular diastolic function. J Am Coll Cardiol. 1997, 30: 474-480. 10.1016/S0735-1097(97)88335-0.
Edvardsen T, Rosen BD, Pan L, Jerosch-Herold M, Lai S, Hundley WG, Sinha S, Kronmal RA, Bluemke DA, Lima JA: Regional diastolic dysfunction in individuals with left ventricular hypertrophy measured by tagged magnetic resonance imaging—the Multi-Ethnic Study of Atherosclerosis (MESA). Am Heart J. 2006, 151: 109-114. 10.1016/j.ahj.2005.02.018.
Nagueh SF, Middleton KJ, Kopelen HA, Zoghbi WA, Quiñones MA: Doppler tissue imaging: a noninvasive technique for evaluation of left ventricular relaxation and estimation of filling pressures. J Am Coll Cardiol. 1997, 30: 1527-1533. 10.1016/S0735-1097(97)00344-6.
Choudhury L, Gheorghiade M, Bonow RO: Coronary artery disease in patients with heart failure and preserved systolic function. Am J Cardiol. 2002, 89: 719-722. 10.1016/S0002-9149(01)02345-1.
Verma A, Solomon SD: Diastolic dysfunction as a link between hypertension and heart failure. Med Clin North Am. 2009, 93: 647-664. 10.1016/j.mcna.2009.02.013.
Aeschbacher BC, Hutter D, Fuhrer J, Weidmann P, Delacrétaz E, Allemann Y: Diastolic dysfunction precedes myocardial hypertrophy in the development of hypertension. Am J Hypertens. 2001, 14: 106-113. 10.1016/S0895-7061(00)01245-0.
Galiuto L, Ignone G, DeMaria AN: Contraction and relaxation velocities of the normal left ventricle using pulsed-wave tissue Doppler echocardiography. Am J Cardiol. 1998, 81: 609-614. 10.1016/S0002-9149(97)00990-9.
Fujii T, Yoshioka K, Nakano M, Nakazawa G, Amino M, Masuda N, Shinozaki N, Kanda S, Ogata N, Deguchi Y: Regional wall motion abnormality at the lateral wall disturbs correlations between tissue Doppler E/e′ ratios and left ventricular diastolic performance parameters measured by invasive methods. J Echocardiogr. 2013, 11: 138-146. 10.1007/s12574-013-0193-x.
Kasner M, Westermann D, Steendijk P, Gaub R, Wilkenshoff U, Weitmann K, Hoffmann W, Poller W, Schultheiss H-P, Pauschinger M: Utility of Doppler Echocardiography and Tissue Doppler Imaging in the Estimation of Diastolic Function in Heart Failure With Normal Ejection Fraction A Comparative Doppler-Conductance Catheterization Study. Circulation. 2007, 116: 637-647. 10.1161/CIRCULATIONAHA.106.661983.
Carlhäl C, Kindberg K, Wigström L, Daughters GT, Miller DC, Karlsson M, Ingels NB: Contribution of mitral annular dynamics to LV diastolic filling with alteration in preload and inotropic state. Am J Physiol Heart Circ Physiol. 2007, 293: H1473-H1479. 10.1152/ajpheart.00208.2007.
Henein MY, Gibson DG: Normal long axis function. Heart. 1999, 81: 111-113.
This study was supported by the National Institutes of Health grant R21-HL108218.
We thank James S. Babb, PhD (Associate Professor, New York University School of Medicine), who provided statistical advice.
The authors declare that they have no competing interests.
Authors’ original submitted files for images
About this article
Cite this article
Wu, V., Chyou, J.Y., Chung, S. et al. Evaluation of diastolic function by three-dimensional volume tracking of the mitral annulus with cardiovascular magnetic resonance: comparison with tissue Doppler imaging. J Cardiovasc Magn Reson 16, 71 (2014). https://doi.org/10.1186/s12968-014-0071-3
- Diastolic function
- Mitral annulus
- Cardiovascular magnetic resonance