Introduction

The pathophysiological mechanisms related to the onset of dyspnea are not well known especially in elderly patients with subtle cardiac abnormalities [1]. Self-reported dyspnea has been associated with incident heart failure (HF) and mortality in the general population, representing the transition phase to overt HF [2, 3]. While dyspnea may be caused by a combination of comorbidities such as obesity, pulmonary disease, muscular disorders and cognitive impairment, dyspnea mainly represents a critical symptom of HF [3] in which left ventricular (LV) systolic and diastolic dysfunction are thought to be key contributors [4]. However, community-based studies have shown conflicting results regarding the association of LV dysfunction and symptoms, leading to underdiagnosis of HF especially among patients with preserved ejection fraction [5,6,7].

Emerging data show that left atrial (LA) structure and function may be independently associated with or may even precede HF onset in asymptomatic patients [8,9,10]. A “stiff LA syndrome” has been described as a consequence of a less compliant LA on exercise capacity despite preserved LV systolic and diastolic function [11]. Since LA plays a key role in the modulation of LV filling and function, a reduced atrial reserve may represent the first sign of a failing heart and thus mark the onset of dyspnea. Until now, majority of the studies have assessed the relationship between LA dysfunction and symptoms in patients with known cardiovascular diseases [12]. However, it remains unknown to what extent LA dysfunction may be associated with the onset of dyspnea in stage A and B HF patients.

Therefore, the aims of this study were: (1) characterizing non-specific symptoms, not referable to HF, in ambulatory patients, who were therefore classified as HF stages A or B and (2) testing the hypothesis that an impairment in LA structure and function may contribute to the presence of non-specific HF symptoms in such patients adding incremental information to the conventional non-invasive measures of cardiac structure and function and comorbidities.

Methods

Patients population

Consecutive patients referred to the echocardiographic laboratory of our Institute for cardiologic evaluation with clinical data available were screened. A recent ambulatory visit performed by a cardiologist was required during the prior 3 months or immediately after echocardiography. Patients included had to fulfil the following criteria: (1) HF Stage A or B defined as presence of conditions that are strongly associated with the development of HF with [stage B] or without [stage A] structural heart disease; (2) no previous history of clear signs or symptoms of HF; (3) absence of major criteria for HF at the time of the cardiologic evaluation according to Framingham criteria [13]; (4) absence of clinical diagnosis or suspicion of unstable angina; (5) echocardiographic evidence of preserved LV systolic function defined as LV ejection fraction (LVEF) > 50%; 6) sinus rhythm at the time of echocardiography. Patients with the following conditions were excluded: (1) unavailable measures of LA structure and function; (2) history of cardiac surgery or myocardial infarction in the previous 6 months; (3) any degree of mitral stenosis, rheumatic involvement of mitral valve, mitral valve prosthesis, moderate to severe aortic disease; (4) known pre-capillary pulmonary hypertension; (5) severe or active known lung disease (asthma exacerbation or advanced stage of chronic obstructive pulmonary disease [COPD]), who have required supplemental oxygen or previous hospitalization. The work has been carried out in accordance with the Declaration of Helsinki. Institutional review board approval was obtained for this project and written informed consent was obtained according to local regulation.

Clinical history and physical examination were recorded for each patient including age, sex, history of coronary artery disease (CAD), percutaneous or surgical revascularization, paroxysmal or persistent atrial fibrillation (AF), hypertension, hyperlipidemia, smoking status, COPD, asthma, chronic kidney disease (defined as GFRMDRD ≤ 60 mL/min/1.73 m2), diabetes. When available, a venous blood samples acquired during the previous 3 months was collected.

The symptomatic status was verified by reviewing clinical reports with particular attention to the presence of perceived breathlessness in daily activities referred by the patient in the absence of any sign of pulmonary congestion, which was considered as non-specific symptom by the cardiologist and justified by age, comorbidities or deconditioning and therefore unsufficient to classify the patient as evolved to the symptomatic phase of HF. Indeed, dyspnea on exertion is a minor criterion for HF according to Framingham criteria and, alone, it doesn’t classify the patient as symptomatic for HF [13].

2D echocardiography

