Abstract
OBJECTIVES:
To evaluate the feasibility of using coronary computed tomography angiography (CCTA)-derived strain to detect regional myocardial dysfunction in coronary artery disease (CAD) patients with normal left ventricular ejection fraction (LVEF).
METHODS:
A total of 1,580 segments from 101 patients who underwent stressed CT myocardial perfusion imaging (CT-MPI) and CCTA were retrospectively enrolled in this study. The CT-derived global and segmental strain values were evaluated using the feature tracking technique. Segments with myocardial blood flow (MBF) < 125 ml/min/100 ml and 95 ml/min/100 ml were categorized as ischemic and infarcted, respectively.
RESULTS:
Segmental radial strain (SRS) and segmental circumferential strain (SCS) in the abnormal segments (including all segments with MBF < 125 ml/min/100 ml) were significantly lower than those in the normal segments (14.81±8.65% vs 17.17±9.13%, p < 0.001; –10.21±5.79% vs –11.86±4.52%, p < 0.001, respectively). SRS and SCS values in infarcted segments were significantly impaired compared with the ischemic segments (12.43±8.03% vs. 15.32±8.71%, p = 0.038; –7.72±5.91% vs. –10.67±5.66%, p = 0.010, respectively). The AUCs for SRS and SCS in detecting infarcted segments were 0.622 and 0.698, respectively (p < 0.05).
CONCLUSIONS:
It is feasible for using CCTA-derived strain parameters to detect regional myocardial dysfunction in CAD patients with preserved LVEF. Segmental radial and circumferential strain have the potential ability to distinguish myocardial ischemia from infarction, and normal from ischemic myocardium.
Keywords
Introduction
Coronary artery disease (CAD), which is characterized by atherosclerotic plaque formation in the coronary arteries, is considered the cause of myocardial ischemia and myocardial infarction [1]. Fractional flow reserve (FFR), single photon emission computed tomography (SPECT) and cardiac magnetic resonance (CMR) are used in current clinical practice to detect myocardial ischemia with different limitations, such as high price, increased radiation and additional requirements for invasive or semi-invasive procedures etc. Coronary computed tomography angiography (CCTA) has been widely used as an accurate method for assessing anatomic stenosis [2]. However, anatomic stenosis and myocardial ischemia may mismatch in CAD patients [3]. Stress dynamic computed tomography myocardial perfusion imaging (CT-MPI) combined with CCTA has been recognized as an ideal one-stop shop imaging modality and can assess both myocardial blood flow (MBF) and stenostic degree of coronary artery [4]. However, stress CT-MPI is subject to several limitations, such as increase of radiation dose, and the use of vasodilator stress agent, increasing the risk of examination [5].
Myocardial ischemia leads to regional myocardial dysfunction in the specific territory supplied by a stenotic coronary artery [6]. Strain analysis is a quantitative tool used in the early evaluation of regional myocardial dysfunction. Strain measured by CMR feature tracking has been explored in the literature in ischemic heart disease. Previous studies showed that myocardial strain parameters based on CMR were reduced in ischemic areas, in comparison to non-ischemic areas [7, 8]. Retrospective electrocardiogram (ECG) gating of CCTA with multi-phase reconstruction allows for three-dimensional (3D) strain analysis by the feature tracking technique [9]. CCTA-derived global and segmental strain has shown good agreement with strain using feature tracking CMR in cardiomyopathy diseases [10, 11]. However, to the best of the knowledge of the authors, few studies have been conducted on the capability of CT-derived strain to distinguish myocardial ischemia from infarction, especially in CAD patients with preserved left ventricular ejection fraction (LVEF). Therefore, the purpose of our study was to primarily investigate the feasibility of using strain parameters derived from CCTA to detect the regional myocardial dysfunction in CAD patients.
Materials and methods
Patients
The study enrolled retrospectively 109 patients with clinically suspected CAD from June 2018 to November 2019 according to inclusion criteria of aged 18 years or older, preserved EF (LVEF ≥ 50% measured by echocardiography) and successfully undergoing stress dynamic CT-MPI combined with CCTA. The exclusion criteria were as follows: 1) inadequate quality of CT-MPI or CCTA images (n = 4); 2) missing MBF information(n = 2); 3) failing to assess strain parameters (n = 2). In the enrolled 101 patients, thirty-six left ventricular myocardial segments were excluded due to poor tracking. Finally, 1580 left ventricular myocardial segments were included in the final analysis. The radiation dose was calculated from the dose-length product in a dose report.
This study was approved by our institutional review board and the informed consent requirement was waived due to its retrospective nature.
CT protocol and post processing
CT examinations were performed using a 96-row detector third-generation dual-source CT (SOMATOM Force, Siemens Healthineers, Forchheim, Germany). Double intravenous accesses in the forearms were implanted. Intravenous adenosine-triphosphate (ATP) (160μg/kg/min) was then administered 3 min prior to stress CT-MPI scan. During the stress dynamic MPI scan, a bolus of contrast medium (50 ml; iopromide (Bayer Healthcare, Berlin, Germany)) was injected into the antecubital vein at a rate of 6 mL/s, and 40 mL of normal saline were flushed to wash out the contrast medium from the superior vena cava and right ventricle. The entire LV was covered under the shuttle mode (i.e., the table was moving back and forth between acquisitions) [12]. Fixed tube voltage of 70 kV and attenuation-based tube current modulation were enabled during CT-MPI. A prospective ECG gating with an acquisition window of 250 ms, was applied to ensure myocardial perfusion captured at the end-systolic phase to reduce the beam-hardening artifact. Depending on the heart rate, a scan was performed every second or third heart cycle, resulting in a series of 10–15 acquired phases in total. Slice thickness of the reconstructed data was 3 mm.
An interval of 5 min between the MPI and CCTA provides optimal contrast and ATP washout. Retrospective ECG-gating CCTA with tube current modulation during the cardiac cycle was performed, using a bolus tracking technique and the region of interest (ROI) located in the ascending aorta [13]. The CCTA scanning parameters were as follows: the pitch, which adapted to the heart rate, tube voltage, and current selected automatically by CARE Dose 4D mode with reference to tube current–time product of 400 mAs, was 0.2 –0.39. Twenty phases were reconstructed in 5% steps of the R-R interval within the full window. Finally, the images were reconstructed using an iterative algorithm (ADMIRE, Siemens, Healthcare, strength set at a level of 3) with section thickness of 0.75mm for coronary stenosis analysis.
Image analysis
Perfusion images were reconstructed using a dedicated kernel for reduction of iodine beam-hardening artifacts and analyzed using a CT-MPI software package (SyngoVia, version VB10A, VPCT, Siemens Healthineers, Germany). Quantification of the MBF was performed using a hybrid deconvolution model [14]. According to the 17-segment model [15] and excluding the apical segment, the ROI was manually placed on the short-axis (SAX) view, and the MBF was collected on a segment base. The ROI was drawn in order to cover the entire area of the suspected perfusion defects within the segment. According to previous studies [16–18], segments were recognized to be myocardial ischemia (MBF < 125 ml/min/100 ml) and myocardial infarction (MBF < 95 ml/min/100 ml), respectively, in this study. Subsequently, the patients were divided into normal, ischemic and infarcted group, according to the segment with lowest MBF. The end-diastolic volume index (EDVi), end-systolic volume index (ESVi), stroke volume index (SVi), LVEF, left ventricular mass index (LVMi), and cardiac index (CI) were obtained by the cardiac function analysis software (SyngoVia, version VB10A, Cardiac Function, Siemens Healthineers, Germany).
A commercial software package (Circle CVI 42, Circle Cardiovascular Imaging Inc., Canada) was used to analyze the global and segmental strain parameters. Twenty phases of CCTA images were imported into CVI 42 and reconstructed for LV SAX slices, and long axis (LAX) slices (two-, three-, and four-chamber views) based on maximal intensity projection (MIP) images. Endocardial and epicardial borders and interventricular septum insertion points were semi-automatically defined at the end-diastolic phase. The software then automatically propagated the delineations over the entire heart cycle. Finally, global and segmental strain parameters were calculated by feature tracking in the radial, circumferential and longitudinal directions, based on the 3D results (Fig. 1).

