Skip to main content

Left atrial morpho-functional remodeling in atrial fibrillation assessed by three dimensional speckle tracking echocardiography and its value in atrial fibrillation screening



Three dimensional speckle tracking echocardiography (3D STE) is a novel technique combining 3D echocardiography and speckle tracking analysis. 3D STE software dedicated to the left atrium (LA) was recently available. Our study aimed to assess (1) atrial fibrillation (AF) related LA morpho-functional remodeling using 3D STE and (2) value of LA function parameters in identifying paroxysmal AF (PAF).


One hundred thirty-nine PAF, 109 persistent AF (Per-AF) and 59 non-AF subjects underwent 3D STE. LA phasic volumes and total LA emptying fraction (LAEF) were obtained and used to calculate passive (pLAEF) and active LA emptying fraction (aLAEF) based on atrial contraction. LA longitudinal and circumferential strain representing reservoir (LASr/LASrc), conduit (LAScd/LAScdc) and pump (LASct/LASctc) function were also assessed.


3D STE was found to have good reproducibility. Increase of LA volumes and decrease of parameters representing LA reservoir and pump function were independently associated with AF as well as AF burden. The correlations between LA emptying fraction and LA circumferential strain representing the same function were always stronger than those with LA longitudinal strain (p < 0.001). Minimal LA volume, LAEF, aLAEF, LASrc and LASctc can be used to accurately differentiate PAF from non-AF subjects (AUC > 0.8) with great sensitivity and specificity.


Assessing LA remodeling in AF using 3D STE was feasible. AF and AF burden were independently associated with LA enlargement and impairment of reservoir and pump function but not conduit function. LA function parameters can indicate underlying PAF and thus can guide AF screening strategy.

Peer Review reports


Atrial fibrillation (AF) is one of the most common arrhythmias, with an estimated prevalence between 2 and 4% worldwide [1]. The burden of the disease is rising due to aging of the population and increasing prevalence of cardiovascular disease [2]. However, the pathogenesis of AF hasn’t been fully clarified until now. Current studies suggested it may begin with hemodynamic overload and/or atrial remodeling mainly including left atrial (LA) enlargement and fibrosis which interferes with normal electric conduction and cause reentry circuits [3,4,5]. On the other hand, the presence of AF in turn attributes to LA dilation and function impairment [3].

The LA plays a phasic role during one whole cardiac circle, which importantly contributes to roughly 30% of the left ventricular (LV) stroke volume. During ventricular systole, the LA receives blood entering LA via pulmonary veins and thus serves as a “reservoir”. During early diastole, blood passively enters LV due to the pressure gradient between LA and LV, reflecting the “conduit” function of LA. Lastly, during late diastole, the LA acts as a booster pump, contracts actively and squeeze more blood into LV [6, 7]. To date, only LA maximal volume index is officially recommended in guidelines as an independent predictor for adverse cardiovascular outcomes [8], however, other LA function parameters have also been proven valuable and might even be more sensitive, including other LA volume parameters, LA emptying fraction and strain [9,10,11,12]. Strain, defined as the fractional change of length/width/thickness of one object in comparison to its original length/width/thickness, describes the deformation capacity of myocardium and thus reflect the mechanic function of heart chambers [13]. LA strain is a relatively new index to characterize LA function compared to other “conventional” echocardiographic measurements, however, some scholars think it might be superior to LA volumetric function parameters as it can distinguish LA active deformation from passive motion [14].

LA strain can be obtained by Tissue Doppler Imagine (TDI) or speckle tracking echocardiography (STE), with the latter being more widely used [14]. “Speckles” are acoustic markers of myocardium that move together with tissue during the cardiac circle. By tracking these speckles through the cardiac circle with certain software, we are able to analyse myocardial deformation function [15]. Three dimensional STE (3D STE) is a novel technique in the field of echocardiography, combining real-time 3D echocardiography with speckle tracking analysis. Theoretically, 3D STE overcomes several inherent limits of TDI and 2D STE such as angle dependency, “out of plane phenomenon”, geometric assumption, etc. [7, 16]. However, till now studies using this technique in LA function assessment is quite scarce and most current studies used softwares originally designed for LV instead of LA. Good news is that 3D STE software package dedicated to the LA was recently launched (4D LAQ, EchoPAC, GE Healthcare) which allows simultaneous analysis of LA volume and strain.

Therefore, in this study, we aimed to investigate if 3D STE can be used to comprehensively assess LA function and detect LA structural and functional changes in patients with AF. Besides, we also investigated the capacity of 3D STE derived LA parameters to identify PAF as well as their value in guiding AF screening strategy in clinical practice.


Study design and population

This was a single center cross-sectional study approved by the Huashan Institutional Review Board. Patients with clearly diagnosed AF visiting Huashan Hospital from 2019.4.1 to 2020.1.15 and 2020.6.1 to 2021.1.10Footnote 1 were continuously enrolled in this study after consent form assigned, meanwhile individuals without AF or major cardiovascular conditions other than a short history of mild hypertension visiting our hospital for health checkup were set as control group. Exclusion criteria included valvular AF, moderate to severe mitral/ tricuspid/ aortic stenosis or regurgitation, recent myocardial infarction within 3 months, acute pulmonary embolism, newly onset AF during or shortly after cardiac surgery, specific type of cardiomyopathy such as cardiac amyloidosis, obstructive hypertrophic cardiomyopathy, etc., decompensation stage of organ failure, and subjects unable to finish a successful echocardiography examination due to any reason.

Classification of AF pattern

Patients with AF were divided into paroxysmal AF group (PAF) and persistent AF group (Per-AF) according to the 2020 ESC guideline by comprehensively considering their symptoms, 12 lead electrocardiogram (ECG), 24 h Holter or other data available. The Per-AF group actually included persistent AF, long standing persistent AF and permanent AF.

Collection of clinical data

