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

Ophthalmic artery reversal predicts contralateral body weakness symptoms better than carotid Doppler velocity-preliminary results

Background

The measurement of blood velocity in the carotid artery has been the most popular noninvasive method of identifying and classifying carotid stenosis for half a century. Carotid stenosis is an indicator of elevated risk of stroke; anatomic revascularization reduces the chance of stroke by more than half. Controversy persists on how patients with severe carotid stenosis should be selected for anatomic revascularization. Patients with a connected circle of Willis (coW) might not benefit from anatomic revascularization; patients with two segments missing in the coW are most likely to benefit from revascularization.

Methods

Based on this analysis of data from carotid duplex examinations and transcranial Doppler examinations including ophthalmic artery (OA) direction in 28 patients, a refined carotid examination protocol is proposed. This refinement includes Doppler measurement of OA flow direction and documentation of internal carotid artery (ICA) bruit in addition to the adoption of an ICA peak systolic velocity (PSV) criterion exceeding 350 cm/s for identification of the patient most likely to benefit from carotid stenosis treatment.

Results

Sensitivity and specificity of OA direction or carotid bruit are 84.6%±5.4%, 71.4%±2.1% and for PSV >350 cm/s are 84.6%±5.4%, 59.5%±2.3% for predicting contralateral body weakness.

Conclusions

The proposed examination can be performed with the same duplex scanner and scan head currently used for carotid examinations with little additional time.

Arteriovenous Fistula Flow Dysfunction Surveillance: Early Detection Using Pulse Radar Sensor and Machine Learning Classification

Vascular Access (VA) is often referred to as the “Achilles heel” for a Hemodialysis (HD)-dependent patient. Both the patent and sufficient VA provide adequacy for performing dialysis and reducing dialysis-related complications, while on the contrary, insufficient VA is the main reason for recurrent hospitalizations, high morbidity, and high mortality in HD patients. A non-invasive Vascular Wall Motion (VWM) monitoring system, made up of a pulse radar sensor and Support Vector Machine (SVM) classification algorithm, has been developed to detect access flow dysfunction in Arteriovenous Fistula (AVF). The harmonic ratios derived from the Fast Fourier Transform (FFT) spectrum-based signal processing technique were employed as the input features for the SVM classifier. The result of a pilot clinical trial showed that a more accurate prediction of AVF flow dysfunction could be achieved by the VWM monitor as compared with the Ultrasound Dilution (UD) flow monitor. Receiver Operating Characteristic (ROC) curve analysis showed that the SVM classification algorithm achieved a detection specificity of 100% at detection thresholds in the range from 500 to 750 mL/min and a maximum sensitivity of 95.2% at a detection threshold of 750 mL/min.

Quest for the Vulnerable Atheroma: Carotid Stenosis and Diametric Strain--A Feasibility Study

The Bernoulli effect may result in eruption of a vulnerable carotid atheroma, causing a stroke. We measured electrocardiography (ECG)-registered QRS intra-stenotic blood velocity and atheroma strain dynamics in carotid artery walls using ultrasonic tissue Doppler methods, providing displacement and time resolutions of 0.1 μm and 3.7 ms. Of 22 arteries, 1 had a peak systolic velocity (PSV) >280 cm/s, 4 had PSVs between 165 and 280 cm/s and 17 had PSVs <165 cm/s. Eight arteries with PSVs <65 cm/s and 4 of 9 with PSVs between 65 and 165 cm/s had normal systolic diametric expansion (0% and 7%) and corresponding systolic wall thinning. The remaining 10 arteries had abnormal systolic strain dynamics, 2 with diametric reduction (>-0.05 mm), 2 with extreme wall expansion (>0.1 mm), 2 with extreme wall thinning (>-0.1 mm) and 4 with combinations. Decreases in systolic diameter and/or extreme systolic arterial wall thickening may indicate imminent atheroma rupture.

Ultrasound speckle tracking strain estimation of in vivo carotid artery plaque with in vitro sonomicrometry validation

