Tissue Doppler imaging of carotid plaque wall motion: a pilot study

Background Studies suggest the physical and mechanical properties of vessel walls and plaque may be of clinical value in the diagnosis and treatment of cardiovascular atherosclerotic disease. The purpose of this pilot study was to investigate the potential clinical application of ultrasound Tissue Doppler Imaging (TDI) of Arterial Wall Motion (AWM) and to quantify simple wall motion indices in normal and diseased carotid arteries. Methods 224 normal and diseased carotid arteries (0–100% stenoses) were imaged in 126 patients (age 25–88 years, mean 68 ± 11). Longitudinal sections of the carotid bifurcation were imaged using a Philips HDI5000 scanner and L12-5 probe under optimized TDI settings. Temporal and spatial AWMs were analyzed to evaluate the vessel wall displacements and spatial gradients at peak systole averaged over 5 cardiac cycles. Results AWM data were successfully extracted in 91% of cases. Within the carotid bifurcation/plaque region, the maximum wall dilation at peak systole ranged from -100 to 750 microns, mean 335 ± 138 microns. Maximum wall dilation spatial gradients ranged 0–0.49, mean 0.14 ± 0.08. The AWM parameters showed a wide variation and had poor correlation with stenoses severity. Case studies illustrated a variety of pertinent qualitative and quantitative wall motion features related to the biophysics of arterial disease. Conclusion Our clinical experience, using a challenging but realistic imaging protocol, suggests the use of simple quantitative AWM measures may have limitations due to high variability. Despite this, pertinent features of AWM in normal and diseased arteries demonstrate the potential clinical benefit of the biomechanical information provided by TDI.


Background
Cardiovascular disease accounts for around half the deaths in the Western world, with stroke being the third leading cause of death and disability worldwide [1]. The predictive ability to identify which patients will have a stroke is poor. Current clinical practice for selecting which patients will have a carotid endarterectomy, to reduce the risk of stroke, relies on assessment of the degree of arterial lumen narrowing. Large European [2] and North American [3] clinical trials have demonstrated the benefit of carotid endarterectomy on groups of patients with severe stenosis. However the degree of stenosis is a poor predictor of individual stroke risk, and treatment of asymptomatic patients remains controversial [4]. Recent studies on improved predictors of stroke risk have focused on the assessment of plaque morphology, composition, and the mechanical properties of the vessel wall and plaque deposit [5][6][7][8][9][10][11][12].
Researchers have long hoped that the physical and mechanical properties of arteries may be of clinical value in the diagnosis and treatment of cardiovascular disease [13,14]. The carotid artery is a useful window for cardiovascular risk, with ultrasound able to provide unique information on morphological, hemodynamic and elastic properties [14]. Numerous studies have shown that the elastic properties of arteries may be useful in the study of cardiovascular disease. The localized mechanical properties of carotid plaque and the vessel wall may be of particular value in identifying plaque vulnerable to rupture. The majority of clinical wall motion studies have derived elasticity parameters from measurements of the wall motion at a single location in vessels. Vessels have included the common carotid artery [15,16], femoral artery [15], brachial artery [17,18], cerebral arteries [19] and aorta [20,21].
Various ultrasound techniques have been used to detect and track the vessel wall motion. Computational techniques have mainly been based on analysis of the B-mode greyscale images [18,22,23], M-mode [16], analyses of the raw RF ultrasound data [24] and Doppler techniques [20,25]. Tissue Doppler Imaging (TDI) is a relatively new commercial technique, originally developed for imaging the myocardium, which has found increasing applications in echocardiography. TDI is basically a colour Doppler technique that has been optimised to provide images of tissue motion rather than blood flow, and is capable of high spatial and temporal resolution. The signal processing techniques employed to extract the tissue velocity information from the RF ultrasound data are typically based on time domain cross-correlation techniques [26] or autocorrelation techniques [27]. Clinical applications of TDI of the myocardium, including theoretical, pathophysiological and methodological considerations are reviewed elsewhere [28]. TDI has been applied to image arterial wall motion patterns of the common carotid artery (CCA) and to assess CCA stiffness by calculating distensibility as the change in normalised mean systolic distension relative to the change in blood pressure, and compliance as the change in cross-sectional area relative to pulse pressure [29].
Atherosclerotic plaque typically develops around the region of the internal and common carotid artery bifurcation. Knowledge of the differential wall motion displacements across the carotid bifurcation in normal and in diseased arteries may thus be useful in the study of atherosclerotic disease. An important clinical application is the identification of high-risk plaques that are vulnerable to rupture, leading to stroke. It is hoped the vulnerable plaque may be predicted by identifying pertinent features of the carotid plaque wall motion behaviour throughout the cardiac cycle. Initial clinical studies, using a TDI technique that was optimised to study carotid arterial wall motion (AWM) highlighted a number of case examples of carotid plaque disease suggesting the potential clinical use of AWM imaging [30].
The purpose of this study was to investigate the possible clinical value of these recent TDI techniques for AWM imaging in the vascular examination. In particular, this study assessed the dynamic wall motion deformation behavior of the carotid plaque and vessel wall. Arterial wall motion indices were quantified in a wide spectrum of normal and diseased carotid arteries to quantify baseline parameters, assess variability and help assess future clinical potential.

