An assessment of intra-patient variability on observed relationships between wall shear stress and plaque progression in coronary arteries
© Molony et al.; licensee BioMed Central Ltd. 2015
Published: 9 January 2015
Wall shear stress (WSS) has been associated with sites of plaque localization and with changes in plaque composition in human coronary arteries. Different values have been suggested for categorizing WSS as low, physiologic or high; however, uncertainties in flow rates, both across subjects and within a given individual, can affect the classification of WSS and thus influence the observed relationships between local hemodynamics and plaque changes over time. This study examines the effects of uncertainties in flow rate boundary conditions upon WSS values and investigates the influence of this variability on the observed associations of WSS with changes in VH-IVUS derived plaque components.
Three patients with coronary artery disease underwent baseline and 12 month follow-up angiography and virtual histology-intravascular ultrasound (VH-IVUS) measurements. Coronary artery models were reconstructed from the data and models with and without side-branches were created. Patient-specific Doppler ultrasound (DUS) data were employed as inflow boundary conditions and computational fluid dynamics was used to calculate the WSS in each model. Further, the influence of representative coronary artery flow waveforms upon WSS values was investigated and the concept of treating WSS using relative, rather than actual, values was explored.
Models that included side-branch outflows and subject-specific DUS velocities were considered to be the reference cases. Hemodynamic differences were caused by the exclusion of side-branches and by imposing alternative velocity waveforms. One patient with fewer side-branches and a scaled generic waveform had little deviation from the reference case, while another patient with several side-branches excluded showed much larger departures from the reference situation. Differences between models and the respective reference cases were reduced when data were analyzed using relative, rather than actual, WSS.
When considering individual subjects, large variations in patient-specific flow rates and exclusion of multiple side-branches in computational models can cause significant differences in observed associations between plaque evolution and ranges of computed WSS. These differences may contribute to the large variability typically found among subjects in pooled populations. Relative WSS may be more useful than actual WSS as a correlative variable when there is a large degree of uncertainty in flow rate data.
Wall shear stress (WSS) derived from computational fluid dynamics (CFD) has emerged as a potential tool to predict coronary artery plaque progression in humans. The PREDICTION study demonstrated that segments exposed to low WSS resulted in an increase in plaque burden  and a study by Samady et al. also found an increase in plaque burden in low WSS regions . Further, a decrease in plaque area accompanied with an increase in necrotic core area was observed in regions exposed to high WSS - suggestive of a transformation to a vulnerable plaque phenotype . As well as these prospective studies, there have been several cross-sectional studies that have shown similar results. Early plaques have been found to have a greater degree of necrotic core exposed to low WSS while advanced plaques (plaque burden > 40%) have a higher degree of necrotic core exposed to high WSS . After adjusting for plaque burden Eshtehardi et al. also found necrotic core to be greater in low WSS segments .
In order to predict WSS precisely in human coronary arteries accurate measurements of patient geometry and flow are necessary. Both of these measurements can be made through interventional techniques; patient geometry by a combination of intravascular ultrasound (IVUS) and angiography, and flow through the use of Doppler ultrasound (DUS). It has previously been shown that coronary artery geometry is the most important determinant of WSS patterns ; and although cardiac motion alters coronary artery geometry, this has been found to have little effect on time averaged WSS calculations [6, 7]. However, the presence of side-branches has not been fully addressed. Local WSS patterns have been shown to be altered by the presence of side-branches, and differences of up to 12 Pa have been reported in arteries with multiple side-branches [8, 9]. The impact of these hemodynamic differences on observed relationships of WSS and plaque progression has not been investigated to date. Further, it has been shown that higher inlet flow rates elevate WSS magnitudes . As patients are at rest and typically sedated during catheterization, velocity measurements represent only a snapshot of the patient's coronary flow, raising the question of relying upon this one-time measurement when coronary flow rates can be expected to vary throughout a normal day. Given the large variability than can be seen in an individual patient, these factors may have an important impact on actual WSS values.
This study investigates the sensitivity of observed relationships between computed values of WSS and changes in plaque composition over time by altering two parameters. Firstly, we look at the impact of excluding coronary side-branches on WSS calculations, an assumption made in some studies [1, 3]; and secondly, we investigate the effect of using physiologically representative velocity waveform inlet data rather than patient-specific velocities. Because of expected variability in velocity/flow boundary conditions, we also introduce the concept of relative WSS in examining WSS/plaque changes relationships.
Patient geometry characteristics.
Left main + LAD
No. side branches
Vessel length (mm)
Mean inlet velocity (m/s)
Baseline plaque burden (%)
Follow-up plaque burden (%)
Baseline lumen area (mm2)
Follow-up lumen area (mm2)
Three-dimensional reconstruction of coronary arteries was achieved by combining biplane angiography data with VH-IVUS images. This process has previously been described extensively [11–13]. Briefly, the 3D spatial location of the IVUS catheter prior to pullback is determined via back-projection (IC-PRO, Paeion Inc., Ha'ayin, Israel), and this serves as a backbone for stacking the IVUS images, which are placed perpendicular to the catheter with the distance between consecutive images determined by the pullback speed (0.5 mm/s). Coronary side-branches are then identified in the images and added to the main vessel reconstruction. Initially, we attempted to extract this information from biplane angiography data, but it was often not possible to view a side-branch simultaneously in both angiographic images. In cases where side-branches were visible in both views and 3D reconstruction was possible, the take-off point was frequently deemed to be anatomically incorrect. For these reasons side-branches were modelled as cylindrical extensions normal to the centreline of the main vessel, and the side-branch take-off location and orientation were obtained from the IVUS images.
Given that atherosclerosis is a focal disease we investigated WSS and plaque composition changes in 45 degree sectors of each IVUS image . Firstly, we co-registered baseline VH-IVUS images with follow-up VH-IVUS images. Using fiduciary anatomical landmarks such as side-branches, images were co-registered in the axial direction. Next, images were circumferentially co-registered through normalized cross-correlation by taking into account plaque thickness, plaque composition and perivascular tissue . After circumferential co-registration was complete, plaque composition changes from baseline to follow-up were quantified in 45 degree sectors. VH-IVUS identifies four different components of plaque, namely Fibrotic (FB), Fibro-fatty (FF), Necrotic core (NC) and Dense Calcium (DC).
Baseline time-averaged wall shear stress data were then associated with baseline plaque components in the 45 degree sectors, resulting in a baseline WSS value, baseline plaque component area, follow-up plaque component area, and change in plaque component area for each sector. Mean plaque component area changes were then placed into low (WSS < 1 Pa), intermediate (WSS 1-2.5 Pa) and high (WSS > 2.5 Pa) WSS categories depending on the WSS value (2). Due to the potential variation in patient data we repeated the analysis using relative WSS, which we defined by dividing all WSS values into 20 equally sized bins. These values and their associated plaque component area changes were then classified into low, intermediate and high WSS categories. Mean plaque component area change in each category was then reported.
Statistical analysis was performed using Microsoft Excel. Agreement in WSS between models with and without side-branches or using alternative waveforms was assessed by Bland-Altman analysis. This method is used to assess whether there is a consistent bias in WSS calculation from a given method. Continuous variables are reported as mean ± 95% confidence interval. An unpaired Student's t-test was used to investigate the association of continuous (e.g., change in plaque area) to categorical (e.g., WSS category) variables. Comparisons were made between changes in VH-IVUS defined constituents in WSS categories within the two analysis methods (e.g., side-branches included or excluded) for each patient. All statistical tests were two-tailed, and P < 0.05 was established as the level of statistical significance.
Time averaged WSS values averaged over model surface for each patient.
WSS without branches (Pa)
3.2 ± 1.5
11.7 ± 6.4
8.5 ± 5.3
2.95 ± 2.3
3.1 ± 2.4
0.15 ± 0.75
5.5 ± 5.1
10.7 ± 11.3
5.2 ± 7.5
Changes in plaque components and wall shear stress
Somewhat in contrast to patient 1, patient 2 experienced plaque progression in low WSS sectors and plaque regression in high WSS sectors (Figure 7), regardless of how WSS was computed. No significant differences were observed whether side-branches were excluded or not in the relationship between plaque progression and all WSS categories. These findings agree very well with the Bland-Altman plot of patient 2 which also showed very little difference in WSS between these models. Characteristics of plaque change using relative WSS as a variable were similar to those using WSS, though there were differences in the intermediate WSS category.
All VH-IVUS components saw an increase in area in Patient 3. As with patient 2, there were similar characteristics for plaque changes in patient 3 between WSS models with and without side-branches and to a lesser extent between WSS and relative WSS (Figure 8). However, there were significant differences seen between models with and without side-branches in the low WSS category for necrotic core and dense calcium progression and also in the intermediate WSS category for fibrotic plaque. There were no significant differences in plaque progression association with relative WSS between models with and without side-branches. Several significant differences could be seen between plaque progression prediction using WSS and relative WSS, although this was almost completely confined to models without side-branches.
Large differences were seen in the association of plaque component changes with WSS when the generic waveform (See Figure 1) was used in patient 2 (Figure 9). This waveform results in a greater flow rate which creates higher WSS values, and as a result of this the number of values in the low WSS category is much smaller (n = 16) compared to the patient-specific (n = 49) or the scaled generic (n = 56) waveform. Use of the scaled generic waveform comes closer to replicating that of the patient-specific case, as there were no significant differences seen between both methods. Interestingly, when analyzed using relative WSS, both waveforms predict similar progression as compared to the patient-specific case (Figure 7). There were, however, several significant differences in plaque progression descriptions when analyzed using WSS versus relative WSS, particularly in the intermediate WSS category as was the case in the patient-specific waveform.
Many factors affect plaque localization and evolution, and mechanical factors are among these. There is considerable interest in employing CFD to compute WSS in coronary arteries and to relate this variable to changes in plaque composition over time, hopefully as one element of a clinical milieu that can be used to improve patient care. Understanding both the potential and limitations of WSS as a predictive variable is essential for addressing clinical needs. This study examined intra-patient variability in hemodynamics by assessing the focal association of WSS and atherosclerotic plaque changes over a twelve month period in three patients with coronary artery disease. Concentrating on flow-related boundary conditions for CFD, we explored the effects of altered patient coronary artery flow waveforms (inflow) and the presence or absence of side-branches (outflow). Given the potential differences these modelling strategies have on WSS values, we also investigated plaque changes using a relative WSS method. Since data acquisition typically occurs at one time point but there are temporal variations in flow rates and geometry over the course of a heartbeat and, more so, over longer time periods, caution may need to be exercised when relating plaque progression to actual, as opposed to relative, WSS values.
When comparing models with and without side-branches, greater differences in WSS values were seen in patients who had a greater number of side-branches. This is due to flow leaving the main vessel, resulting in a reduced distal flow rate and hence reduced WSS. Without these side-branches WSS remains artificially high in distal segments. For example, in patient 3 the presence of 5 side-branches resulted in a mean WSS difference of 5.2 Pa (Figure 3). Though patient 1 had only one more side-branch than patient 2, their respective differences in WSS for models with and without side-branches were 0.16 and 8.5 Pa. This can be attributed to the presence of the LCX branch in patient 1 which, due to its large size, has a large outflow. Similarly to the exclusion of side-branches, the use of a non patient-specific waveform with a high flow rate as the inflow condition resulted in a large mean difference (6.6 Pa). This difference was effectively eliminated when the same waveform was scaled to the same mean velocity as the patient-specific case (Figure 4), suggesting that the specific waveform shape is less important than the imposed mean flow. Qualitatively, the spatial patterns of relative WSS were similar regardless of the presence or absence of side branches or of the level of flow rate used as an input (Figure 5).
We then explored whether these hemodynamic differences would result in differences in interpreting plaque changes as related to WSS in individual patients. No significant differences were seen between inclusion/exclusion of side branches in patient 2 which had the least number of side branches and a shorter LAD length. On the other hand, more noticeable differences could be seen in patient 3 in low and intermediate WSS categories. This is due to the artificially higher distal flow rate caused by neglecting outflows in the five side-branches, thus increasing computed WSS and resulting in a greater number of sectors being categorized as high. For patient 1, despite the large hemodynamic differences seen, only fibrotic plaque change in the intermediate WSS category was statistically different between the side-branch/no side-branch cases. Large differences were seen in several of the components in the low WSS category, though as some of these (FF, NC, DC) contained only 1 sector they could not be tested for significance. Notably, neglecting branch outflow distorted the WSS field to the point that low WSS regions were essentially eliminated in patient 1.
The use of a generic inlet velocity waveform in patient 2 resulted in significant changes in plaque progression association with WSS compared to the use of the patient-specific waveform. This patient had an overall trend of increasing plaque and VH-IVUS component area in low WSS categories, changing to regression in high WSS sectors (Figure 7). While this trend was still evident in the scaled generic waveform model, the same could not be said for the basic generic waveform case. Scaling the generic waveform to match that of the mean velocity of the patient-specific waveform was seen to eliminate hemodynamic differences and as a result maintained the reference case relationships between changes in plaque component area and WSS.
Because of the sensitivity of computed values of WSS to inflow and outflow conditions, we also compared plaque changes in each scenario using a relative WSS metric under the hypothesis that areas of high WSS will remain relatively high and areas of low WSS will remain relatively low (this is more likely the case where the inlet waveform is altered as opposed to the side-branch exclusion cases, as there may be local changes in high and low WSS at these branch locations). When relative WSS is the hemodynamic variable being considered for plaque progression prediction, the generic waveform did not result in any significant differences compared to the patient-specific case (Figures 7 and 9). This is in contrast to using actual WSS where there were differences between the generic and patient-specific cases in observed plaque progression/WSS relationships. Using relative WSS also suppressed differences in observed plaque progression relationships when the models excluded side-branches (Figures 6, 7, 8, right panels). These findings suggest that relative WSS could be used as a surrogate for actual WSS in predicting plaque progression when patient-specific data are unavailable or unreliable.
A limitation of this study relates to how the outflow boundary conditions are treated. Ideally, patient pressure and flow rate should be measured in each branch, but this is not realistic in routine clinical practice. Previously, other investigators have specified outlet boundary conditions based on Murray's Law or empirical data from a study by Doriot et al. [15, 16]. They reported that average WSS was higher by 8% when using Murray's Law, though the size of low WSS regions varied by up to 68% in bifurcation regions. Accurate knowledge of branch diameters is required in order to use these power laws, and these measurements were not available due to the limited penetration depth of IVUS. Variation of the flow rate under pressure free boundary conditions has previously been found to have no significant changes on predictive values of WSS . Though we have not performed a sensitivity analysis of the boundary conditions in this study, we do not expect the uncertainty in outflow boundary conditions to have an important effect given that the WSS is ultimately parsed into low, intermediate and high categories.
This study investigated the LAD coronary artery of patients with non-obstructive coronary artery disease and assessed effects on computed WSS arising from intra-patient variability in flow boundary conditions that may arise due to experimental measurement errors, flow changes to be expected under various physiological conditions, and assumptions made relating to the presence of side-branch outflow. The reference case for each individual employed patient-specific DUS velocity measurements from which LAD flow waveforms were computed, and the CFD model geometry included side-branches. Other cases examined were models without side-branches, and two alternative inlet flow waveforms.
Neglecting branch outflows can result in large differences in computed WSS (patients 1 and 3; Figures 2 and 3). Differences are exacerbated if the arterial segment being studied has multiple branches, as can be expected in longer segments.
The relationships between plaque component changes and computed WSS within a given patient can be significantly different depending upon the flow conditions. Specifically, neglecting side branches can lead to significantly different predictions of plaque progression than the reference case (patients 1 and 3; Figures 6 and 8, left panels) in individual cases.
High flow rates can also affect observed plaque progression relationships by comparison with the patient-specific flow rate case (patient 2; Figure 7, upper left panel, and Figure 9, upper left panel). This implies that use of a "snapshot" flow measurement at a baseline study may not fully capture this dynamic relationship experienced over periods of days or weeks.
These studies should be considered exploratory, since we investigated a small number of patients, and more extensive work is needed in order to develop definitive, optimum approaches that can give greater confidence in WSS/plaque prediction methods for individuals. Nonetheless, the present work serves to illustrate important issues that should be considered and demonstrates approaches that offer improvements in interpretation of data within the clinical setting.
All human data used in this study were obtained from consenting patients who presented to the Cardiac Catheterization Laboratory at Emory University Hospital and was approved by the hospital's Ethical Research Committee.
We acknowledge funding from the American Heart Association Greater Southeast Affiliate - Postdoctoral Fellowships (DSM), the Georgia Research Alliance (DPG, LHT), Toshiba America Medical Systems (HS, DPG), Volcano Corp. (HS) and Forest Pharmaceuticals (HS).
Publication of this article was paid with funding from the Georgia Research Alliance Eminent Scholar Program
This article has been published as part of BioMedical Engineering OnLine Volume 14 Supplement 1, 2015: Cardiovascular Disease and Vulnerable Plaque Biomechanics. The full contents of the supplement are available online at http://www.biomedical-engineering-online.com/supplements/14/S1
- Stone PH, Saito S, Takahashi S, Makita Y, Nakamura S, Kawasaki T, Takahashi A, Katsuki T, Nakamura S, Namiki A, Hirohata A, Matsumara T, Yamazaki S, Yokoi H, Tanaka S, Otsuji S, Yoshimachi F, Honye J, Harwood D, Reitman M, Coskun AU, Papafaklis MI, Feldman CL, PREDICTION Investigators: Prediction of progression of coronary artery disease and clinical outcomes using vascular profiling of endothelial shear stress and arterial plaque characteristics: the PREDICTION study. Circ 2012, 126: 172–181. 10.1161/CIRCULATIONAHA.112.096438View ArticleGoogle Scholar
- Samady H, Eshtehardi P, McDaniel MC, Suo J, Dhawan SS, Maynard C, Timmins LH, Quyyumi AA, Giddens DP: Coronary artery wall shear stress is associated with progression and transformation of atherosclerotic plaque and arterial remodelling in patients with coronary artery disease. Circ 2011, 124: 770–788.View ArticleGoogle Scholar
- Wentzel JJ, Schuurbiers JCH, Lopez N, Gijsen FJH, van der Giessen AG, Groen HC, Dijkstra J, Garcia-Garcia HM, Serruys P: In vivo assessment of the relationship between shear stress and necrotic core in early and advanced coronary artery disease. EuroIntervention 2013, 9: 989–995.Google Scholar
- Eshtehardi P, McDaniel MC, Suo J, Dhawan SS, Timmins LH, Binongo JNG, Golub LJ, Corban MT, Finn AV, Oshinski JN, Quyyumi AA, Giddens DP, Samady H: Association of coronary wall shear stress with atherosclerotic plaque burden, composition, and distribution in patients with coronary artery disease. J Am Heart Assoc 2012, 1: e002634.View ArticleGoogle Scholar
- Myers JG, Moore JA, Johnston KW, Ethier CR: Factors influencing blood flow patterns in the human right coronary artery. Ann Biomed Eng 29: 109–120.Google Scholar
- Torii R, Heegan J, Wood NB, Dowsey AW, Hughes AD, Yang GZ, Firmin DN, Thom SA, Xu XY: The effect of dynamic vessel motion on haemodynamic parameters in the right coronary artery: a combined MR and CFD study. Brit J Radiol 2009, 82: S24-S32. 10.1259/bjr/62450556View ArticleGoogle Scholar
- Zeng D, Ding Z, Friedman MH, Ethier CR: Effect of cardiac motion on right coronary artery hemodynamics. Ann Biomed Eng 2004, 31: 420–429.View ArticleGoogle Scholar
- Gijsen FJH, Wentzel JJ, Thury A, Lamers B, Schuurbiers JCH, Serruys PW: A new imaging technique to study 3-D plaque and shear stress distribution in human coronary artery bifurcations in vivo. J Biomech 2007, 40: 2349–2357. 10.1016/j.jbiomech.2006.12.007View ArticleGoogle Scholar
- Wellnhofer E, Osman J, Kertzcher U, Affeld K, Fleck E, Goubergrits L: Flow simulation studies in coronary arteries - Impact of side-branches. Atherosclerosis 2010, 213: 475–481. 10.1016/j.atherosclerosis.2010.09.007View ArticleGoogle Scholar
- Liu B, Tang D: Computer simulations of atherosclerotic plaque growth in coronary arteries. Mol Cell Biomech 2010, 7: 193–202.MathSciNetGoogle Scholar
- Slager CJ, Wentzel JJ, Schuurbiers JCH, Oomen JA, Kloet J, Krams R, von Birgelen C, van der Giessen WJ, Serruys PW, de Feyter PJ: True 3-Dimensional reconstruction of coronary arteries in patients by fusion of angiography and IVUS (ANGUS) and its quantitative validation. Circulation 2000, 102: 511–516. 10.1161/01.CIR.102.5.511View ArticleGoogle Scholar
- Prause GPM, DeJong SC, McKay CR, Sonka M: Towards a geometrically correct 3-D reconstruction of tortuous coronary arteries based on biplane angiography and intravascular ultrasound. Int J Card Imag 1997, (13):451–462.Google Scholar
- Timmins LH, Molony DS, Eshtehardi P, McDaniel MC, Oshinski JN, Samady H, Giddens DP: Examination of the circumferential versus focal association between wall shear stress and clinical coronary artery disease progression. Ann Biomed Eng, in press.Google Scholar
- Timmins LH, Suever JS, Eshtehardi P, McDaniel MC, Oshinski JN, Samady H, Giddens DP: Framework to co-register longitudinal virtual histology-intravascular ultrasound in the circumferential direction. IEEE Trans Med Imaging 2013, (11):1989–1996.Google Scholar
- Van Der Giessen AG, Groen HC, Doriot P, de Feyter P, van der Steen A, van de Vosse FN, Wentzel JJ, Gijsen FJH: The influence of boundary conditions on wall shear stress distribution in patients specific coronary trees. J Biomech 2011, 44: 1089–1095. 10.1016/j.jbiomech.2011.01.036View ArticleGoogle Scholar
- Doriot PA, Dorsaz PA, Dorsaz L, De Benedetti E, Chatelain P, Delafontaine P: In-vivo measurements of wall shear stress in human coronary arteries. Coron Artery Dis 2000, (11):495–502.Google Scholar
- Rikhtegar F, Knight JA, Olgac U, Saur SC, Poulikakos D, Marshall W Jr, Cattin PC, Alkadhi H, Kurtcuoglu V: Choosing the optimal wall shear parameter for the prediction of plaque location - A patient specific computational study in human coronary arteries. Atherosclerosis 2012, 221: 432–437. 10.1016/j.atherosclerosis.2012.01.018View ArticleGoogle Scholar
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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.