Our objective was to validate a previously developed speckle tracking (ST) algorithm to assess strain in common carotid artery plaques. Radial and longitudinal strain was measured in common carotid artery gel phantoms with a plaque-mimicking inclusion using an in-house ST algorithm and sonomicrometry. Moreover, plaque strain by ST for seven patients (77 ± 6 y) with carotid atherosclerosis was compared with a quantitative visual assessment by two experienced physicians. In vitro, good correlation existed between ST and sonomicrometry peak strains, both radially (r = 0.96, p < 0.001) and longitudinally (r = 0.75, p < 0.01). In vivo, greater pulse pressure-adjusted radial and longitudinal strains were found in echolucent plaques than in echogenic plaques. This illustrates the feasibility of ultrasound ST strain estimation in plaques and the possibility of characterizing plaques using ST strain in vivo.

Effects of parameters on the accuracy and precision of ultrasound-based local pulse wave velocity measurement: a simulation study

Quantification of arterial stiffness, such as pulse wave velocity (PWV), is increasingly used in the risk assessment of cardiovascular disease. Pulse wave imaging (PWI) is an emerging ultrasound-based technique to noninvasively measure the local PWV instead of the global PWV, as in conventional methods. In PWI, several key parameters, including the frame rate of ultrasound imaging, motion estimation rate (MER), number of scan lines, image width, PWV value, and sonographic signal-to-noise ratio (SNRs), play an important but still unclear role in the accuracy and precision of PWV measurement. In this study, computer simulations were performed to investigate the fundamental effects of these parameters on the PWV measurement. The pulse waveform was estimated by speckle tracking on ultrasound RF signals acquired at a frame rate of 2083 Hz from a location on the common carotid artery of a healthy subject. By applying different time delays on the estimated waveform based on specific PWI parameters, the pulse waveforms at others locations were simulated. Ultrasound RF signals of the artery during the pulse wave propagation were generated from a 2-D convolutional image formation model. The PWI technique was applied to estimate the PWV at different values of frame rate, MER, number of scan lines, image width, PWV, and SNRs. The performance of the PWV estimation was evaluated by measuring the relative error, coefficient of variation (CV) and coefficient of determination (R(2)). The results showed that PWVs could be correctly measured when the frame rate was higher than a certain value (i.e., minimum frame rate), below which the estimated error increased rapidly. The minimum frame rate required for PWV estimation was found to increase with the value of PWV. An optimal MER was found (i.e., about 200 Hz) and allowed better performance of PWV measurement. The CV of PWV estimation decreased and R(2) increased with number of scan lines and image width, indicating that the performance of the PWV estimation could be improved with a larger number of scan lines and image width. For a given sufficiently high frame rate, a higher PWV value was found to deteriorate the PWV estimation, as indicated by an increasing CV and decreasing R(2). The simulation results were in good agreement with the theoretical analysis. Finally, high-quality PWV estimation could be obtained as long as the SNRs was higher than about 30 dB. The quantitative effects of the key parameters obtained from this study might provide important guidelines for parameter optimization in ultrasound-based local PWV measurement in vivo.

Brachial artery vasomotion and transducer pressure effect on measurements by active contour segmentation on ultrasound

PURPOSE

To use feed-forward active contours (snakes) to track and measure brachial artery vasomotion on ultrasound images recorded in both transverse and longitudinal views; and to compare the algorithm's performance in each view.

METHODS

Longitudinal and transverse view ultrasound image sequences of 45 brachial arteries were segmented by feed-forward active contour (FFAC). The segmented regions were used to measure vasomotion artery diameter, cross-sectional area, and distention both as peak-to-peak diameter and as area. ECG waveforms were also simultaneously extracted frame-by-frame by thresholding a running finite-difference image between consecutive images. The arterial and ECG waveforms were compared as they traced each phase of the cardiac cycle.

RESULTS

FFAC successfully segmented arteries in longitudinal and transverse views in all 45 cases. The automated analysis took significantly less time than manual tracing, but produced superior, well-behaved arterial waveforms. Automated arterial measurements also had lower interobserver variability as measured by correlation, difference in mean values, and coefficient of variation. Although FFAC successfully segmented both the longitudinal and transverse images, transverse measurements were less variable. The cross-sectional area computed from the longitudinal images was 27% lower than the area measured from transverse images, possibly due to the compression of the artery along the image depth by transducer pressure.

CONCLUSIONS

FFAC is a robust and sensitive vasomotion segmentation algorithm in both transverse and longitudinal views. Transverse imaging may offer advantages over longitudinal imaging: transverse measurements are more consistent, possibly because the method is less sensitive to variations in transducer pressure during imaging.

