|
|
||||||||
Original Research |
1 Department of Radiology, Division of Medical Intelligence and Informatics,
Programs for Applied Biomedicine, Graduate School of Biomedical Sciences,
Hiroshima University, Hiroshima, Japan.
2 Department of Clinical Radiology, Hiroshima University Hospital, 1-2-3,
Kasumi-cho, Minami-ku, Hiroshima 734-8551, Japan.
3 Department of Molecular and Internal Medicine, Division of Clinical Medical
Science, Programs for Applied Biomedicine, Graduate School of Biomedical
Sciences, Hiroshima University, Hiroshima, Japan.
4 Department of Environmetrics and Biometrics, Research Institute for Radiation
Biology and Medicine, Hiroshima University, Hiroshima, Japan.
Received September 7, 2007;
accepted after revision January 2, 2008.
Address correspondence to J. Horiguchi
(horiguch{at}hiroshima-u.ac.jp).
Abstract
|
|
|---|
SUBJECTS AND METHODS. Ninety-one patients with suspected coronary
artery disease were scanned twice on retrospective ECG-gated helical scans.
Images with 2.5-mm thickness and 1.25-mm interval at nine cardiac phases
(center of cardiac phase: 40-80% in 5% increments) were reconstructed. The
interscan variability of coronary artery scores (Agatston, volume, and mass)
per patient and motion artifact scores per branch, subjectively assigned by
motion artifact grading (1, none; 2, minor; and 3, major), were compared
between cardiac phases for all patients, low (< 65 beats per minute [bpm])
and high (
65 bpm) heart rate patient groups.
RESULTS. For all patients, two-factor factorial analysis of variance revealed that the interscan variability was different between cardiac cycles (p < 0.01); however, this was not statistically significant between scoring algorithms (p = 0.46). The least variability was obtained at 70% on Agatston (8%) and volume (7%) and at 75% on mass (7%). Adjacent categories logit model analysis revealed that the motion artifact score was the least at 75% (left anterior descending coronary artery, 1.3; left circumflex coronary artery, 1.4; and right coronary artery, 1.9 in all patients) and that a smaller difference in calcium scores between the scans led to a smaller motion artifact score (p < 0.05).
CONCLUSION. Middiastole reconstruction (center of cardiac phase: 70-75%), with the least interscan variability and the least motion artifacts, is recommended on 64-MDCT.
Keywords: calcification calcium score cardiac imaging coronary artery CT
|
|
|---|
Among many factors influencing interscan variability in coronary artery calcium, such as the measuring method, motion artifact, signal-to-noise ratio, partial volume effects, triggering errors, point of image acquisition, and varying analysis tools, the main factor is suggested to be the motion artifact of the coronary artery on the image [11]. For images free of cardiac motion artifacts, acquisition times shorter than 19.1 milliseconds are necessary [12, 13]. By reducing motion artifacts in the optimization of ECG triggering, however, cardiac images with the least motion can be obtained with an acquisition time of 50-70 milliseconds for a patient baseline heart rate of 50-100 beats per minute (bpm) [13]. The e-Speed scanner (GE Healthcare) and dual-source CT scanners with high temporal resolution may permit freezing the heart in most patients; however, no single-source 64-MDCT scanner with half reconstruction has a temporal resolution high enough to achieve this.
There have been many single-source MDCT studies [14-25] investigating the optimal cardiac window for coronary artery calcium scoring and coronary CT angiography (Table 1). Thus far, to the best of our knowledge, no 64-MDCT studies have been reported that document the optimal cardiac phase for coronary artery calcium scoring. In a study using 16-MDCT, Schlosser et al. [14] concluded that "To assess the accurate score for the follow-up examinations, it seems mandatory to perform reconstructions that cover the entire cardiac cycle. Further studies will be needed to show whether the mean values of diastolic or systolic measurements are better suited for follow-up examinations than are reconstructions performed at a fixed time point within the cardiac cycle." On coronary CT angiography using 64-MDCT with a gantry rotation time of 0.37 second, Leschka et al. [20] stated that the best image quality was obtained in systole at 80 bpm or higher. Wintersperger et al. [21], using 64-MDCT with a gantry rotation time of 0.33 second, concluded that the time point of best image quality shifts from middiastole to systole with increasing heart rate. In contrast to this, Leschka et al. [19], using 64-MDCT with a gantry rotation time of 0.33 second, stated that the best overall image quality was obtained in middiastole, although images were rated as nondiagnostic in a small number of patients.
|
The main purpose of this study is to investigate the cardiac phase with the least interscan variability and motion artifacts on coronary artery calcium studies using single-source 0.35-second rotation speed 64-MDCT.
|
|
|---|
For patients with a positive coronary artery calcium measurement, the mean
heart rate, the change in heart rate during scanning, and number of patients
with the change in heart rate during scanning
5 bpm were obtained in all
patients, both the low (< 65 bpm) and the high (
65 bpm) heart rate
groups. These values were compared between the low and high heart rate patient
groups.
|
|
|
75 bpm. The temporal
resolution ranged from 90 to 175 milliseconds, according to the heart rate and
the number of cardiac cycles used for image reconstruction. The matrix size
was 512 x 512 pixels and the display field of view was 26 cm. The
reconstruction kernel Standard was used. Nine cardiac phase image data sets
were reconstructed by positioning the center of the temporal window at regular
5% increments from 40% to 80% of the R-R interval.
Coronary Artery Calcium Scoring
The Agatston [1], volume
[7], and mass
[8] scores were determined by a
radiologist (coronary artery calcium scoring experience of 1 year) on a
commercially available external workstation (Advantage Windows version 4.2, GE
Healthcare) and coronary artery calcium-scoring software (Smartscore version
3.5, GE Healthcare) according to the following equations:
Agatston score = slice increment/slice thickness x
(area
x cofactor)
Volume =
(area x slice increment)
Mass =
(area x slice increment x mean CT density) x
calibration factor [27].
For patients with a positive coronary artery calcium measurement, each of the Agatston, volume, and mass scores was compared among cardiac phases and two repeated scans (repeated measures analysis of variance). This is because, in cases where coronary artery calcium scores in repeated scans are both negative, the variability value is not obtained. If this occurs in specific cardiac phases or coronary artery calcium measurement algorithms, cross comparison between cardiac phases and algorithms becomes impossible. The coronary artery calcium scores were transformed to logarithmic scale to reduce skewness.
For each coronary artery calcium score in each patient, the interscan variability for the first (scan1) and the second scan (scan2) was evaluated using the percentage difference:
Interscan variability = [absolute (scan1 - scan2)/0.5 x (scan1 + scan2)] x 100.
The interscan variability was compared between cardiac phases and scoring
algorithms (two-factor factorial analysis of variance) for all patients, both
the low-(< 65 bpm) and the high-(
65 bpm) heart-rate patient groups. We
performed monthly scanning of a calibration phantom (Anthropomorphic Cardio
Phantom, Institute of Medical Physics and QRM GmbH).
Quantitative Assessment of Motion Artifacts
The motion artifact score was assigned for calcium-score-positive main
coronary artery branches (left anterior descending, left circumflex, and right
coronary arteries) using a 3-point grading scale (grade 1, none; grade 2,
minor; and grade 3, major) (Figs.
1A,
1B and
1C) by two independent
radiologists with 1 and 8 years of coronary artery calcium scoring experience,
respectively. In cases in which two or more calcified plaques were present in
the coronary artery branch, the score was assigned to the largest plaque. If
the assigned grades were different between observers, the determination was
achieved by consensus. The larger motion artifact scores between the
scans—that is, the score assigned to the most severe motion
artifact—was used for further investigation.
For the three heart rate groups and for the three coronary artery branches, the cardiac phase that provided the least motion artifact and the relationship between the degree of motion artifact and the difference in coronary artery calcium scores of two scans and cardiac phases were investigated by the adjacent categories logit model [28], which is a statistical model for describing the regression relationship between ordinal categoric response (e.g., motion artifact score) and explanatory variables (e.g., difference of coronary artery calcium score of two scans and cardiac phases). Furthermore, whether the motion artifact score was different between the low and high heart rate groups was tested using the Mann-Whitney U test.
|
|
|
|
|
|---|
5 bpm tended to be greater compared with
the low-heart-rate patient group; however, it did not reach the level of
significant difference (chi-square test, scan 1: p = 0.42, and scan
2: p = 0.06).
|
Coronary Artery Calcium Scores
The Agatston, volume, and mass scores on the nine cardiac phases in the two
scans are summarized in Table
3. On volume and mass scores, repeated measures analysis of
variance revealed that there were no statistical differences between the two
scans (p = 0.99 and 0.98); however, there were statistically
significant differences between cardiac phases (p < 0.01). On the
Agatston score, repeated measures analysis of variance revealed that there
were no statistical differences between the two scans (p = 0.99) or
between cardiac phases (p = 0.1).
|
Interscan Variability
The interscan variability in Agatston, volume, and mass scores on the nine
cardiac phases are shown in Figures
2A,
2B and
2C. For the analysis of all
patients, two-factor factorial analysis of variance revealed that there were
significant differences between cardiac phases (p < 0.01);
however, there were none between scoring algorithms (p = 0.46). For
the all-heart-rate patient group, the least variability was obtained in
diastole—that is, at 70% on Agatston (8% ± 8%) and volume (7%
± 10%) and at 75% on mass (7% ± 8%). For the low-heart-rate
patient group, images in diastole showed low variability; however, images in
end-systole, especially at 45%, showed high variability. In contrast to this,
for the high-heart-rate patient group, images at 45% showed the lowest vari
ability, although images in diastole showed comparable results to those in
end-systole (variability
of 10%).
Quantitative Assessment of Motion Artifacts
Of a total of 2,574 assignments of motion artifact score (143 coronary
arteries x 9 cardiac cycles x 2 scans), the same grading was
observed in 2,446 (95%). The remaining 128 (5%) assignments differed between
observers 1 and 2, and this was settled by consensus. Motion artifact scores
(larger scores between the scans) per branch in the three heart rate groups
are shown in Table 4. The
differences between cardiac phases were much greater than the interobserver
variability of 5%. The adjacent categories logit model revealed that in the
all-heart-rate groups the least motion artifacts were shown at 75-80% on the
right coronary artery, 70-75% on the left anterior descending coronary artery,
and 75-80% on the left circumflex coronary artery (p < 0.05).
|
In the adjacent categories logit model, the reason the cardiac phase with the least motion artifact included two neighboring cardiac phases (e.g., 75-80%) was binding of cardiac phases without any statistical difference in adjacent categories. The test also revealed that a smaller difference of calcium mass scores between the scans was associated with a smaller motion artifact score (p < 0.01). Motion artifact scores for the low-heart-rate group in the three coronary artery branches were lower than those for the high-heart-rate group on images at 75% (Mann-Whitney U test, p < 0.01), whereas there was no significant difference on images at 40-60% (p = 0.09-0.93). Motion artifact scores in the low-heart-rate group were lower than those in the high-heart-rate group on images at 65%, 70%, and 80%, with partially statistical differences depending on the branch—that is, at 65%, left anterior descending coronary artery: p < 0.01; left circumflex coronary artery: p = 0.32; and right coronary artery: p = 0.09.
|
|
|---|
Coronary Artery Calcium Score
The volume and mass scores were significantly different between cardiac
phases; however, they were not significantly different between the two scans.
The Agatston scores tended to be different between cardiac phases, although
not reaching statistical significance (p = 0.1). The results are in
line with previous studies showing the difference of Agatston and volume
scores between cardiac phases on 4-MDCT
[17,
30] and 16-MDCT
[14]. The facts indicate that
cardiac images on 64-MDCT, despite accelerated gantry rotation speed, are
profoundly under the influence of cardiac motion. Therefore, reconstruction of
the cardiac phase still needs optimization. In contrast to this, the mass
scores were not significantly different between cardiac phases or between the
two scans. The reason for this is not clear; however, it might be that the
mass is the most reproducible algorithm in coronary artery calcium
scoring.
Optimal Cardiac Phase for Coronary Artery Calcium Scoring
Although there is likely to be a consensus that cardiac images with the
least motion artifacts can be obtained in diastole on a low-heart-rate
patient, the optimal cardiac phase on a high-heart-rate patient is still under
debate. We have found no MDCT studies investigating the optimal cardiac phase
for coronary artery calcium scoring on high-heart-rate patients. Different
conclusions for the optimal cardiac phase were found in 64-MDCT angiography
studies with 0.33-second and 0.37-second gantry rotation speeds. In the
current study, the least variability on high-heart-rate patients was observed
at 45%; however, diastolic images also showed favorable results. This
indicates that reconstruction in middiastole in all patients (or
reconstruction at diastole in low-heart-rate patients while switching to
end-systole in high-heart-rate patients) is reasonable.
The adjacent categories logit model revealed that a smaller difference of calcium scores between the scans was associated with fewer motion artifacts. Motion artifacts in all coronary artery branches were mostly reduced at 75%, which corresponded to the phase with the least variability in the mass algorithm. These facts indicate that motion artifacts are a dominant factor increasing variability and that the improvement of temporal resolution by acceleration of gantry rotation speed will contribute to further reducing variability.
Study Limitations
This study has limitations. We used retrospective ECG-gated CT to search
for the optimal cardiac phase for 64-MDCT. To reduce partial volume effect
increasing the variability, we used overlapping reconstruction. We have shown
low variability values in the current study; however, these levels of
variability may not be preserved for prospective ECC-triggered coronary artery
calcium scoring. Next, we used multisector reconstruction for patients with a
heart rate
75 bpm. In a CT angiography study on 0.33-second rotation
64-MDCT, Wintersperger et al.
[21] stated that multisector
reconstruction did not significantly improve image quality in any heart rate
group (
65, > 65 to
75, and > 75 bpm). We, however, are not
able to fully understand the effect of multisector reconstruction. Ideally,
cross-comparisons between two groups with single-sector and multisector
reconstructions should have been performed for evaluating the effect of
multisector reconstruction on variability and motion artifacts.
In conclusion, middiastole reconstruction (center of cardiac phase, 70-75%), with the least interscan variability and the least motion artifacts, is recommended on a 64-MDCT scanner.
Acknowledgments
We thank Megu Ohtaki for statistical analysis using the adjacent categories
logit model.
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |