1. Introduction
Left ventricular (LV) dysfunction is a common consequence of coronary artery disease (CAD), particularly in individuals who have experienced a myocardial infarction or have chronic critical coronary artery narrowing. Assessing LV systolic function is crucial for patients’ diagnosis, treatment, and prognosis.
The left ventricular ejection fraction (LVEF) is the most commonly used single parameter for evaluating LV function, but it has several limitations, such as reproducibility, image quality, and endocardial border definition. LVEF is a measure of global LV systolic function that can be assessed using various noninvasive methods, including two- (2D) and three-dimensional echocardiography, cardiac magnetic resonance imaging, computed tomography, and nuclear medicine imaging, such as single-photon emission computer tomography [
1]. LVEF can also be estimated invasively during cardiac catheterization. However, LVEF is typically measured using 2D echocardiography in clinical practice due to its low cost and wide availability, including at the bedside.
Different methods may produce slightly different LVEF values, and the choice of method depends on the individual patient and the specific clinical context [
1]. While the description of LV systolic function using LVEF is approximate, many studies have demonstrated its practicality and clinical usefulness, particularly in patients with severely compromised LV function.
A diminished LVEF < 40% defines heart failure with reduced LVEF (HFrEF) [
2]. Based on LVEF > 50% as an indicator of preserved global LV systolic function, many patients with CAD are at risk of developing heart failure (HF).
Global longitudinal peak systolic strain (GLPSS) measured using speckle tracking echocardiography can detect subtle impairment of LV systolic function [
3]. Like LVEF, GLPSS reflects global LV systolic function, but it is more sensitive to LV dysfunction than LVEF and provides additional prognostic information. However, GLPSS may be afterload-dependent [
4]. Recently, some new indices of LV systolic function, including the first-phase ejection fraction (LVEF-1) and myocardial work (MW), may be less dependent on afterload than LVEF and GLPSS [
5,
6,
7].
Gu et al. [
8] introduced LVEF-1 as a measure of the earliest part of left ventricular (LV) mechanical contraction related to blood ejection into the aorta. Unlike LVEF and GLPSS, LVEF-1 is not a global marker of LV systolic function. Instead, LVEF-1 is an index of the most dynamic part of LV ejection, occurring immediately after the opening of the aortic valve until the peak flow through the aortic valve. LVEF-1 represents the rapid myocardium shortening deactivation and myocardial wall stress reduction. It has been shown to predict preclinical HF in hypertension. Additionally, LVEF-1 may approximate an increased risk of adverse events in asymptomatic patients with aortic stenosis and preserved LVEF, and premature mortality in patients with COVID-19 [
6,
8,
9,
10,
11,
12,
13,
14,
15,
16].
Russel et al. [
17] explored left ventricular pressure-strain loops and proposed several indices to approximate the global work of LV myocardium. MW generally refers to all the LV’s mechanical work performed to pump blood to the aorta. However, only some MW relates to effective blood ejection to the aorta, representing the constructive MW. The ratio of the constructive MW to the total MW measures global work efficiency (GWE) [
18]. In contrast, the proportion of mechanical systolic LV work not contributing to ejection is the global wasted work (GWW) [
19]. However, most wasted work appears during the late systole and early diastole, right before and after the aortic valve closure when diastole begins. The LV myocytes start isometric relaxation, followed by their lengthening, but some segments of the LV continue to shorten in the early diastole, causing a partial waste of mechanical LV contraction. Indices of MW have been studied in patients with cardiac dyssynchrony, HF, CAD, or survivors of myocardial infarction [
18,
20,
21]. Usually, if the GWW increases, values of GWE decline.
To date, no study has directly compared LVEF and GLPSS with GWW, GWE, and LVEF-1 in CAD patients with no LV contractile abnormalities and LVEF > 50%. As such, the present investigation aims to explore this issue. For this reason, the diagnostic features of standard and new echocardiographic measures of the LV systolic function in healthy control subjects were compared to patients with established CAD and preserved LVEF.
2. Materials and Methods
2.1. Sample Size Estimation
We used the following assumptions for sample size computation:
2.2. Participants
2.2.1. Healthy Volunteers
In the past (November 2018 and May 2019), we gathered healthy volunteers for the Department’s Echo Lab internal reference database for a newer echo machine (Vivid E95, GE Vingmed Ultrasound AS, Horten, Norway) purchased in October 2018. Routinely, all people enrolled in our reference database undergo a clinical evaluation, which includes brachial blood pressure measurements and standard 12-lead ECG recording before transthoracic 2D echocardiography. We used the following criteria to define participants as healthy: no known acute or chronic illness, no signs and symptoms of any disease, and not taking any chronic medication except for oral hormonal contraception in women of reproductive age. Smokers were allowed to participate. Occasional use of non-steroid anti-inflammatory drugs for occasional pain (e.g., headache) was allowed but not in the 48 h preceding ECG and echo studies. Sinus rhythm and normal findings on ECG, brachial blood pressure (<130/80 mmHg) at supine rest, and echocardiographic measurements were obligatory for all healthy people.
For this study, we used anonymized data from our reference database from all consecutive volunteers who had available cine loops to measure LVEF-1 and MW indices. No other data selection criteria were applied, and the data were used retrospectively.
2.2.2. Patients with a Stable Coronary Artery Disease
From our clinical database of patients who have undergone clinically indicated coronary angiography (between 1800 and 2000 procedures per year), we selected data from 50 consecutive CAD individuals who also had performed routine 2D echocardiography on the Vivid E95 and fulfilled the following criteria:
- −
Good quality images and cine loops allowing the measurement of LVEF-1 and MW in the post-processing (directly on Vivid E95 or EchoPAC platform by GE)
- −
No contractile abnormalities of the LV and right ventricle
- −
Presence of a significant (over 50%) narrowing of the lumen of at least one coronary artery. All coronary lesions were assessed by invasive angiography according to the current guidelines of the European Society of Cardiology [
23,
24]. Visual angiographic assessment and computer-assisted quantitative coronary analysis (QCA) using the Philips Allura Xper FD 10 software version R1.3.2 (10.2.0.10090) (Philips Medical Systems Nederland B.V., Best, The Netherlands) were performed to evaluate coronary artery stenosis.
Individuals with any clinically significant valvular disease, defined as more than mild regurgitation or stenosis, were excluded from the study.
Because the study had no experimental features and was retrospective, the local Ethics Committee waived the requirement to obtain patients’ agreement and permitted the use of data from our department’s databases. This approach is consistent with the Declaration of Helsinki for retrospective studies [
25].
2.3. Clinical Examination
Peripheral arterial blood pressure was measured using an Omron 705 IT device (Omron Healthcare Co., Ltd., Kyoto, Japan). Echocardiographic analyses were performed using a Vivid E95 equipped with an M5S 3.5 MHz transducer. Standard techniques outlined by the American Society of Echocardiography were used to obtain 2D, M-mode, and Doppler measures [
26].
Digital images were captured and stored for further offline analysis. LV end-diastolic volume (LVEDV), LV end-systolic volume (LVESV), and left ventricular ejection fraction (LVEF) were calculated in four- and two-chamber views using the biplane modified Simpson’s method. A pulsed-wave Doppler (PW) was used to measure the peak velocity of the E wave, and Doppler tissue imaging for estimation of the mitral annular early diastolic velocity (E′) in both the septal and lateral annulus, and these were averaged to calculate the E/E′ ratio.
2.4. First-Phase Ejection Fraction
LVEF-1 is calculated as the relative LV volume change from end-diastole to the time of maximal velocity of aortic flow (T1V):
T1V is measured using a biplane-modified Simpson’s method. The time difference between the peak of the ECG R-wave and a maximal aortic flow velocity recorded by PW gate in the LV outflow tract, close to the aortic valve leaflets, from an apical five-chamber view, were used to estimate T1V.
2.5. Myocardial Work Quantification
Echocardiographic continuous dynamic images of the LV apical four-chamber, three-chamber, and two-chamber views at a frame rate of at least 60 frames per second were obtained. Depth, sector width, and gain were optimized for adequate LV myocardium visualization. Automated Function Imaging software of the Vivid E95 ultrasound system (GE Vingmed Ultrasound AS, Horten, Norway) was used to perform speckle-tracking analysis. GLPSS and brachial blood pressure recordings measured before the echocardiography were used to quantify global MW indices through the noninvasive LV pressure-strain analysis [
27]. In this method, the LV pressure-strain loop curve is constructed and then adjusted to the duration of isovolumic and ejection phases, which are defined by the timing of mitral and aortic valvular events evaluated through 2D echocardiography. Myocardial work and related indices were calculated using a specific module within the Automated Function Imaging software version 202.
2.6. Statistical Analysis
Qualitative data were coded as either 0 (for women or individuals from the control group) or 1 (for men or CAD patients) and then represented as numbers and percentages. The gender distribution between the control group and CAD patients was compared using the Fisher exact test. As most continuous data did not have a normal distribution (according to the D’Agostino–Pearson normality test), results are summarized as medians and the 25th and 75th percentiles (IQR).
Comparisons between healthy individuals and CAD patients were made using the unpaired Mann–Whitney test. For exploratory purposes, only LV systolic function parameters that differed significantly were compared using analysis of covariance (ANCOVA), a priori adjusted for subjects’ age, gender, and body mass index (BMI). Only ANCOVA models with a normal distribution of residuals (as determined by the Cook’s and Mahalanobis distances) were accepted and presented, even if the original covariates did not have a normal distribution.
Nonparametric Spearman tests were performed to establish correlations between LVEF or GLPSS and newer LV systolic function indices. Multivariate linear regression was analyzed to study the association of LV systolic function indices with sex, gender, blood pressure, heart rate, and LV mass index. Receiver operating characteristic (ROC) characteristics with the area under the curve (AUC) were computed for various measures of LV systolic function to determine their potential for differentiating control subjects from CAD patients.
Additionally, the obtained ROC curves were compared using the method of DeLong et al. [
28] to test whether newer indices of LV systolic function have significantly larger AUCs than LVEF and GLPSS. For all parameters with AUCs significantly departing from 0.5, the optimum cut-off points with corresponding sensitivity, specificity, positive predictive value, negative predictive value, and accuracy were calculated using the maximum value of Youden’s index. Finally, the cut-offs of all LV systolic function indices were tested first in the unadjusted and then in a priori adjusted models for participants’ age, gender, and BMI using univariate logistic regression; the purpose of this was to obtain the odds ratio (OR) with the corresponding 95% confidence intervals (CI).
As some critical parameters, such as gender distribution and BMI, were not comparable between the control and CAD groups, we used sex for the case-control matching. The results of comparisons of continuous data for the 31 case-control matched pairs by the paired Wilcoxon test are shown in the
supplementary material.
A p-value < 0.05 was considered statistically significant. MedCalc® Statistical Software version 20.110 (MedCalc Software Ltd., Ostend, Belgium; 2022) and PQStat Software version 1.8.4.136 (PQStat Software, Poznań, Poland; 2022) were used for statistical analyses.
4. Discussion
The newer measures of LV systolic function significantly differ between patients with stable CAD and normal LVEF and healthy individuals. Regardless of the comparable LVEF in both groups, GLPSS, GWW, GWE, and LVEF-1 indicate poorer LV systolic function in CAD patients. Larger LV mass correlates with worse systolic function, which also declines (except for LVEF-1) if heart rate and blood pressure increase. LVEF-1 outperforms the diagnostic properties of other LV systolic function measures—people with reduced LVEF-1 < 30% are at the highest risk for CAD regardless of age, gender, and BMI.
An impaired LV systolic function with LVEF > 50% was demonstrated with GLPSS in studies on patients with aortic stenosis, hypertension complicated by HF with diastolic dysfunction, or those treated with certain chemotherapeutics [
29,
30,
31,
32,
33]. Here, we present that GLPSS is also worse in patients with CAD and LVEF > 50% than in healthy subjects. This observation deserves some comments.
GLPSS is a direct measure of the LV walls’ deformation during systole. On the other hand, LVEF quantifies the relative contribution of the blood pumped from the LV (i.e., stroke volume) to the overall LV end-diastolic volume. Although LVEF does not reflect myocardial contractions and deformation, it is an indirect measure of LV systolic function.
A normal LVEF provides approximate information that the LV global pumping function is preserved. Nevertheless, an increased GLPSS detects some early abnormalities in myocardial deformation. In CAD patients with no visible contractile impairment and normal LVEF, deformation abnormalities could be caused by active myocardial ischemia. However, our patients had stable CAD with no angina during the examination. Long-term consequences of chronic CAD (e.g., myocardial fibrosis) might contribute to this finding. We also show that GLPSS increases with higher DBP, HR, and LVMI. Our CAD patients had significantly larger LVMI, whereas their HR and blood pressure were comparable to healthy people.
LVEF-1 measures the early systolic shortening of cardiomyocytes during the most dynamic phase of LV contraction, which is associated with blood ejection. At this time, the highest myocardial wall stress is generated. As Gu et al. [
8] suggested, a reduced LVEF-1 may result from delayed shortening or deactivation of the LV myocardium. To preserve global LVEF, the myocardial contraction (i.e., the shortening of cardiomyocytes) may last longer and sometimes exceed the systole duration. As a result, some LV segments continue their contraction early in the diastole, leading to wasted and less effective myocardial work. Our findings confirm this reasoning, as we have observed that, compared to healthy controls, GWW (global wall work) increases while GWE (global wall efficiency) declines in CAD patients.
Our study found that LVMI negatively impacts all newer LV systolic function indices. Wasted MW increased, while GWE and LVEF-1 decreased with a larger LV mass, which alters LV geometry, raises myocardial wall stress, and changes LV systolic performance. Higher GWW appears to be a consequence of LVEF-1 reduction. At normal LVEF, the LV pumping function is preserved at the expense of diastolic function when some myocardial segments continue their contraction. In our study, LVEF-1 and GWW were negatively correlated (rho = −0.34;
p = 0.0001—data not presented in Results). As shown by Gu et al. [
8] and Chirinos et al. [
34], people with hypertension and reduced LVEF-1 or increased late systolic myocardial loading gradually develop diastolic dysfunction—isovolumetric relaxation is impaired, LV diastolic pressure increases, and the left atrium dilates and shows abnormal function [
8,
34]. In CAD, diastolic dysfunction is a common finding, usually preceding LV systolic impairment [
35].
Compared to other studies on LVEF-1 and MW indices, we demonstrate that LVEF-1 is also a sensitive marker of subclinical LV dysfunction in CAD patients with normal myocardial contractility and preserved LVEF. We also demonstrate that LVEF-1 < 30% is significantly associated with a substantial risk of CAD, regardless of the patient’s age, gender, and BMI.
4.1. Study Limitations
This was a retrospective study using two distinct sets of data. The control group consisted of clinical data from consecutive healthy volunteers who underwent a 2D transthoracic echocardiography for a reference database for a newer echocardiograph (Vivid E95). The CAD group consisted of consecutive patients who met certain criteria, including having no contractile abnormalities, LVEF > 50%, and available 2D echocardiography recordings on the same echo machine. Although the minimum sample size was estimated to be 38 people per group, data from 45 healthy individuals and 50 patients with CAD were used.
Despite the expected clinical differences between healthy individuals and those with CAD, the sex distribution and BMI significantly varied between the two groups. This resulted from using consecutive data from healthy people and CAD patients. Therefore, an additional case-control matching based on sex was performed (results shown in the
Supplementary Material). Paired comparisons (Wilcoxon test) confirmed that GLPSS, GWW, GWE, and LVEF-1 were significantly worse in CAD patients than in healthy individuals (
Table S1). Moreover, the thresholds defined for the unmatched participants were also applied after sex-based case-control matching. GLPSS > −17%, GWE < 93%, GWW > 123 mmHg%, and LVEF-1 < 30% were all associated with an increased odds ratio for CAD, even in the analysis adjusted for age and BMI (
Table S2).
The sample size for this study was relatively small, but the statistical power for the Mann–Whitney test was strong (0.8848) for GLPSS. Additionally, for LVEF-1, the statistical power was even stronger (0.9998).
Our CAD patients share several common features, such as having stable CAD, no LV contractility impairment, and normal LVEF > 50%. However, this group consists of quite heterogeneous patients. Some of them have survived MI, TIA, or stroke, undergone PCI (one person had CABG) in the past, and have diabetes. Among them are single-, two-, and three-vessel CAD patients. All of these factors might be considered bias factors. Unfortunately, it is impossible to control and adjust for all of these factors in the multivariate analyses with a group of 50 people and obtain reliable results. However, we acknowledge the presence of potential biases and that their effects should be explored in more detail in a larger CAD group.
It is important to note that our findings may not be applicable to other cardiac patients, as the study specifically focused on the diagnostic properties of newer LV systolic function indices for a specific subpopulation of CAD patients (stable with normal contractility and LVEF > 50%). Additionally, as the study population was predominantly Polish and of Caucasian race, the results may not be generalizable to other racial groups.
Finally, this study utilized 2D echocardiography to measure cardiac structure and function. While other imaging modalities such as 3D echocardiography, cardiac magnetic resonance, or cardiac computed tomography may be more accurate, they also have limitations (specific clinical indications, longer examination times, and higher costs). Additionally, LVEF-1 and MW indices analyzed in this study were initially developed for use with 2D echocardiography.
4.2. Perspectives and Conclusions
Overall, we demonstrate that GLPSS, GWE, GWW, and particularly LVEF-1 help distinguish subjects with CAD and normal LVEF from healthy people. LVEF-1 focuses mainly on the early part of the LV’s blood ejection, while GWW partially covers the early diastole. The newer echocardiographic indices emphasize the importance of going beyond global LV systolic function analysis and using other parameters that better explore other phases of LV contraction and possibly the early diastole.
It is unclear which pharmaceutical agents may target the early ejection or early diastole, but those interacting with ATP metabolism and mitochondria are the most promising. Drugs acting at this time might and targeting cellular and mitochondrial energy production might significantly improve the treatment of CAD and possibly prevent its progression to HF.
In summary, contemporary markers of LV systolic function aid in assessing subjects with CAD and normal LVEF; indeed, GWW, GWE, and LVEF-1 demonstrate better functional phenotyping of CAD patients with normal LVEF but early signs of systolic dysfunction.