Methods
Subjects 224 carotid arteries were scanned in 126 patients and volunteers (age 25-88 years, mean 68 ± 11) with a wide spectrum of carotid disease. Patients presenting through the Leicester Royal Infirmary vascular studies unit for carotid artery studies were recruited. Thus, the majority of patients presented with medical conditions including cardiovascular disease (over 90% of subjects) and hypertension (blood pressure greater than 140 mmHg systole or 90 mmHg diastole in approximately 65% of subjects, as measured following ultrasound scan, using an Omron 705CP monitor and cuff over the brachial artery). The 4 volunteers were young and healthy, with no history of cardiovascular disease and no ultrasound evidence of atherosclerotic disease. 70 subjects were male and 56 were female with mean systolic and diastolic blood pressures of 152 ± 24 mmHg and 83 ± 13 mmHg (mean ± 1 standard deviation) respectively. In accordance with the Helsinki declaration, informed consent and Ethical approval for the study was obtained. In total, 58 carotid arteries were imaged from patients with no ultrasound evidence of atherosclerotic disease, and 166 arteries imaged with ultrasound evidence of atherosclerotic disease.

Ultrasound equipment
This study used an Arterial Wall Motion (AWM) imaging facility developed by Philips Research Laboratories using novel Tissue Doppler Imaging signal processing techniques [26]. The AWM software was able to assess, along a 3-4 cm length, the arterial wall velocities and displacements, and to derive vessel compliance, tissue elasticity and stress/strain images. This additional information on the biomechanical nature of the artery and plaque deposit is available throughout the cardiac cycle.
An ATL HDI5000 ultrasound scanner and an L12-5 linear array probe (Philips Medical Systems, Andover, US) with Tissue Doppler Imaging (TDI) mode was used to scan a longitudinal section of the vessel. Ultrasound scanning was performed by experienced vascular sonographers. Ultrasound data was collected from diseased carotid arteries at the location of stenotic plaque at the bifurcation of the common carotid artery (CCA) and the proximal internal carotid artery (ICA). In patients with no ultrasound evidence of stenoses, data was collected from the common carotid bifurcation including longitudinal segments of the ICA and the CCA. TDI data was collected over 6 cardiac cycles using optimal and standardized machine settings as follows: image depth adjusted to display the carotid image in longitudinal view; 2D Opt Parameter set to Pen.; Color and TDI imaging set to ON; Tissue Opt Parameter set to Gen.; Colour Line Density set to A; colour sensitivity set to Medium; colour smoothing set to Off; colour DMD set to Off; colour Persistence set to Off; no beam steering of TDI colour box; Velocity scale set to 1.4 cm/sec to avoid aliasing of wall motion data and maintain high frame rate; colour gain set to 100%.

