- Research article
- Open Access
Shear rate specific blood viscosity and shear stress of carotid artery duplex ultrasonography in patients with lacunar infarction
BMC Neurology volume 13, Article number: 36 (2013)
This study describes a new method for determining site-specific vascular shear stress using dynamic measures of shear rate and blood viscosity (BV) in the carotid arteries, and examines characteristics of carotid arterial shear stress among patients with lacunar infarction.
Vascular shear stress measurements were conducted in 37 patients (17 lacunar infarction patients and 20 control subjects) using duplex ultrasonography. Vessel wall diameters and velocities were measured in each arterial segment at peak-systolic (PS) and end-diastolic (ED) phases, for calculation of PS/ED shear rates. PS/ED shear stresses [dyne/cm2] were determined with PS/ED shear rates and shear-rate dependent BV values. For comparison, both values of hematocrit-derived BV and BV measurements at 300 s-1 were used for calculation of shear stress.
All cardiovascular disease (CVD) risk factors including BV values were similar between the two groups. In both common carotid arteries, PS and ED shear stresses were significantly lower in the patients with lacunar infarction than in controls in multivariate models that included age, sex, and other major CVD risk factors. PS and ED shear stresses using the shear rate specific BV were 4.5% lower and 7.3% higher than those using the two other BVs, respectively.
Lacunar infarction was associated with reduced carotid arterial shear stress. The use of estimated BV for calculating carotid arterial shear stress provides more accurate assessment of the hemodynamic contribution of shear stress than previous models that have arbitrarily assigned a constant value to this dynamic flow property.
The vascular system can generally be viewed as a closed system, constraining a series of biochemical processes and mechanical stresses that are regulated by various homeostatic mechanisms . From this theoretical perspective, if mechanical stresses such as wall shear stresses and tensile stresses are not properly dispersed, then an artery may sustain injury. This mechanical injury ultimately results in atherothromboembolic disease and its clinical manifestations of stroke, myocardial infarction and claudication .
Blood is a fluid suspension of plasma and cells such as erythrocytes, leukocytes, and platelets, and it demonstrates non-Newtonian fluid mechanics. A non-Newtonian fluid has a variable, non-linear relation between blood viscosity and blood flow . Specifically, blood viscosity is higher at low shear rates () and is reduced as shear rate increases . This shear rate dependent aspect of blood viscosity (BV, μ) has presented challenges to the accurate calculation of vascular wall shear stress, the frictional force per unit area acting tangentially to the arterial wall. To date, nearly all previous studies on arterial wall shear stress have neglected the flow-dependency of blood  and improperly presented BV as a constant Newtonian fluid [6–8]. Although the non-Newtonian characteristics of blood have been assumed in some previous studies, BV over ranges of shear rates have not been measured directly, but the laws of viscosity were applied empirically [9, 10].
This study was designed to evaluate a new method for calculating vascular shear stresses along carotid arterial segments using duplex ultrasonography and dynamic BV measurements, and to examine whether the vascular shear stress of carotid artery was different between patients with lacunar infarction and control subjects using the new method (shear rate-specific BV and shear stress). This hypothesis is based on the observation that lacunar infarction is accompanied by endothelial dysfunction  and low cerebral arterial blood flow velocities . Finally, we compared the shear stresses of the new method with those of the conventional methods using other BV determinations such as hematocrit-derived BV or constant values at shear rate of 300 s-1.
Seventeen patients with lacunar infarction and 20 control subjects were examined for carotid duplex ultrasonography and measured for BV at the Department of Neurology, Chonbuk National University Hospital, Jeonju, South Korea. A diagnosis of lacunar infarction was made based upon acute onset of focal neurological deficit and a relevant lesion was seen on diffusion-weighted magnetic resonance imaging (MRI) of the brain . A lacunar infarct was diagnosed when the brain MRI identified a lesion less than 15 mm in diameter with characteristics of high signal intensity on diffusion-weighted imaging and low signal intensity on the apparent diffusion coefficient map. The location of the lacunar infarction included either the basal ganglia, corona radiata, thalamus, or brainstem, but not cortex . The control group was defined as patients with cardiovascular risk factors but no previous history of ischemic vascular diseases (e.g. stroke, myocardial ischemia, and lower extremity arterial disease). The institutional ethics committee at Chonbuk National University Hospital approved the present study. All subjects provided written informed consent.
Assessments and measurements
The medical history included data on prior CVD, type 2 diabetes mellitus, hypertension, dyslipidemia, and medication usage. Smoking status was recorded in pack-years. The patients were classified as hypertensive based on a persistent elevation of blood pressure (≥140/90 mmHg) or treatment with antihypertensive medications. Type 2 diabetes mellitus was defined by fasting blood glucose >7.0 mmol/L or previous use of glucose-lowering medications. Hypercholesterolemia was defined by a total cholesterol >6.2 mmol/L or low-density lipoprotein cholesterol >4.1 mmol/L, or current use of cholesterol-lowering medication. Plasma concentration of total homocysteine (tHcy) was measured by fluorescence polarization immunoassay (AxSYM, Abbott Laboratories, Abbott Park, IL).
In the present study, internal diameter, maximum centerline velocity, and intima-media thickness (IMT) were measured bilaterally along the common carotid artery (CCA), carotid bulb, and internal carotid artery (ICA), as shown in Figure 1. To calculate site-specific shear rates, diameters and velocities of both CCAs were used . In the patients with lacunar infarction, carotid ultrasonography was performed at least one week after the vascular attack. A higher-frequency 5-to-12 MHz (12 L5) linear transducer (Terason t3000, Teratech, Inc, Burlington, MA) was used with ECG triggering, and all the examinations were performed by a certified neurosonologist (SKJ) who has performed more than 5,000 cases (10,000 carotid arteries) during the past 8 years. Intra-luminal diameters were measured separately at peak-systole (PS) and end-diastole (ED) using two-dimensionally guided continuous M-mode tracings of the intimal-luminal interface of the near and far walls of carotid arterial segments, as reported previously . In all segments, video images of the interfaces between lumen and intima over 5 cardiac cycles were captured, stored, and the diameters were measured from fixed images. The axial resolution of the M-mode system was 0.1 mm.
Velocities were measured at both PS and ED with the sample volume reduced to the smallest possible size (1 mm) and placed in the center of flow, thereby enabling the maximum centerline velocity to be determined. The Doppler angle was generally maintained between 45 and 55 degrees. The sample volume box was placed in the mid-lumen parallel to the vessel wall in all of the studies, as recommended .
The degree of ICA stenosis at both gray-scale and Doppler ultrasonography was stratified into the following 6 categories: (1) normal (no stenosis), (2) less than 50% stenosis, (3) 50~69% stenosis, (4) ≥70% stenosis to near occlusion, (5) near occlusion, and (6) total occlusion, as recommended .
The patients with the following characteristics were excluded: 1) patients who showed ICA stenosis more than 50% in at least one side, 2) patients who have been performed for carotid artery revascularization including carotid endarterectomy or stenting. Finally, patients with lacunar infarction and control subjects, who had ICA of normal or less than 50% stenosis, were included in the present study.
According to the PS/ED velocities and diameters along the carotid arteries, PS/ED shear rates were determined as follows :
where V and D are the maximum centerline velocity and local lumen diameter, respectively. The PS and ED shear rates were calculated with corresponding PS/ED velocities and diameters, which were designated as PS shear rate ( PS) and ED shear rate ( ED).
For a measurement of BV, blood was sampled at antecubital vein and stored in ethylenediaminetetraacetic acid (EDTA) tube. BV was measured using a computerized scanning capillary viscometer (Hemathix®, Health Onvector, NJ) immediately after the carotid ultrasonography. The scanning capillary viscometer computes fluid viscosity using the Casson equation , which provides a relationship between the shear stress (τ) and shear rate () [18, 19]. The wall shear stress and shear rate in the Casson equation were calculated from flow velocity and pressure drop measurements in a U-shaped capillary tube. Qualitatively, the fluid dynamic theory used in the scanning capillary viscometer is similar to the Poiseuille flow relation, which describes BV (μ) as:
where Q is the volume flow rate of blood, d is the inside diameter of a capillary tube, △P is the pressure drop along the capillary tube length, μ is BV, and L is the length of the tube.
Shear rate () specific blood viscosity
Shear rate specific BV values were determined at the corresponding PS/ED shear rates (i.e., μ PS and μ ED) in each carotid arterial segment using the two Casson model constants (k and τy):
Blood viscosity at shear rate of 300 s-1and of hematocrit-derived
Hematocrit-derived BV (μ Hct) was calculated using the formula :
where H meant hematocrit (%)/100.
Values of BV at 300 s-1 (μ 300 s-1) were chosen from the measurements using the computerized scanning capillary viscometer.
The wall shear stresses specific to PS/ED BV values (e.g., μ 300 s-1, μ Hct, and μ PS and μ ED) and PS/ED shear rates were calculated as follows:
Descriptive data for the major characteristics were expressed as means ± standard deviations (SD) or percentage as appropriate. An independent t-test was used to determine the statistical differences in the continuous variables, whereas a chi-square test was used for categorical variables. For comparison of adjusted PS/ED shear stresses, analysis of covariance was performed separately in the carotid segments and expressed as adjusted means ± standard errors. The mean values of PS/ED shear stresses using 3 methods of BV determination were compared with paired t-test. All statistical analyses were conducted using PASW Statistics 18 (SPSS, Chicago, IL).
Seventeen patients with lacunar infarction and 20 control subjects were enrolled, as shown in Table 1. All the demographic characteristics and major CVD risk factors were not different between cases and controls. Triglyceride and HDL cholesterol showed differences of borderline significance. Plasma tHcy, hematocrit, and blood viscosities of both calculated and measured were not significantly different.
Carotid IMT showed significantly higher values in patients with lacunar infarction than the control subjects only at the right bulb region, as shown in Table 2. The left CCA also showed higher IMT values in the patients, which reached significance. On the contrary, all the PS and ED shear rates at the both-sided CCAs were significantly lower in the patients with lacunar infarction than the controls. The PS and ED shear stresses at the both CCAs using shear rate specific BVs were significantly lower in the patients with lacunar infarction than the controls, even adjusted for age, sex, and cardiovascular risk factors, as shown in Table 3. Shear stresses of the both CCAs using blood viscosities at 300 s-1 or hematocrit-derived BV also showed significant differences between the two groups. There was no significant difference of the PS and ED shear stresses at the both CCAs according to the laterality of lacunar infarction (data not shown).
PS and ED shear stresses of CCAs and blood viscosities (shear rate specific BV, BV at 300 s-1, and Hct-derived BV) were measured, as shown in Table 4. Compared with the shear stresses using shear rate specific BV, overall 4.5% higher and 7.3% lower PS and ED shear stresses were calculated using the two other BVs (Hct-derived BV or BV at 300 s-1). Blood viscosities showed 0.9% and 13.0% differences for PS and ED phases, respectively. There was no significant difference between BV of hematocrit-derived and BV at 300 s-1, and also for their calculated shear stresses.
The present study introduces a novel method for calculating site-specific vascular shear at the CCA using dynamic measures of shear rate and BV. For the calculation, Casson model constants (k and τ y ) were derived from BV data obtained using a computerized scanning capillary viscometer. In this study, we demonstrate that shear stress is a dynamic parameter that varies widely along the carotid artery. With the novel method, vascular shear stress at the both CCAs was significantly lower in patients with lacunar infarction than the control subjects. Although CVD risk factors were not significantly different between the two groups, and even carotid IMT showed variable findings, vascular shear stresses at the CCAs showed discriminative features. Low shear stress along carotid artery was reported to be associated with ischemic stroke of large artery atherosclerosis, especially in the same side to the affected hemisphere . To our knowledge, this study is the first to compare vascular shear stress in the patients with lacunar infarction with control subjects, on a shear rate specific basis using the individual subject’s own BV profile.
The vascular shear stress of the CCAs, irrespective of the methods of shear stress calculation, showed similar patterns of difference between the patients with lacunar infarction and controls, as shown in Table 3. But, there were significant differences for the mean values of CCA shear stresses and blood viscosities according to the application methods of BV. Using the shear rate specific BV, the lowest values of BV were observed at the PS phase and the highest values were detected at the ED phase. As such, if blood is assumed to behave as a Newtonian fluid and a certain representative value is used, vascular shear stresses will be calculated as higher than their true values at the PS phases, and lower in the ED phases. In the present study, when shear rate specific BVs were used, PS and ED shear stresses were 4.5% higher and 7.3% lower than shear stresses using constant values of BV, as shown in Table 4. This discrepancy would become more pronounced in arterial segments with elevated pulsatility indexes . While the diagnostic methods described herein require further validation, this rationale underscores the non-Newtonian behavior of blood flow and makes a meritorious case for further research.
Our study results are consistent with the report that lacunar infarction is associated with endothelial dysfunction , which is controlled by vascular shear stress . Previously, we reported patients with lacunar infarction had lower blood flow velocities in cerebral arteries than patients with other ischemic stroke subtypes . The both low vascular shear stress in the CCAs and low blood flow velocities in cerebral arteries of patients with lacunar infarction may foster understanding of the pathophysiology of the disease.
Low carotid arterial shear stress in the patients with lacunar infarction who have no obvious carotid stenosis suggests increased resistance in small cerebral arteries  that may potentially increase arterial pressure in the more proximate carotid arteries. The transmission of increased pressure drop of cerebral arteries would cause an increase of carotid arterial pressure, a reduction of flow velocity, a compensatory increase in arterial diameter, and a decrease of vascular shear stress. When the increase of arterial pressure and the decrease of shear stress reinforce themselves through a feedback loop toward greater instability, a vicious cycle might cause lacunar infarction. Although thrombosis due to endothelial denudation under high shear could occur [24, 25], thrombosis in the milieu of low blood flow velocities and low shear stress may also increase the propensity for ischemic events, especially at vulnerable regions . Under low shear, increased transudation of blood or its components in small arteries of the deep cerebral tissue might also occur through the increased porosity of gap junctions, the morphological changes of endothelial cells , to cause lacunar (small-sized) hemorrhage or potentially lacunar infarction .
The present study has several limitations. First, we used a local centerline velocity to determine a local shear rate, which is expressed as 4 V/D. This equation is derived from a parabolic velocity profile for a fully developed laminar flow in a straight circular tube ; so, the assumptions in this model have limitations for estimation of shear stress in regions of complex flow even in the relatively straight CCAs. However, the use of the centerline velocity together with 4 V/D expression allows the estimation of local changes in the velocity gradient if there is no significant hemodynamic compromise . Further investigations are required to determine if centerline arterial shear stress is a parameter that has some critical threshold values for vascular diseases or pathologic conditions. Our study showed 76 carotid arteries in patients with lacunar infarction at CCA peak systolic shear stress of 12~13 dyne/cm2, which is slightly lower than the previous report of carotid arterial shear stress in the subjects with high risk vascular profiles  or with carotid plaque . Second, there were differences of shear rate and shear stress between the right and left CCAs. That might be largely due to the differences of cerebral arterial resistance, and in part, the proximal arterial geometry: the right CCA comes from a celiac trunk and might have more tortuosity than the left CCA which comes from aorta directly. Third, the present cross-sectional study could not address the causal relationship between low carotid artery shear stress and lacunar infarction.
We describe a new method to calculate dynamic vascular shear stress in human carotid arteries using flow rates obtained from ultrasonography and shear rate-specific BV. With the innovative method, the carotid arterial shear stress was significantly lower in patients with lacunar infarction than the control subjects. The prognostic contribution of decreased vascular shear stress to incident lacunar infarction requires further study.
Cunningham KS, Gotlieb AI: The role of shear stress in the pathogenesis of atherosclerosis. Lab Invest. 2005, 85: 9-23.
Cecchi E, Giglioli C, Valente S, Lazzeri C, Gensini GF, Abbate R, Mannini L: Role of hemodynamic shear stress in cardiovascular disease. Atherosclerosis. 2010, 214: 249-256.
Soulis JV, Farmakis TM, Giannoglou GD, Chatzizisis YS, Giannakoulas GA, Parcharidis GE, Louridas GE: Molecular viscosity in the normal left coronary arterial tree. Is it related to atherosclerosis?. Angiology. 2006, 57: 33-40. 10.1177/000331970605700105.
Fahraeus R, Lindqvist T: The viscosity of blood in narrow capillary tubes. Am J Physiol. 1931, 96: 562-568.
Cho YI, Kensey KR: Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: Steady flows. Biorheology. 1991, 28: 241-262.
Efstathopoulos EP, Patatoukas G, Pantos I, Benekos O, Katritsis D, Kelekis NL: Wall shear stress calculation in ascending aorta using phase contrast magnetic resonance imaging. Investigating effective ways to calculate it in clinical practice. Phys Med. 2008, 24: 175-181. 10.1016/j.ejmp.2008.01.004.
Box FMA, van der Grond J, de Craen AJM, Palm-Meinders IH, van der Geest RJ, Jukema JW, Reiber JHC, van Buchem MA, Blauw GJ, for the PSG: Pravastatin decreases wall shear stress and blood velocity in the internal carotid artery without affecting flow volume: results from the PROSPER MRI Study. Stroke. 2007, 38: 1374-1376. 10.1161/01.STR.0000260206.56774.aa.
Groen HC, Gijsen FJ, van der Lugt A, Ferguson MS, Hatsukami TS, van der Steen AF, Yuan C, Wentzel JJ: Plaque rupture in the carotid artery is localized at the high shear stress region: a case report. Stroke. 2007, 38: 2379-2381. 10.1161/STROKEAHA.107.484766.
Soulis JV, Giannoglou GD, Chatzizisis YS, Seralidou KV, Parcharidis GE, Louridas GE: Non-Newtonian models for molecular viscosity and wall shear stress in a 3D reconstructed human left coronary artery. Med Eng Phys. 2008, 30: 9-19. 10.1016/j.medengphy.2007.02.001.
Gijsen FJ, Wentzel JJ, Thury A, Mastik F, Schaar JA, Schuurbiers JC, Slager CJ, van der Giessen WJ, de Feyter PJ, van der Steen AF: Strain distribution over plaques in human coronary arteries relates to shear stress. Am J Physiol Heart Circ Physiol. 2008, 295: H1608-H1614. 10.1152/ajpheart.01081.2007.
Kim JS, Lee HS, Park HY, Kim SS, Kang HG, Kim NH, Park JS, Kim Y: Endothelial function in lacunar infarction: a comparison of lacunar infarction, cerebral atherosclerosis and control group. Cerebrovasc Dis. 2009, 28: 166-170. 10.1159/000226115.
Kim JT, Lee SH, Hur N, Jeong SK: Blood flow velocities of cerebral arteries in lacunar infarction and other ischemic strokes. J Neurol Sci. 2011, 308: 57-61. 10.1016/j.jns.2011.06.013.
Norrving B: Lacunar infarcts: no black holes in the brain are benign. Pract Neurol. 2008, 8: 222-228. 10.1136/jnnp.2008.153601.
Jiang Y, Kohara K, Hiwada K: Association between risk factors for atherosclerosis and mechanical forces in carotid artery. Stroke. 2000, 31: 2319-2324. 10.1161/01.STR.31.10.2319.
Tahmasebpour HR, Buckley AR, Cooperberg PL, Fix CH: Sonographic examination of the carotid arteries. Radiographics. 2005, 25: 1561-1575. 10.1148/rg.256045013.
Grant EG, Benson CB, Moneta GL, Alexandrov AV, Baker JD, Bluth EI, Carroll BA, Eliasziw M, Gocke J, Hertzberg BS: Carotid artery stenosis: gray-scale and Doppler US diagnosis–Society of Radiologists in Ultrasound Consensus Conference. Radiology. 2003, 229: 340-346. 10.1148/radiol.2292030516.
Kim S, Cho YI, Hogenauer WN, Kensey KR: A method of isolating surface tension and yield stress effects in a U-shaped scanning capillary-tube viscometer using a Casson model. J Non-Newtonian Fluid Mech. 2002, 103: 205-219. 10.1016/S0377-0257(02)00004-6.
Rosenson RS, Wolff D, Green D, Boss AH, Kensey KR: Aspirin. Aspirin does not alter native blood viscosity. J Thromb Haemost. 2004, 2: 340-341. 10.1111/j.1538-79333.2004.0615f.x.
Jeong SK, Cho YI, Duey M, Rosenson RS: Cardiovascular risks of anemia correction with erythrocyte stimulating agents: should blood viscosity be monitored for risk assessment?. Cardiovasc Drugs Ther. 2010, 24: 151-160. 10.1007/s10557-010-6239-7.
Guyton AC, Hall JE: Textbook of Medical Physiology. 2000, Philadelphia: W. B. Saunders, 10
Carallo C, Lucca LF, Ciamei M, Tucci S, de Franceschi MS: Wall shear stress is lower in the carotid artery responsible for a unilateral ischemic stroke. Atherosclerosis. 2006, 185: 108-113. 10.1016/j.atherosclerosis.2005.05.019.
Lim MH, Cho YI, Jeong SK: Homocysteine and pulsatility index of cerebral arteries. Stroke. 2009, 40: 3216-3220. 10.1161/STROKEAHA.109.558403.
Kobayashi S, Okada K, Yamashita K: Incidence of silent lacunar lesion in normal adults and its relation to cerebral blood flow and risk factors. Stroke. 1991, 22: 1379-1383. 10.1161/01.STR.22.11.1379.
Fry DL: Acute vascular endothelial changes associated with Increased blood velocity gradients. Circ Res. 1968, 22: 165-197. 10.1161/01.RES.22.2.165.
Li ZY, Howarth SP, Tang T, Gillard JH: How critical is fibrous cap thickness to carotid plaque stability? A flow-plaque interaction model. Stroke. 2006, 37: 1195-1199. 10.1161/01.STR.0000217331.61083.3b.
Nesbitt WS, Westein E, Tovar-Lopez FJ, Tolouei E, Mitchell A, Fu J, Carberry J, Fouras A, Jackson SP: A shear gradient-dependent platelet aggregation mechanism drives thrombus formation. Nat Med. 2009, 15: 665-673. 10.1038/nm.1955.
Levesque MJ, Liepsch D, Moravec S, Nerem RM: Correlation of endothelial cell shape and wall shear stress in a stenosed dog aorta. Arterioscler Thromb Vasc Biol. 1986, 6: 220-229. 10.1161/01.ATV.6.2.220.
Labovitz DL, Boden-Albala B, Hauser WA, Sacco RL: Lacunar infarct or deep intracerebral hemorrhage. Neurology. 2007, 68: 606-608. 10.1212/01.wnl.0000254619.98089.43.
Sutera SP, Skalak R: The history of Poiseuille’s law. Annu Rev Fluid Mech. 1993, 25: 1-19. 10.1146/annurev.fl.25.010193.000245.
Irace C, Cortese C, Fiaschi E, Carallo C, Farinaro E, Gnasso A: Wall shear stress is associated with intima-media thickness and carotid atherosclerosis in subjects at low coronary heart disease risk. Stroke. 2004, 35: 464-468. 10.1161/01.STR.0000111597.34179.47.
Gnasso A, Irace C, Carallo C, De Franceschi MS, Motti C, Mattioli PL, Pujia A: In vivo association between low wall shear stress and plaque in subjects with asymmetrical carotid atherosclerosis. Stroke. 1997, 28: 993-998. 10.1161/01.STR.28.5.993.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2377/13/36/prepub
This article was supported by Fund of Biomedical Research Institute of Chonbuk National University Hospital and the measurement of blood viscosity was supported by Medical Device Clinical Trial Center (MDCTC), Chonbuk National University Hospital.
The authors declare that they have no competing interest.
SKJ participated in data collection, conceived the study, drafted and reviewed the manuscript. RSR conceived the present study and reviewed the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Jeong, S., Rosenson, R.S. Shear rate specific blood viscosity and shear stress of carotid artery duplex ultrasonography in patients with lacunar infarction. BMC Neurol 13, 36 (2013). https://doi.org/10.1186/1471-2377-13-36
- Vascular shear stress
- Blood viscosity
- Shear rate