Subjects’ baseline clinical data including general health information and comorbidity were collected after enrolment mainly by reviewing their official medical records and comprehensive interviews with individual participants. The definitions for all clinical conditions were in accordance with the related guidelines [17,18,19,20,21,22].

Conventional echocardiography examination

GE Vivid E95 System and a M5Sc-D transducer (frequency: 3.5-5 MHz) were used for the transthoracic echocardiography (TTE) examination, color doppler imaging, pulsed-wave Doppler echocardiography and tissue doppler imaging. Patients were positioned in left lateral decubitus position and connected to ECG monitor. All procedures were implemented by a sophisticated echocardiography technician with over 10 year experience using standard imaging plane following the 2019 and 2015 ASE guideline [8, 23]. The measurements we took included left ventricular end-diastolic diameter (LVEDD), left ventricular end-diastolic volume (LVEDV), left ventricular end-systolic volume (LVESV), left ventricular ejection fraction (LVEF) calculated by the Biplane Simpson method. Early (E) and late (A) filling velocity at the tip of mitral valve leaflet, E to A ratio, peak early (e’) and late (a’) diastolic mitral annular velocities at both septal and lateral sides were also obtained and E/e’ were calculated to evaluate LV filling pressure with e’ being the average of septal and lateral early diastolic velocities. LA diameter in 3 directions (suprainferior, mediolateral and anteroposterior) were measured at end systole just before the mitral valve opening. We also measured LA maximal volume (2DLAVmax) using Biplane Simpson method. LV diastolic function of each patient were comprehensively evaluated based on septal and lateral e’, E/e’, 2DLAVmax index to body surface area (2DLAVImax) as well as maximal tricuspid regurgitation velocity (TRMax) and graded as normal diastolic function, indeterminate diastolic function and diastolic dysfunction according to the 2016 EACVI guideline [24]. The above criteria was not applicable for patients with Per-AF, in which case E/(septal e’) ratio ≥ 11 was used to determine the existence of LV diastolic dysfunction. Effort was made to select matched RR intervals for E and e’ that were within 20% of the average heart rate.

Speckle tracking echocardiography of LA

3D speckle tracking echocardiography was performed right after the 2D echo exams using the same GE E95 system and a 4V-D (frequency: 1.5-4 MHz) transducer. Subjects were asked to maintain their position and breathe effortlessly. With good ECG monitor and clear visualization of the entire LA cavity at the apical 4-chamber view, 4D Large mode were switched into. Single-beat (if patients were under AF rhythm) or multi-beat (3 beats, if patients were under sinus rhythm) 3D LA datasets were acquired and stored with one breath hold at end expiration for optimal image reconstruction. It’s worth mentioning that all control and PAF subjects underwent acquisition under sinus rhythm while patients with Per-AF were under AF rhythm. Care was taken to avoid artifacts and volume rate was set at least above 12 volumes per second to satisfy speckle tracking between consecutive frames. Offline analysis was then performed using a commercially available semiautomatic software package (4D LAQ, EchoPAC, GE Healthcare). Three long-axis planes including apical 4-chamber, apical 2-chamber and apical 3-chamber together with 12 short-axis slices from mitral annular plane to LA roof were displayed on the screen (Fig. 1A). The user then needed to set 3 landmarks on the 3 long-axis planes by placing them at the centre of MV at the annular level in each plane (Fig. 1B). Then the endocardial border of LA as well as its movement through the entire cardiac cycle was traced and tracked by the system automatically and could be adjusted by the operator if necessary (Fig. 1C). After the motion-tracking analysis, the system output the following parameters: LA phasic volumes (maximal LA volume (LAVmax), minimal LA volume (LAVmin), LA volume at the onset of atrial contraction (LAVpreA) and LAVmax index to body surface area (LAVImax)), volumetric function parameters (total LA emptying volume (LAEV) and emptying fraction (LAEF)) and strain parameters including peak atrial longitudinal/circumferential strain (LASr/LASrc), early diastolic longitudinal/circumferential strain (LAScd/LAScdc), late diastolic longitudinal/circumferential strain (LASct/LASctc) along with an LA volume-time curve or LA strain-time curve (Fig. 1D&E). Note that the LA volume-time curve and LA strain-time curve followed the same shape. Another example of the (longitudinal/circumferential) LA strain-time curve was as presented in Fig. 2. The zero point refers to the end ventricular diastole. As the blood enters LA, the LA strain increases and achieves its maximal value at end ventricular systole on mitral valve opening. The peak strain value measures the reservoir function of LA and thus is also known as LA strain reservoir (LASr). LA serves as a booster pump and contracts actively during late ventricular diastole. A plateau can be seen on the strain-time curve, which corresponds to the beginning of LA contraction. The plateau value reflects the contraction function of LA, and thus is known as LA strain contraction (LASct). During early ventricular diastole, LA serves as a conduit, allowing blood passively flows into LV. LA strain that measures the LA conduit function, namely LAScd, can be calculated as LASct minus LASr. LASrc, LAScdc and LASctc hold the same meaning as LASr, LAScd and LASct, respectively, but describe LA deformation in the circumferential direction. It’s worth mentioning that LAScd/LAScdc and LASct/LASctc are negative values since the atrial muscle fibres shorten during ventricular diastole and only the absolute value reflects the function condition. The three components of LA function can also be indicated by volumetric parameters. LAEF, representing the LA reservoir function, can be further divided into passive LAEF (pLAEF) and active LAEF (aLAEF) based on atrial contraction, reflecting conduit function and pump function respectively.

Fig. 1
figure 1

An example of motion tracking analysis of LA using 4D Auto LAQ, Echopac, GE Health. A Display of acquired LA 3D images showing apical four-chamber, two-chamber, three-chamber plane on the left and 9 slices of short-axis views from mitral annular plane to the atrial roof. B Landmark setting stage with three landmarks set at proper places according to the instruction. C Review stage with the LA endocardial contour traced. D A simulated model of reconstructed LA and output parameters. E LA strain/volume-time curve

Fig. 2
figure 2