All patients underwent a comprehensive echocardiographic evaluation according to the last Echocardiographic Guidelines of American Society of Echocardiography (ASE)/European Society of Cardiology (ESC) [14]. LV end-diastolic volume (EDV) and LV end-systolic volume (ESV) were measured using the biplane Simpson method and EF was calculated. LV mass was calculated according to the ASE convention using two-dimensional echocardiography from the parasternal long-axis view: LV wall mass (gr) = 0.8 × (1 0.04((LVIDd + IVDs + PWd) 3 − (LVIDd) 3)) + 0.6 [14]. LA volume was measured at end-systole in apical-four-chamber view, by tracing the LA inner border excluding the area under the mitral valve annulus and the inlet of the pulmonary veins. Measures were indexed by body surface area (EDVi and LAVi). Lateral and septal Doppler tissue imaging-derived systolic (S-TDI) and diastolic velocities were measured, and mean E/E' ratio was calculated. The presence of diastolic dysfunction was evaluated and graded according to recommended algorithms. In presence of mitral regurgitation (MR), effective regurgitant orifice area (ERO) was calculated by means of proximal iso-velocity surface area method [15]. Pulmonary artery systolic pressure (PASP) was obtained by the following formula as recommended: 4*(peak tricuspid regurgitation velocity)2 + right atrial pressure.

2D‑speckle tracking echocardiography

For LA echocardiographic speckle-tracking analysis, 2D grayscale images were acquired in the standard apical four-chamber view at a frame rate of 50 frames/s. Three consecutive heart cycles were recorded during a quiet breath hold. The analysis was performed offline using dedicated semi-automated acoustic tracking software (QLAB9; Philips). To calculate LA longitudinal strain, the QRS onset was used as reference point. After identifying the endocardial inner region of interest composed of six segments, tracking quality analysis and eventual manual adjustment were performed. Peak atrial longitudinal strain (PALS), which corresponds to the end of the reservoir phase, was calculated by averaging the peak values observed in all LA segments. Peak atrial contraction strain (PACS) corresponds to the positive peak at the onset of the p wave, at the end of the conduit phase in late diastole. Conduit strain (CS) was derived by the difference between PALS and PACS.

In patients in whom a maximum of two out of six segments were excluded because of inadequate tracking, strains values were calculated by averaging the values measured in the remaining segments. The reproducibility and feasibility of LA speckle-tracking measurement has been previously reported by studies conducted in our echocardiographic laboratory [16]. Patients were divided according to LA enlargement (≥ 35 mL/m2) and/or dysfunction. LA reservoir function (PALS) with a cut-off of 38% and LA booster function (PACS) with a cut-off of 16% were used as representative of LA dysfunction [17]. Four stages of atrial myopathy were defined: 0) normal LAVi, PALS and PACS; (1) normal LAVi, reduced PALS or PACS; (2) dilated LAVi, normal PALS and PACS or one of the two reduced; (3) dilated LAVi, reduced PALS and PACS.

Statistical analyses

Data are expressed as mean ± standard deviation, median (25th–75th percentiles) or frequencies (%). The Shapiro–Wilk test was performed to assess data normal distribution. Baseline characteristics were compared between asymptomatic and symptomatic patients, regardless of discrimination between HF stage A and B, using t-test or nonparametric test for continuous variables [18] and χ2 test for binary variables. Associations between symptomatic status and measures of LA structure (LAVi) and function (PALS, PACS and CS) were assessed with logistic regression analysis. After testing the unadjusted association, two models were constructed. Model 1 included age, sex, body mass index, hypertension, eGFR, history of CAD, total cholesterol, diabetes and COPD. Model 2 included S-TDI, E/E', LV mass-i and MR-ERO. A third model was constructed (Model 3), which included all the covariates considered in Model 1 and in Model 2. The covariates included in the models were a selection of variables significantly associated with symptomatic status (Tables 1 and 2) and plausible demographic and clinical confounders with symptoms. Receiver operating characteristic (ROC) curves were generated to assess the performance of individual LA structural (LAVi) and functional parameters (PALS, PACS and CS) to predict the symptomatic status. Furthermore, the presence of enlarged LA and/or reduced LA function were considered as expression of atrial myopathy and patients were divided in the previously described four stages. To assess the incremental value of atrial myopathy in symptoms prediction in comparison to clinical or echocardiographic parameters alone, ROC curves analyses and DeLong test were performed. Three covariates were selected from each of the two previously described models (Model 1 and Model 2) based on statistical significance at t-test and χ2 test. Two ROC curves analyses were performed, which included in addition to atrial myopathy, age, hypertension and COPD; LV mass-i, E/E′ and mitral ERO, respectively. All tests were two-tailed. A p value < 0.05 was considered to indicate statistical significance. Analyses were performed using SPSS version 20.0 (SPSS, Chicago, IL).