Arterial mechanical motion estimation based on a semi-rigid body deformation approach

Arterial motion estimation in ultrasound (US) sequences is a hard task due to noise and discontinuities in the signal derived from US artifacts. Characterizing the mechanical properties of the artery is a promising novel imaging technique to diagnose various cardiovascular pathologies and a new way of obtaining relevant clinical information, such as determining the absence of dicrotic peak, estimating the Augmentation Index (AIx), the arterial pressure or the arterial stiffness. One of the advantages of using US imaging is the non-invasive nature of the technique unlike Intra Vascular Ultra Sound (IVUS) or angiography invasive techniques, plus the relative low cost of the US units. In this paper, we propose a semi rigid deformable method based on Soft Bodies dynamics realized by a hybrid motion approach based on cross-correlation and optical flow methods to quantify the elasticity of the artery. We evaluate and compare different techniques (for instance optical flow methods) on which our approach is based. The goal of this comparative study is to identify the best model to be used and the impact of the accuracy of these different stages in the proposed method. To this end, an exhaustive assessment has been conducted in order to decide which model is the most appropriate for registering the variation of the arterial diameter over time. Our experiments involved a total of 1620 evaluations within nine simulated sequences of 84 frames each and the estimation of four error metrics. We conclude that our proposed approach obtains approximately 2.5 times higher accuracy than conventional state-of-the-art techniques.

Shear wave elastography assessment of carotid plaque stiffness: in vitro reproducibility study

This study assessed inter- and intra-observer reproducibility of shear wave elastography (SWE) measurements in vessel phantoms simulating soft and hard carotid plaque under steady and pulsatile flow conditions. Supersonic SWE was used to acquire cine-loop data and quantify Young's modulus in cryogel vessel phantoms. Data were acquired by two observers, each performing three repeat measurements. Mean Young's modulus was quantified within 2-mm regions of interest averaged across five frames and, depending on vessel model and observer, ranged from 28 to 240 kPa. The mean inter-frame coefficient of variation (CV) was 0.13 (range: 0.07-0.18) for observer 1 and 0.14 (range: 0.12-0.16) for observer 2, with mean intra-class correlation coefficients (ICCs) of 0.84 and 0.83, respectively. The mean inter-operator CV was 0.13 (range: 0.08-0.20), with a mean ICC of 0.76 (range: 0.69-0.82). Our findings indicate that SWE can quantify Young's modulus of carotid plaque phantoms with good reproducibility, even in the presence of pulsatile flow.

Wall motion in the stenotic carotid artery: association with greyscale plaque characteristics, the degree of stenosis and cerebrovascular symptoms

Background

Systolic dilation of the atherosclerotic carotid artery depends on several factors including arterial compliance and the haemodynamic environment. The purpose of this study was to quantify wall motion in stenotic carotid arteries and investigate any associations with the ultrasound greyscale plaque characteristics, the degree of stenosis, and the presence of cerebrovascular symptoms.

Methods

Variations in the lumen diameters of 61 stenotic carotid arteries (stenosis range 10%-95%) from 47 patients were measured before the proximal shoulder of the atherosclerotic plaque using ultrasound image sequences over several cardiac cycles. Absolute and percentage diameter changes from diastole to systole were calculated and their relationship to the degree of stenosis, greyscale plaque characteristics, and the presence of ipsilateral hemispheric symptoms were studied.

Results

The mean absolute diameter change from diastole to systole was 0.45 mm (s.d. 0.17), and the mean percentage diameter change was 6.9% (s.d. 3.1%). Absolute and percentage diameter changes did not have a statistically significant relationship to the degree of stenosis, greyscale plaque characteristics, or the presence of ipsilateral hemispheric symptoms (p > 0.05). Parameters significantly correlated with the presence of symptoms were the degree of stenosis (p = 0.01), plaque greyscale median (p = 0.02) and the plaque surface irregularity index (p = 0.02).

Conclusions

Our study confirmed the degree of stenosis, plaque greyscale median and our surface irregularity index were significant predictors of symptoms, but found no significant correlation between diameter changes of stenosed carotid arteries and the presence of ipsilateral hemispheric symptoms.

BIOMECHANICAL INVESTIGATION OF PULSATILE FLOW IN A THREE-DIMENSIONAL ATHEROSCLEROTIC CAROTID BIFURCATION MODEL