An example of the (longitudinal/circumferential) LA strain-time curve. ED = end ventricular diastole, ES = end ventricular systole, PreA = pre atrial contraction

For patients under AF rhythm, the LA lost active contraction and thus only pLAEF existed and was equalled to LAEF, which was also the case for LAScd/LAScdc.

Statistical analysis

Continuous variables are expressed as mean ± standard deviation (SD) or median (quartiles) depending on the results of normality test. Kruskal-Wallis H test or one way ANOVA was used to assess difference among 3 groups when appropriate, followed by post hoc test with adjusted p value for multiple comparisons. Differences between two groups were compared using t test or wilcoxon rank-sum test. Categorical variables are presented as amount (percentage) and compared using Chi-square test or Fisher’s exact test where appropriate.

We randomly selected 30 3D volume datasets from control group, 60 from PAF group, 30 from Per-AF group and had the same operator as well as a different observer repeat the offline analysis at least one month after the initial measurements to assess intraobserver and interobserver variability. The intraclass correlation coefficient (ICC) were used to assess reproducibility.

Multi-linear regression was performed to determine whether LA function parameters were independently corelated with the presence and type of AF. Variables found significant in univariant analysis (p < 0.1) or considered clinically relevant were adjusted including age, body mass index (BMI), heart rate (HR), diastolic blood pressure (DBP), LVEF, E/e’, LV diastolic function grade, comorbidity, etc..

Correlation between two variables was analysed using Spearman rank correlation and the difference between 2 correlation coefficients were tested using z test.

Receiver operating characteristics (ROC) curves of different LA function parameters to identify PAF from non-AF subjects were constructed. Cutoff value was determined with the greatest Youden index (sensitivity + specificity − 1) and we calculated the positive predictive value (PPV) as well as negative predictive value (NPV) of different LA parameters in populations with PAF prevalence of 1 and 20%.

SPSS (IBM, Armonk, New York, USA) was used for statistical analysis. P < 0.05 were considered statistically significant.


Feasibility of 3D STE

Three hundred thirty four subjects including 59 non-AF subjects, 144 PAF and 131 Per-AF patients were initially enrolled in our study and completed conventional transthoracic echocardiographic exam as well as 3D STE. Five PAF as well as 22 Per-AF were then excluded during offline analysis process due to poor image quality, existence of artifacts or tachycardia leading to unsatisfying endocardial tracking. Eventually, 59 non-AF subjects, 139 PAF and 109 Per-AF patients were included in further analysis. The overall success rate of 3D STE was 91.9%.

Description of study population

Baseline characteristics of the 3 groups were summarized in Table 1. The disease course of patients with Per-AF was longer than that of patients with PAF. There was no difference in terms of gender, height, weight, BMI, body surface area (BSA) or blood pressure at examination among three groups. However, subjects with Per-AF were older and had higher heart rate than the other two groups. Both AF groups had higher prevalence of hypertension, coronary artery disease (CAD) and stroke compared with control group meanwhile more subjects in Per-AF group were comorbid with congestive heart failure (CHF) and diabetic mellitus (T2DM). In addition, patients in PAF group had relatively lower CHA2DS2-VASc score, lower concentration of N terminal pro B type natriuretic peptide (NT pro-BNP), higher percentage of using type I and type III antiarrhythmic medication as well as ablation history compared to Per-AF group.

Table 1 Baseline characteristics among three groups

Conventional TTE derived LV and LA parameters

Results of conventional 2D echocardiography exam were shown in Table 2. The LV systolic function of subjects in Per-AF group were relatively compromised demonstrated by lower LVEF than subjects in the other two groups. The LV diastolic function were sequentially decreased among 3 groups. All 3 LA diameters significantly increased along with the progression of AF.

Table 2 Conventional TTE results among three groups

Reproducibility of 3D STE derived parameters under sinus and AF rhythm

Ninety 3D volume data sets acquired under sinus rhythm and 30 under AF rhythm were randomly selected for reproducibility test. ICC as well as its 95% confidence interval (95%CI) for intraobserver and interobserver variability were summarized in Table 3. For data sets acquired under sinus rhythm, all 3D STE derived parameters except LAScd and LAScdc had excellent intraobserver and interobserver reproducibility demonstrated by an ICC greater than 0.75. On the other hand, 3D STE derived LA function parameters acquired under AF rhythm had only modest interobserver reproducibility in despite of great intraobserver reproducibility.

Table 3 Intraclass correlation coefficient for 3D STE derived parameters testing intraobserver and interobserver reproducibility under sinus or AF rhythm

Association between AF and LA function demonstrated by 3D STE derived parameters

LA structural and functional changes in subjects with PAF and per-AF

Comparison of 3D STE derived LA parameters among 3 groups were shown in Table 4. Note that LAVpreA, pLAEF, aLAEF, LASct and LASctc weren’t available in subjects under AF rhythm due to loss of atrial contraction. All four LA volume parameters were increased in PAF group compared with control group, and further increased in Per-AF group. LAEF, LASr and LASrc representing LA reservoir function decreased sequentially within the three groups. However, pLAEF reflecting LA conduit function seemed paradoxically higher in Per-AF group than in PAF group though no significant difference were shown by LAScd or LAScdc. LA pump function was notably impaired in patients with PAF even under sinus rhythm indicated by reduced aLAEF, LASct and LASctc.

Table 4 Changes of 3D STE derived LA function parameters in PAF and Per-AF

AF and higher burden of AF were independently associated with impairment of LA reservoir and pump function

Since there were some baseline imbalances among 3 groups and LA function were largely influenced by LV function, multivariant analysis was performed to illustrate LA morpho-functional changes independently associated with AF as well as AF type, and the partial regression coefficient(β) and adjusted p value were summarized in Table 5. It’s notable that pLAEF and LAScd no longer showed correlation with AF except that LAScdc was reduced in PAF group. Other than that, LA enlargement and impairment of the reservoir and pump function were independently associated with AF, and the remodeling was more significant in Per-AF group.