Table 1 Clinical characteristics of the study population stratified by symptomatic status
Table 2 Echocardiographic characteristics of the study population stratified by symptomatic status

Results

The cohort consisted of 185 patients [mean age 63 ± 16 years, 47% women] of whom 64 (35%) were classified as stage A HF and 121 (65%) as stage B HF. One-hundred-nine patients (59%) had hypertension, 48 (26%) diabetes and 20 (11%) history of AF. Forty patients (22%) were routinely followed for CAD, 28 (15%) for heart valve disease and 20 (11%) for exclusion of drug-induced cardiotoxicity. Ninety-six patients (57%) were evaluated due to a specific murmur, palpitation, or cardiovascular risk factors.

Overall, mean EDVi was 63 ± 15 mL, EF 62 ± 6% and S-TDI 8.0 ± 1.7 cm/s. LA was enlarged in 86 patients (46%). Mean E/E′ was 9.7 ± 3.8; E/E′ was > 14 in 16 patients (9%). Diastolic dysfunction was evaluated and graded according to recommended algorithms [19] and shown in Supplementary Table 1. Mean PALS was 33 ± 11%, PACS 15 ± 6% and CS 18 ± 7%. MR was present in 84 patients (45%) with a mean ERO of 0.07 ± 0.12 cm2.

Tables 1 and 2 show clinical and echocardiographic characteristics of patients according to presence of symptoms. One-hundred-thirty-three patients (72%) were asymptomatic and 52 (28%) referred breathlessness, classified as non-specific HF symptom. Symptomatic patients were more likely to be older, to suffer from hypertension, chronic kidney disease, to have higher levels of total and low-density lipoprotein cholesterol and to have history of non-severe COPD or asthma (Table 1). Overall measures of cardiac structure and function were impaired in the symptomatic group (Table 2). In particular, patients with symptoms had lower S-TDI, higher LV-mass, E/E′, MR-ERO, PASP but not lower TAPSE.

Measures of LA structure, LAVi, and function, PALS, PACS and CS, were significantly different between the two groups (Table 2; Fig. 1). A significant association between each measure of LA structure and function and symptomatic status was observed (Table 3). LAVi, PALS and PACS, but not CS, remained significantly associated with symptomatic status after adjustment for clinical and echocardiographic confounders (Table 3). Moreover, only LAVi, PALS and PACS, among all echocardiographic variables, had independent association with symptoms after comprehensive adjustment for clinical variables included in Model 1 (Table 2). LAVi was the most accurate in predicting symptoms (AUC:0.81) followed by PALS (AUC:0.75), E/E′ (AUC: 0.73) and PACS (AUC: 0.71).

Fig. 1
figure 1

Measures of left atrial (LA) structure (left atrial volume index [LAVi]) and function (peak atrial-longitudinal strain [PALS], -contraction strain [PACS] and—conduit strain [CS]) in the asymptomatic and the symptomatic group (between groups: p < 0.0001 for all). Light grey boxes represent asymptomatic patients, dark grey boxes represent symptomatic patients

Table 3 Association between Symptomatic Status and LA structure and function