It is a well-established fact that atherosclerosis in carotid bifurcation depends on flow parameters such as wall shear stress, flow pulsatility, and blood pressure. However, it is still not clearly verified how atherosclerosis can become aggravated when plaque experiences a high level of shear stress during advance stages of this disease. In this paper, fluid and structural properties in idealistic geometries are analyzed by using fluid-structure interaction (FSI). From our results, the relationship among blood pressure, stenotic compression, and deformation was established. We show that a high level of compression occurs at the stenotic apex, and can potentially be responsible for plaque progression. Moreover, wall shear stress and deformation are significantly affected by the degree of stenosis. Finally, through analysis of the FSI-based simulation results, we can better understand the parameters that influence flow through a stenotic artery and plaque aggravation, and apply the knowledge for the enhancement of clinical research and prediction of treatment outcomes.

Identification of high-risk carotid artery stenosis: motion of intraplaque contents detected using B-mode ultrasonography

Object

Identification of the risk of rupture and vulnerability of arterial plaque is not yet clearly understood. The aim of this study was to assess the clinical features of the motion of intraplaque contents (MIC) detected by B-mode ultrasonography. The MIC is characterized by the peculiar movement of the intraplaque contents that is not synchronized with the heartbeat; however, the movement of the carotid artery (CA) wall depends on the heartbeat.

Methods

From January 2008 to November 2010, 1798 consecutive patients with transient ischemic attacks (TIAs) or acute ischemic stroke underwent CA ultrasonography for the examination of the MIC. Patients with CA stenosis greater than 50% were followed up until they underwent carotid endarterectomy or CA angioplasty and stent placement. If neither of these procedures were used, the patients were followed up at 90 days. Chi-square and Mann-Whitney tests were performed to compare the categorical and continuous demographic data and risk factors. The effect of the MIC on the rate of recurrent cerebral ischemia was examined using Kaplan-Meier and univariate Cox regression analyses.

Results

One hundred and fifteen patients had CA stenosis greater than 50%. Among these 115 patients, 58 with a total of 59 CA stenoses had MIC. Twenty-four recurrent ischemic events were associated with MIC, whereas only 6 such events occurred in the absence of MIC. The MIC decreased event-free survival (log-rank test = 15.8, p < 0.001); univariate Cox analysis confirmed that MIC increased the risk of a recurrent ischemic event (HR 5.12, 95% CI 2.08–12.58; p < 0.001).

Conclusions

The MIC is one of the findings of vulnerable plaques. The MIC is more useful in predicting the recurrence of TIAs or ischemic events in patients with symptomatic CA stenosis.

Dependence of cyclic stretch-induced stress fiber reorientation on stretch waveform

Cyclic uniaxial stretching of adherent nonmuscle cells induces the gradual reorientation of their actin stress fibers perpendicular to the stretch direction to an extent dependent on stretch frequency. By subjecting cells to various temporal waveforms of cyclic stretch, we revealed that stress fibers are much more sensitive to strain rate than strain frequency. By applying asymmetric waveforms, stress fibers were clearly much more responsive to the rate of lengthening than the rate of shortening during the stretch cycle. These observations were interpreted using a theoretical model of networks of stress fibers with sarcomeric structure. The model predicts that stretch waveforms with fast lengthening rates generate greater average stress fiber tension than that generated by fast shortening. This integrated approach of experiment and theory provides new insight into the mechanisms by which cells respond to matrix stretching to maintain tensional homeostasis.

Dynamic assessment of carotid plaque motion

Assessment of dynamic plaque behaviour may help identify vulnerable carotid plaque before rupture and hence has potential clinical value for screening patients at risk of stroke. The aim of this study was to develop non-invasive ultrasound methods for quantifying dynamic plaque and vessel wall behaviour and assess their potential clinical utility. Ultrasound data from the carotid arteries of one normal subject and four patients with atherosclerotic disease were acquired using a 10 MHz linear array transducer recording raw RF/IQ data at a frame rate up to 80 Hz for 3–6 seconds. Image reconstruction and processing was performed using Matlab. Speckle tracking techniques were developed to characterize: (1) intraplaque deformation; and (2) plaque surface and vessel wall motion. Speckle tracking techniques were able to measure the range of intraplaque tissue deformation (–1.3 to 1.7 mm), plaque surface displacement (0.2–0.7 mm) and vessel wall radial strain (0.02–0.13) throughout the cardiac cycle. The feasibility of using an intraplaque deformation parameter, based on the deformation of a square template, is demonstrated. Speckle tracking techniques can be used to assess dynamic carotid plaque behaviour. Further work is required to evaluate how best to quantify biomechanical behaviour to help predict plaque rupture and hence improve risk stratification models for stroke.