Table 5 Multi-linear regression of AF type and 3D STE derived LA parameters

3D STE detected LA function impairments in subjects without LA enlargement

Fifty-five(93.2%) control, 113(74.8%) PAF and 38(25.2%) Per-AF in our study had normal LA sizes defined as LAVImax≤34 ml/m2 measured by 2D Biplane Simpson method according to the 2015 ASE guideline. We further analyzed the 3D STE derived LA parameters in this sub-group and results were shown in Table 6. Significant increases in LAVmin, LAVmax and LAVImax were observed among 3 groups along with the progression of AF. LAEF, LASr, LASrc were reduced in PAF group compared with the control group, and the reductions were more pronounced in Per-AF group. LAScd and LAScdc were decreased in both PAF and Per-AF groups. LA pump function was also significantly declined demonstrated by lower aLAEF and smaller LASct/LASctc in PAF subjects than control subjects.

Table 6 3D STE derived LA function parameters in subjects with normal LA size (defined as 2DLAVImax ≤34 ml/m2)

Notable correlations between volumetric function parameters and strain parameters representing the same LA function

As the three components of LA function could be derived from volumetric parameters or strain, we analyzed the relationship between LA emptying fraction and LA strain reflecting the same function and the correlation coefficients were shown in Table 7. It was notable that the correlation between LASrc and LAEF was significantly stronger than that between LASr and LAEF. Similarly, LAScdc and LASctc also correlated better with pLAEF and aLAEF respectively, compared with LAScd and LASct.

Table 7 Correlations between LA volumetric function parameters and strains indicating the same function

3D STE derived parameters can be used to differentiate subjects with PAF from non-AF subjects

We further sought to assess whether 3D STE derived function parameters could be used to differentiate PAF patients from subjects with no recorded history of AF. ROC curves were constructed and related parameters were demonstrated in Table 8. LASrc yielded the largest AUC for identification of PAF, with 75.5% sensitivity and 81.4% specificity when the cutoff value was set at 24%. PPV and NPV of the parameters with an AUC over 0.8 when applied to populations with AF prevalence of 1% or 20% were calculated and summarized in Table 9.

Table 8 ROC analysis for differentiating PAF patients from non-AF subjects with 3D STE derived parameters
Table 9 PPV and NPV of LA parameters when applied to populations with AF prevalence of 1% or 20%


In this study, we aimed to comprehensively assess LA morphofunctional remodeling using 3D STE which was shown to be feasible and reproducible. Our main finding is that AF is independently associated with LA enlargement as well as impairment of LA reservoir and pump function but not conduit function. Besides, 3D STE derived function parameters, especially LAVmin, LAEF, aLAEF, LASrc and LASctc can be used to differentiate PAF from non-AF subjects with good sensitivity and specificity.

Feasibility of 3D STE in assessing LA function

In our study, the success rate of 3D STE analysis in each groups was 100, 96 and 83% respectively. In terms of reproducibility, both intraobserver and interobserver reproducibility were excellent when conducted under sinus rhythm, however, in patients with Per-AF, the interobserver reproducibility was less satisfying. The reasons were as follows: 1) LA enlargement was more significant in Per-AF and thus wider sector angle was needed during acquisition to include the entire LA which lowered the frame rates and image quality, 2) multi-beat acquisition might be feasible under AF rhythm, but the reconstruction effect would be largely compromised, 3) only single-beat acquisition can be used with obvious heart rate variation which significantly lowered the frame rates, 4) patients with Per-AF tended to have higher heart rates as shown by our study, and thus higher frame rates were needed in order to improve temporal resolution which again compromised image quality. To sum up, practically it’s more difficult to obtain 3D images with good quality in subjects with Per-AF, which lowers the accuracy of automatic LA contour recognition and tracing by the software, and therefore more manual correction is needed which is not always easy and introduces more error in measurements since there still lack widely accepted standard operational guide of this new technique, especially in subjects with Per-AF.

Echopac 4D Auto LAQ by GE Health, as a 3D STE software package dedicated to LA, was not available until recent years. Experience in using this software is scarce. One study demonstrated LA functional changes in end stage renal dysfunction in a total of 71 subjects, but didn’t discuss the usage of the software at the same time [25]. To our knowledge, our study is the first to analyze the feasibility and reproducibility of this software as well as its utility to assess LA morpho-functional remodeling in patients with both PAF and Per-AF, at the meantime with a decent sample size and thus filled certain gap in this field.

LA structural and functional changes associated with AF and AF pattern

As LA function is affected by many factors such as age, gender, BMI, comorbidity, LV function status, etc., therefore, in order to demonstrate LA remodeling independently associated with AF, confounders were adjusted based on univariant linear regression and clinical judgement. Our study showed that impairment of LA reservoir and contraction function was independently associated with PAF, additionally, LA reservoir function was further decreased in Per-AF patients. However, pLAEF and LAScd were not independently associated with PAF or Per-AF, suggesting that AF has relatively small direct impact on LA conduit function. Interestingly, we noticed that pLAEF was even higher in Per-AF group compared to PAF group, which we believe was caused by different calculating methods of pLAEF between 2 groups as it was considered equal to LAEF under AF rhythm. Therefore, comparison of LA conduit function between subjects under different rhythm is not suitable.

