- Open Access
Validation of a rapid semi-automated method to assess left atrial longitudinal phasic strains on cine cardiovascular magnetic resonance imaging
Journal of Cardiovascular Magnetic Resonancevolume 20, Article number: 71 (2018)
Abnormal left atrial (LA) function is a marker of cardiac dysfunction and adverse cardiovascular outcome, but is difficult to assess, and hence not, routinely quantified. We aimed to determine the feasibility and effectiveness of a fast method to measure long-axis LA strain and strain rate (SR) with standard cardiovascular magnetic resonance (CMR) compared to conventional feature tracking (FT) derived longitudinal strain.
We studied 50 normal controls, 30 patients with hypertrophic cardiomyopathy, and 100 heart failure (HF) patients, including 40 with reduced ejection fraction (HFrEF), 30 mid-range ejection fraction (HFmrEF) and 30 preserved ejection fraction (HFpEF). LA longitudinal strain and SR parameters were derived by tracking the distance between the left atrioventricular junction and a user-defined point at the mid posterior LA wall on standard cine CMR two- and four-chamber views. LA performance was analyzed at three distinct cardiac phases: reservoir function (reservoir strain εs and strain rate SRs), conduit function (conduit strain εe and strain rate SRe) and booster pump function (booster strain εa and strain rate SRa).
There was good agreement between LA longitudinal strain and SR assessed using the fast and conventional FT-CMR approaches (r = 0.89 to 0.99, p < 0.001). The fast strain and SRs showed a better intra- and inter-observer reproducibility and a 55% reduction in evaluation time (85 ± 10 vs. 190 ± 12 s, p < 0.001) compared to FT-CMR. Fast LA measurements in normal controls were 35.3 ± 5.2% for εs, 18.1 ± 4.3% for εe, 17.2 ± 3.5% for εa, and 1.8 ± 0.4, − 2.0 ± 0.5, − 2.3 ± 0.6 s− 1 for the respective phasic SRs. Significantly reduced LA strains and SRs were observed in all patient groups compared to normal controls. Patients with HFpEF and HFmrEF had significantly smaller εs, SRs, εe and SRe than hypertrophic cardiomyopathy, and HFmrEF had significantly impaired LA reservoir and booster function compared to HFpEF. The fast LA strains and SRs were similar to FT-CMR for discriminating patients from controls (area under the curve (AUC) = 0.79 to 0.96 vs. 0.76 to 0.93, p = NS).
Novel quantitative LA strain and SR derived from conventional cine CMR images are fast assessable parameters for LA phasic function analysis.
Left atrial (LA) function has been increasingly recognized as an important determinant of cardiovascular morbidity and mortality . The key role of LA is to modulate left ventricular (LV) filling through three phases : 1) Reservoir phase – LA collects pulmonary venous return during LV contraction and isovolumetric relaxation; 2) Conduit phase – LA conducts blood passively into the LV; 3) Booster pump phase – atrial contraction that actively forces blood to the LV.
LA function can be assessed using LA pressure-volume loops , which require invasive measurements. Non-invasive assessment of dynamic changes in LA size, including diameters, areas and volumes, might be inadequate for describing complex LA phasic function . Tissue Doppler imaging (TDI) has been used to measure late diastolic mitral annular velocity during atrial contraction for assessing LA function . Errors arising from angle dependency mitigate its diagnostic accuracy and negate its utility . Strain analyses have been performed using speckle-tracking echocardiography  or cine cardiovascular magnetic resonance (CMR) feature tracking (FT) . However, the former is challenging for the LA due to low signal-to-noise ratio and the thin atrial wall; while the latter is hindered by the complex LA anatomy comprising the LA appendage and pulmonary veins . Furthermore, lack of standardized methodology in strain analysis affects reproducibility of LA strain measurements . Herein, we aimed to evaluate a novel and rapid semi-automatic post-processing method for assessing long-axis strain and strain rate (SR) for the determination of LA phasic longitudinal function from standard cine CMR images. The new method does not require total delineation of the LA contours, but only requires annotations of three distinct anatomical reference points, thus it is less affected by the presence of the LA appendage and pulmonary vein. We hypothesized that LA strain and SR derived with this novel method would require less processing time without compromising accuracy and reproducibility.
The study population consisted of 50 normal controls, and 130 patients (40 with heart failure (HF) with reduced ejection fraction (HFrEF, LV ejection fraction (EF) < 40%); 30 HF with mid-range ejection fraction (HFmrEF, LVEF 40–49%); 30 HF with preserved ejection fraction (HFpEF, LVEF ≥ 50%); and 30 patients with hypertrophic cardiomyopathy (HCM)). Inclusion criteria for HF required the presence of signs or symptoms of HF based on modified Framingham criteria  and prior hospitalization with primary diagnosis of HF. Exclusion criteria for HF included specific subgroups of HF (e.g., amyloidosis, eosinophilic myocarditis, etc.) and isolated right heart disease. HCM patients were recruited from specialized cardiomyopathy clinics. All patients had sinus rhythm during examination. The protocol was approved by the SingHealth Centralised Institutional Review Board and informed consent was obtained from all participants.
CMR acquisitions were performed using a 3 T system (Ingenia, Philips Healthcare, Best, The Netherlands). Balanced steady steady state free precession (bSSFP) end-expiratory breath hold cine images were acquired in multi-planar long-axis views including the 2- and 4-chamber views. Typical parameters were as follows: TR/TE, 3/1 ms; matrix, 240 × 240; flip angle, 45°; field of view, 300 × 300 mm2; pixel bandwidth, 1776 Hz; pixel spacing, 1.25 × 1.25 mm; slice thickness, 8 mm; number of cardiac frames, 30/40 per cardiac cycle.
Feature tracking (FT)
With dedicated QStrain software (Version 2.0, Medis BV, Leiden, The Netherlands), FT was used to track tissue voxel motion on cine CMR images in deriving LA longitudinal strain and SR in both 2- and 4-chamber views. In each view, the LA endocardial contour was manually drawn when the atrium was at its minimum volume after atrial contraction. The length of the contour was then partitioned into a series of 48 evenly spaced points. Automatic contour tracking was then performed by tracking each single point based on a hierarchical algorithm that combined 1- (1D) and 2-dimensional (2D) tracking . Manual contour adjustments were made at the frame when the LA volume was maximal. The LV end-diastolic and end-systolic phases coincide with the minimum and maximum LA volume phases with 30/40 cardiac frames per cycle that was employed in the cine CMR acquisition in this study. Tracking was repeated for three times and global LA longitudinal strain and SR measurements were averaged across all three repetitions in both 2- and 4-chamber views. Figure 1 (left column) shows a representative example of the LA tracking in the 2- and 4-chamber views and the strain data (i.e. the endocardial global longitudinal strain (GLS)). SR results were then obtained by taking the first-order derivative of the strain curve.
Fast long-axis strain
Long-axis strain was assessed by automatically tracking the distance (D) between the left atrioventricular junction and a user-defined point at the mid posterior LA wall on standard CMR 2- and 4-chamber views (Fig. 1 right column). The atrioventricular junctions were selected as the mitral valve insertion points at the septal and lateral borders of the annulus on the 4-chamber view, and the anterior and inferior annular insertion points on the 2-chamber view [12,13,14,15]. The mid-point of posterior LA wall was defined as the intersection point of the LA posterior wall and the LA long-axis (the length of which is used for estimating the LA volume based on area-length formula) . The tracking system in this study used the method of template matching . Briefly, a small square (called mask) was manually drawn in the LV end-diastole frame containing the anatomical point of interest (white squares in Fig. 1 right column), and automatically tracked within a search region that was coincident with the mask and automatically selected by the program in the target frame, typically the next temporal frame in the cardiac cycle. A correlation map was produced using template matching algorithm based on normalized cross correlation, in which the point with highest correlation coefficient indicated the location of the best match. This point was used to update the mask in the target frame and the same procedure was automatically executed iteratively for all subsequent frames over the cardiac cycle. The point tracking was semi-automatic with mask selection in the initial frame as the only user input. Technical details including the processing parameters can be found in our prior study .
Longitudinal strain (ε) at any time point (t) in the cardiac cycle from LV end-diastole (time 0) was calculated as [18, 19]: ε(t) = (D(t) − D0) × 100/D0. Similar to previous definitions from FT-CMR [20, 21], LA reservoir strain (εs), conduit strain (εe) and booster strain (εa) were calculated at t equal to LV end-systole, diastasis and pre-LA systole, respectively, and the corresponding peak SRs (SRs, SRe, SRa) were derived (Fig. 1). Strain and SR parameters were derived for each of the four measured walls (septal and lateral in 4-chamber view and anterior and inferior in 2-chamber view) and the average values were used for analysis.
Left atrium volumetric analysis
LA volumetric analysis was performed using QMass software (Version 8.0, Medis BV). Apical 2- and 4-chamber views were reviewed offline to measure the maximal LA volume at LV end-systole (LAVmax), diastasis LA volume at LV diastole before LA contraction (LAVpreA), and minimal LA volume at LV end-diastole (LAVmin) using biplane area-length method : LA volume (ml) = 0.85*A2C* A4C/L, where A2C and A4C are the LA areas on the 2- and 4-chamber views, and L is the shorter long-axis length of LA in both 2- and 4-chamber views determined as the distance of the perpendicular line bisecting mitral annular plane that intersects with the posterior wall of the LA. LA phasic function, viz. reservoir, conduit and booster pump function were characterized by total LA emptying fraction (LAEF), passive LAEF and active LAEF, respectively, according to the following equations :
Statistical analysis was conducted using SPSS (version 17.0, International Business Machines, Armonk, New York, USA). Continuous data were summarized as mean ± SD. Comparisons of characteristics and strain data between patient groups and control subjects were performed using independent t tests for normally distributed data, Mann-Whitney U tests for non-normally distributed data, and chi-square tests for categorical data. Pearson’s r correlation, Passing-Bablok non-parametric regression, and Bland-Altman plots were used to assess the agreement of fast long-axis strain data and FT derived longitudinal strain. Correlation r larger than 0.7 was interpreted as a strong relationship between two variables. Pearson’s correlation coefficients were calculated to determine the linear association between LA strain and volumetric measurements. A p value less than 0.05 was considered statistically significant.
To evaluate the reproducibility, intra-observer and inter-observer variability were studied on a randomly selected subgroup of 20 cases (10 normal controls and 10 patients) using Bland-Altman analysis and coefficient of variation. For inter-observer variability, measurements were repeated by a second-independent observer, blinded to the first observer’s results. For intra-observer variability study, the analysis was repeated by the same observer who re-analyzed the same 20 cases after 1 week.
Baseline characteristics of study population
Baseline demographics and clinical characteristics of controls (n = 50), all HF groups (HFrEF, n = 40; HFmrEF, n = 30; HFpEF, n = 30) and HCM patients (n = 30) are summarized in Table 1. All HF groups were highly symptomatic (69% New York Heart Association [NYHA] II-IV) with 83% of them on diuretics. There were no significant differences among the groups in age, gender and body surface area (BSA). Patients with HFrEF had significantly lower LVEF, larger LV end-diastolic volume and end-systolic volume indices, and smaller LAEF compared to the other groups (controls, HCM, HFpEF, HFmrEF). Patients with HCM, HFpEF and HFmrEF had significantly higher LA volumes and lower LAEF in all phases compared with normal controls. There were no significant differences between HCM and HFpEF with respect to LA volumes and LAEF.
Fast LA strain measurements
Fast LA strain analysis was successfully performed in all subjects. Averaged strain and SR profiles were calculated from both the 2- and 4-chamber views. Figure 2 shows representative examples of fast LA strain and SR data in normal control, HCM, HFpEF, HFmrEF and HFrEF. Mean values in all normal controls were 35.3 ± 5.2% for reservoir strain, 18.1 ± 4.3% for conduit strain, 17.2 ± 3.5% for booster strain, and 1.8 ± 0.4, − 2.0 ± 0.5, − 2.3 ± 0.6 s− 1 for the respective phasic SRs. In comparison to controls, patients with HCM, HFpEF, HFmrEF and HFrEF had significantly decreased LA strain and SR measurements. Among the patients, HFrEF had lowest LA strain and SR compared to other patient groups (HCM, HFpEF, HFmrEF). Patients with HFpEF and HFmrEF had significantly smaller LA reservoir strain and SR, and conduit strain and SR than HCM, and HFmrEF had significantly impaired LA reservoir and booster function compared to HFpEF (Table 2). The strain and SR results derived by conventional FT-CMR are presented in Additional file 1: Table S1.
Strong correlations were exhibited between volumetric indices and fast LA strain/SR measurements for reservoir (εs vs. Total LAEF: r = 0.92; SRs vs. Total LAEF: r = 0.83), conduit (εe vs. Passive LAEF: r = 0.83; SRe vs. Passive LAEF: r = − 0.76) and booster pump functions (εa vs. Active LAEF: r = 0.88; SRa vs. Active LAEF: r = − 0.84) (Table 3).
LA strain was negatively associated with NYHA class (Fig. 3) in the patient groups (HCM, HFpEF, HFmrEF, HFrEF). Significantly reduced fast LA strain measurements were also observed in patients with mitral regurgitation (n = 54) compared to patients without mitral regurgitation (n = 76) in the patient groups (HCM, HFpEF, HFmrEF, HFrEF) (Fig. 4).
Significant negative correlations were found between Ln NT-proBNP and fast εs (r = − 0.67, p < 0.001), εe (r = − 0.62, p < 0.001), εa (r = − 0.55, p < 0.001), SRs (r = − 0.64, p < 0.001), SRe (r = − 0.60, p < 0.001, absolute values were used for correlation analysis here) and SRa (r = − 0.57, p < 0.001, absolute values were used for correlation analysis here) (Fig. 5) in HF patients.
Comparison with feature tracking strain
Fast LA strain and SR measurements exhibited strong correlation and agreement with FT based strain derived from dedicated software (reservoir strain and SR, r = 0.99 and 0.91; conduit strain and SR, r = 0.96 and 0.91; booster strain and SR, r = 0.96 and 0.89) with low bias and narrow 95% limits of agreement (Figs. 6 and 7).
Prediction of disease states
Table 4 reports results of the receiver operating characteristic (ROC) analysis performed, on a per subject basis, to assess the utility of LA volumetric and functional parameters for characterizing LA functional alterations in the patient groups (HCM, HFpEF, HFmrEF) compared with controls. Compared to LA volumetric and FT based strain parameters, fast LA strain and SR measurements characterizing LA alterations achieved higher levels of accuracy and greater areas under ROC curves (AUCs). The best variable for discriminating patients with HCM, HFpEF or HFmrEF from normal controls was fast LA SR for conduit phase (SRe) with AUC 0.96, sensitivity 93%, and specificity 90%.
Reproducibility and time requirement
Table 5 shows the intra- and inter-observer variability. Intra-observer CV was 3.6–6.2% for fast LA strain and 6.4–7.9% for the SR measurements, with corresponding inter-observer CVs of 4.7–8.9% and 8.1–9.7%, respectively. The mean ± SD time per subject required for fast LA strain measurement in 2- and 4-chamber views was 85 ± 10 s, which was significantly shorter than the mean measurement time of 190 ± 12 s using the FT derived strain analysis which was also evaluated in the 2- and 4-chamber views.
The LA has previously been considered a neglected chamber. The interest in LA function evaluation has resurged over the recent years, with more focus on FT methods based on whole LA delineation that is, in our opinion, more suitable for ventricular function assessment. The LA wall is made up of circumferential and longitudinal muscular bundles. The former is arranged at the base of the atria, while the latter predominates at the parietal walls . Barbier et al. revealed that LA reservoir function is primarily determined by longitudinal descent of the cardiac base and LA chamber stiffness . In this study, we demonstrated the feasibility and effectiveness of a novel and rapid assessable parameter for the analysis of global longitudinal LA function using standard CMR cine images. Our results showed that the fast LA longitudinal strain correlated strongly with FT derived strain measurements, while sensitivity and specificity data were similar. Fast LA strain and SR performed better at differentiating diseased groups from normal controls than LA volumetric measurements with respect to LA reservoir, conduit and contractile booster pump function. Therefore, our approach could supplant conventional FT when the LA longitudinal function is being assessed.
The clinical significance of global longitudinal LA strain has been demonstrated in echocardiographic studies of patients with HF and HCM [5, 25, 26]. LA strain in HFpEF was associated with worsening NYHA functional class . Echocardiographic speckle tracking-assessed LA longitudinal strain correlated with pulmonary capillary wedge pressure – and therefore LV filling pressure – in subjects with advanced systolic HF . LA longitudinal strain and SR were shown to be sensitive discriminators of HCM, non-HCM LV hypertrophy and healthy controls . Other than HF and HCM, impaired echocardiography-assessed LA longitudinal strain is also seen in other conditions – valvular heart diseases [27, 28], atrial stunning , and hypertension and diabetes . Of note, LA longitudinal strain demonstrated prognostic utility in HF , aortic stenosis , and after myocardial infarct . In a study of HFpEF patients, LA reservoir strain correlated negatively with adverse events, and was an independent predictor of the composite outcome of cardiovascular hospitalization or death .
Fast LA strain vs. feature tracking strain
Instead of tracking 48 points on the contour line by the FT analysis, we presented a fast strain parameter that required the automatic tracking of only 3 anatomically discrete points, and demonstrated its feasibility and capability to characterize the phasic longitudinal LA function compared to the conventional FT-based strain and SR measurements. Our technique is analogous to , in which Riffel et al. assessed LV longitudinal function by using the long-axis strain derived from the distance of the epicardial apical border to the midpoint of the line connecting the mitral valve insertion points (LAS-epi/mid). The LAS-epi/mid showed a high correlation with conventional FT analysis results, and was non-inferior for discriminating patients with cardiomyopathies from healthy controls. Another similar approach  was applied to analyze long-axis strain in the LV by measuring the distances from the mitral valve insertions to the epicardial apex. They observed good correlation between long-axis strain – assessed with echocardiography and CMR – and infarct mass in patients with prior myocardial infarction. The authors in these two studies [18, 19] used manual tracking to measure the displacement of the mitral annulus between LV end-systolic and end-diastolic phases of the cardiac cycle, which yield only LV systolic strain parameters. Ours was the first study to apply semi-automatic tracking through all temporal phases of the cardiac cycle, which allowed for the dynamic evaluation of the complex phasic LA function.
FT in LA is more challenging and time-consuming than in the LV. The thin atrial wall, complex LA anatomy, presence of the LA appendage and pulmonary veins all influence tracking quality, and are potential sources of error . In , the tracking quality was inadequate in 10.8% of total segments; while in another study , CMR images were not interpretable in 13% of the study population due to the failure of LA wall tracking. In our study, we observed that FT was not accurate in determining the valve annulus position in about 15% of cardiac frames, which vitiated the accuracy of FT-derived phasic LA strain and SR measurements. In contrast, fast LA strain analysis was successful in all study subjects with good reproducibility. The target feature was not correctly located automatically in less than 2% of total discrete points tracked. This was due to temporal blurring, and necessitated manual correction.
LA function in patients with HCM, HFpEF, HFmrEF and HFrEF
This study examined LA structure and function in HF by volumetric and strain analyses in the three HF phenotypes. Compared with controls, all patients with HFpEF, HFmrEF and HFrEF displayed abnormal LA size and function characterized by increased LA maximal and minimal volume indices, decreased LAEF, and decreased phasic LA strains and SRs.
In general, LA function was less impaired in HFpEF and HFmrEF than in HFrEF. HFmrEF is a new category of HF, intermediate HFrEF and HFpEF . Patients with HFmrEF had clinical characteristics that are more similar to those of HFpEF than HFrEF . However, studies of LA structure and function in the HFmrEF cohort are lacking. We found progressively decreased fast LA strain and SR measurements in HFpEF, HFmrEF and HFrEF. In all HF groups, LA function was significantly associated with right heart function (i.e. right ventricular ejection fraction (RVEF), r = 0.73 for RVEF vs. reservoir strain εs; r = 0.59 for RVEF vs. conduit strain εe; and r = 0.64 for RVEF vs. booster strain εa, Additional file 2: Figure S1); and negatively with NYHA class. These data revealed the important role of LA dysfunction in HF. LA function can be used to monitor HF progression, and plausibly, strategies to maintain or restore normal LA function can help to improve NYHA scores and mitigate progression of right heart dysfunction . Significant negative associations between NT-proBNP levels and LA phasic functions were found in our study. On multivariate analysis, only LVEF (β = − 0.037, p = 0.001) and fast LA reservoir strain εs (β = − 0.073, p < 0.0001) persisted as significant predictors of NT-proBNP, indicating the incremental contribution of fast LA strain measurements in HF.
Decreases in reservoir and conduit function and increases in booster pump function in HCM patients have been reported previously [1, 38, 39]. In our study, fast LA strain analysis also revealed significantly decreased reservoir and conduit strains and SRs in HCM compared to controls. Booster strain and SR, however, were also decreased in our HCM patients, which differed from the findings in . In a recent study, Kowallick et al.  characterized LA strain in HCM according to the extent of LV hypertrophy and fibrosis, and found that LA booster pump function was impaired in HCM with severe late gadolinium enhancement (LGE) only (LGE ≥ 20%). Indeed, our HCM subjects may represent those with more advanced disease. In our data, 67% (20/30) of the HCM patients had moderate to severe patchy fibrosis. In addition, 43% (13/30) of our HCM patients had significant left ventricular outflow tract (LVOT) obstruction, which has also been reported to negatively influence LA mechanics .
The current study found significant deteriorations in LA function with presence of mitral regurgitation. Among the subjects with NYHA class II-IV in the patient group (HCM and HF), patients with moderate and severe mitral regurgitation had LA reservoir strain εs reduced by 35% and 49%, respectively, compared to those with mild or no mitral regurgitation. The corresponding reductions were 25% and 40%, and 44% and 56% for LA conduit strain εe and booster strain εa, respectively. Studies have shown that global peak LA longitudinal strain is inversely correlated with the degree of mitral regurgitation , which may possibly be explained by LA abnormalities such as myocyte hypertrophy, interstitial fibrosis and decreased metalloproteinase expression .
All fast strain and SR measurements demonstrated better intra- and inter-observer reproducibility than FT-CMR (Table 5). The tracking system in this study used the method of template matching, which is a semi-automatic algorithm for searching and finding the location of a template image within a larger image. The template selection in the initial frame is the only user input. Hence, the observed intra- and inter-observer variability are expected to be insignificant since operator dependence is relegated exclusively to variation in site selection of the atrioventricular junction points and the user-defined point at the mid posterior LA wall. Fast strain and SR were evaluated with the tracking of only three discrete spatially separated points, and were thus less affected by the presence of the LA appendage and pulmonary veins, which are the primary factors hindering effective endocardial tracking in conventional FT strain analysis.
Prior studies have investigated the impact of repeated measures on reproducibility of FT-CMR [43, 44] where the differences in intra- and inter-observer variability were assessed based on single and averaged measurements (two and three repetitions with subsequent averaging of results, respectively). It was found that averaging of the results of repeated analyses improves the reproducibility of LA strain measurements. Although the benefit of repeated measures is relatively low, considering that doubling or tripling of analysis times would be required [43, 45], FT-CMR analysis was performed three times in the current study in order to further maximize its reproducibility . The evaluation time presented herein and the ~ 55% in time savings were based on single tracking (i.e. the average time of three repetitions). In addition, inter-vendor comparisons have been performed in prior studies [43, 44, 46] to determine differences in strain measurements between commercially available FT-CMR software packages (TomTec, Medis QStrain and Circle Cardiovascular Imaging). The inter-vendor agreement was reasonably good for LV global circumferential strain and longitudinal strain, but was lower in right ventricular global longitudinal and radial strains [44, 46]. The time for post-processing of a given case did not vary between the different types of software . No prior studies could be found that assessed the inter-vendor agreement for LA strain measurements, therefore, further studies are warranted to determine the interchangeability of LA longitudinal strain and strain rate between different FT-CMR software solutions.
Limitation and long-term prospects
The fast LA strain is an index of global long axis LA function in the longitudinal direction, and does not provide information on segmental and regional deformation, circumferential and radial strain, or myocardial strain outside the LA. In the assessment of LA phasic function, our proposed fast LA strain/SR measurement method using standard CMR cine images yields similar number and type of metrics as FT-based assessment using proprietary software , i.e. LA longitudinal reservoir, conduit and active strains and the corresponding SRs.
The reported experience of using FT-CMR LA strain parameter assessment for clinical diagnosis is quite recent, and encompasses cohorts of patients undergoing atrial fibrillation catheter ablation [8, 47, 48], implantable cardioverter-defibrillator implant [47, 49], community volunteers without cardiovascular disease at baseline [50, 51], patients with acute myocarditis , obesity , and Ebstein’s anomaly . These studies established the feasibility and validity of LA strain parameter measurement using CMR; and in  found negative correlation between LA strain and LA fibrosis. Notably, none of the published studies investigated the role of FT-derived LA strain parameters for HF diagnosis.
We focused on the use of CMR for HF diagnosis and have demonstrated that the fast method for assessing phasic LA longitudinal strains and SRs was useful in discriminating subjects with HFrEF and HFpEF – including a recently described and guideline-adopted classification of HFmrEF  – as well as HCM, which shares morphological similarities with HFpEF (concentric LV hypertrophy). Our findings corroborated the results of prior echocardiographic studies [5, 25, 26], and underscore the clinical diagnostic utility of and rationale for measuring LA longitudinal strain using either echocardiography or CMR.
We believe our fast method is interchangeable with, and can eventually supplant FT for several reasons. First, in this study we performed both FT and fast method for LA strain and SR measurements among HF, HCM and normal subjects, and demonstrated high degrees of correlation (Fig. 6) and agreement (Fig. 7) between the two methods. Importantly, the fast method had higher intra- and inter-observer reproducibility with tighter coefficient of variation (Table 5). Second, the FT method based on whole chamber contour delineation is, in our opinion, more suited for LV function assessment rather than the LA, where LA appendage and pulmonary vein anatomy interrupt the LA wall outline and complicate contour tracing and frame-to-frame phasic FT. Third, the LA wall is made up of circumferential muscle fibres at the base of the atria and predominantly longitudinal muscular bundles at the parietal walls . Whereas FT includes the basal segment at the posterior aspect of the LA in the chamber contour for calculation, our fast method excludes it – arguably more directly mimicking the function of the longitudinal muscle bundles, which along with LA chamber stiffness, determines LA reservoir function . Fourth, our automatic algorithm dispenses with manual LA contour tracing, requires minimal operator input and significantly shorter processing time, which should garner wider access and acceptance than FT.
The long-term prospect is to make the analysis fully automated, by automatically detecting the three anatomical reference points based on deep learning based technique in a large dataset. Second, the present study focused primarily on assessment of LA longitudinal function. A similarly comprehensive approach for the right atrial (RA) long-axis strain represents a potentially researchable goal in the future for evaluation of RA function. Third, this approach can be applied to derive papillary muscle longitudinal strain by automatically tracking the distance between papillary muscle tips and the left atrioventricular junction points throughout the cardiac cycle. The papillary muscle tethering distance and strain have been demonstrated to have significant association with ischemic mitral regurgitation . Fourth, prior study has revealed distinct associations between metabolic perturbations with LA phasic function in healthy aging population . Future investigations using both imaging (e.g. our fast LA strain analysis) and molecular approach in the patient groups may help to identify mechanisms involved in cardiovascular diseases in specific cohorts. Further studies are warranted on the association between booster LA strain and LA ejection force , as the latter is a measure of atrial contractile function and plays an important role in LV diastolic filling.
The presented fast LA strain and SR measurements represent reliable and expeditiously calculable parameters for quantifying longitudinal LA deformation using routine clinical cine CMR images without specific acquisition protocol or special software tools. Fast LA strain measurements exhibited high correlation with conventional LA FT strain analysis and were comparable in discriminating patients with HF and HCM from normal controls, were less affected by the LA appendage and pulmonary veins, and required significantly less processing time.
- ε a :
- ε e :
- ε s :
Area under ROC curve
Body surface area
Balanced steady state free precession
Cardiovascular magnetic resonance
Global longitudinal strain
Heart failure with mid-range ejection fraction
Heart failure with preserved ejection fraction
Heart failure with reduced ejection fraction
Left atrium/left atrial
Left atrial emptying fraction
- LAVmax :
Maximal LA volume
- LAVmin :
Minimal LA volume
- LAVpreA :
LA volume before LA contraction
Late gadolinium enhancement
Left ventricle/left ventricular
Left ventricular outflow tract
New York Heart Association
Right atrium/right atrial
Receiver operating characteristic
- SRa :
Booster strain rate
- SRe :
Conduit strain rate
- SRs :
Reservoir strain rate
Tissue Doppler imaging
Hoit BD. Left atrial size and function: role in prognosis. J Am Coll Cardiol. 2014;63:493–505.
Blume GG, Mcleod CJ, Barnes ME, Seward JB, Pellikka PA, Bastiansen PM, et al. Left atrial function: physiology, assessment, and clinical implications. Eur J Echocardiogr. 2011;12:421–30.
Stefanadis C, Dernellis J, Toutouzas P. A clinical appraisal of left atrial function. Eur Heart J. 2001;22:22–36.
Rossi A, Gheorghiade M, Triposkiadis F, Solomon SD, Pieske B, Butler J. Left atrium in heart failure with preserved ejection fraction: structure, function, and significance. Circ Heart Fail. 2014;7:1042–9.
Paraskevaidis IA, Panou F, Papadopoulos C, Farmakis D, Parissis J, Ikonomidis I, et al. Evaluation of left atrial longitudinal function in patients with hypertrophic cardiomyopathy: a tissue Doppler imaging and two-dimensional strain study. Heart. 2009;95:483–9.
Kadappu KK, Thomas L. Tissue Doppler imaging in echocardiography: value and limitations. Heart Lung Circ. 2015;24:224–33.
Pathan F, D'Elia N, Nolan MT, Marwick TH, Negishi K. Normal ranges of left atrial strain by speckle-tracking echocardiography: a systematic review and meta-analysis. J Am Soc Echocardiogr. 2017;30:59–70.
Inoue YY, Alissa A, Khurram IM, Fukumoto K, Habibi M, Venkatesh BA, et al. Quantitative tissue-tracking cardiac magnetic resonance (CMR) of left atrial deformation and the risk of stroke in patients with atrial fibrillation. J Am Heart Assoc. 2015;4:e001844.
Smiseth OA, Torp H, Opdahl A, Haugaa KH, Urheim S. Myocardial strain imaging: how useful is it in clinical decision making? Eur Heart J. 2016;37:1196–207.
Ho KK, Anderson KM, Kannel WB, Grossman W, Levy D. Survival after the onset of congestive heart failure in Framingham heart study subjects. Circulation. 1993;88:107–15.
Schuster A, Hor KN, Kowallick JT, Beerbaum P, Kutty S. Cardiovascular magnetic resonance myocardial feature tracking: concepts and clinical applications. Circ Cardiovasc Imaging. 2016;9:e004077.
Leng S, Zhao XD, Huang FQ, Wong JI, Su BY, Allen JC, et al. Automated quantitative assessment of cardiovascular magnetic resonance-derived atrioventricular junction velocities. Am J Physiol Heart Circ Physiol. 2015;309:H1923–35.
Leng S, Jiang M, Zhao XD, Allen JC, Kassab GS, Ouyang RZ, et al. Three-dimensional tricuspid annular motion analysis from cardiac magnetic resonance feature-tracking. Ann Biomed Eng. 2016;44:3522–38.
Leng S, Zhang S, Jiang M, Zhao XD, Wu R, Allen J, et al. Imaging 4D morphology and dynamics of mitral annulus in humans using cardiac cine MR feature tracking. Sci Rep. 2018;8:81.
Xiao PD, Zhao XD, Leng S, Tan RS, Wong P, Zhong L. A software tool for heart AVJ motion tracking using cine cardiovascular magnetic resonance images. IEEE J Transl Eng Health Med. 2017;5:1900412.
Lang RM, Badano LP, Mor-Avi V, Afilalo J, Armstrong A, Ernande L, et al. Recommendations for cardiac chamber quantification by echocardiography in adults: an update from the American Society of Echocardiography and the European Association of Cardiovascular Imaging. J Am Soc Echocardiogr. 2015;28:1–39.
Gonzalez RC, Woods RE. Digital image processing. 4th ed. NY: Pearson/Prentice Hall; 2018.
Riffel JH, Andre F, Maertens M, Rost F, Keller MG, Giusca S, et al. Fast assessment of long axis strain with standard cardiovascular magnetic resonance: a validation study of a novel parameter with reference values. J Cardiovasc Magn Reson. 2015;17:69.
Gjesdal O, Almeida AL, Hopp E, Beitnes JO, Lunde K, Smith HJ, et al. Long axis strain by MRI and echocardiography in a postmyocardial infarct population. J Magn Reson Imaging. 2014;40:1247–51.
Kowallick JT, Kutty S, Edelmann F, Chiribiri A, Villa A, Steinmetz M, et al. Quantification of left atrial strain and strain rate using cardiovascular magnetic resonance myocardial feature tracking: a feasibility study. J Cardiovasc Magn Reson. 2014;16:60.
Evin M, Cluzel P, Lamy J, Rosenbaum D, Kusmia S, Defrance C, et al. Assessment of left atrial function by MRI myocardial feature tracking. J Magn Reson Imaging. 2015;42:379–89.
Sievers B, Kirchberg S, Addo M, Bakan A, Brandts B, Trappe HJ. Assessment of left atrial volumes in sinus rhythm and atrial fibrillation using the biplane area-length method and cardiovascular magnetic resonance imaging with TrueFISP. J Cardiovasc Magn Reson. 2004;6:855–63.
Wang K, Ho SY, Gibson DG, Anderson RH. Architecture of atrial musculature in humans. Br Heart J. 1995;73:559–65.
Barbier P, Solomon SB, Schiller NB, Glantz SA. Left atrial relaxation and left ventricular systolic function determine left atrial reservoir function. Circulation. 1999;100:427–36.
Morris DA, Gailani M, Vaz Pérez A, Blaschke F, Dietz R, Haverkamp W, et al. Left atrial systolic and diastolic dysfunction in heart failure with normal left ventricular ejection fraction. J Am Soc Echocardiogr. 2011;24:651–62.
Cameli M, Lisi M, Mondillo S, Padeletti M, Ballo P, Tsioulpas C, et al. Left atrial longitudinal strain by speckle tracking echocardiography correlates well with left ventricular filling pressures in patients with heart failure. Cardiovasc Ultrasound. 2010;8:14.
Cameli M, Lisi M, Giacomin E, Caputo M, Navarri R, Malandrino A, et al. Chronic mitral regurgitation: left atrial deformation analysis by two-dimensional speckle tracking echocardiography. Echocardiography. 2011;28:327–34.
O'Connor K, Magne J, Rosca M, Piérard LA, Lancellotti P. Left atrial function and remodelling in aortic stenosis. Eur J Echocardiogr. 2011;12:299–305.
Kaya EB, Tokgözoglu L, Aytemir K, Kocabas U, Tülümen E, Deveci OS, et al. Atrial myocardial deformation properties are temporarily reduced after cardioversion for atrial fibrillation and correlate well with left atrial appendage function. Eur J Echocardiogr. 2008;9:472–7.
Mondillo S, Cameli M, Caputo ML, Lisi M, Palmerini E, Padeletti M, et al. Early detection of left atrial strain abnormalities by speckle-tracking in hypertensive and diabetic patients with normal left atrial size. J Am Soc Echocardiogr. 2011;24:898–908.
Freed BH, Daruwalla V, Cheng JY, Aguilar FG, Beussink L, Choi A, et al. Prognostic utility and clinical significance of cardiac mechanics in heart failure with preserved ejection fraction: importance of left atrial strain. Circ Cardiovasc Imaging. 2016;9:e003754.
Galli E, Fournet M, Chabanne C, Lelong B, Leguerrier A, Flecher E, et al. Prognostic value of left atrial reservoir function in patients with severe aortic stenosis: a 2D speckle-tracking echocardiographic study. Eur Heart J Cardiovasc Imaging. 2016;17:533–41.
Antoni ML, ten Brinke EA, Atary JZ, Marsan NA, Holman ER, Schalij MJ, et al. Left atrial strain is related to adverse events in patients after acute myocardial infarction treated with primary percutaneous coronary intervention. Heart. 2011;97:1332–7.
Habibi M, Chahal H, Opdahl A, Gjesdal O, Helle-Valle TM, Heckbert SR, et al. Association of CMR-measured LA function with heart failure development: results from the MESA study. JACC Cardiovasc Imaging. 2014;7:570–9.
Ponikowski P, Voors AA, Anker SD, Bueno H, Cleland JG, Coats AJ, et al. 2016 ESC Guidelines for the diagnosis and treatment of acute and chronic heart failure: The Task Force for the diagnosis and treatment of acute and chronic heart failure of the European Society of Cardiology (ESC). Developed with the special contribution of the Heart Failure Association (HFA) of the ESC. Eur J Heart Fail. 2016;20:592.
Hsu JJ, Ziaeian B, Fonarow GC. Heart failure with mid-range (borderline) ejection fraction: clinical implications and future directions. JACC Heart Fail. 2017;5:763–71.
Melenovsky V, Hwang SJ, Redfield MM, Zakeri R, Lin G, Borlaug BA. Left atrial remodeling and function in advanced heart failure with preserved or reduced ejection fraction. Circ Heart Fail. 2015;8:295–303.
Anwar AM, Soliman OI, Nemes A, Geleijnse ML, ten Cate FJ. An integrated approach to determine left atrial volume, mass and function in hypertrophic cardiomyopathy by two-dimensional echocardiography. Int J Cardiovasc Imaging. 2008;24:45–52.
Anwar AM, Soliman OI, Geleijnse ML, Michels M, Vletter WB, Nemes A, et al. Assessment of left atrial ejection force in hypertrophic cardiomyopathy using real-time three-dimensional echocardiography. J Am Soc Echocardiogr. 2007;20:744–8.
Kowallick JT, Silva Vieira M, Kutty S, Lotz J, Hasenfu G, Chiribiri A, et al. Left atrial performance in the course of hypertrophic cardiomyopathy: relation to left ventricular hypertrophy and fibrosis. Investig Radiol. 2017;52:177–85.
Williams LK, Chan RH, Carasso S, Durand M, Misurka J, Crean AM, Ralph-Edwards A, et al. Effect of left ventricular outflow tract obstruction on left atrial mechanics in hypertrophic cardiomyopathy. Biomed Res Int. 2015;2015:481245.
Cameli M, Incampo E, Mondillo S. Left atrial deformation: useful index for early detection of cardiac damage in chronic mitral regurgitation. Int J Cardiol Heart Vasc. 2017;17:17–22.
Schuster A, Stahnke VC, Unterberg-Buchwald C, Kowallick JT, Lamata P, Steinmetz M, et al. Cardiovascular magnetic resonance feature-tracking assessment of myocardial mechanics: Intervendor agreement and considerations regarding reproducibility. Clin Radiol. 2015;70:989–98.
Gertz RJ, Lange T, Kowallick JT, Backhaus SJ, Steinmetz M, Staab W, et al. Inter-vendor reproducibility of left and right ventricular cardiovascular magnetic resonance myocardial feature-tracking. PLoS One. 2018;13:e0193746.
Dick A, Schmidt B, Michels G, Bunck AC, Maintz D, Baeßler B. Left and right atrial feature tracking in acute myocarditis: a feasibility study. Eur J Radiol. 2017;89:72–80.
Bourfiss M, Vigneault DM, Aliyari Ghasebeh M, Murray B, James CA, Tichnell C, et al. Feature tracking CMR reveals abnormal strain in preclinical arrhythmogenic right ventricular dysplasia/cardiomyopathy: a multisoftware feasibility and clinical implementation study. J Cardiovasc Magn Reson. 2017;19:66.
Tao S, Ciuffo LA, Lima JAC, Wu KC, Ashikaga H. Quantifying left atrial structure and function using single-plane tissue-tracking cardiac magnetic resonance. Magn Reson Imaging. 2017;42:130–8.
Gucuk Ipek E, Marine JE, Habibi M, Chrispin J, Lima J, Rickard J, et al. Association of left atrial function with incident atypical atrial flutter after atrial fibrillation ablation. Heart Rhythm. 2016;13:391–8.
Tao S, Ashikaga H, Ciuffo LA, Yoneyama K, Lima JAC, Frank TF, et al. Impaired left atrial function predicts inappropriate shocks in primary prevention implantable cardioverter-defibrillator candidates. J Cardiovasc Electrophysiol. 2017;28:796–805.
Yoneyama K, Venkatesh BA, Bluemke DA, McClelland RL, Lima JAC. Cardiovascular magnetic resonance in an adult human population: serial observations from the multi-ethnic study of atherosclerosis. J Cardiovasc Magn Reson. 2017;19:52.
Imai M, Ambale Venkatesh B, Samiei S, Donekal S, Habibi M, Armstrong AC, et al. Multi-ethnic study of atherosclerosis: association between left atrial function using tissue tracking from cine MR imaging and myocardial fibrosis. Radiology. 2014;273:703–13.
Evin M, Broadhouse KM, Callaghan FM, McGrath RT, Glastras S, Kozor R, et al. Impact of obesity and epicardial fat on early left atrial dysfunction assessed by cardiac MRI strain analysis. Cardiovasc Diabetol. 2016;15:164.
Steinmetz M, Broder M, Hösch O, Lamata P, Kutty S, Kowallick JT, et al. Atrio-ventricular deformation and heart failure in Ebstein’s anomaly – a cardiovascular magnetic resonance study. Int J Cardiol. 2018;257:54–61.
Uemura T, Otsuji Y, Nakashiki K, Yoshifuku S, Maki Y, Yu B, et al. Papillary muscle dysfunction attenuates ischemic mitral regurgitation in patients with localized basal inferior left ventricular remodeling: insights from tissue Doppler strain imaging. J Am Coll Cardiol. 2005;46:113–9.
Koh AS, Gao F, Leng S, Kovalik JP, Zhao XD, Tan RS, et al. Dissecting clinical and metabolomics associations of left atrial phasic function by cardiac magnetic resonance feature tracking. Sci Rep. 2018;8:8138.
Zhong L, Tan LK, Finn CJ, Ghista D, Liew R, Ding ZP. Effects of age and gender on left atrial ejection force and volume from real-time three-dimensional echocardiography. Ann Acad Med Singapore. 2012;41:161–9.
The authors appreciate the support and medical editing assistance from Duke-NUS/SingHealth Academic Medicine Research Institute. The authors thank Dr. Heng Ge, MD, Dr. Jie He, MD (Shanghai Renji Hospital, China), and Ms. Han Yiying for their excellent technical support.
This study received funding support from the National Medical Research Council of Singapore (NMRC/EDG/1037/2011; NMRC/TA/0031/2015; NMRC/OFIRG/0018/2016; NMRC/BnB/0017/2015), Hong Leong Foundation and Edwards Lifesciences. The funder had no role in the design and conduct of the study; collection; management, analysis, and interpretation of the data; and preparation, review, or approval of the manuscript.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The protocol was approved by the SingHealth Centralised Institutional Review Board and informed consent was obtained from all participants.
Consent for publication
Written informed consent was obtained from all participants for inclusion of their data in publications.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Comparison of feature tracking derived left atrial strain (ε) and strain rate (SR) measurements among subject groups. (PDF 193 kb)
Figure S1. Linear relation between right ventricular ejection fraction (RVEF) and fast left atrial (A) reservoir strain, (B) conduit strain and (C) booster strain in patients with heart failure. (PDF 98 kb)