Assessment of tissue Doppler imaging measurements of arterial wall motion using a tissue mimicking test rig

The aim of this in vitro study is to assess the accuracy of the tissue Doppler imaging arterial wall motion (TDI AWM) technique in measuring dilation over a range of distances and velocities. A test rig, consisting of two parallel blocks of tissue mimicking material (TMM), has been developed to generate known wall motion. One block remains stationary while the other moves in a cyclical motion. A calibrated laser range finder was used to measure the TMM motion. The TDI AWM measurements were found to underestimate the dilation by 21% +/- 4.7% when using the recommended scanner parameters. The size of the error was found to increase with a decrease in ultrasound output power. Results suggested that errors in the TDI AWM dilation measurements relate to underestimates in the velocity measured by the TDI technique. The error demonstrated in this study indicates a limitation in the value of TDI AWM result obtained in vivo. (E-mail: abigail.thrush@bartsandthelondon.nhs.uk).

Derivation of the distensibility coefficient using tissue Doppler as a marker of arterial function

To date, the main cardiovascular application of TDI (tissue Doppler imaging) has been in myocardial evaluation. In the present study, we investigated the feasibility and reproducibility of assessing arterial elasticity using the DC (distensibility coefficient) measured by TDI, the correlation of this with the DC obtained by other methods and the DC in patients with various degrees of cardiovascular risk. We studied 450 subjects (256 men; age, 51±10 years) with and without risk factors of cardiovascular disease. Arterial displacement was measured from TDI, and B-mode and M-mode images of the common carotid artery in the longitudinal plane, and the DC with each method was compared. Linear regression showed a good correlation between all three methods. The results for TDI and B-mode were comparable [(21±10) compared with (21±10)×10−3/kPa respectively; P=not significant], but there were significant differences between TDI and M-mode [(21±10) compared with (31±13)×10−3/kPa respectively; P<0.0001] and between B-mode and M-mode [(21±10) compared with (31±13)×10−3/kPa respectively; P<0.0001]. Similarly, Bland–Altman analysis showed the least variability in the DC between TDI and B-mode, and there were no significant differences between the average measurements. The TDI DC also had the lowest paired difference for inter-observer variability [(−0.1±1.1)×10−3/kPa; P=not significant]. In conclusion, the results of the present study suggest that TDI of the carotid arteries is feasible, comparable with B-mode measurements, more robust than M-mode and less variable than the other methods.

Noninvasive Young's modulus evaluation of tissues surrounding pulsatile vessels using ultrasound Doppler measurement

This paper presents an indirect approach to estimating the mechanical properties of tissues surrounding the arterial vessels using ultrasound (US) Doppler measurements combined with an inverse problem-solving method. The geometry of the structure and the dynamic behavior of the inner fluid are first evaluated using a novel dual-beam US system. A numerical phantom associated with a parametric finite element simulator that calculates the hydrodynamic pressure and the displacement on the walls' boundaries is then built. The simulation results are iteratively compared to the US measurement results to deduce the value of the unknown parameters, i.e., the Young's modulus and the pressure resulting from the downstream load. The feasibility of the proposed approach was experimentally tested in vitro using a phantom composed of a latex tube surrounded by a cryogel tissue-mimicking material.

A new technique for assessing arterial pressure wave forms and central pressure with tissue Doppler

Background

Non-invasive assessment of arterial pressure wave forms using applanation tonometry of the radial or carotid arteries can be technically challenging and has not found wide clinical application. 2D imaging of the common carotid arteries is routinely used and we sought to determine whether arterial waveform measurements could be derived from tissue Doppler imaging (TDI) of the carotid artery.

Methods

We studied 91 subjects (52 men, age 52 ± 14 years) with and without cardiovascular disease. Tonometry was performed on the carotid artery simultaneously with pulsed wave Doppler of the LVOT and acquired digitally. Longitudinal 2D images of the common carotid artery with and without TDI were also acquired digitally and both TDI and tonometry were calibrated using mean and diastolic cuff pressure and analysed off line.