Our results are consistent with previous studies, demonstrating association between AF and enlarged LA volumes as well as impaired LA function indicated by conventional 2D echocardiography derived volume measurements [10, 26] or longitudinal LA strain and strain rate derived from 2D STE [27] or DTI [28]. However, 3D STE, as a novel technique combining real-time echocardiography with speckle tracking analysis, is a lot less applied in LA function assessment compared to LV. Mochizuki et al. used 2D STE and 3D STE to compare atrial deformation and synchrony in healthy subjects and patients with AF and found that 3D STE was less time consuming and more sensitive [16]. Significant decreases in longitudinal, circumferential and area strain across AF types were detected, however, confounders were not adjusted in their study in despite of some significant imbalances at baseline and thus couldn’t yet manifest AF-induced changes in LA strain as they tried to address in their article. The results of our study, with a much greater sample size, is a valuable supplement to theirs and more conclusive as explicating LA structural and functional changes independently associated with AF. A sub-study of the MYGYAR-path study conducted by Chadaide et al. compared 3D STE derived LA strain parameters in 20 AF patients and 11 healthy controls. Decrease of peak circumferential and radial strain was shown in all LA segments while reduction of peak longitudinal strain was only detected in mid and superior LA segments, and no difference was found in terms of LAEF [29]. However, their results may hold certain bias due to small sample size. Furukawa et al. sought to assess the effect of PAF on top of hypertension on LA function using 3D STE (44 PAF + HTN vs. 50 non-AF HTN), and demonstrated decreased peak atrial longitudinal strain, LAEF and synchrony compared to hypertensive patients without PAF which indicated diminished LA compliance and electro-mechanical abnormality [30]. The above studies shared one common limitation that they used 3D STE software originally designed for LV which was prone to introducing errors. Good news is that GE Health has recently launched 3D STE software package dedicated to LA, and to our knowledge, our study is the first to explore its application in AF population including subjects with Per-AF whom were usually excluded in many studies, and thus provides valuable experience.

Relationship between LA volumetric function parameters and strain

Echopac 4D Auto LAQ allows simultaneous analysis of LA volume and strain, and therefore we were able to compare the relationship between different LA function parameters. The three components of LA function can be either indicated by volumetric function parameters (LAEF, pLAEF, and aLAEF) or strain parameters (LASr/LASrc, LAScd/LAScdc, and LASct/LASctc). It’s not surprising to find strong correlation between LA emptying fraction and its corresponding strain parameter that represents the same LA function, however, on top of that, we found that the correlation between LA emptying fraction and its corresponding circumferential strain was always significantly stronger than that with the longitudinal strain, which we believe suggest that changes of LA deformation in the circumferential direction might have greater impact on LA global hemodynamic function. To our knowledge, our study is the first to compare the relationship between LA volumetric function parameters and longitudinal strain with the relationship between LA volumetric function parameters and circumferential strain since few studies were able to provide all these parameters at the same time using other techniques. At present, many studies are focusing on 2D STE derived LA longitudinal strain and its prognostic value for various adverse cardiovascular outcome, however, we assume that LA circumferential strain might possibly hold greater value, and therefore we call for more studies with this index in the future.

AF related fibrosis and LA remodeling

LA enlargement and fibrosis are hallmark of AF which not only initiate but also contribute to the maintenance and progression of AF [4, 31]. The degree of LA fibrosis can be evaluated via delayed-enhancement cardiac magnetic resonance imaging (DE-CMR), however its use in daily clinical practice is limited due to its high cost, operational complexity and significant time consumption [32]. Till now, echocardiography remains the most useful and convenient tool to evaluate cardiac structure and function. Since it’s now believed that fibrosis is associated with LA remodeling and decrease of LA compliance, LA function parameters assessed by echocardiography can be considered as a surrogate for LA fibrosis. There have been some studies supporting this idea. Cameli et al. found that peak atrial longitudinal strain (PALS) derived from 2D STE had strong inverse relationship with the extent of LA fibrosis (r = − 0.82) based on histological sample obtained from valve replacement surgery in 46 AF patients with severe mitral regurgitation while LAVImax and LAEF showed medium association with LA fibrosis (r = 0.51 and − 0.61, respectively) [33]. Another study by Kuppahally et al. demonstrated that strain value of LA mid-lateral wall was inversely related to the extent of LA fibrosis detected by DE-CMR as well as AF burden [34]. Therefore, comprehensive assessment of LA function via 3D STE can provide doctors with more disease information of different AF patients.

Identification of PAF using LA function parameters

Another important finding of our study is that LA volume and function parameters can be used to distinguish PAF from non-AF subjects with great accuracy, sensitivity and specificity.

According to current guideline, diagnosis of AF should be made based on a standard 12-lead ECG or a single-lead ECG [35]. The current prevalence of AF is estimated to be 2–4% [1], however, since some AF episodes are completely asymptomatic, or instantly self-terminate before ECG examinations are conducted, many AF patients are undetected. AF can lead to some fatal and highly-disabling complications such as stroke and congestive heart failure, which can be prevented or managed through medication. Therefore, it’s of great importance to identify those PAF patients who visited the clinics under sinus rhythm and make timely diagnosis. In the REVEAL AF study, insertable cardiac monitor (ICM) was used to screen for AF in patients with high risk for a mean observation period of 22.5 months. The detection rate of AF at 18 months was 29.3% which further increased to 40.0% if the monitoring period was prolonged to 30 months, revealing a surprising amount of undetected AF via conventional 12-lead ECG or regular 24 h ambulatory ECG monitoring [36]. The CRYSTAL-AF study aimed to detect underlying AF in patients with cryptogenic stroke using ICM and the detection rate at 6, 12, and 36 months was 8.9, 12.4 and 30%, respectively [37]. It’s easy to understand that the more intense the screening strategy is, the more undiagnosed AF can be detected, yet with higher cost. Therefore, the major issue is to determine the high risk population that truly needs prolonged rhythm monitoring. In our opinion, LA function parameters are useful under this scenario because unlike ECG abnormalities of PAF patients which is “paroxysmal”, changes of LA function are stable and can be detected by echocardiography under sinus rhythm, even before morphology changes take place, and thus can help to identify subjects with potentially underlying AF. In this study, we calculated the PPV and NPV of LA function parameters to identify PAF patients in populations with different AF prevalence. As for the general public with an estimated PAF prevalence of 1% [1, 38], the NPV of various LA volume/function parameters was excellent (> 99%), and therefore can be used to effectively rule out PAF and so we know who don’t need intensified ECG monitoring. On the other hand, in populations with a stroke history where the prevalence of AF is approximately 20% [39], the PPV of LAEF and LASrc for identification of PAF exceeded 50%, which is of great value because current commonly used tools for investigating cardiogenic stroke are echocardiography and 24 h ambulatory electrocardiogram with the purpose to find relevant structure abnormalities (such as patent foramen ovale, LA thrombus, endocardial vegetations, etc.) and AF. However, it’s highly likely that no AF episodes occur during monitoring period after the stroke event, which may lead to misclassification of the stroke type as cryptogenic stroke and delay the initiation of appropriate anticoagulant therapy. What we propose is that comprehensive assessment of LA volume and function using 3D STE can further suggest whether the patients still hold great chance (> 50%) of underlying AF even when Holter results are negative, and therefore require prolonged ECG monitoring. On the contrary, intensified AF screening can be reasonably saved when LA volume and function are within normal ranges. It needs to be point out that our aim is not to use LA function parameters derived from echocardiography to directly diagnose AF, but to help determine the optimal screening strategy for different individuals. In this way, valuable medical resources can be reasonably allocated.