Example of CT-MPI and strain analysis. Upper row shows images of a 66-year-old male patient with normal myocardial perfusion in CT-MPI, (A) whose lowest MBF was 134 ml/min/100 ml. No significant reduction in SRS) (B), SCS(C) and SLS(D) was detected. Lower row shows images of a 77-year-old male ischemic patient with infarcted myocardium, (E) whose lowest MBF was 95ml/min/100 ml in the inferior-septal segment at basal level. There was reduction of SRS(F), SCS(G) and SLS(H) in the inferior-septal myocardium compared to the normal patients. SRS, segmental radial strain; SCS, segmental circumferential strain; SLS, segmental longitudinal strain.
All image analysis was performed by Z.M.M. (with 2 years of cardiac imaging experience). The extracted features were shown in the Table 1. Then twenty patients were randomly selected to be evaluated the inter (Z.M.M. and G.Y.J.) and intra (4-week interval by Z.M.M.) reproducibility analysis.
Extracted features
CCTA, coronary computed tomography angiography; CT-MPI, computed tomography myocardial perfusion imaging.
All continuous data were checked for normality using the Kolmogorov-Smirnov test. Continuous variables are expressed as mean±standard deviation, while categorical variables as frequency with percentage. Global strain and function parameters were compared among normal, ischemic and infarcted patients by one-way analysis of variance (ANOVA). Segment strain parameters were compared using ANOVA for comparing difference of regional function. Two-by-two comparison of continuous variables was used for the LSD (L) method. Independent sample Student’s t-test was used to compare segmental strain values between normal and abnormal segments, including all MBF < 125 ml/min/100 ml segments. Receiver operating characteristic (ROC) curve analysis and logistic regression analysis were performed at segment-level to determine the association between the strain values and infarcted segments, with the former expressed as the area under the curve (AUC) and the latter as odds ratio (OR). Reproducibility between observers in strain analysis was evaluated by intra-class correlation coefficient (ICC). Statistical significance was set at p≤0.05 in all tests. Statistical analyses were carried out using SPSS software (SPSS 21.0 for Windows, IBM, USA).
Results
Study population and segments baseline characteristics
Stress CT-MPI detected 51 infarcted, 236 ischemic and 1293 normal segments (Fig. 2). The patient characteristics are displayed in Table 2. Fifteen patients had severe stenosis, with at least one lesion lumen area stenosis > 70%. Eighteen patients had moderate stenosis, with at least one lesion lumen area stenosis between 50 –70%. Twenty-five patients underwent coronary stent implantation and four underwent coronary artery bypass grafting. The radiation dose was 3.48±0.52 mSv and 4.56±1.99 mSv for MPI and CCTA scan, respectively.

