Prognostic value of temporal patterns of left atrial reservoir strain in patients with heart failure with reduced ejection fraction

Background We investigated whether repeatedly measured left atrial reservoir strain (LASr) in heart failure with reduced ejection fraction (HFrEF) patients provides incremental prognostic value over a single baseline LASr value, and whether temporal patterns of LASr provide incremental prognostic value over temporal patterns of other echocardiographic markers and NT-proBNP. Methods In this prospective observational study, 153 patients underwent 6-monthly echocardiography, during a median follow-up of 2.5 years. Speckle tracking echocardiography was used to measure LASr. Hazard ratios (HRs) were calculated for LASr from Cox models (baseline) and joint models (repeated measurements). The primary endpoint (PEP) comprised HF hospitalization, left ventricular assist device, heart transplantation, and cardiovascular death. Results Mean age was 58 ± 11 years, 76% were men, 82% were in NYHA class I/II, mean LASr was 20.9% ± 11.3%, and mean LVEF was 29% ± 10%. PEP was reached by 50 patients. Baseline and repeated measurements of LASr (HR per SD change (95% CI) 0.20 (0.10–0.41) and (0.13 (0.10–0.29), respectively) were both significantly associated with the PEP, independent of both baseline and repeated measurements of other echo-parameters and NT-proBNP. Although LASr was persistently lower over time in patients with PEP, temporal trajectories did not diverge in patients with versus without the PEP as the PEP approached. Conclusion LASr was associated with adverse events in HFrEF patients, independent of baseline and repeated other echo-parameters and NT-proBNP. Temporal trajectories of LASr showed decreased but stable values in patients with the PEP, and do not provide incremental prognostic value for clinical practice compared to single measurements of LASr. Graphical abstract Supplementary Information The online version contains supplementary material available at 10.1007/s00392-023-02244-x.


Conclusions
Temporal paƩerns of LA reservoir strain (LASr) are associated with the primary endpoint, independent of GLS, E/A, E/e', LAVI and NT-proBNP.A single measurement of LASr provides sufficient informaƟon for determining prognosis in clinical pracƟce compared to repeated measurements.

Introduction
Most of the contemporary risk scores for heart failure with reduced ejection fraction (HFrEF) focus on systolic echocardiographic determinants, while the influence of diastolic determinants on prognosis has been studied less extensively [1].Categorization of HFrEF patients based on diastolic determinants is mainly used to non-invasively estimate left atrial pressure (LAP), which can be useful to guide medical treatment and provide information on prognosis [2,3].However, the algorithm currently in use for estimating LAP carries an important limitation; it requires multiple parameters that are often affected by cardiac rhythm and/or mitral valve disease, with the consequence that a substantial part of HFrEF patients remain uncategorized [4].
Recently, there has been an emerging interest in the use of left atrial reservoir strain (LASr) as a measure of left atrial (LA) function and as a derived measure for LAP in HFrEF patients [5][6][7][8].LASr is predominantly determined by LV GLS and filling pressure [7].Although LA strain could have limitations like load dependency, Doppler echocardiography may carry more limitations.It is angle-dependent, which is not an issue in strain analysis [9].In addition, LA strain is preload dependent, but to a lesser degree than LA volume [10].Studies have demonstrated that LASr is decreased in HFrEF patients and that an abnormal LASr is associated with increased LAP [7,8].Additionally, studies have shown that LASr carries prognostic information in patients with atrial fibrillation and mitral valve disease.This implies that LASr is less affected by these conditions as compared to conventional echocardiographic parameters, to estimate LAP [11,12].Only a few studies have demonstrated that LASr may have prognostic value in HFrEF [13][14][15].These studies only examined single measurements of LASr, which merely represent a snapshot of a patient's condition, and related these measurements to clinical endpoints occurring over several years thereafter.The prognostic value of repeatedly measured LASr has never been examined before and has never been compared to that of other echocardiographic parameters in chronic HFrEF patients.
Therefore, we investigated whether repeatedly measured LASr provides incremental prognostic value over a single baseline LASr value in stable chronic HFrEF patients.In addition, we hypothesized that temporal patterns of LASr are associated with adverse clinical events, and that temporal patterns of LASr may provide incremental prognostic value to temporal patterns of other prognostic markers.

Study design
The design of the Serial Biomarker Measurements and New Echocardiographic Techniques in Chronic Heart Failure Patients Result in Tailored Prediction of Prognosis (Bio-SHiFT) study has previously been described [16].Bio-SHiFT is a prospective, observational study of stable patients with chronic heart failure (CHF), conducted at the Erasmus MC, Rotterdam, and Northwest clinics, Alkmaar, The Netherlands.Recruitment was conducted during the patient's regular outpatient visits while in clinically stable condition (i.e., they had not been hospitalized for HF in the 3 months prior to inclusion).The main inclusion criteria were diagnosis of HF according to the then prevailing guidelines of the European Society of Cardiology 3 or more months before inclusion and age ≥ 18 years [17].Patients with an atrial pacemaker were excluded from the current analysis.Patients were observed for a maximum of 30 months, with follow-up visits scheduled every 3 months.A brief medical examination and blood samples were taken at each visit.All patients' usual outpatient follow-up with their treating physician continued throughout the study, independently of the study visits.This study was approved by the medical ethics committees, conducted in accordance with the Declaration of Helsinki, and registered in Clinical-Trials.gov(NCT01851538).Informed consent was obtained from all patients.In total, 398 patients were included in Bio-SHiFT.The repeated echo study that we currently report was performed at the Erasmus MC only and consisted of 175 HFrEF patients with echocardiographic assessment every 6 months during follow-up [18].Two patients had insufficient image quality, and therefore the remaining 173 patients were included in the current study.

Echocardiography measurements and evaluation
Two-dimensional gray-scale harmonic images were obtained in the left lateral decubitus position.Standard apical four-, three-, and two-chamber views were recorded.A commercially available ultrasound system was used (iE33, Philips, Best, The Netherlands), equipped with a broadband (1-5 MHz) S5-1 transducer (frequency transmitted 1.7 MHz, received 3.4 MHz).Images were stored in the echo core lab of Erasmus MC.Using specialized software (2D Cardiac Performance Analysis version 4.5; TomTec Imaging Systems, Unterschleissheim, Germany), LVEF, tricuspid regurgitation (TR) velocity, and the function of the mitral valves were assessed according to the then prevailing guidelines [4].The diastolic parameters were evaluated using Philips Excellera version R4.1 (Philips Medical Systems, The Netherlands) or TomTec Imaging Systems.To assess diastolic function indexed left atrial volume (LAVI), the peak early filling velocity (E)/late filling velocity (A) ratio (E/A ratio) and the ratio of the E and early diastolic mitral annular velocity (e') (E/é ratio) were calculated [4].For the e', we used the mean of the lateral and medial e' when available; however, if only one of the two was available, this value was used.All echocardiographic measurements were performed blinded to biomarker and clinical event data [16].
Strain parameters were measured with speckle tracking echocardiography (also using TomTec Imaging Systems) by a single operator, according to the prevailing consensus document of the EACVI/ASE/Industry Task Force to standardize deformation imaging [17].The apical 4-chamber view was used preferably for the analysis.LA endocardial borders were automatically traced.The LA is contoured extrapolating across the pulmonary veins and LA appendage orifice [19].We used end-diastole as a reference.Fine-tuning was performed manually if the tracking was suboptimal.If the quality of the 4-chamber view was poor, the 2-chamber view was used.Patients with insufficient image quality to perform LA strain analysis or patients with an atrial pacemaker were excluded.LA strain was assessed according to the three phases of the LA cycle: LA reservoir strain (LASr) which starts at the end of ventricular diastole (mitral valve closure) and continues until mitral valve opening (20).An example of a LA strain curve is provided in Fig. 1.Global longitudinal strain (GLS) was assessed in 18 LV segments on the standard apical four-, three-, and two-chamber views, where the endocardial border was traced manually at end-systole.The Fig. 1 Example of LA strain measurement.LASr left atrial reservoir strain, LAScd left atrial conduit strain, LASct left atrial contractile strain mean GLS from the three apical views was considered the LV GLS.If a patient had AF during the echocardiography, the index beat method was used.This is a validated method to measure echocardiographic parameters during AF [21].
Patients underwent echocardiographic assessment at baseline and every 6 months during follow-up.Due to logistic reasons, 58% of the first available echoes were performed at baseline (follow-up time zero), 15% of the first available echoes were performed within 3 months after the start of the study, 18% within 6 months after the start of the study, and the remaining 9% thereafter (Fig. 2).Missing echocardiograms occurred due to logistic circumstances (e.g., the unavailability of an ultrasound technician during the study visit).Intra-observer reproducibility was assessed by remeasuring GLS in 20 echocardiograms and calculating the intraclass correlation coefficient.The intraclass correlation coefficient was 0.93 for LASr.

Clinical study endpoints
The primary endpoint (PEP) comprised the composite of hospitalization for the management of acute or worsened HF, left ventricular assist device (LVAD) implantation, cardiac transplantation, and cardiovascular death, whichever occurred first in time.All events were adjudicated by a clinical event committee blinded to the echocardiographic assessments and biomarker measurements, after reviewing corresponding hospital records and discharge letters [18].

Statistical analyses
Distributions of continuous variables were tested for normality using the Shapiro-Wilk test.Normally distributed continuous variables are presented as mean ± standard deviation (SD), and non-normally distributed variables as a median and 25th-75th percentile.Categorical variables are presented as numbers and percentages.Differences in baseline characteristics between patients who experienced the PEP and those who did not were tested using the t-test and Mann-Whitney test, according to variable distributions, for continuous variables.For categorical variables, χ 2 -tests and Fisher's exact tests were used.
First, we examined single measurements of LASr and other echo parameters of interest in relation to the PEP using Cox models (only the first available echo was used), and we adjusted for age, sex, duration of HF and N-terminal pro-b-type natriuretic peptide (NT-proBNP).In addition, we calculated the Pearson correlation coefficients for the echo variables of interest to assess the correlation with LASr.
Then, we assessed the value of repeated strain measurements for prediction of the PEP, as well as their incremental value to sole, baseline measurements.For this purpose, joint models for longitudinal and survival data were used [22].In these joint models, a linear mixed effects (longitudinal) model provided estimates of the individual temporal trajectories of the echo parameters, and in combination with a relative risk model, the association of the trajectories with the risk of the PEP was estimated.The associations between the temporal evolutions of LASr and the PEP, resulting from the relative risk model, were first only adjusted for age, sex and duration of HF (model 1).Thereafter, baseline NT-proBNP (model 2), GLS and LVEF (model 3) were added consecutively.Furthermore, we adjusted for the diastolic parameters (E/A ratio, E/e ratio, LAVI) in model 4. Lastly, all variables with significant differences between those with and without the PEP were added (model 5).To investigate the incremental value of repeatedly measured LASr to repeatedly measured echo parameters and NT-proBNP, multivariable joint models were used.
We performed the above-described analyses on the full cohort.Subsequently, because LASr could be affected by the presence of AF and mitral regurgitation (MR), we performed subgroup analyses in patients with and without AF and MR.
To enable comparisons of effect sizes of different variables, we calculated the Z-scores for all investigated echo parameters, and NT-proBNP.Hazard ratios were obtained from both the Cox as the joint models.Thus, the results of the regression analyses of the Cox and joint models can be directly compared and are presented as HRs, which represent Fig. 2 Study design: first available and follow-up echocardiograms.This figure provides 3 example patients to illustrate which echocardiograms were the first available echocardiograms, considered as 'baseline' in the analysis (red circles), and at which time-points follow-up echocardiograms were scheduled (white circles).55% of the first available echocardiograms were performed at baseline (followup time zero), 12.8% were performed during the first study follow-up visit (target follow-up time 3 months) and 18% were performed during the second follow-up visit (target 6 months).Subsequently, echoes were performed every six months a risk per Z-score unit, along with the corresponding 95% confidence interval (CI).
Missing values in LASr and the other echo parameters (besides the A wave) were due to poor image quality and were therefore considered missing completely at random.Accordingly, we chose to perform a complete case analysis.Missing values for the A wave were due to atrial fibrillation during the echo or due to mitral valve replacement or clipping.In these patients the A wave can never be measured, thus imputation of missing values is inappropriate.Therefore, we again chose a complete case analysis here.
All analyses were performed with R Statistical Software using packages nlme [23] and JMbayes [22].All tests were two-tailed, and p values < 0.05 were considered statistically significant.

Baseline characteristics
Between October 2011 to January 2018, 173 patients were included in the Bio-SHiFT echo study.Twenty patients had an atrial pacemaker and were therefore excluded from the current analysis.Of the remaining 153 patients, 76% of the patients were male, mean age was 58 ± 11 years, and mean BMI was 27.5 ± 4.6 kg/m 2 .A total of 27% were in NYHA class I, and 55% were in NYHA class II.Ischemic heart disease was the most prevalent HF etiology (44%).The median time between diagnosis of HF and inclusion in the study was 6.5 (6.1-7.3)years.During a median follow-up time of 2.5 (2.3-2.6)years, a total of 50 patients (33%) reached the PEP, out of whom 37 were re-hospitalized for acute or worsened HF, 6 patients received a heart transplantation, 4 patients received an LVAD, and 3 patients died from a cardiovascular cause.Patients who reached the composite PEP had lower systolic and diastolic blood pressure, had a higher NT-proBNP (303 [180-540] vs. 71  pmol/L, p < 0.001), and comorbidities such as atrial fibrillation and renal failure were more prevalent in this group (resp.46% vs. 21%, p = 0.009; 58% vs. 30%, p = 0.003) (Table 1).

Echocardiographic characteristics
During a median follow-time of 2.5 years, 410 echocardiograms were performed with a median of 3 (2-4) per patient.Patients had up to 8 consecutive echocardiographic evaluations performed, with 65% having at least 3 evaluations.An overview of the characteristics of the first available echocardiogram for each patient is presented in Table 2.

Subgroup analyses in patients with atrial fibrillation and in patients with mitral regurgitation
The subgroup analyses showed that the repeated measurements of LASr were still significantly associated with the endpoint both in patients with SR as patients with AF, with an HR (95% CI) of 0.15 (0.06-0.34) and 0.03 (0.01-0.27), respectively (corrected for age, sex, duration of HF).This was also the case in patients with MR compared to patients with no MR, with HRs (95% CI) of 0.09 (0.03-0.22) and 0.14 (0.06-0.33), respectively.The full results are presented in Supplementary Table 1.

Temporal evolution of LASr
In the total population, there was a decrease in LASr over time as the PEP or censoring approached (Beta [95% CI] − 1.72 [− 2.46 to − 0.98]) per LASr (%) change per year).Supplementary Fig. 2 and Fig. 3 show the temporal evolution of patients who experienced the PEP and those who did not.Although, as described above, repeatedly measured LASr was associated with the occurrence of the PEP, and the average LASr was lower in patients who experienced the PEP compared to those who did not, this difference remained stable over time.LASr did not diverge further between patients with vs. without the PEP, as the PEP or censoring approached.

Discussion
We demonstrated that during a median follow-up of 2.5 years, repeated measurements of LASr were significantly associated with adverse cardiovascular events in HFrEF patients, independent of repeatedly measured GLS, LAVI, E/A ratio, E/e' ratio, and NT-proBNP.LASr was a stronger predictor than GLS, LAVI, E/A ratio, and E/e' ratio.Although repeated measurements of LASr were associated with the primary outcome, the difference in LASr remained stable over time, and temporal LASr evolutions did not further diverge in patients with events versus those without events.Therefore, for clinical purposes, repeated measurements of LASr do not seem to provide additional value over single measurements over a time frame of 2.5 years.To our knowledge, this is the first study that investigated the prognostic value of repeated LASr measurements in HFrEF.

Echocardiographic determinants of survival in patients with HFrEF
Although significant improvements in HF therapy have been made in the last two decades, the mortality and morbidity due to HF remain substantial [24].Numerous multivariable risk models have been proposed to identify patients with a poor prognosis, but the usefulness of these models in clinical practice has been limited [24].Most of these risk scores incorporate parameters of systolic function.Yet, the assessment of diastolic function may be equally important, as diastolic parameters provide a non-invasive estimation of LAP [4].However, in a substantial part of HFrEF patients, guideline-based estimation of LAP is not possible, as crucial parameters are often affected by mitral regurgitation (MR) and/or atrial fibrillation (AF) [4].In addition, some of the conventional parameters that are used for diastolic function in HFrEF have several limitations.For instance, LAVI, which is widely used as an indicator of LAP, does not always provide an accurate estimation of LAP, as LAVI can be increased in the presence of normal filling pressures (e.g.healthy athletes).A study by Benfari et.al. demonstrated that E/e' ratio outperformed other diastolic parameters as a prognosticator in HFrEF patients, but LA strain was not included in their study [25].The use of LASr as a non-invasive estimate for LAP has recently gained more interest, as LASr has been shown to be correlated with invasively measured LV filling pressures and LASr showed better prognostic performance than LAVI and E/e' ratio [7,[13][14][15].Our results are in line with these previous studies, as we show that LASr is a strong predictor of adverse cardiovascular events; and a stronger predictor than other echocardiographic parameters, such as E/e'.In contrast, in our study LAVI was not associated with an increased risk of PEP, which is similar to a study by Modin et al. [26].This could in part be explained by the fact that LASr is a more sensitive parameter than a volumetric parameter such as LAVI, and that an impairment in LA function is detected earlier than changes in LA volume [27].
Our study confirms and extends previous evidence on the added prognostic value of LASr in HFrEF.A few previous studies have investigated the prognostic value of single measurements of LASr in HFrEF patients [13][14][15].In a study consisting of 405 patients with an LVEF < 40%, LASr strongly predicted adverse outcomes, independent of other clinical and echocardiographic predictors of prognosis [13].A study by Malagoli et al. showed similar results; patients with lower LASr showed worse event-free survival than those with higher LASr [14].In acute HF, LASr was also shown to be a significant prognosticator [15].These studies only examined baseline measurements of LASr.Our study is the first to investigate the prognostic value of repeatedly measured LASr, and its added value over a single baseline LASr assessment, and over repeated measurements of other echocardiographic variables.We showed that repeated measurements of LASr were associated with the PEP, and the  Each dot represents a single measurement association persisted after consecutively adding repeated measurements of GLS, LAVI, E/A ratio, E/e' ratio, and NTpro-BNP.However, the difference in LASr between patients with events versus those without events remained stable over time, and temporal LASr evolutions did not further diverge as the PEP or censoring approached.A factor that may in part have contributed to these findings, is that in advanced stages of HF further reduction of already low values of LASr is less likely.

Results in the context of LA physiology and function in HFrEF
The LA plays a pivotal role in the filling of the LV and contributes to the cardiac output as the LA interacts with the LV and the pulmonary veins.The LA cycle is composed of three phases, which reflect the three main LA functions, reservoir, conduit, and contractile function [20].A recent meta-analysis found a normal value of > 39% for LASr in healthy individuals [28,29].Mean LASr in our population of HFrEF patients was 20.9%.Therefore, profound LA dysfunction exists in our cohort of HFrEF patients, which is in line with previous literature [28].
In the cardiac cycle, LASr and GLS are tightly coupled, as a maximal expansion of the LA takes place during LV systole.This is supported by the observation that LASr and GLS are significantly correlated in HFrEF [7].Our results confirm this, as a more advanced impairment of GLS was significantly correlated with an impairment in LASr.Previously, we have demonstrated that baseline and repeated measurements of GLS provide incremental prognostic value over LVEF [30].In the current investigation, we observed that LASr outperforms GLS as a prognostic marker in chronic HFrEF patients.This finding is in line with previous studies that have shown that LASr was superior in predicting outcomes compared to GLS [31].A potential explanation is that LASr might be affected by atrial inflammation and atrial fibrosis, which restricts atrial stretching, independent of LV longitudinal contraction and subsequent impairment of GLS.Our study is the first to report that repeated measurements of LASr were associated with clinical outcomes, independent of repeated measurements of GLS.Our results extend and add to previous studies and underline that LASr has more value as a prognostic marker in clinical practice than GLS, as well as other known prognostic markers, in stable patients with chronic HFrEF.

Study limitations
Several limitations of our study should be noted.First, treating physicians were not blinded to the conventional parameters assessed by echocardiography and therefore echocardiographic characteristics might have influenced treatment.
However, LASr values were not available for the treating physicians as these were assessed retrospectively.Secondly, the sample size of the study was modest and so was the number of endpoints.To prevent overfitting, we fitted multiple multivariable models containing different confounders.In addition, we adjusted for the duration of HF at baseline, to control for possible lead-time or length-time bias.Lastly, our cohort consisted of patients who were relatively young and in NYHA class I and II.Our results can therefore not be extrapolated to older patients in a more advanced stage of HF.

Conclusion
Repeatedly measured LASr was significantly associated with adverse cardiovascular events in patients with HFrEF.However, although the temporal trajectories of LASr were different in patients who reached the PEP compared to those who did not, they did not diverge as the PEP or censoring approached, and therefore repeatedly measuring LASr does not seem to provide additional incremental prognostic information over a single baseline measurement over a median follow-up time of 2.5 years.A single measurement of LASr showed a stronger prognostic value than conventional echocardiographic parameters.Therefore, LASr should be considered for routine use in clinical practice in patients with HFrEF, for prognostication and potentially for guiding treatment.
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Fig. 3
Fig. 3 Mean temporal patterns of LA reservoir strain until occurrence of the primary endpoint or censoring.LASr left atrial reservoir strain; Continuous lines represent mean temporal patterns for patients with the PEP (red) and patients who remained PEP-free (blue), as extracted from the joint model.Time-point zero represents the occurrence of an event in the PEP patients and censoring in patients who remained PEP-free.Dotted lines represent 95% confidence intervals.Each dot represents a single measurement

Table 1
Baseline patient characteristics of the total study populationNormally distributed data are presented as mean ± SD, non-normally distributed data are presented as median (25th-75th percentile).p values represent an overall comparison between PEP and no PEP PEP primary endpoint, HF heart failure, NYHA new york heart association, PCI percutaneous coronary intervention, CABG coronary artery bypass graft, COPD chronic obstructive pulmonary disease, ACE angiotensine converting enzyme

Table 2
Echocardiographic characteristics of the first available echo in relation to the composite endpoint Normally distributed data are presented as mean ± sd, non-normally distributed data are presented as median (25th-75th percentile).p values represent overall comparison between PEP and no PEP PEP primary endpoint, LASr left atrial reservoir strain, LAScd left atrial conduit strain, LASct left atrial contractile strain, LV GLS left ventricular global longitudinal strain, LVEF left ventricular ejection fraction, E/A ratio the ratio of the peak early left ventricular filling velocity over the late filling velocity, E/e' ratio E to early diastolic mitral annular tissue velocity, TR tricuspid regurgitation

Table 3
Associations of baseline and repeatedly measured LASr with the primary endpointLASr left atrial reservoir strain, GLS global longitudinal strain, LAVI left atrial volume indexed, E/A ratio the ratio of the peak early left ventricular filling velocity over the late filling velocity E/e' ratio E to early diastolic mitral annular tissue velocity *Corrected for age, sex, duration of HF, baseline NT-proBNP † Multivariable Joint Models: Corrected for age, sex, duration of HF, atrial fibrillation, renal failure, systolic and diastolic blood pressure