Advantages and limitations of 3D STE

STE is the current preferred method for strain measurement compared to TDI, as the latter is angle dependent, time-consuming due to requirement of wall-by-wall sampling, relatively more dependent on pre-load, and poor in reproducibility. Up till now, 2D STE is the most commonly used method to obtain strain, however, it faces problems such as lacking consensus on acquisition and offline analysis, foreshortening of the 2D images, the out-of-plane phenomenon, etc.. The Echopac 4D Auto LAQ we used in our study is a 3D STE software package dedicated to the LA, which technically overcomes the inherent disadvantages of 2D STE. Its high degree of automation saves the need to manually trace LA endocardium, which not only reduces analysis time, but also greatly improves reproducibility. In addition, the information it provides is more elaborate, including reliable LA phasic volumes which is not based on geometric assumption, LA volumetric function parameters as well as LA strain in multiple directions. On the other hand, 3D STE faces certain challenges as well. As it’s a rather new technique, it lacks operation consensus and experience. Besides it has a relatively higher requirement for image quality that might lead to unsuccessful analysis or unreliable results if not fulfilled. Specific software and transducer are needed to perform 3D STE and values vary between venders which limit its application in clinical practice. However, since 3D echocardiography is a promising tool in the future, and the operation of 3D STE is actually easier compared to 2D STE, yet providing a lot more information, we believe this technique hold great value in both clinical practice and scientific research.


Our study is a single centre, cross-sectional study took place in a public teaching hospital and though we tried our best to enrol patients with AF continuously during the study period, some might still be missed due to various reasons, and thus our study may carry certain bias. The classification of PAF and Per-AF was made based on patients’ description of their symptoms combining with their previous ECG and Holter reports, and therefore misclassification between groups might exist. We detected stronger correlations between LA circumferential strains and volumetric function parameters, but the underlying mechanism couldn’t yet be drawn from this study. LA circumferential strains might possibly hold greater value than longitudinal strains, but more studies will be needed in the future.


3D STE allows us to comprehensively assess LA function with good feasibility and reproducibility, but more experience will be needed in the future especially in patients with Per-AF. Changes of LA deformation in the circumferential direction might have greater impact on LA global hemodynamic function. LA enlargement as well as impairment of LA reservoir and pump function are independently associated with AF, and are more significant with higher AF burden. LA volume and function parameters can be used to identify PAF with good accuracy and thus can guide AF screening strategy in clinical practice. As LA pump function is lost in patients under AF rhythm and conduit function is relatively less affected by AF, we believe parameters representing LA reservoir function, especially LASrc, hold the greatest clinical value. Therefore, we call for more studies to investigate the prognostic value of this index in the future.


  1. The study was temporarily suspended from 2020.1.15 to 2020.5.31 due to COVID-19.



Two dimensional speckle tracking echocardiography


Left atrial maximal volume measured by Biplane Simpson method index to body surface area


Left atrial maximal volume measured by Biplane Simpson method


Three dimensional speckle tracking echocardiography


Late filling velocity of trans-mitral Doppler flow


Atrial fibrillation


Active left atrial emptying fraction


Left atrial anteroposterior diameter


Area under the curve


Body mass index


Body surface area


Coronary artery disease


Congestive heart failure


Confidence interval


Diastolic blood pressure


Delayed-enhancement cardiac magnetic resonance imaging


Early filling velocity of trans-mitral Doppler flow


E to A ratio


E to e’ ratio (e’ being the average value of septal e’ and lateral e’)




European Society of Cardiology


Heart rate


Intraclass correlation coefficient


Insertable cardiac monitor


Left atrium


Left atrial emptying fraction


Left atrial emptying volume


Longitudinal left atrial strain conduit


Circumferential left atrial strain conduit


Longitudinal left atrial strain contraction


Circumferential left atrial strain contraction


Longitudinal left atrial strain reservoir


Circumferential left atrial strain reservoir

Lateral a’:

Peak late diastolic mitral annular velocity at lateral wall

Lateral e’:

Peak early diastolic mitral annular velocity at lateral wall


Maximal left atrial volume index to body surface area


Maximal left atrial volume


Minimal left atrial volume


Left atrial volume at the onset of atrial contraction


Left ventricle


Left ventricle diastolic function


Left ventricular end-diastolic diameter


Left ventricular end-diastolic volume


Left ventricular ejection fraction


Left ventricular end-systolic volume


Left atrial mediolateral diameter


Negative predictive value

NT pro-BNP:

N terminal pro B type natriuretic peptide


Paroxysmal atrial fibrillation


Peak atrial longitudinal strain


Persistent atrial fibrillation


Passive left atrial emptying fraction


Positive predictive value


Peripheral vascular disease


Receiver operating characteristics


Systolic blood pressure


Standard deviation

Septal a’:

Peak late diastolic mitral annular velocity at septal wall

Septal e’:

Peak early diastolic mitral annular velocity at septal wall