Results

Correlation between central pressure by TDI and tonometry was excellent for maximum pressure (r = 0.97, p < 0.0001). The mean differences between central pressures derived by TDI and tonometry were minimal (systolic 5.36 ± 5.5 mmHg; diastolic 1.2 ± 1.2 mmHg).

Conclusion

Imaging of the common carotid artery motion with tissue Doppler may permit acquisition of a waveform analogous to that from tonometry. This method may simplify estimation of central arterial pressure and calculation of total arterial compliance.

Noninvasive simultaneous assessment of wall shear rate and wall distension in carotid arteries

A novel technique has been developed for the noninvasive real-time simultaneous assessment of both blood velocity profile and wall displacements in human arteries. The novel technique is based on the use of two ultrasound beams, one set at optimal angle for wall motion measurements and the other for blood velocity profile measurements. The technique was implemented on a linear array probe divided into two subapertures. A modified commercial ultrasound machine and a custom PC board based on a high-speed digital signal processor was used to process the quadrature demodulated echo signals and display results in realtime. Flow phantom experiments demonstrated the validity of the technique, providing wall shear rate (WSR) estimates within 10% of the theoretical values. The system was also tested in the common carotid arteries of 16 healthy volunteers (age 30 to 53 y). Results of simultaneous diameter distension and WSR measurements were in agreement with published data.

A non-invasive study of alterations of the carotid artery with age using ultrasound images

A region-based method for measurement of arterial diameter to find out the elasticity of the vessel is proposed in this paper. Arterial segments are studied by using images obtained through ultrasound scanning in B-mode. Pulsatile changes of the common carotid artery during diastole and systole are computed. To achieve this, thinned segmentation is done by suitably adjusting the contrast of the image. The diameter changes of the artery wall from the centre of artery are calculated. Fifty-three normal subjects with age group 20-40 years are taken for measurement. Measured diameter is plotted as a graph and pulsatile changes of the artery are obtained. Since no atherosclerotic lesions are detected in the studied subjects, it is suggested that the common carotid artery is a highly compliant artery with a strong alteration of viscoelastic properties with age.

Oxygen Mass Transport in a Compliant Carotid Bifurcation Model

The purpose of the present study was to investigate oxygen mass transfer in the human carotid bifurcation, focusing on the effects of the wall compliance and flow field on the temporal variation and spatial distribution of the oxygen wall flux. Details of unsteady convective-diffusive oxygen transport were examined numerically using a compliant model of the human carotid bifurcation and realistic blood flow waveforms. Results reveal that axial flow separation at the outer common-internal carotid wall can significantly alter the flow field, oxygen tension field, and oxygen wall flux distribution. At the outer wall of the sinus, the Sherwood number, Sh (non-dimensional oxygen wall flux), takes on significantly lower values than at other sites due to the attenuation of transport rates by convective flow away from wall. More specifically, the lowest value of Sh was Sh approximately 6 (in the sinus), which is much lower than the value of the non-dimensional oxygen consumption rate (Damkohler number, Da) in the reactive wall tissue (Da=29-39). At the inner wall of the sinus, Sh approximately 170 is far above the expected value of Da. This implies that flow separation on the outer wall of the sinus provides a very strong fluid mechanical barrier to oxygen transport; whereas at the inner wall of the sinus, the mechanism of transport is controlled by the wall consumption rate.

A Computational Study of Flow in a Compliant Carotid Bifurcation–Stress Phase Angle Correlation with Shear Stress

The present study presents a three-dimensional, unsteady supercomputer simulation of the coupled fluid-solid interaction problem associated with flow through a compliant model of the bifurcation of the common carotid artery into the internal and external carotid arteries. The fluid wall shear stress (WSS) and solid circumferential stress/strain (CS) are computed and analyzed for the first time using the complex ratio of CS to WSS (CS/WSS). This analysis reveals a large negative phase angle between CS and WSS (stress phase angle--SPA) on the outer wall of the carotid sinus where atherosclerotic plaques are localized. This finding is consistent with other measurements and computations of the SPA in coronary arteries and the aortic bifurcation that show large negative SPA correlating with sites of plaque location and in vitro studies of endothelial cells showing that large negative SPA induces pro-atherogenic gene expression and metabolite release profiles.