- Open Access
Influence of vascular geometry on local hemodynamic parameters: phantom and small rodent study
BioMedical Engineering OnLine volume 17, Article number: 30 (2018)
Many studies have demonstrated that the geometry of the carotid bifurcation enables prediction of blood flow variation associated with atherosclerotic plaque formation. The phase angle between the arterial wall circumferential strain and its instantaneous wall shear stress is known as stress phase angle (SPA). This parameter is used to evaluate hemodynamic factors of atherogenesis. Note that SPA can be numerically computed for the purpose of locating atherosclerosis in different artery geometries. However, there is no experimental data to verify its role in the location of atherosclerosis in different artery geometries. In this study, we use an ultrasonic biomechanical method to experimentally evaluate the role of SPA for locating atherosclerosis in carotid bifurcation.
For carotid anthropomorphic vascular phantom experiments, the SPAs of common carotid arteries (CCAs), external carotid arteries (ECAs) and internal carotid arteries (ICAs) are − 148.53 ± 6.92°, − 153.95 ± 5.11°, and − 238.69 ± 1.72°, respectively. The corresponding SPAs are − 173.47 ± 0.065°, − 115.57 ± 4.83° and − 233.9 ± 8.12° for the polyvinyl alcohol (PVA-c) phantoms. In vivo mouse experiments indicated that the wall shear stress and circumferential strain were out of phase in the ICAs (− 280.08 ± 13.12°) to a greater extent as compared to CCAs (− 141.97 ± 8.03°) and ECAs (− 170.07 ± 9.24°).
The results suggested that SPA may be a useful indicator to locate the atherosclerosis position in carotid bifurcation.
Atherosclerotic plaques are localized and they are most likely to occur at arterial bifurcations and bends [1, 2]. Many studies have demonstrated that the mechanical forces exposed to arterial endothelial cells play an important role in the development and progression of atherosclerosis [3,4,5,6]. These forces mainly include wall shear stress (WSS) as a result of blood flow and the circumferential strain (CS) caused by periodical wall motion and pulsatile pressure. Since these forces are primarily determined by the arterial geometry, some regions of the artery may have an increased risk of developing atherosclerosis , particularly for arterial bifurcations. Recent research has shown that the plaque formation is most obvious at the outer wall of the carotid sinus where the flow separates [8,9,10]. Furthermore, arterial remodeling and the development of atherosclerotic disease need to be evaluated by the circumferential wall stress, an essential parameter for determining atherogenesis [7, 11,12,13]. Therefore, the CS and WSS can be used as indicators for carotid bifurcation atherosclerosis. The temporal phase angle between CS and WSS, which we have termed as the stress phase angle (SPA), determines the influencing effect of arterial wall shear stress and blood shear stress [13, 14]. Note that SPA has been proven as a potential important parameter due to its role in the pathophysiology of vascular disorders . Furthermore, several numerical investigations on SPA have been carried out to locate the positions of atherosclerosis in different geometric arteries. Carotid bifurcation modelling using fluid–structure interaction shows that large negative SPA values are located at the atherosclerotic plaques attached to the outer wall of the carotid sinus . In addition, some studies also demonstrated that SPA had significant value for predicting stenosis severity of artery [17, 18]. However, the role of SPA for locating atherosclerosis in different geometric arteries has not been verified experimentally yet.
Recently, an ultrasonic biomechanical (UBM) method was established to characterize the displacement and flow pattern of arterial wall simultaneously [19,20,21,22]. This method is employed to quantitatively determine and evaluate the function of SPA in predicting early atherosclerosis based on an experimental platform. Results have shown that highly negative SPA correlates to arterial wall stiffening . The UBM method has been utilized in this experimental study to analyze the role of SPA for location of atherosclerosis in different geometric arteries. Firstly, the changes of SPA in different locations of the carotid artery are studied by UBM. High levels of negative SPA could be found in the internal carotid artery (ICA) as compared to the common carotid artery (CCA) and external carotid artery (ECA). Secondly, the results are verified by in vivo mouse experiments. Experimental results demonstrated that a greater negative SPA value was indicative of arterial disease risk.
The UBM method is capable of measuring the displacement of vessel wall and blood flow simultaneously, which has been demonstrated in . Further detail on this method for calculating the SPA is given in , and this method has been used to analyze SPA temporal distribution with arterial stiffening in an experimental setup. The data file acquired contains frames over several cardiac cycles. Contrast-enhanced images are computed between consecutive frames using the UBM method. Firstly, ROIs are selected including the vessel wall and lumen. Each ROI is divided into a grid of small sections known as interrogation windows. Then, the UBM method is used to calculate the displacement of each interrogation window. This process is repeated at each interrogation window within the image, resulting in a map of displacement vectors to describe the arterial wall movement and flow. The CS can be obtained from the displacement of vessel wall, and WSS can be determined based on gradient of the flow velocities. Therefore, the SPA can be acquired by the temporal phase angle between CS and WSS.
A closed-loop compression system can be used to pressurize an arterial phantom (Fig. 1) and an ultrasound imaging system to acquire the image sequence. Pulsatile flow, mimicking the ventricular action of the heart, was generated by a Harvard Apparatus pulsatile pump (Model 55-3305; Harvard Apparatus, Holliston, MA, USA). Ultrasound contrast microbubbles containing cores of octafluoropropane (C3F8) gas surrounded by a phospholipid shell are developed. A 0.025 ml volume of microbubbles was injected into the upstream chamber and mixed to ensure uniform bubble concentration during each measurement. Then, we performed microbubble seeding as described in the previous section . The transducer was placed perpendicular to the flow direction to obtain the ultrasound contrast images.
A carotid anthropomorphic vascular phantom manufactured by Shelley Medical Imaging Technologies (London, ON, Canada) for mimicking a normal carotid bifurcation was used to study the change of SPA at different locations of the carotid bifurcation (Fig. 2). An open-architecture ultrasound system (Sonix RP, Vancouver, Canada) with a 10 MHz linear array transducer was used to image the vascular phantom in the longitudinal section. The frame rate was 125 Hz and the field of view was 40 mm (depth) × 20.8 mm (width).
In addition, the arterial phantom was made of PVA cryogel (PVA-c). The PVA aqueous solution was filled into a steel mold, and then frozen and subsequently thawed to form a cryogel with rubber-like properties . The composition (by weight) of the solution was 10% PVA powder, 87% deionized water, and 3% scattering particles composed of sigmacell cellulose with diameter of 20 μm (Sigma-Aldrich, USA). The mold consisted of two polished stainless-steel shells, a bifurcate inner rod (6 mm CCA diameter, 4.5 mm ICA diameter, 4 mm ECA diameter), and two hexagonal screws, as shown in Fig. 3a, b. The solution was injected into the gap between the rod and the shells (Fig. 3c). The whole mold underwent a number of freeze–thaw cycles to generate higher stiffness. Each freeze–thaw cycle was composed of 12-h freezing at − 20 °C and 12-h thawing at 20 °C. The optical photo and the ultrasound image of the carotid phantom are shown in Fig. 3e, f. The stiffness of the PVA phantom was also tested on an electronic universal material testing machine (CMT6104; New Sans Machinery Co., Ltd., Shenzhen, Guangdong, China). The phantom with three freeze–thaw cycles (Young’s modulus, 162.45 kPa) was used in flow experiments. The carotid phantom was imaged using a high frequency small-animal imaging system with a 40-MHz linear transducer (Visualsonics). The image acquisition frame rate was 125 Hz and the number of the image sequence was 800.
In vivo mouse carotid artery experiments
The change of SPA in carotid bifurcation was further investigated in vivo. Our animal experiments were conducted in strict accordance to protocols that are approved by the Institutional Animal Care and Use Committee of Shenzhen Institutes of Advanced Technology, Chinese Academy of Sciences. During the injection and image acquisition process, mice were anesthetized with 1.0–2.0% isoflurane in oxygen delivered at a flow rate of 1.0 min/l and monitored with electrocardiogram. Using a rectal temperature probe, body temperature was carefully maintained between 36.7 and 37.3 °C throughout the study. Ultrasound contrast agents composed of octafluoropropane (C3F8) gas encapsulated by a phospholipid shell fabricated in house were injected through the tail vein into each mouse with a 2.0 × 105 microbubbles/ml concentration. Ultrasound biomicroscopy (Vevo 2100, Visualsonics, Toronto, Canada) with a transducer frequency of 40 MHz was used for vascular imaging in anesthetized mice. The frame rate was 200 Hz.
The analysis of covariance (ANOVA) was used to examine the role of SPA for locating atherosclerosis in carotid bifurcation. A p value less than 0.05 was accepted as indicating statistical significance. All statistical analyses were performed using the Statistical Package for Social Sciences statistical software package, version 17.0 (SPSS Inc.).
Carotid anthropomorphic vascular phantom experiments
The UBM method was used to measure the measure hemodynamic parameters in the carotid anthropomorphic vascular phantom experiments. The distribution of arterial diameter and wall shear rate for vascular phantom in the CCA (a), ECA (b) and ICA (c) was shown in Fig. 4. It is worthwhile noting that the diameter and WSR curves exhibit a periodical variation and phase shift. The SPA values located at the CCA, ECA and ICA are − 148.53 ± 6.92°, − 153.95 ± 5.11°, and − 238.69 ± 1.72°, respectively. It can be observed that SPA is more negative for the ICA where the atherosclerotic plaque develops.
PVA-c phantom experiments
The PVA phantoms are performed to verify the accuracy of the results in carotid anthropomorphic vascular phantom. The temporal diameter and WSR curves of the CCA, ECA and ICA in the PVA phantom are shown in Fig. 5a–c, respectively. The diameter and WSR curves show a periodical variation and phase shift. Note that the SPA, WSR and the strain at different locations of the PVA phantoms are shown in Table 1. There is a 35% increment in negative SPA, a 27% reduction in WSR, and a 26% reduction in strain pertaining to the ICA in comparison to that of the CCA.
In vivo mouse carotid artery experiments
The UBM method is utilized in determining the flow velocity distributions and the arterial diameter changes of the CCAs, ECAs and ICAs present in mice. The change of diameter and WSR profiles with time is depicted in Fig. 6. It is notable that there are periodic patterns between the diameter and WSR and the phase shift at various times. The SPA, WSR and the strain at various locations of mice are shown in Table 2. Here, more negative SPA is achieved in the ICAs (− 280.08 ± 13.12°) compared with CCAs (− 141.97 ± 8.03°) and ECAs (− 170.07 ± 9.24°), while correspondingly, WSR and strain are lower in the ICAs. A comparison of different locations was made by analysis of variance. Significant differences were found among locations for SPA (p = 0.0001), WSR (p = 0.045) and strain (p = 0.009).
This study showed that SPA was highly dependent on geometries and it might be an essential predictor in early atherosclerosis. Note that CS and WSS were not synchronous in the ICA as verified by in vitro carotid anthropomorphic vascular phantoms. Moreover, the SPA is large and negative for the ICA where atherosclerotic plaques located. The results were further demonstrated by PVA-c phantom experiments and in vivo mouse carotid artery experiments. Previous studies have indicated that SPA appears to be negative at higher degrees in positions that depict flow separation having low shear stress in comparison to the non-separated positions that have high shear stress. Lee and Tarbell performed experimental studies in a compliant aortic bifurcation model and found that the SPA was − 40° on the flow divider (high shear) and − 100° on the outer wall (low shear) . Tada and Tarbell showed that the SPA became highly negative (approaching − 180°) along the entire length of the carotid sinus outer wall . The results are consistent with our measurements and computations of the SPA in the carotid bifurcation and large negative SPA in precisely those regions of arteries where atherosclerotic disease usually develops.
Many studies have shown that the development of atherosclerosis in the naturally bulbic ICA may be due to the local hemodynamic conditions, such as reduced and oscillating WSS [7, 28, 29]. Although cardiovascular risk factors typically cause the thickening and stiffening of the CCA wall, the development of atherosclerosis at the naturally dilated ICA bulb is mainly related to the bifurcation geometry . Our findings indicate that the WSR is found to be significantly reduced in the ICA, as shown in Tables 1, 2. Our results offer further evidence that the area of atherosclerotic plaque development has a higher tendency for low WSS values and indicates that the geometry of the artery may contribute to subclinical atherosclerosis. The regions at relatively low WSS are also characterized by more negative values of SPA in the phantom and mice experiments as described in this study. Therefore, the geometry and unique hemodynamic characteristics (SPA = − 238.69 ± 1.72° for vascular phantom; SPA = − 233.9 ± 8.12° for PVA-c phantom; SPA = − 280.08 ± 13.12° for mice) of ICAs may contribute to predicting high-risk atherosclerosis regions in these vessels.
The SPA combines the information of both WSS, which was resulted from blood flow, and the CS caused by periodical wall motion and pulsatile pressure. The results suggest that flow separation and low WSS areas comprise the most negative SPA (Figs. 4, 5, 6), which is consistent with previous studies [2, 29, 30]. Biological studies have demonstrated that negative SPA can inhibit anti-atherogenic gene expression and release, but at the same time, they also increase pro-atherogenic gene expression and metabolite release [31, 32]. These results proved that ECs driven by hemodynamics and wall mechanics can potentially cause a pro-atherogenic effect on regions of the local circulation.
There were several limitations of this study: (1) The ultrasound data for CCA, ICA and ECA were collected independently at different initial phases. The SPA is independent of the initial phase, thus the results are unaffected. (2) Previous studies addressed that the acoustic properties of the experimental phantoms (artificial tissue, vessels, and blood) should match to those of actual human tissue . Future studies will take into account these effects. (3) Given the small number of animals, our findings must be regarded as preliminary. Future research with larger sample size will be required to further verify these findings .
The determination and analysis of SPA is proposed for locating the positions of atherosclerosis in different geometric arteries based on the UBM method. The SPA of CCA, ECA and ICA is − 148.53 ± 6.92°, − 153.95 ± 5.11°, and − 238.69 ± 1.72° for carotid anthropomorphic vascular phantom experiments, respectively. The corresponding SPA is − 173.47 ± 0.065°, − 115.57 ± 4.83° and − 233.9 ± 8.12° for the PVA phantoms, respectively. The in vitro experimental results indicated that in the ICAs, WSS and CS were of different phases. The SPA was more negative in the ICAs as verified and indicated by our in vivo mouse experiments (− 280.08 ± 13.12°) compared with CCAs (− 141.97 ± 8.03°) and ECAs (− 170.07 ± 9.24°). The results indicate that SPA may have a key role in prediction of the atherosclerosis-prone regions, which now motivates future study on larger sample sizes for further verification.
Malek AM, Alper SL, Izumo S. Hemodynamic shear stress and its role in atherosclerosis. J Am Med Assoc. 1999;282(21):2035–42.
Lee S-W, Antiga L, Spence JD, Steinman DA. Geometry of the carotid bifurcation predicts its exposure to disturbed flow. Stroke. 2008;39(8):2341–7.
Jiang Y, Kohara K, Hiwada K. Association between risk factors for atherosclerosis and mechanical forces in carotid artery. Stroke. 2000;31(10):2319–24.
Papafaklis MI, Takahashi S, Antoniadis AP, Coskun AU, Tsuda M, Mizuno S, Andreou I, Nakamura S, Makita Y, Hirohata A. Effect of the local hemodynamic environment on the de novo development and progression of eccentric coronary atherosclerosis in humans: insights from prediction. Atherosclerosis. 2015;240(1):205–11.
Dumont DM, Doherty JR, Trahey GE. Noninvasive assessment of wall-shear rate and vascular elasticity using combined ARFI/SWEI/spectral Doppler imaging system. Ultrason Imaging. 2011;33(3):165–88.
Wong KK, Kelso RM, Worthley S, Sanders P, Mazumdar J, Abbott D. Medical imaging and processing methods for cardiac flow reconstruction. J Mech Med Biol. 2009;9(01):1–20.
Liu G, Wu J, Huang W, Wu W, Zhang H, Wong KK, Ghista DN. Numerical simulation of flow in curved coronary arteries with progressive amounts of stenosis using fluid-structure interaction modelling. J Med Imaging Health Inform. 2014;4(4):605–11.
Zarins CK, Giddens DP, Bharadvaj B, Sottiurai VS, Mabon RF, Glagov S. Carotid bifurcation atherosclerosis. Quantitative correlation of plaque localization with flow velocity profiles and wall shear stress. Circ Res. 1983;53(4):502–14.
Fisher M. Stroke: geometry is destiny for carotid atherosclerotic plaques. Nat Rev Neurosci. 2012;8(3):127–9.
Wong KK, Thavornpattanapong P, Cheung SC, Tu J. Biomechanical investigation of pulsatile flow in a three-dimensional atherosclerotic carotid bifurcation model. J Mech Med Biol. 2013;13(01):1350001.
Pries AR, Reglin B, Secomb TW. Remodeling of blood vessels: responses of diameter and wall thickness to hemodynamic and metabolic stimuli. Hypertension. 2005;46(4):725.
Laurent S, Boutouyrie P. Dose-dependent arterial destiffening and inward remodeling after olmesartan in hypertensives with metabolic syndrome. Hypertension. 2014;64(4):709–16.
Wong KKL, Kelso RM, Worthley SG, Sanders P, Mazumdar J, Abbott D. Cardiac flow analysis applied to phase contrast magnetic resonance imaging of the heart. Ann Biomed Eng. 2009;37(8):1495–515.
Qiu Y, Tarbell JM. Interaction between wall shear stress and circumferential strain affects endothelial cell biochemical production. J Vasc Res. 2000;37(3):147–57.
Förstermann U, Sessa WC. Nitric oxide synthases: regulation and function. Eur Heart J. 2012;33(7):829–37.
Tada S, Tarbell J. A computational study of flow in a compliant carotid bifurcation–stress phase angle correlation with shear stress. Ann Biomed Eng. 2005;33(9):1202–12.
Torii R, Wood NB, Hadjiloizou N, Dowsey AW, Wright AR, Hughes AD, Davies J, Francis DP, Mayet J, Yang G-Z. Stress phase angle depicts differences in coronary artery hemodynamics due to changes in flow and geometry after percutaneous coronary intervention. Am J Physiol Heart Circ Physiol. 2009;296(3):H765.
Sadeghi M, Shirani E, Tafazzoli-Shadpour M, Samaee M. The effects of stenosis severity on the hemodynamic parameters-assessment of the correlation between stress phase angle and wall shear stress. J Biomech. 2011;44(15):2614–26.
Niu L, Qian M, Song R, Meng L, Liu X, Zheng H. A 2D non-invasive ultrasonic method for simultaneous measurement of arterial strain and flow pattern. Clin Physiol Funct Imaging. 2012;32(4):323–9.
Tortoli P, Morganti T, Bambi G, Palombo C, Ramnarine KV. Noninvasive simultaneous assessment of wall shear rate and wall distension in carotid arteries. Ultrasound Med Biol. 2006;32:1661–70.
Hasegawa H, Kanai H. Simultaneous imaging of artery-wall strain and blood flow by high frame rate acquisition of RF signals. IEEE Trans Ultrason Ferroelectr Freq Control. 2008;55:2626–39.
Luo J, Konofagou EE. Imaging of wall motion coupled with blood flow velocity in the HEART and vessels in vivo: a feasibility study. Ultrasound Med Biol. 2011;37:980–95.
Niu L, Meng L, Xu L, Liu J, Wang Q, Xiao Y, Qian M, Zheng H. Stress phase angle depicts differences in arterial stiffness: phantom and in vivo study. Phys Med Biol. 2015;60(11):4281.
Niu L, Qian M, Song R, Meng L, Liu X, Zheng H. A texture matching method considering geometric transformations in noninvasive ultrasonic measurement of arterial elasticity. Ultrasound Med Biol. 2012;38(3):524–33.
Niu L, Qian M, Yan L, Yu W, Jiang B, Jin Q, Wang Y, Shandas R, Liu X, Zheng H. Real-time texture analysis for identifying optimum microbubble concentration in 2-D ultrasonic particle image velocimetry. Ultrasound Med Biol. 2011;37(8):1280–91.
Dineley J, Meagher S, Poepping T, McDicken W, Hoskins P. Design and characterisation of a wall motion phantom. Ultrasound Med Biol. 2006;32(9):1349–57.
Lee C, Tarbell JM. Influence of vasoactive drugs on wall shear stress distribution in the abdominal aortic bifurcation: an in vitro study. Ann Biomed Eng. 1998;26(2):200–12.
Sitzer M, Puac D, Buehler A, Steckel DA, von Kegler S, Markus HS, Steinmetz H. Internal carotid artery angle of origin a novel risk factor for early carotid atherosclerosis. Stroke. 2003;34(4):950–5.
Markl M, Wegent F, Zech T, Bauer S, Strecker C, Schumacher M, Weiller C, Hennig J, Harloff A. In vivo wall shear stress distribution in the carotid artery effect of bifurcation geometry, internal carotid artery stenosis, and recanalization therapy. Circ Cardiovasc Imaging. 2010;3(6):647–55.
Wong KK, Tu J, Kelso RM, Worthley SG, Sanders P, Mazumdar J, Abbott D. Cardiac flow component analysis. Med Eng Phys. 2010;32(2):174–88.
Dancu MB, Berardi DE, Vanden Heuvel JP, Tarbell JM. Asynchronous shear stress and circumferential strain reduces endothelial NO synthase and cyclooxygenase-2 but induces endothelin-1 gene expression in endothelial cells. Arterioscler Thromb Vasc Biol. 2004;24(11):2088.
Wong KK, Tu J, Sun Z. Methods in research and development of biomedical devices. Singapore: World Scientific Publishing Co.; 2013.
Hoskins PR. Simulation and validation of arterial ultrasound imaging and blood flow. Ultrasound Med Biol. 2008;34(5):693–717.
Wong KK, Hui SC. Ethical principles and standards for the conduct of biomedical research and publication. Australas Phys Eng Sci Med. 2015;38(3):377–80.
LN, LM and XZ designed the project. LN, LX and LM performed the experiments. MP, DA, HZ and XZ discussed the results. LN and LM wrote the manuscript. All authors reviewed the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article.
Ethics approval and consent to participate
Our animal experiments were conducted in strict accordance to protocols that are approved by the Institutional Animal Care and Use Committee of Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences.
The work was supported by National Key R&D Program of China (2016YFC0104700), National Science Foundation Grants (NSFC Grant Nos. 11574341, 11774371, 11674347, 81527901, 61773110, 61374015), Natural Science Foundation of Guangdong Province (2016TQ03X716, 2016TQ03R472) and Shenzhen Basic Science Research (JCYJ20160429190550139, JCYJ20160429184552717).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Niu, L., Zhu, X., Pan, M. et al. Influence of vascular geometry on local hemodynamic parameters: phantom and small rodent study. BioMed Eng OnLine 17, 30 (2018) doi:10.1186/s12938-018-0458-8
- Ultrasound imaging
- Wall shear stress
- Carotid bifurcation