Left atrial suprainferior diameter


Speckle tracking echocardiography


Type 2 diabetic mellitus


Tissue doppler imaging


Peak tricuspid regurgitant velocity


Transthoracic echocardiography


  1. Benjamin EJ, Muntner P, Alonso A, et al. Heart disease and stroke statistics — 2019 update: a report from the American Heart Association. Circulation. 2019;139(10):e56–e528.

    Article  PubMed  Google Scholar 

  2. Staerk L, Sherer JA, Ko D, et al. Atrial fibrillation: epidemiology, pathophysiology, and clinical outcomes. Circ Res. 2017;120(9):1501–17.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  3. Delgado V, Di Biase L, Leung M, et al. Structure and function of the left atrium and left atrial appendage: AF and stroke implications. J Am Coll Cardiol. 2017;70(25):3157–72.

    Article  PubMed  Google Scholar 

  4. Melissa L, Rachid A, Van RPJ, et al. Relation of echocardiographic markers of left atrial fibrosis to atrial fibrillation burden. Am J Cardiol. 2018;122(4):584-591.

  5. European Heart Rhythm A, European Association for Cardio-Thoracic S, Camm AJ, et al. Guidelines for the management of atrial fibrillation: the task force for the management of atrial fibrillation of the European Society of Cardiology (ESC). Eur Heart J. 2010;31(19):2369–429.

    Article  Google Scholar 

  6. Blume GG, McLeod CJ, Barnes ME, et al. Left atrial function: physiology, assessment, and clinical implications. Eur J Echocardiogr. 2011;12(6):421–30.

    Article  PubMed  Google Scholar 

  7. Olsen FJ, Bertelsen L, de Knegt MC, et al. Multimodality cardiac imaging for the assessment of left atrial function and the Association with atrial arrhythmias. Circ Cardiovasc Imaging. 2016;9(10).

  8. Lang RM, Badano LP, Mor-Avi V, 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. Eur Heart J Cardiovasc Imaging. 2015;16(3):233–70.

    Article  PubMed  Google Scholar 

  9. Pagola J, Gonzalez-Alujas T, Flores A, et al. Left atria strain is a surrogate marker for detection of atrial fibrillation in cryptogenic strokes. Stroke. 2014;45(8):e164–6.

    Article  PubMed  Google Scholar 

  10. Gupta DK, Shah AM, Giugliano RP, et al. Left atrial structure and function in atrial fibrillation: ENGAGE AF-TIMI 48. Eur Heart J. 2014;35(22):1457–65.

    Article  PubMed  Google Scholar 

  11. Biering-SøRensen T, Christensen LM, Krieger DW, et al. LA emptying fraction improves diagnosis of paroxysmal AF after cryptogenic ischemic stroke: results from the SURPRISE study. JACC Cardiovasc Imaging. 2014;7(9):962–3.

    Article  Google Scholar 

  12. Wu VC, Takeuchi M, Kuwaki H, et al. Prognostic value of LA volumes assessed by transthoracic 3D echocardiography: comparison with 2D echocardiography. JACC Cardiovasc Imaging. 2013;6(10):1025–35.

    Article  PubMed  Google Scholar 

  13. Saraiva RM, Demirkol S, Buakhamsri A, et al. Left atrial strain measured by two-dimensional speckle tracking represents a new tool to evaluate left atrial function. J Am Soc Echocardiogr. 2010;23(2):172–80.

    Article  PubMed  Google Scholar 

  14. Dandel M, Lehmkuhl H, Knosalla C, et al. Strain and strain rate imaging by echocardiography - basic concepts and clinical applicability. Curr Cardiol Rev. 2009;5(2):133–48.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Geyer H, Caracciolo G, Abe H, et al. Assessment of myocardial mechanics using speckle tracking echocardiography: fundamentals and clinical applications. J Am Soc Echocardiogr. 2010;23(4):351-; quiz 453-355.

    Article  PubMed  Google Scholar 

  16. Mochizuki A, Yuda S, Oi Y, et al. Assessment of left atrial deformation and synchrony by three-dimensional speckle-tracking echocardiography: comparative studies in healthy subjects and patients with atrial fibrillation. J Am Soc Echocardiogr. 2013;26(2):165–74.

    Article  PubMed  Google Scholar 

  17. Hindricks G, Potpara T, Dagres N, et al. 2020 ESC guidelines for the diagnosis and management of atrial fibrillation developed in collaboration with the European Association of Cardio-Thoracic Surgery (EACTS): the task force for the diagnosis and management of atrial fibrillation of the European Society of Cardiology (ESC) developed with the special contribution of the European Heart Rhythm Association (EHRA) of the ESC. Eur Heart J. 2020.

  18. Unger T, Borghi C, Charchar F, et al. 2020 International Society of Hypertension Global Hypertension Practice Guidelines. Hypertension. 2020;75(6):1334–57.

    CAS  Article  PubMed  Google Scholar 

  19. Ponikowski P, Voors AA, Anker SD, 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 Heart J. 2016;37(27):2129–200.

    Article  PubMed  Google Scholar 

  20. Knuuti J, Wijns W, Saraste A, et al. 2019 ESC guidelines for the diagnosis and management of chronic coronary syndromes: the task force for the diagnosis and management of chronic coronary syndromes of the European Society of Cardiology (ESC). Eur Heart J. 2019;41(3):407–77.

    Article  Google Scholar 

  21. Association AD. 2. Classification and diagnosis of diabetes: standards of medical Care in Diabetes-2019. Diabetes Care. 2019;42(Suppl 1):S13–s28.

    Article  Google Scholar 

  22. Kronzon I, Tunick PA. Aortic atherosclerotic disease and stroke. Circulation. 2006;114(1):63–75.

    Article  PubMed  Google Scholar 

  23. Mitchell C, Rahko PS, Blauwet LA, et al. Guidelines for performing a comprehensive transthoracic echocardiographic examination in adults: recommendations from the American Society of Echocardiography. J Am Soc Echocardiogr. 2019;32(1):1–64.

    Article  PubMed  Google Scholar 

  24. Nagueh SF, Smiseth OA, Appleton CP, et al. Recommendations for the evaluation of left ventricular diastolic function by echocardiography: an update from the American Society of Echocardiography and the European Association of Cardiovascular Imaging. Eur Heart J Cardiovasc Imaging. 2016;17(12):1321–60.

    Article  PubMed  Google Scholar 

  25. Minm S, Yumeng X, Yao G, et al. Impact of indoxyl sulfate on left atrial strains in patients with maintenance hemodialysis: a three-dimensional speckle tracking study. Chin J Ultrasonogr. 2020;29(11):934–40.

    Google Scholar 

  26. Shin SY, Lim HE, Choi UJ, et al. Impaired transport function of the left atrium in patients with lone paroxysmal atrial fibrillation. Echocardiography. 2011;28(1):44–51.

    Article  PubMed  Google Scholar 

  27. Shang Z, Su D, Cong T, et al. Assessment of left atrial mechanical function and synchrony in paroxysmal atrial fibrillation with two-dimensional speckle tracking echocardiography. Echocardiography. 2017;34(2):176–83.

    Article  PubMed  Google Scholar 

  28. Zhang Q, Yip GW, Yu CM. Approaching regional left atrial function by tissue Doppler velocity and strain imaging. Europace. 2008;10(Suppl 3):iii62–9.

    Article  PubMed  Google Scholar 

  29. Chadaide S, Domsik P, Kalapos A, et al. Three-dimensional speckle tracking echocardiography-derived left atrial strain parameters are reduced in patients with atrial fibrillation (results from the MAGYAR-path study). Echocardiography. 2013;30(9):1078–83.

    Article  PubMed  Google Scholar 

  30. Atsuko F, Katsuhisa I, Eiichi H, et al. Three-dimensional speckle tracking imaging for assessing left atrial function in hypertensive patients with paroxysmal atrial fibrillation. Int Heart J. 2016;57(6):705.

    Article  Google Scholar 

  31. Wijesurendra RS, Casadei B. Mechanisms of atrial fibrillation. Heart. 2019;105(24):1860–7.

    CAS  Article  PubMed  Google Scholar 

  32. Lee DK, Shim J, Choi JI, et al. Left atrial fibrosis assessed with cardiac MRI in patients with paroxysmal and those with persistent atrial fibrillation. Radiology. 2019;292(3):575–82.

    Article  PubMed  Google Scholar 

  33. Cameli M, Lisi M, Righini FM, et al. Usefulness of atrial deformation analysis to predict left atrial fibrosis and endocardial thickness in patients undergoing mitral valve operations for severe mitral regurgitation secondary to mitral valve prolapse. Am J Cardiol. 2013;111(4):595–601.

    Article  PubMed  Google Scholar 

  34. Kuppahally SS, Akoum N, Burgon NS, et al. Left atrial strain and strain rate in patients with paroxysmal and persistent atrial fibrillation: relationship to left atrial structural remodeling detected by delayed-enhancement MRI. Circ Cardiovasc Imaging. 2010;3(3):231–9.

    Article  PubMed  Google Scholar 

  35. Hindricks G, Potpara T, Dagres N, et al. 2020 ESC guidelines for the diagnosis and management of atrial fibrillation developed in collaboration with the European Association for Cardio-Thoracic Surgery (EACTS): the task force for the diagnosis and management of atrial fibrillation of the European Society of Cardiology (ESC) developed with the special contribution of the European Heart Rhythm Association (EHRA) of the ESC. Eur Heart J. 2020;42(5):373–498.

    Article  Google Scholar 

  36. Reiffel JA, Verma A, Kowey PR, et al. Incidence of previously undiagnosed atrial fibrillation using Insertable cardiac monitors in a high-risk population: the REVEAL AF study. JAMA Cardiol. 2017;2(10):1120–7.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Brachmann J, Morillo CA, Sanna T, et al. Uncovering atrial fibrillation beyond short-term monitoring in cryptogenic stroke patients: three-year results from the cryptogenic stroke and underlying atrial fibrillation trial. Circ Arrhythm Electrophysiol. 2016;9(1):e003333.

    Article  PubMed  Google Scholar 

  38. Lip GY, Hee FL. Paroxysmal atrial fibrillation. Qjm. 2001;94(12):665–78.

    CAS  Article  PubMed  Google Scholar 

  39. Reiffel JA. Atrial fibrillation and stroke: epidemiology. Am J Med. 2014;127(4):e15–6.

    CAS  Article  PubMed  Google Scholar 

