Dynamic variations in the ultrasound greyscale median of carotid artery plaques
© Kanber et al.; licensee BioMed Central Ltd. 2013
Received: 4 May 2013
Accepted: 6 June 2013
Published: 14 June 2013
Several studies have found that the ultrasound greyscale median (GSM) of carotid artery plaques may be useful for predicting the risk of cerebrovascular events. However, measurements of GSM are typically performed on still ultrasound images ignoring any variations that may be observed on a frame-by-frame basis. The aim of this study was to establish the existence and investigate the nature and extent of these variations.
Employing a novel method that enabled plaque boundaries to be tracked semi-automatically, variations in the plaque GSM and observed cross-sectional area were measured for 27 carotid artery plaques (19 consecutive patients, stenosis range 10%-80%) over image sequences of up to 10 seconds in length acquired with a mean frame rate of 32 frames per second.
Our results showed a mean inter-frame coefficient of variation (CV) of 5.2% (s.d. 2.5%) for GSM and 4.2% (s.d. 2.9%) for the plaque area. Thirteen of the 27 plaques (48%) exhibited CV in GSM greater than 5% whereas only six plaques (22%) had CV in plaque area of greater than 5%. There was no significant correlation between the CV of GSM and plaque area.
Inter-frame variations in the plaque GSM such as those found in this study have implications on the reproducibility of GSM measurements and their clinical utility. Studies assessing the GSM of carotid artery plaques should consider these variations.
KeywordsUltrasound Carotid plaque GSM Frame-by-frame Inter-frame Variations
The North American Symptomatic Carotid Endarterectomy Trial (NASCET) and the European Carotid Surgery Trial (ECST) have shown that surgery in symptomatic patients with severe internal carotid artery stenosis results in a six-fold reduction of stroke risk [1, 2]. However, patients who do not have severe stenoses and patients who are asymptomatic can also go on to develop stroke. It is, therefore, important to be able to determine whether any of these patients have carotid plaques which are high-risk or unstable. Ultrasound greyscale median (GSM) is commonly used to quantify the ultrasound appearance of carotid plaques, and several studies have found that it may be valuable for predicting the risk of cerebrovascular events. In particular, statistically significant associations have been reported between plaque GSM and the presence of cerebrovascular symptoms [3, 4], cerebral infarction in symptomatic and asymptomatic patients [5–7], recurrent cerebrovascular events before undergoing carotid endarterectomy , and the overall risk of stroke in symptomatic patients , asymptomatic patients , and during and after carotid artery stenting .
GSM measurements currently have poor reproducibility across studies. This can be partly attributed to the differences in the acquisition settings used in separate studies. In order to reduce this variability, investigators have attempted to standardise ultrasonic images of carotid plaques by specifying certain acquisition settings to be used for carotid artery scanning and normalizing the resultant ultrasound images . However, existing studies typically measured GSM on still ultrasound images, and thus ignored any variations that may have been observed on a frame-by-frame basis. The purpose of this study was to establish the existence of and investigate the nature and extent of any frame-by-frame variations in the plaque GSM using a novel technique for tracking plaque boundaries in ultrasound image sequences. We hypothesized that variations in the GSM of carotid artery plaques may occur due to deformation of the plaque during the cardiac cycle, and other confounding factors such as out-of-plane plaque, patient or probe motion. Furthermore, it was hypothesized that plaques of different composition and morphology may exhibit different inter-frame variations in GSM in otherwise equivalent hemodynamic circumstances and hence the measurement of these variations may give useful insight into the dynamic behaviour of plaques and help identify vulnerable plaques.
Frame-by-frame variations in the plaque GSM and area of 27 carotid artery plaques (19 consecutive patients, 11 males, mean age 76, stenosis range 10%-80%) were studied by measuring the GSM and area of plaques on each image frame separately and computing the mean, the standard deviation (s.d.) and the coefficient of variation (s.d./mean) across the frames. The image sequences used were of up to 10 seconds in length (average 4.4 seconds) and were acquired with a mean frame rate of 32 frames per second. The degrees of stenosis of the corresponding arteries were measured using criteria consistent with the NASCET methodology utilizing blood flow velocities in conjunction with the B-Mode and colour flow imaging [2, 13, 14]. Eleven of the plaques studied were found to be asymptomatic and the remaining sixteen symptomatic after assessment at the University Hospitals of Leicester NHS Trust’s Rapid Access Transient Ischemic Attack (TIA) Clinic. The use of the clinical data for our research had been approved by the National Research Ethics Service (NRES) Committee East Midlands - Northampton (reference 11/EM/0249), and each patient gave informed consent before participating in the study. The ultrasound data were obtained as longitudinal cross-sections using a Philips iU22 ultrasound scanner (Philips Healthcare, Eindhoven, The Netherlands) with an L9-3 probe and included B-Mode (i.e. greyscale) and Colour Doppler image sequences. The vascular carotid preset on the scanner was used (Vasc Car preset, persistence low, XRES and SONOCT on) and the gain was optimized by the operator (TCH) who is an experienced vascular sonographer. In the case of B-Mode acquisitions, the greyscale transfer curve was kept set to Gray Map 2, as this was reported to be the most linear transfer curve on this scanner . Colour Doppler cine-loops were used as a qualitative aid to identifying the location and extent of the plaques, while the B-Mode data were used for the quantitative analyses of the plaque GSM and cross-sectional area.
The arterial lumen tracking technique based on this algorithm was previously found to have good arterial wall tracking performance, comparable to that of Tissue Doppler Imaging .
Speckle tracking, which was used to track the plaque-underlying tissue boundaries, is a standard image analysis technique that involves measuring the similarity between a template and a search image . Given a point to speckle track, a region is defined around the point and used as a template. The process is then essentially to find the position in the search image that has the largest similarity to this template. There are many different measures of similarity between a template and a search image; in this study the normalized correlation coefficient was used since it is invariant to changes in image amplitude . Square regions of approximate area 1.4 × 1.4 cm2 were employed. This template size was found to produce optimum speckle tracking quality in our study as was verified by observing the plaque tracking results (Additional file 1).
Regions of plaques that could not be distinguished from the arterial lumen (e.g. echo-free regions) and regions of plaques in areas of acoustic shadowing were excluded from analysis. Plaques for which anechoic regions and regions of shadowing exceeded more than 70% of the total plaque cross-sectional area as observed on Colour Doppler sequences were not included in the study.
Image normalisation was carried out using two different methods in order to observe their effects on the frame-by-frame variations observed. The first normalisation (NORM1) was performed by linearly scaling the ultrasound image intensities so that the GSM of a user-selected blood region inside the vessel lumen was mapped to 0 and the brightest region of the adventitia was mapped to 190. Both of these regions were 5 × 5 pixels in size, corresponding to an approximate area of 0.4 × 0.4 mm2. The reference regions were selected on the first image of the sequence and the reference GSM values calculated on the first frame were applied to that and all subsequent images. The second normalisation (NORM2) involved selecting blood and adventitia regions on each image separately, thus applying separate reference values to individual images.
A semi-qualitative assessment of whether physiologically reasonable (e.g. of the order of 60/min) periodical variations were visually apparent on the GSM and plaque area waveforms was also carried out. This involved measuring the frequency of any periodical variations seen on the GSM and cross-sectional area waveforms and considering frequencies in the range 50/min - 160/min to be potentially attributable to cardiac variations. Conversely, variations with frequencies lower than 50/min or higher than 160/min were not considered to be due to physiological sources and such plaques were placed in the same category as those not showing any apparent periodical variations in the GSM and cross-sectional area.
Statistical analyses were carried out using MATLAB version 7.14, release 2012a (MathWorks, Natick, Massachusetts, USA) and SPSS version 20 (IBM Corporation, Armonk, New York, USA). Spearman’s test was used to study the correlation between the inter-frame variations in GSM and area, since neither of these parameters was expected to follow a Gaussian distribution and any correlation between the two was likely to be non-linear. Multi-variable linear regression was used to study the contribution of other plaque GSM and area parameters to the differences observed in the magnitude of the GSM variations for each plaque. An unpaired, non-parametric Mann–Whitney U-test was used to investigate whether the GSM values averaged across all frames, as well as their standard deviations and the coefficients of variation, differed significantly between the asymptomatic and symptomatic plaque groups. Two-tailed values of significance were used in each case.
Intra-observer coefficients of variation for eight selected plaque samples of varying echogenicities were studied by measuring the frame-by-frame variations in the plaque GSM and cross-sectional area five times for each plaque. The measurements were made by the same operator (BK) and in sequential order. The same ultrasound acquisition sequences were used for each plaque respectively. The eight plaques concerned were selected from the available dataset to give a wide spectrum of plaque echogenicities for reproducibility analysis.
Comparison against manual measurements
In order to compare the plaque GSM and cross-sectional areas obtained using our method with those obtained using manual delineation, plaque GSM and cross-sectional area were measured by the same operator (BK) using manual delineation for every 5th frame, for each of the same eight plaque samples used for our study of reproducibility. This enabled the magnitude of and variation in the plaque GSM and cross-sectional areas to be compared between the two techniques. A Bland-Altman analysis was also carried out to assess the agreement between the GSM measurements made using our method and manual delineation on matching image frames.
Variations observed in the plaque GSM and area
Periodical variations observed?
Periodic variations with frequencies of the order of 60/min in either or both of the plaque GSM or area waveforms were observed for 12 plaques (50%) but not observed for 12 other plaques (50%, Table 1). Three plaques were excluded from this analysis as they had short acquisition sequences.
Results of multi-variable linear regression, testing for the influences of (a) mean frame-by-frame GSM values, (b) mean frame-by-frame plaque areas, and (c) the standard deviations of the frame-by-frame plaque areas on the standard deviations of the frame-by-frame GSM values
Standardized coefficient (β)
Intra-observer coefficients of variation (standard errors) for the measurement of the inter-frame mean GSM (un-normalized and NORM1 normalized) and mean area, for eight plaque samples
Comparison with manual delineation for eight selected plaque samples
Our investigation highlighted variations in the GSM and area of plaques when measured on a frame-by-frame basis throughout ultrasound image sequences. Image normalisation did not reduce the extent of the GSM variations and in some cases resulted in greater variation. These results demonstrate that frame-by-frame variations in the plaque GSM cannot be offset by applying normalisation factors based on the selection of blood and adventitia regions in one of the frames. Furthermore, selecting separate reference regions in all images introduced an additional source of variability due to the subjective nature of the process. The reference regions were user-selected and not computerized as they are easily identified by the operator and it would have been difficult to ensure the accuracy of a computerized selection. In NORM1 normalisation, the coefficients of variation for GSM changed after normalization for some plaques. This occurred when the blood reference regions had non-zero GSM, causing an intercept to be introduced into the linear relationship between the normalized and un-normalized greyscale values. Furthermore, as GSM values are limited to the range 0 to 255, normalization could result in the clipping of the GSM values outside this range, thus affecting the coefficients of variation. The increased coefficients of variation for GSM, in the case of the NORM2 normalisation, provided evidence that the frame-by-frame variations seen in GSM were unlikely to be, at least significantly, due to a general temporal variability in the overall image brightness.
The extent of the GSM variations seen in the study were similar in magnitude to those reported by Elatrozy et al. , who found that, after normalisation, the coefficient of variation among 4 different observers was 4.7% for the GSM of plaques. However, the variations captured by that study did not include variations that may have been observed to the selection of different still images as each of the four observers appeared to have used the same image to assess the GSM. Therefore, the true inter-observer variabilities may have been greater than that suggested by the results of that study.
The findings of our study have two implications. Firstly, in the case of studies which have already considered intra/inter-observer variabilities, the variabilities found may have been under-estimated unless the same image frames were not used in multiple assessments of the GSM. Secondly and conversely, in the case of inter/intra-sonographer or across-study variabilities of GSM measurements, some of these variabilities may have been due to the selection of different image frames corresponding to the differences in the exact cross-sections being imaged.
The variations seen in the plaque GSM and area may be due to a number of different factors. While changes to the acquisition settings during a single acquisition would not be expected, changes in the plaque GSM could occur, for example, if the distance between the plaque and the transducer face changed during an acquisition. Patient or probe motion may also change the location and orientation of the scan plane with respect to the plaque being imaged, affecting both the measured GSM and the observed cross-sectional area. These are likely to be significant contributors to the variations seen in the GSM and area of plaques in this study. Deformation or compression of the plaque under the pulse pressures may also cause changes in the measured plaque GSM and cross-sectional area and the observation of periodical variations with physiologically reasonable frequencies in the plaque GSM and area for several of the plaques provided evidence to support this hypothesis. However, it should be noted that such cyclic variations could also have been caused by periodic variations in the scan plane location and orientation due to out-of-plane plaque, patient, or probe motion. The poor correlation between the inter-frame coefficients of variation of GSM and plaque area suggested that at least some of the variations seen in the plaque GSM were likely to have occurred due to factors other than changes in the observed plaque area. This was also supported by the results of the regression analysis, which did not highlight the parameters based on the plaque area as being statistically significant contributors to the variabilities seen, across plaque samples, in the extent of the inter-frame GSM variations.
Other factors that may have caused apparent changes in the plaque GSM and cross-sectional area included unclear plaque boundaries (e.g. poor quality image or substantial image noise), which may have caused fluctuations in the detected plaque boundaries. However, the images used in this study were of sufficiently good quality that any variations due to such fluctuations were not thought to be a major contributor to the GSM variations observed.
The statistically significant difference found in the mean GSM of plaques in the symptomatic and asymptomatic groups is in accord with previous findings that have shown symptomatic plaques, in general, had lower GSM [3, 4]. The differences between the two groups in the case of the parameters describing the inter-frame variations in the GSM were not statistically significant which was plausible as out-of-plane plaque, patient, and probe motion appeared to be a significant sources of variation for GSM measurements.
Manual delineation of the plaque boundaries separately for each image frame was found to increase the extent of the frame-by-frame variations observed in the plaque GSM and area (7.7% and 8.0%, respectively, compared with 5.4% and 4.0% for our method) due to the greater subjectivity of the manual delineation process.
The main limitations of our study were the use of two dimensional ultrasound and the absence of any attempts to fix the scan plane location and orientation with respect to the plaque being imaged, other than those measures normally taken in the clinics (e.g. holding the probe fixed and asking the patient to remain still). It should be remembered that the method of ultrasound acquisition commonly used in carotid plaque GSM studies, namely two dimensional ultrasound, provides only a cross section of the whole plaque volume. Since the plaque GSM measured using two dimensional techniques reflects only a cross-section, these measurements are susceptible to variations due to out-of-plane plaque, probe and patient motion. Studies incorporating three dimensional techniques may overcome this limitation and enable further investigation of the nature of any intrinsic frame-by-frame variations in the plaque GSM. Such follow-up studies may also identify whether any inter-frame variations seen in the GSM and volume of plaques can provide additional insight into the dynamic behaviour of carotid plaques, thus improving clinical utility.
Another limitation of our assessment of the plaque GSM and cross-sectional area was with regard to anechoic regions of plaques and regions of plaques in areas of acoustic shadowing. These regions were excluded from analysis. The cross-sectional areas of plaques that had such regions were, thus, under-estimated and neither the plaque area nor the GSM reflected the true values. The excluded regions also had an effect on the magnitude of the frame-by-frame variations that were measured for the affected plaques. In the case of anechoic regions of plaques, the inclusion of the anechoic regions would have reduced the magnitude of the variations observed. However, this would have been only because of the absence of echogenicity in these regions. In the case of the regions of acoustic shadowing, these regions need to be excluded from analysis due to the absence or deterioration of plaque texture information resulting from acoustic shadowing. Although Colour Doppler is useful for subjectively defining the plaque-lumen boundary, it is not suitable for quantifying the plaque area throughout the cardiac cycle, since colour filling of the lumen is dependent on blood flow velocity . Nevertheless, our results demonstrated that variations in the plaque GSM and cross-sectional area were observed, in the visible parts of the plaques that were not in regions of acoustic shadowing. Such variations are important as they could lead to an error in a potential diagnostic test that uses the GSM as the selection criterion, particularly for plaques of intermediate echogenicity where a coefficient of variation of 5% may provide enough bias to move the plaque between the high-risk and low-risk groups. As the plaque GSM is not generally used to inform clinical decision making, these variations do not currently have a clinical impact. Nevertheless, the variations should be appreciated for research studies which increasingly utilize the plaque GSM.
Our results did not find the cross-sectional plaque area to affect the extent of the frame-by-frame variations observed in the plaque GSM significantly. A major source of variation in the inter-frame plaque GSM may in fact be the movement of the plaque cross-section with respect to the scan plane and this may be a bigger problem for smaller plaques. However, our results showed that the observed GSM could vary on a frame-by-frame basis substantially for large plaques as well the minor stenoses.
Since previous studies typically quantified GSM in single frames of ultrasound images, the variations found in this study have been previously neglected. Improved attempts to standardise GSM measurements and reduce variability between centres should account for these findings, for example, by performing GSM measurements at peak systolic/diastolic frames or by carrying out an assessment of the average GSM throughout the cardiac cycle. Techniques such as GSM assessment using multiple cross-sectional views of plaques and plaque texture analysis can also be used to improve diagnostic precision compared to a single cross-sectional assessment. The best option would be to carry out these processes in three dimensions, however, three dimensional ultrasound techniques are still under development and are not widely available. It should be noted that we do not propose the technique we have used in our study as a replacement for other methods but we highlight the variations in plaque GSM and cross-sectional area that may be observed on a frame-by-frame basis using single-view, two dimensional ultrasound.
In conclusion, this investigation found that the GSM of carotid artery plaques can vary when measured on a frame-by-frame basis throughout ultrasound image sequences. These variations affect the reproducibility of studies and have implications for the use of GSM as a predictor of cerebrovascular events. Future studies looking at the GSM of carotid artery plaques may need to take these variations into consideration.
The research was funded by and took place at the National Institute for Health Research (NIHR) Collaboration for Leadership in Applied Health Research and Care based at the University Hospitals of Leicester NHS Trust. The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health.
- European Carotid Surgery Trialists’ Collaborative Group: MRC European Carotid Surgery Trial: interim results for symptomatic patients with severe (70–99%) or with mild (0–29%) carotid stenosis. Lancet. 1991, 337: 1235-1243. 10.1016/0140-6736(91)92916-P.View Article
- North American Symptomatic Carotid Endarterectomy Trial Collaborators: Beneficial effect of carotid endarterectomy in symptomatic patients with high-grade stenosis. N Engl J Med. 1991, 325: 445-453.View Article
- Elatrozy T, Nicolaides A, Tegos T, Griffin M: The objective characterisation of ultrasonic carotid plaque features. Eur J Vasc Endovasc Surg. 1998, 16 (3): 223-230. 10.1016/S1078-5884(98)80224-4.View ArticlePubMed
- Tegos TJ, Stavropoulos P, Sabetai MM, Khodabakhsh P, Sassano A, Nicolaides AN: Determinants of carotid plaque instability: echoicity versus heterogeneity. Eur J Vasc Endovasc Surg. 2001, 22 (1): 22-30. 10.1053/ejvs.2001.1412.View ArticlePubMed
- Biasi GM, Sampaolo A, Mingazzini P, DeAmicis P, El-Barghouty N, Nicolaides AN: Computer analysis of ultrasonic plaque echolucency in identifying high risk carotid bifurcation lesions. Eur J Vasc Endovasc Surg. 1999, 17 (6): 476-479. 10.1053/ejvs.1999.0789.View ArticlePubMed
- El-Barghouty N, Geroulakos G, Nicolaides A, Androulakis A, Bahal V: Computer-assisted carotid plaque characterisation. Eur J Vasc Endovasc Surg. 1995, 9: 389-393. 10.1016/S1078-5884(05)80005-X.View ArticlePubMed
- Sabetai MM, Tegos TJ, Nicolaides AN, Dhanjil S, Pare GJ, Stevens JM: Reproducibility of computer-quantified carotid plaque echogenicity: can we overcome the subjectivity?. Stroke. 2002, 31: 2189-2196.View Article
- Salem MK, Sayers RD, Bown MJ, West K, Moore D, Nicolaides A, Robinson TG, Naylor AR: Patients with recurrent ischaemic events from carotid artery disease have a large lipid core and Low GSM. Eur J Vasc Endovasc Surg. 2012, 43: 147-153. 10.1016/j.ejvs.2011.11.008.View ArticlePubMed
- Grønholdt MLM, Nordestgaard BG, Schroeder TV, Vorstrup S, Sillesen H: Ultrasonic echolucent carotid plaques predict future strokes. Circulation. 2001, 104: 68-73. 10.1161/hc2601.091704.View ArticlePubMed
- Polak JF, Shemanski L, O’Leary DH, Lefkowitz D, Price TR, Savage PJ, Brant WE, Reid C: Hypoechoic plaque at US of the carotid artery: an independent risk factor for incident stroke in adults aged 65 years or older. Radiology. 1998, 208: 649-654.View ArticlePubMed
- Biasi GM, Froio A, Diethrich EB, Deleo G, Galimberti S, Mingazzini P, Nicolaides AN, Griffin M, Raithel D, Reid DB, Valsecchi MG: Carotid plaque echolucency increases the risk of stroke in carotid stenting: the imaging in carotid angioplasty and risk of stroke (ICAROS) study. Circulation. 2004, 110: 756-762. 10.1161/01.CIR.0000138103.91187.E3.View ArticlePubMed
- Elatrozy T, Nicolaides A, Tegos T, Zarka AZ, Griffin M, Sabetai M: The effect of B-mode image standardisation on the echodensity of symptomatic and asymptomatic carotid bifurcation plaques. Int J Angiol. 1998, 17: 179-186.
- Oates CP, Naylor AR, Hartshorne T, Charles SM, Fail T, Humphries H, Aslam M, Khodabaksh P: Joint recommendations for reporting carotid ultrasound investigations in the united kingdom. Eur J Vasc Endovasc Surg. 2009, 37 (3): 251-261. 10.1016/j.ejvs.2008.10.015.View ArticlePubMed
- Grant EG, Benson CB, Moneta GL, Alexandrov AV, Baker JD, Bluth EI, Carroll BA, Eliasziw M, Gocke J, Hertzberg BS, Katanick S, Needleman L, Pellerito J, Polak JF, Rholl KS, Wooster DL, Zierler RE: Carotid artery stenosis: gray-scale and Doppler US diagnosis—society of radiologists in ultrasound consensus conference. Radiology. 2003, 229 (2): 340-346. 10.1148/radiol.2292030516.View ArticlePubMed
- Jansen M, van Alfen N, van der Sanden MWGN, van Dijk JP, Pillen S, de Groot IJM: Quantitative muscle ultrasound is a promising longitudinal follow-up tool in Duchenne muscular dystrophy. Neuromuscul Disord. 2012, 22: 306-317. 10.1016/j.nmd.2011.10.020.View ArticlePubMed
- Kanber B, Ramnarine KV: A probabilistic approach to computerized tracking of arterial walls in ultrasound image sequences. ISRN Signal Processing. 2012, 10.5402/2012/179087.
- Ramnarine KV, Kanber B, Panerai RB: Assessing the performance of vessel wall tracking algorithms: the importance of the test phantom. J Phys Conf Ser. 2004, 1: 199-204.View Article
- Lewis JP: Fast Template Matching. 1995, Quebec City, Canada: In Vision Interface 95, Canadian Image Processing and Pattern Recognition Society, 120-123.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.