Data analysis
The ultrasound data were stored as cineloop digital data files. Data were analyzed off-line using Philips Research Link HDILAB analysis software and the Arterial Wall Motion proprietary software developed by Philips Research Laboratories [26,30]. The software provided values for the axial wall dilations (defined as the difference between the anterior and posterior wall diameters relative to the reference diastole values) along a longitudinal segment of the artery for each image frame. Scan line spatial resolution was 0.148 mm, temporal resolution between 30-45 Hz and wall displacement resolution of a few microns. These data were exported to MATLAB™ (The MathWorks, Inc., Natick, MA, US) for data analysis. The wall displacements at peak systole were averaged over 5 cardiac cycles. The data were then divided into 3 segments corresponding to the ICA segment, the plaque region (or the carotid bulb in normal arteries) and the CCA segment. The location of the vessel boundaries were based on the ultrasound images. An additional quantitative criteria was used to precisely define the boundary based on the first derivative (gradient) of the spatial wall dilation having a value of 0. In each of the 3 regions, the maximum and minimum wall displacements were quantified and also the maximum and minimum spatial gradients along the longitudinal artery segments. Statistical analysis was performed using the unpaired Student's t-Test to determine if there is a significant difference between the means of two sample populations (p < 0.05). Figure 1 shows an example of synchronized, spatially homogeneous arterial wall motion of the common carotid artery proximal to the carotid bifurcation of a normal vessel with no ultrasound evidence of atherosclerotic disease. Figure 2 shows an example of synchronized, spatially homogeneous arterial wall motion at the carotid bifurcation of a normal vessel with no ultrasound evidence of atherosclerotic disease. See Additional file 1: AVI movie file of Figure 2. Figure 3 shows an example of inhomogeneous arterial wall motion at the carotid bifurcation of a normal vessel with no ultrasound evidence of atherosclerotic disease. See Additional file 2: AVI movie file of Figure 3. Figure 4 shows an example of arterial wall motion across a plaque with a mild 30% stenosis (diameter reduction). See Additional file 3: AVI movie file of Figure 4. Figure 5 shows an example of arterial wall motion across a plaque with 70% stenosis (diameter reduction). High dilation gradients at the plaque vessel interface and localized stiffening with constricted wall motion at the plaque are demonstrated. See Additional file 4: AVI movie file of Figure 5. Figure 6 shows an example of arterial wall motion across a plaque with a severe 80% stenosis (diameter reduction). See Additional file 5: AVI movie file of Figure 6. Figure 7 shows an interesting example of an ulcerated plaque and highlights the uncertainty over the wall tracking across the lumen in these tuning fork bifurcations when the ICA and external carotid artery (ECA) are in the same plane. The lower ICA has been selected for wall tracking by the user selecting a point within the lumen of the ICA. At the bifurcation, the wall tracking has a discontinuity across the lumen of the ECA. The example highlights the requirement for careful interpretation of the AWM data. Figure 8 shows an example of a severe stenosis in which the vessel wall is poorly defined and the reliability of the AWM tracking algorithm is uncertain. Figures 9 and 10 show the respective quantification results for the minimum and maximum wall dilations, and the wall dilation spatial gradients in each of the 3 arterial segments. Maximum and minimum wall dilations at peak systole were plotted as a function of stenosis severity. The significant difference between the maximum and minimum values of the wall dilation within the plaque highlights the inhomogeneous wall motion within the plaque, but it is also due to the criteria used to identify the vessel boundaries. The greater wall motion inhomogeniety in the plaque region is also indicated by the significantly Homogeneous arterial wall motion of a normal CCA segment higher spatial gradient values (p < 0.0001) compared with the ICA and the CCA. The variability of the wall motion indices are wide with no significant dependence on plaque severity (% stenosis). Even in arteries with no plaque, the 0% stenoses cases, the wall dilation measures vary widely. The segments of the ICA and the CCA also display wide variability in the wall motion indices. Table 1 summarises the AWM dilation measures within segments of the ICA, the plaque (carotid bulb) and the CCA. Table 2 summarises the arterial wall motion spatial gradients within segments of the ICA, the plaque (carotid bulb) and the CCA.