Download references


Not applicable.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.


This study was funded by China international medical foundation, special award for cardiovascular disease.

Author information

Authors and Affiliations



LB was in charge of data collection, offline speckle tracking analysis, statistical analysis and was a major contributor in writing the manuscript. LC and HF performed echocardiographic acquisition, and LC also performed offline speckle tracking analysis as a second observer. XG, FY, YS, LwB and HS were in charge of subjects enrolment. LB, GH, YL and LwB were in charge of study design. LwB was also a major contributor in writing the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Liwen Bao.

Ethics declarations

Ethics approval and consent to participate

This study was approved by Huashan Institutional Review Board (February 6, 2020, Ethics approval number: 2020 No. 005) and performed in accordance with the Declaration of Helsinki. All subjects enrolled in our study had signed a consent form. Details related to subjects’ identity were omitted.

Consent for publication

Subjects signed informed consents regarding publishing their data.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not 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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Bao, L., Cheng, L., Gao, X. et al. Left atrial morpho-functional remodeling in atrial fibrillation assessed by three dimensional speckle tracking echocardiography and its value in atrial fibrillation screening. Cardiovasc Ultrasound 20, 13 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Three dimensional speckle tracking echocardiography
  • Atrial fibrillation
  • Left atrial function
  • Left atrial strain
  • Remodeling