Mean values of LA strain and LAVi according to normal/abnormal LAVi or PALS are reported in Supplementary Table 2. In the overall population, 46 patients (25%) had normal LAVi and normal LA function (PALS and PACS) (stage 0); 53 (29%) had reduced PALS or PACS but normal LAVi (stage 1); 28 (15%) had dilated LAVi with preserved function (PALS and PACS) or dilated LAVi with impairment in PALS or PACS (stage 2); 58 (31%) showed both structural and functional (PALS and PACS) impairment (stage 3). On comparison among the four groups, symptoms were more likely to be present in advanced stage of atrial myopathy (p for trend < 0.0001) (Supplementary Table 3; Fig. 2; Supplementary Fig. 1). The incremental value of the presence of atrial myopathy in symptoms prediction on top of clinical and echocardiographic confounders was demonstrated using ROC curves analyses (Fig. 3a, b). The clinical model included age, hypertension and COPD; the echocardiographic model included E/E′, LV mass-i and mitral ERO. The further addition of atrial myopathy added predictive value to the models with AUC increase from 0.80 to 0.88 (p = 0.004) and from 0.79 to 0.84 (p = 0.06), respectively (Fig. 3a, b).

Fig. 2
figure 2

Prevalence of non-specific HF symptoms in the four stages of atrial myopathy

Fig. 3
figure 3

Incremental predictive value of LA myopathy. Atrial myopathy was graded as abnormal left atrial volume (LAVi) and/or peak longitudinal strain (PALS) and/or contraction strain (PACS). Notably, the addition of atrial myopathy to clinical (age, Hypertension, COPD) (a) or echocardiographic confounders (E/E′, LV mass-i, mitral effective regurgitant orifice area) (b) significantly improved the predictive value for symptoms (AUC increases from 0.80 to 0.88, p = 0.004, and from 0.79 to 0.84, p = 0.06, respectively)

Results remained similar in the subgroup of patients without history of AF (n = 165) [Supplementary Table 4].

Discussion

We report several key findings in this study. First, in patients with HF stages A–B and preserved LVEF, measures of LA structure and function are the only echocardiographic parameters associated with non-specific referred symptoms independently of non-cardiac comorbidities. Second, LA measures are also significantly associated with symptoms after accounting for echocardiographic parameters that are strongly linked with the development of HF such as LV mass, indices of systolic and diastolic function, quantitatively graded mitral regurgitation. Third, LA myopathy evaluation through structural (LAVi) and functional parameters (PALS and PACS) provides incremental information in identifying symptomatic status of patients after assessment of comorbidities and structural pre-requisites for the development of HF (Fig. 4).

Fig. 4
figure 4

Central role of left atrial myopathy in the transition from the asymptomatic (stages A–B) to the symptomatic phase (stage C) of heart failure with preserved left ventricular ejection fraction (HF-pEF)