Discussion
This is the first study to assess the clinical feasibility of TDI of carotid plaque wall motion using subjective and quantitative measures in a wide spectrum of normal and diseased arteries. Currently, clinically useful endpoints are based on morphological and hemodynamic properties of arteries. Extension of the ultrasound vascular examination to provide clinically useful vascular endpoints associated with the mechanical behaviour of carotid plaque may also be useful for characterising plaque.
Wall motion behaviour within a straight segment of the CCA typically displays synchronized, spatially homogeneous wall motion (Figure 1). This type of homogeneous wall motion is typically observed in normal straight Homogeneous arterial wall motion at the carotid bifurcation of normal arteries  segments of individual vessels and is in accord with previous studies on the normal CCA [31]. An important potential application of measuring carotid plaque wall motion is to help identify the unstable plaque at high risk of rupture and consequent stroke. As plaque typically develops in the region of the carotid bifurcation, this study used a challenging imaging protocol that included the plaque (carotid bulb in normal arteries), the proximal segment of the ICA and the CCA. A variety of wall motion features were observed in a wide spectrum of normal and diseased arteries. The wall motion across the carotid bifurcation, even in vessels with no plaque may exhibit more complex behaviour than single vessels segments. Figure 2 shows an example of homogeneous wall motion across the bifurcation in a normal subject with no atherosclerotic disease. This is in contrast to more complex wall motion across the bifurcation (Figure 3) observed in a different subject, also with no evidence of atherosclerotic disease. This differential wall motion across the carotid bifurcation is not uncommon and represents a confounding factor in the use of our simple quantification measures for the differentiation of normal and diseased vessels. The variability may reflect true differences in wall motion between the different regions (ICA, carotid bulb and CCA). However, it is also likely, at least in part, to be caused by limitations and artefacts from various potential sources that we discus later.
The dynamic wall motion deformation behaviour appears to provide potentially useful information on the elastic Inhomogeneous arterial wall motion at the carotid bifurcation of normal torturous arteries Figure 3 Inhomogeneous arterial wall motion at the carotid bifurcation of normal torturous arteries. This shows an example of inhomogeneous arterial wall motion at the carotid bifurcation of normal arteries with no ultrasound evidence of atherosclerotic disease. Note torturous artery with significant compression within the ICA due to posterior wall displacement. See Additional file 2 for the real-time movie of the top right image. The lower right graph displays the dilation data as a surf plot with dilation plotted as a function of the spatial position and the time. These included examples of AWM constriction or stiffening associated with calcified plaque, constriction due to pressure drop through a severe stenosis (Venturi effect), high spatial gradient and intraplaque/vessel-interface inhomogeneities which were related to plaque type and severity. Pertinent features associated with clinical end points such as plaque development, plaque rupture and stroke are important and will require additional follow up studies.
It is also important to appreciate the potential problems and reliability of the technique. The clinical case examples in Figures 3, 7 and 8 showed examples of variable AWM behavior that did not appear to relate to plaque disease. These highlighted potential limitations of the technique and the use of simple quantification indices. Our simple AWM indices, which included the minimum and maximum wall dilations and spatial wall dilation gradients, showed a wide variation and generally had poor correlation with stenosis severity. The sources of error and variability may arise from various sources including the equipment, the operator, the patient and the limitations and design of our pilot study. A brief summary of some important considerations is given below.
Local arterial wall motion stiffening around a mild stenosis calcified plaque Figure 4 Local arterial wall motion stiffening around a mild stenosis calcified plaque. This shows an example of arterial wall motion across a calcified plaque with a mild 30% stenosis (diameter reduction). Note the localised stiffening and high dilation gradients. See Additional file 3 for the real-time movie of the top right image. The lower right graph displays the dilation data as a surf plot with dilation plotted as a function of the spatial position and the time.

Limitations and sources of error Equipment
The AWM system was fundamentally limited to the measurement of wall motion within the beam scan plane and in the direction of the ultrasound pulse. TDI velocity data is angle dependent and provides the velocity component in scan plane only. The transverse component (perpendicular to the beam) of carotid plaque wall motion is also significant and may be of potential clinical importance. Indeed this transverse motion of the artery, synchronized with the arterial blood flow, is often apparent from the grayscale and color Doppler images.
The TDI AWM signal processing algorithm is based on a cross correlation technique on multiple RF signals. The characteristics, performance and limitations of the algorithm from a signal processing perspective are discussed elsewhere [26,32]. The performance of the TDI technique has been evaluated using in-vitro test phantoms developed to assess the performance for measurements of the tissue motion of the myocardium [33]. Studies have also used in-vitro phantoms and flow models to assess accuracy, variability and reproducibility of the TDI AWM technique [34]. One study [34], using a commercial polyurethane flow phantom, showed the intrinsic spatial and temporal stability was 2.46% with accuracy of 5.36%. However, the applicability of these in-vitro measurements in providing a realistic assessment of the in-vivo clinical measurements is likely to depend on the characteristics and properties of the test phantoms. We have shown that use of agar based Arterial wall motion stiffening around a severe stenosis plaque Figure 5 Arterial wall motion stiffening around a severe stenosis plaque. Figure 5 shows an example of arterial wall motion across a plaque with 70% stenosis (diameter reduction). High dilation gradients at the plaque vessel interface and localized stiffening with constricted wall motion at the plaque are demonstrated. See Additional file 4 for the real-time movie of the top right image. The lower right graph displays the dilation data as a surf plot with dilation plotted as a function of the spatial position and the time. tissue mimicking wall-less flow phantoms [35], with Rayleigh scattering characteristics are likely to provide a more challenging assessment of wall motion compared to phantoms with well defined specular wall /interface reflections [36].

Operator
Variations in wall motion quantification may arise due to inter-and intra-operator variability. Experienced ultrasonographers are required to optimise wall tracking performance. Positioning of the transducer to maintain good vessel wall definition across the carotid bifurcation whilst avoiding excessive transducer pressure and potential arte-facts from the neighbouring jugular vein are important considerations.
Important factors include the orientation of the scan plane and beam vessel angle, optimization of image quality, operator steadiness, transducer contact pressure and selection of machine settings (settings were standardized in this study).