Study flow chart. CT-MPI, computed tomography myocardial perfusion imaging; CCTA, coronary computed tomography angiography; MBF, myocardial blood flow.
Patient characteristics
Data are presented as the mean±standard deviation or the number (%) of subjects.
Out of the 101 patients, 50 were characterized as normal, 33 as ischemic and 18 as infarcted. In cardiac function analysis, significant difference of LVMi was found among normal, ischemic and infarcted patients (p = 0.032). LVMi in the infarcted patients was significantly higher than that in normal patients (65.61±11.12 g/m2 vs. 58.46±9.77 g/m2, p = 0.009). Meanwhile, no significant differences were detected in other cardiac function parameters and global strain values in three directions among these three groups (Table 3).
A comparison of cardiac function and global strain values among the normal, ischemic, and infarcted patients
A comparison of cardiac function and global strain values among the normal, ischemic, and infarcted patients
Data are presented as the mean±standard deviation. LVEF, left ventricular ejection fraction; EDVi, end-diastolic volume index; ESVi, end-systolic volume index; SVi, stroke volume index; LVMi, left ventricular mass index; GRS, global radial strain; GCS, global circumferential strain; GLS, global longitudinal strain. *p < 0.05, a p < 0.05 compared with the normal patients.
Segmental radial strain (SRS) and segmental circumferential strain (SCS) in abnormal segments were significantly lower than those in normal segments (14.81±8.65% vs 17.17±9.13%, p < 0.001; –10.21±5.79% vs –11.86±4.52%, p < 0.001). On the other hand, no significant difference was found in the segmental longitudinal strain (SLS) values between normal and abnormal segments (Fig. 3). SRS values were significantly impaired in both ischemic and infarcted segments, compared with the normal segments (15.32±8.71% vs 17.17±9.13%, p = 0.004; 12.43±8.03% vs 17.17±9.13%, p < 0.001). The same results were also found in the SCS values among normal, ischemic and infarcted segments (–10.67±5.66% vs –11.86±4.52%, p = 0.008; –7.72±5.91% vs –11.86±4.52%, p < 0.001). Moreover, infarcted segments had further reduction in SRS and SCS compared with ischemic segments at a significant level (12.43±8.03% vs. 15.32±8.71%, p = 0.038; –7.72±5.91% vs. –10.67±5.66%, p = 0.010). Finally, no significant difference in the SLS values was detected among normal, ischemic and infarcted segments (Table 4).