We found that measures of LA remodeling are significantly associated with non-specific symptoms showing a significant discrimination accuracy compared with other cardiac measures. The LA is commonly considered a buffer chamber between the pulmonary circulation and the LV and its changes are indirectly thought to be related to LV function [10]. Chronic exposure to high LV filling pressure initiates an adaptive process leading to LA enlargement and dysfunction [20]. However, our data showed that in patients with preserved LVEF, LA impairment was predictive of symptoms regardless of age, sex and the common measures of LV structure and systolic/diastolic function (i.e. LV hypertrophy, S-TDI, E/E`). Thus, it is possible that, especially in the context of preserved EF, underlying LA functional impairment per se may contribute to the onset of symptoms and should not be considered just an innocent bystander [21].

LA modulates LV filling and function by the reservoir phase, when LA chamber collects blood coming from the pulmonary veins, the conduit phase, when blood passively transfers to the LV in early diastole and the contraction phase when blood is forced to fill LV at end-diastole. LA also acts as interplay between the LV and the pulmonary circulation, buffering pressure and flow oscillations during the cyclic cardiac work [11, 22]. The loss of global atrial function may directly affect LV output, especially during exercise [10]. Also, reduced LA wall compliance, which influences atrial reservoir function, results in elevated LA pressure which consequently can lead to increased pulmonary venous and arterial pressure and thus symptoms [23]. From this perspective, assessment of LA function by speckle-tracking-echocardiography may enhance symptoms identification and improve diagnostic accuracy in the ambulatory settings, especially in subtle cardiac structural and functional abnormalities. Speckle-tracking-echocardiography provides a direct measurement of intrinsic LA myocardial deformation and, although not load independent, LA strain seems to be less dependent on loading conditions and geometric assumptions than traditional parameters with the further advantage to show good feasibility and reproducibility [24, 25]. Moreover, it correlates well with the burden of myocardial fibrosis, i.e. histological marker of chamber stiffness [26, 27].

A “stiff LA syndrome” has been described [11] as clinical consequence of LA dysfunction in the presence of preserved LV function. The clinical usefulness to assess symptoms has been previously investigated in homogeneous study population affected by HFpEF, hypertrophic cardiomyopathy or heart valve disease [28, 29]. However, the coexistence of other conditions, which are hierarchically considered the source of symptoms, make it challenging to understand the additional role of LA. It is possible that the LA chronic exposure to hemodynamic, neuro-endocrinal and metabolic stressors such as those characterizing hypertension, diabetes and CAD favor structural and functional changes that may determine symptoms onset even before macroscopic cardiac structural changes. This concept is in line with the description of a specific “atrial cardiomyopathy” caused primarily by an atrial remodeling process related to external etiological factors directly affecting the LA [30]. From this perspective, atrial cardiomyopathy per se may involved in the risk of stroke before the onset of AF, and as well it may be responsible of neuro-hormonal activation regardless of the hemodynamic status [31]. The pathological changes related to these mechanisms may reflect a spectrum of atrial tissue structural alterations leading to atrial impairment, including myocyte hypertrophy, necrosis, apoptosis and changes in extracellular matrix composition [10, 20] favoring atrial stiffness.

Our results have important clinical implications since it has been recognized that patients with systolic dysfunction and isolated diastolic dysfunction may remain free of HF symptoms [4, 5] and it is still not clear which are the determinants of the transition from the asymptomatic to the symptomatic phase. Moreover, such discrepancy between detected cardiac abnormalities and unrecognized HF symptoms contributes to misdiagnosis of HF [32,33,34] especially in the context of preserved LVEF. We demonstrated that the presence of atrial myopathy is associated with non-specific symptoms, and it could contribute to their interpretation adding incremental information to the conventional non-invasive measures of cardiac structure and function and comorbidities. The presence of atrial myopathy should raise the suspicion for overt HF and therefore suggest further functional evaluation.

There are several limitations in this study. First, this was a single-center study with a limited number of patients. Second, an erroneous classification of patients as asymptomatic or symptomatic for HF or the presence of non-cardiac sources of symptoms cannot be excluded. Third, measures of NT-pro-BNP plasma levels were not available. Nevertheless, our results suggested the utility of LA measures in an ambulatory setting whereas NT-pro-BNP has been shown to be a marker associated with HF symptoms especially in the acute setting thanks to its high negative predictive value. Of note, it is still not known whether LA functional parameters diagnostic value could apply similarly to natriuretic peptides, i.e. normal LA strain could be used to exclude HF related symptoms. Large cohost studies are needed to investigate LA function thresholds for this purpose.

The echocardiographic estimation of PASP was not included in the analyses since pulmonary hemodynamic alterations were considered as consequence of left chambers remodeling and dysfunction. However, the association between symptoms and LA structural and functional impairment was not affected by introducing PASP values in the analyses (data not shown). Fourth, the diastolic dysfunction grading was not considered in the analyses, since LAVi is one of the four variables considered in diastolic dysfunction stratification [19] and therefore, diastolic dysfunction assessment is strongly dependent on left atrial remodeling. Moreover, the four recommended parameters did not receive an extensive validation in patients with HFpEF and only a modest correlation with invasive filling pressure has been described [35]. Our data showed that LA myopathy was associated with symptomatic status with at least the same predictive value of LA pressure estimated by means of recommended algorithms. Therefore, it could be speculated that the predictability of diastolic dysfunction grading is strongly guided by LA remodeling.

Conclusions

In ambulatory patients with HF stages A−B and preserved LVEF, measures of LA structure and function are significantly associated with non-specific symptoms, even after accounting for risk factors, comorbidities and measures of cardiac structure and systolic/diastolic function. LA myopathy assessment through LA dilatation and/or reservoir and booster function impairment may add incremental value in the recognition of mild symptoms, which represent the transition from the asymptomatic (stage B) to the symptomatic (stage C) phase of HF. Therefore, a comprehensive LA remodeling evaluation may help clinicians in the appropriate identification of overt HF.