- Study population: Initially, 100 patients (58±8 yr, male) CAD with suspected acute coronary syndrome who were candidates for angiography and angioplasty were enrolled. The inclusion criteria for the patients were: chest pain in the last three days, and LVEF≥50% with preserved systolic function. The exclusion criteria were: previous myocardial infarction, previous heart surgery, more than mild heart valve disease, hypertension or diabetes, advanced renal failure (estimated glomerular filtration rate <30 ml/min), left ventricular dysfunction (LVEF<30%), atrial fibrillation with a heart rate of >100 beats per minute or a persistent arrhythmia that could affect image analysis, and images that were unavailable in one or more views. After hospitalization, diagnostic electrocardiograms (ECG) were obtained from the patients. Unstable chest pain was diagnosed by routine laboratory testing with cardiac troponin levels>0.112 mg/l. Echocardiographic examination was performed for all the patients before angiography. The study was approved by the Ethics Committee of Tarbiat Modares University and Shahid Rajaei Hospital (approval No. 1397.095 and 1398.091). The study period was from March 2018 to January 2020, and all of the patients signed an informed consent form.
- Conventional echocardiography: Echocardiography examination was performed 24 h before angiography using a Philips Affiniti 50 system (Philips Healthcare, Andover, Massachusetts, USA) equipped with a cardiac transducer sector S4-2 (2-4 MHz) with 80 elements, 20.3 mm. Blood pressure was recorded in the left brachial artery in the supine position using a semiautomatic device (Riester 0124, Jungingen, Germany). The systolic stress of the LV was estimated non-invasively by systolic blood pressure. All the patients lay down in the position of left lateral decubitus and were attached to the electrocardiography (ECG) device. 2D-echocardiography was superimposed in three consecutive cycles on the images. Each plane was obtained in three consecutive cardiac cycles in breath-hold with 60-100 frames/s. The end-diastolic phase was initiated at the peak R-wave of the ECG. The end-systolic phase was marked as aortic valve closure in the apical view of the long axis, and grayscale images were obtained with B-mode second harmonic imaging and using the transmitted and received frequencies (1.9/4 MHz). 2D-echocardiographic images were taken with standard parasternal short-axis views and two, three, and four-chamber apical views in accordance with the guidelines of the American Society of Echocardiography (ASE) and European Association of Cardiovascular Imaging (EACVI) [10]. To obtain circular parasternal views of the left ventricular wall at three levels, the position of the transducer was changed by different intercostal spaces. The base, mid, and apex levels in the short-axis view were defined respectively at mitral valve level, at the level of papillary muscle, and distal LV without observation of papillary muscle. Circumferential strain and transverse diameter were provided by the parasternal short-axis view. In the long-axis view, basal and apical levels [11] were displayed as the 1/3 highest and lowest point of the LV axis, respectively (Fig. 1).
In the LV long-axis view, the LV end-diastolic dimension and end-systolic dimension were measured and fractional shortening (FS) was calculated. The left ventricular ejection fraction (LVEF) was calculated by Simpson’s biplane method and based on the end-systolic volume and end-diastolic volume in 2D images [10]. The data were digitally stored in DICOM and AVI formats on a memory drive and transferred to a PC for subsequent processing.
- Coronary angiography: All the patients underwent angiography within 24 hr after echocardiography examination according to specific instructions [12]. An interventional cardiologist evaluated angiographic findings. Coronary artery stenosis was visually estimated in two planes perpendicular to each other. The location of the lesion was assessed, and the percent of stenosis diameter for each coronary lesion was determined based on the American Heart Association classification [12]. We selected CAD patients that includes the LAD territory. The patients were divided into two groups according to the results of coronary artery angiography (CAG): one or two-vessel stenosis (stenosis only in one or two vessels; LAD (left anterior descending), LAD+RCA (right coronary artery), or LAD+LCX (left circumflex), and control group (subjects with no CAD). Initially, 100 patients with suspected CAD were enrolled. Echocardiographic examination was performed for all the patients before angiography. Five patients, due to personal reasons, did not consent to CAG. Also, five patients, because of poor imaging quality and motion artifacts, were excluded from the study. Ninety patients were evaluated with 2D-STE before CAG. Sixty-four cases (71.11%), with more than 70% coronary stenosis, were successfully treated with DES (drug-eluting stents). In two cases (2.22%), due to stenosis, 3 vessel disease patients were not included for PCI, they were a candidate for CABG (coronary artery bypass grafting). Twenty-four cases (26.67%) had nonsignificant stenosis and did not undergo angioplasty. Thus, the final sample included 64 patients treated with DES and 24 patients in the control group. Significant CAD patients used medications such as β-blockers, angiotensin II receptor blockers, statins, and nitrates.
- Theory of ellipsoidal thick-walled model: Ghista and Sandler 3D elasticity model [9] is a 3D elasticity of LV in which the shape of the LV varies during systolic and diastolic phases. In the diastolic phase, the oval shape becomes smaller and the internal cavity becomes larger. During the systolic phase, LV is thicker and more ellipsoidal. Varying wall thickness and geometrical cavity of LV represent more accurate assessment and are closer to the real shape of LV and the clinical situation. The myocardium is assumed with an elastic, homogeneous, and isotropic material. The geometrical data employed in the evaluation of the stresses were obtained from 2D echocardiography.
- Regional wall stress: For the distribution of regional wall stress in equatorial, the ellipsoid thick-walled model was considered (Ghista-Sandler). With a cylindrical coordinate system (r, θ, z), we can show the equatorial surface of the closed elliptical shell [13], where L and W are the semi-major and semi-minor axes of the ellipsoid, respectively (Fig. 2).
The Ghista-Sandler model [9] is based on the ellipsoidal thick-shell theory and the 3D left ventricular geometry is simulated with a quasi-ellipsoidal model. For the geometrical model to be matched with the actual model in size, the value of "a" factor (size parameter of the model) is also important. Consequently, in this model, the stress distribution at the equatorial level for radial, longitudinal, and circumferential directions was obtained from the following equations [9].
σrr, σyy, and σww are radial, longitudinal, and circumferential stresses, respectively; A and B are intensity
stress parameters; W and H are width and wall thickness, respectively; and dc is defined as follows:
To estimate the stress based on this model, geometrical parameters, including cavity dimensions and ventricular wall thickness, were measured by echocardiography images. Off-line analysis was performed by the Philips DICOM Viewer analysis software, version 3.0 (Philips Healthcare, the Netherlands). The anterior and inferoseptal wall thicknesses and length of the LV were measured in the frozen views of apical four- and two-chamber in the peak systolic phase. The transverse diameter of the basal level was considered in the standard LV short-axis view (Fig. 3), so longitudinal, radial, and circumferential stresses were calculated by the given formulas.
- Speckle tracking echocardiography for strain analysis: The strain is a dimensionless quantity. Because of the 3D geometry of the heart and the orientation of myofibrils, the strain is defined in three directions (longitudinal, radial, and circumferential). Peak GLS (global longitudinal strain) and GCS (global circumferential strain) describe the relative length change of the LV myocardium between the end-diastolic phase and the end-systolic phase [10]. Positive strain indicates lengthening, and negative strain indicates shortening. Longitudinal and circumferential strain values are negative, and radial strain values are positive. All images were analyzed offline before having angiographic information on the central computer using the QLAB software version 12 (advanced quantification software, cardiology QLAB Philips Healthcare, Andover, the Netherlands). Myocardial function as a strain analysis in the form of frame-frame was evaluated during the cardiac cycle by semi-automatic tracking of acoustic markers. Endocardial borders were automatically traced in all frames of 2D images for each cycle at the end of systole, and poor tracking segments were manually adjusted. The QLAB software divided the left ventricular wall into 18 segments for three levels (base, mid, and apex), and every three levels were divided into six segments (anteroseptal, inferoseptal, anterior, anterolateral, inferolateral, and inferior wall) [10, 14]. Anteroseptal and anterior walls are supplied by the LAD artery, inferolateral and anterolateral by the LCX, and inferior and inferoseptal by the RCA artery (Fig. 4).
Longitudinal strain in the apical view and circumferential strain in the parasternal view of global longitudinal strain segments were calculated between aortic valve opening and closing, in the peak systolic phase. In this study, CAD patients had LAD and RCA stenosis, therefore, base anterior and inferoseptal segments were selected. An example of a strain curve is depicted in a two-chamber view for six segments in Fig. 5.
- Elastic modulus: The elastic modulus of the LV is a ratio of stress to strain and explains the functional and structural system of LV. For determination of the elastic modulus of LV muscle, passive stress in the geometrical model and strain by STE was calculated. Anterior and inferoseptal elastic moduli were estimated in the direction of longitudinal and circumferential in the systolic phase.
- Statistical analysis: All the data were presented as mean±standard deviation (SD). A comparison of differences between the variables in two groups was performed with a t-test at the significance level of 0.05. The data were tested for normality of distribution by the Kolmogorov-Smirnov (K-S) test. Sample size was estimated on 24 samples, in each group with a confidence level of 95% and test power of 85%. P-value less than 0.05 was chosen as the levels of statistical significance. The receiver operating characteristic (ROC) curve (the plot of test sensitivity versus 1-specificity) as a nonparametric analysis, as well as the area under the curve (AUC), was used. The ROC curve was employed to determine the quality of the diagnostic modality to establish cut-off points for optimal sensitivity and specificity. Intraobserver and interobserver variabilities were the differences between the measurements expressed as a percentage of the error of the means. Reproducibility was analyzed with intra-class correlation (ICC). All the statistical analyses were performed using the SPSS software package (SPSS Inc. Chicago, IL, USA).