Comparison of SRS (A), SCS (B) and SLS (C) in the normal and abnormal segments. SRS was significantly impaired in the abnormal segments, compared with the normal segments (A). The abnormal segments had significantly decreased SCS, compared with the normal segments (B). No significant difference was found between these two groups, with regards to SLS (C). SRS, segmental radial strain; SCS, segmental circumferential strain; SLS, segmental longitudinal strain. *p < 0.05.
A comparison of strain values among the normal, ischemic, and infarcted segments
SRS, segmental radial strain; SCS, segmental circumferential strain; SLS, segmental longitudinal strain. *p < 0.05, a p < 0.05 compared with the normal segments, b p < 0.05 compared with the ischemic segments.
The results of the ROC analysis and logistic regression using strain parameters to predict infarcted segments from non-infarcted segments were illustrated in Table 5. Univariate logistic regression analysis showed that SRS and SCS were significantly associated with infarcted segments, while in the multivariate logistic regression, SCS was independently associated with infarcted segments. The AUCs for SRS and SCS in detecting infarcted segments were 0.622 and 0.698, respectively (p < 0.05) (Fig. 4). A cut-off of –10.73% for SCS allowed the differentiation of infarction segments with sensitivity 72.7%, specificity 57.5% and accuracy 57.9%.
Results of the ROC analysis and logistic regression analysis demonstrating the association between strain measurements and infarcted segments
Results of the ROC analysis and logistic regression analysis demonstrating the association between strain measurements and infarcted segments
SRS, segmental radial strain; SCS, segmental circumferential strain; SLS, segmental longitudinal strain. *p < 0.001.