Patient
For the clinical application of identifying the vulnerable plaque, or characterizing plaque type, confounding factors and sources of variability may arise from the patient.
Complex arterial wall motion across a plaque with a severe 80% stenosis Figure 6 Complex arterial wall motion across a plaque with a severe 80% stenosis. This shows an example of arterial wall motion across a plaque with a severe 80% stenosis (diameter reduction). The high velocity jet flow hits against the posterior wall which shows marked wall displacements distal to the plaque. Within the plaque there is significant vessel wall constriction due to the pressure drop within the lumen. This is an example of the Venturi effect, well known in fluid dynamics in which the increased kinetic energy of blood flow through the constriction corresponds to a pressure drop and consequent constriction of the vessel walls. See Additional file 5 for the real-time movie of the top right image. The lower right graph displays the dilation data as a surf plot with dilation plotted as a function of the spatial position and the time. Limitations of arterial wall motion tracking across the vessel lumen. This shows an interesting example of an ulcerated plaque and highlights the uncertainty over the wall tracking across the lumen in these tuning folk bifurcations. At the bifurcation, the wall tracking of the ICA has a discontinuity across the lumen of the ECA. The example highlights the requirement for careful interpretation of the AWM data. S c a n l i n e n u m b e r Wall motion and the elastic properties of arteries have been shown to be related to the patient's age, sex, pathology and physiology [37,38]. Variations in anatomy, and breathing/motion artifacts are also likely to increase the variability of wall motion quantification measures. Geometrical deviations of the vessels from the perpendicular and variation of image quality and vessel wall definition will have a deleterious effect on the wall tracking algorithm.

Pilot Study
Our pilot study of clinical potential used a challenging but realistic imaging protocol that encompassed the carotid bifurcation. In contrast to measurements on a single Limitations of arterial wall motion tracking in severe disease Figure 8 Limitations of arterial wall motion tracking in severe disease. This shows an example of a severe stenosis in which the vessel wall is poorly defined and the reliability of the AWM tracking is uncertain.  Only the simplest of wall motion measures, the dilations and spatial gradients at peak systole were quantified. Additional parameters and more sophisticated characterization of the dynamic wall motion patterns in normal and diseased arteries may provide a more suitable and clinically useful vascular measure. Other measures, for spatial and temporal wall motion homogeneity have been described [31].
The design of our pilot study accentuated the variability as no attempt was made to discard cases due to technical or other artifacts. Only 9% of cases were discarded prior to quantification. These were due to difficulties with data collection at the carotid bifurcation or failure of the AWM software to discriminate between consecutive cardiac cycle wall motion. In order to avoid potential bias by subjective discrimination, all AWM results were pooled, including outliers with apparent artifacts. This approach also helped assess the potential for routine application of the TDI technique as an extension to the vascular examination. A further subjective assessment of the scans could reject cases that suffer from poor image quality, wall tracking problems, anatomical variations etc.
There was no gold standard for plaque characterization and the type of plaque in case examples were based on the ultrasound grayscale appearance. Further work and fol-Arterial wall dilation as a function of stenosis severity for segments of the ICA, the plaque (carotid bulb) and the CCA Figure 9 Arterial wall dilation as a function of stenosis severity for segments of the ICA, the plaque (carotid bulb) and the CCA. low-up studies are required to assess interesting characteristics of the plaque (area, homogeneity, histology, composition, mechanical and physical properties) as well as the carotid vessels (intima-media thickness, kinking, stenosis in series, stenosis of contralateral carotid bifurcation). Finally, there is no gold standard to compare the invivo equipment performance, as the spatial and temporal resolution performance capabilities are unique to the ultrasound TDI technique.

Conclusions
Our pilot study demonstrates that the TDI Arterial Wall Motion imaging technique provides an impressive biomechanical window to atherosclerotic disease. It provides potentially useful information on the dynamic behavior of normal and diseased carotid arteries. Our initial clinical experience, using a challenging but realistic protocol for imaging carotid plaque wall motion, suggests the use of simple quantitative AWM measures may have limitations due to high variability. The technique must be robust and shown to be clinically useful if it is to play a major role in the vascular examination, where it is hoped that assessment of the mechanical properties of arteries, may help to provide a more rational basis for referring patients to carotid surgery.

Competing interests
None declared.

Authors' contributions
KVR conceived of the study, participated in its design and coordination, analyzed the data and drafted the manuscript.
TH, YS, MN and JW recruited patients, collected and reviewed the ultrasound data.

ARN recruited patients.
Arterial wall dilation spatial gradients as a function of stenosis severity for segments of the ICA, the plaque (carotid bulb) and the CCA Figure 10 Arterial wall dilation spatial gradients as a function of stenosis severity for segments of the ICA, the plaque (carotid bulb) and the CCA. Peak systole max spatial gradient