Receiver-operator characteristic curves with corresponding area under the curves of SRS and SCS for detecting infarcted segments. The AUCs for SRS and SCS detecting infarcted segments were 0.622 and 0.698, respectively. SRS, segmental radial strain; SCS, segmental circumferential strain.
The global strain showed higher ICC of variability than the segmental strain. The ICC of radial strain values was slightly low (0.70). SLS had slightly higher ICC compared with SRS and SCS (0.76 –0.87 vs 0.70 –0.82, 0.75 –0.81) (details are shown in the Table 6).
Intra- and inter-observer variability of strain values
Intra- and inter-observer variability of strain values
GRS, global radial strain; GCS, global circumferential strain; GLS, global longitudinal strain; SRS, segmental radial strain; SCS, segmental circumferential strain; SLS, segmental longitudinal strain; ICC, intra-class correlation coefficient.
The major findings of this study were as follows: (1) CCTA-derived segmental strains were reduced in ischemic and infarcted segments compared with normal myocardial segments respectively; (2) Segmental radial and circumferential strains were furtherly reduced in infarcted segments compared with ischemic segments, and (3) the segments with circumferential strain below –10.73% could be infarcted segments.
At patient level comparison, neither cardiac function except LVMi nor global strain parameters showed differences among normal, ischemic and infarcted patients. Previous studies have shown that the global strain values were significantly lower in patients with obstructive coronary heart disease [19]. In the present study, the average global strain value was lower in the ischemic and infarcted group, but not at a statistically significant level. This may be related to the small sample size, as well as to the overall mild degree of ischemia in the ischemic patients or less segments involved in the infarcted patients. Additionally, abnormal LV systolic function occurs further along in CAD pathway, while patients enrolled in this study had preserved LVEF [20].
At segment level comparison, CT-derived SRS and SCS values were lower in segments with reduced perfusion. SCS and SRS derived from CMR were able to differentiate ischemic from non ischemic segments according to previous studies [7]. The conclusion was corroborated in this study. SCS was independently associated with infarcted segments. SCS could be superior to SRS in discriminating infarcted from non-infarcted myocardium. Everaars et al. found that SCS was the most accurate tool for quantifying regional function compared with SRS and wall thickening [21]. Noriko et al. also demonstrated that the SCS assessed by CMR had superior discrimination ability of detecting infarcted segments with the AUC = 0.82 [22]. The proposed SCS cut-off value for diagnosing infarcted segments in this study is –10.73% with AUC = 0.698. A possible reason why AUC for detecting infarcted segments was lower in our study could be that the temporal resolution and soft tissue resolution of CCTA were lower than CMR. On the other hand, CMR is still not as widely used as CCTA due to its technical complexity and long examination time. Stress CT MPI+CCTA can provide coronary anatomy and function assessment, while CCTA-derived strain can provide myocardial motion assessment. We support that a quantitative strain analysis from the rest CCTA is a useful tool in practice to help detecting CAD, especially for the patients who are contraindicated to receive stress inducing agents. Moreover, CCTA-derived SCS may be useful for predicting the necessity of further tests such as CT-MPI, SPECT or coronary angiography [23].
Longitudinal strain is a stable and reproducible parameter in many previous literature studies [24, 25], which is consistent with the findings of our study. However, it opposes studies stating that circumferential strain had higher reproducibility [26, 27]. Firstly, different software programs such as Circle CVI42 and TomTec Arena may influence the reproducibility. Secondly, strain parameters were calculated by feature tracking based on the 3D results in this study. The 3D feature tracking considers the short- and long-axis data simultaneously for strain calculation. Previouse study also showed that 3D feature tracking presented worse intra-observer reproducibility in radial and circumferential, comparable reproducibility in longitudinal tracking compared with 2D feature tracking [28].
This study has several limitations. First, this study was a retrospective design and a small data sample was collected from a single center, limiting the generalization ability of the conclusions. Second, our study showed that normal GLS was –6.58% based on CT images, which was lower than the value based on CMR [29]. Further studies with larger study cohorts compared with CMR are required to confirm the detection value of CT-derived strain analysis. Additionally, the best cut-off value for predicting ischemia or infarcted myocardium remains a controversial subject. Last, radiation exposure for retrospective ECG-gated cardiac CT was relatively high. However, tube-current was adjusted during R-R interval to reduce the radiation dose in this study.
In conclusion, it is feasible for CCTA-derived strain to detect regional myocardial dysfunction in coronary artery disease (CAD) patients with preserved LVEF. Segmental radial and circumferential strain have the potential ability to distinguish myocardial ischemia from infarction, and normal from ischemic myocardium.
Funding information
This study received funding from Research Foundation of Jiangsu Medical Association under grant of SYH-3201150-0006 (Yi Xu).
