Original Contributions |
From the Departments of Physiology (L.K., J.L., R.S.R.) and Biophysics (A.P.G.H.), Cardiovascular Research Institute Maastricht, Maastricht University, Maastricht, The Netherlands.
Correspondence to L. Kornet, PhD, Department of Physiology, Cardiovascular Research Institute Maastricht, Maastricht University, PO Box 616, 6200 MD Maastricht, The Netherlands. E-mail: L.Kornet{at}FYS.Unimaas.nl
| Abstract |
|---|
|
|
|---|
Key Words: carotid artery intima-media thickness shear rate ultrasound
| Introduction |
|---|
|
|
|---|
Recent developments in ultrasonography have made it possible to assess noninvasively IMT and near-wall shear rate (WSR), and hence, near-wall shear stress, in humans at the same site in the common carotid artery.6 This makes it possible to relate areas of low or high WSR to IMT in humans at various sites in 1 artery. In the common carotid artery, reflections originating from the periphery and the flow divider may affect the shape of the time-dependent flow velocity profile, and hence, the time-dependent WSR, differently just before the bifurcation and 20 to 30 mm farther proximally. Therefore, we investigated whether both locations are indeed subjected to a different WSR, and hence, a different near-wall shear stress, and whether this difference, if any, affects IMT and mechanical arterial wall properties differently at these locations. Moreover, the effect of age on the possible differences was assessed because IMT is known to increase with age.7 8 9 10 11 The common carotid artery was chosen for this study because it is easy to assess with ultrasonography and because of the pronounced difference in the echocardiographic density of the intima and media, which facilitates assessment of IMT.
| Methods |
|---|
|
|
|---|
Characterizing the Behavior of Blood
Shear stress is the product of shear rate and viscosity.
Local blood viscosity is difficult to assess. The assumption that
viscosity is the same at 2 locations closely spaced within 1 carotid
artery allows the use of shear rate as an indicator of shear stress for
direct comparison at both locations. WSR is the radial derivative of
velocity (v) at the wall:
WSR=dv(r)/dr||r=R,
where v(r) is the velocity distribution,
determined as a function of radial position r, and
R is the local radius of the artery. WSR can be determined
directly by considering the velocity gradient of blood at the
wall-lumen boundary as determined by ultrasonography (see below).
Protocol
Data gathering started after an acclimation period of 10 to 15
minutes. Subjects were in the supine position with the head slightly
tilted in the contralateral direction. Each subject was connected to an
electrocardiograph to generate a pulse to signal the onset of a cardiac
cycle. In each subject, only 1 carotid artery was investigated, because
in a previous study, no significant differences in the
parameters to be measured were found between the left and
right carotid arteries.12 Measurements were
obtained using the right carotid artery because the coefficient of
variation of WSR was slightly smaller for the right than for the left
common carotid artery in case the investigator was
right-handed.13 First, the wall of the common
carotid artery was visualized just proximal to the carotid artery
bifurcation (Figure 1
, location B), where
the artery was not yet widened, and the carotid artery bifurcation was
checked for the presence of plaque or stenosis. If present,
the subject was excluded. If not, the IMT of the posterior wall was
determined 6 times at location B.14 Next, the IMT
was determined 6 times at a position 20 to 30 mm proximal to the
bifurcation (Figure 1
, location A). Subsequently, at both locations we
determined (18 times) the time-dependent velocity profile and the
arterial diameter over a complete cardiac
cycle.6 13 The time-dependent velocity
distribution (Figure 2A
and 2B
) was used
to calculate the time-dependent WSRs at the posterior
(WSRp) and anterior walls
(WSRa) (Figure 2C
). For each location, maximum
and mean WSRp and WSRa were
derived from the WSR waveforms. Only the WSRs determined at the
posterior wall were used to study its relation with IMT, which was also
determined at the posterior wall. End-diastolic diameter
and distension (systolic minus diastolic diameter)
were extracted from the diameter waveforms obtained
simultaneously with the velocity distribution. The values
averaged over 18 measurements were taken as the subject's reading.
Arterial blood pressure and heart rate were determined
automatically and noninvasively every 3 minutes from blood pressure
measurements obtained from the brachial artery using a cuff
(Dinamap, Critikon).
|
|
Assessment of IMT
The method used to assess IMT is described in detail
elsewhere.14 In brief, the common carotid artery was
visualized in B mode (ATL Mark 9, Advanced Technology Laboratories).
Using the tip of the flow divider as a landmark, the position of the
L10-5 38-mm probe (operating frequency, 5 to 10 MHz) was manipulated
until a suitable line of sight was obtained (Figure 1
, solid vertical
lines). Along this line, starting synchronously with a trigger derived
from the top of the R wave on the electrocardiogram,
radiofrequency (RF) data were collected at a sample frequency of 20 MHz
over a period of 4 seconds. After acquisition of RF data, the first
line was displayed on a personal computer screen, allowing
identification of a window of 3 mm covering the posterior
lumen-wall transition. Data from this window over time were then stored
on hard disk for further off-line processing. Because of
arterial wall motion, averaging the power distribution over
time along a line will reduce speckle artifacts. Wall motion, however,
will complicate processing, because echocardiographic
signals will be displaced over time along the line of view. To avoid
the latter complication, RF signals were aligned in phase before
averaging on the basis of the displacement detected between
observations. Subsequently, an edge detection algorithm was applied to
the average envelope (amplitude distribution) of the processed RF
signals. The intima position was assigned to the point halfway along
the first up slope, where the spatial derivative of the envelope
exceeded for the first time a preselected level. The same procedure was
repeated for the next significant upstroke, which is the
media-adventitia transition. The outer boundary of the adventitia
cannot be distinguished. The difference between the intima position and
the media-adventitia transition was taken as the IMT. In vivo
registrations exhibit a variation in the order of 45
µm.14
Assessment of WSR, Diameter, and Distension
The ultrasonic Shear Rate Estimating System has been described
in detail elsewhere.6 13 The ultrasound echo
system (ATL Mark 9, HDI), in combination with dedicated signal
processing,6 is able to measure the blood flow
velocity distribution along a selected line of observation across the
center of the vessel with a spatial resolution of 300 µm and a
temporal resolution of 10 milliseconds. To obtain a more detailed
velocity distribution in both the spatial resolution and temporal
direction, intermediate values were calculated using half-overlapping
data segments. The common carotid artery was first visualized in B mode
using the C9-5 ICT probe (operating frequency, 5 to 9 MHz). After
positioning the M line (Figure 1
, dashed line), the ultrasound system
was switched to echo M mode with a high pulse-repetition frequency. A
registration started synchronously with a trigger derived from the top
of the R wave on the electrocardiogram, facilitating
the detection of the maximum (systolic) and minimum
(diastolic) velocity as well as the initial diameter. The
captured RF signals (reflected and scattered ultrasound signals) were
digitized at 20 MHz and transferred to the memory of the computer. The
size of the memory allows recording over 1.2 seconds, which is
normally sufficient for capturing 1 complete heart beat. The first
digitized RF line, as a function of depth, was displayed on the
computer screen. From the shape and position in depth of the
reflections, the wall-lumen interfaces on both sides were identified
manually by placing sample volumes, indicated by markers, on the
reflections of the anterior and posterior vessel walls. The distance
between both markers, corrected for the angle of observation (70°),
was considered the initial (end-diastolic) diameter of the
vessel. Processing of the RF data within the sample volumes as a
function of time, with the sample volumes adaptively tracking the
observed displacement, resulted in the time-dependent change of the
arterial wall position (displacement waveform). The
difference between the displacement curves of both walls reflected the
pulsatile change in diameter over time (ie, distension as a function of
time). To obtain the time-dependent blood flow velocity distribution, a
modeled cross-correlation function was applied to the RF data between
the markers to estimate the mean velocity over time segments of 10
milliseconds spaced at 5-milliseconds (50% overlap) time intervals.
The length of the RF segments corresponded to 300 µm in depth,
and the segments were spaced at 150-µm intervals (50%
overlap).6 15 Calculating the mean velocity of
all RF segments resulted in a time-dependent velocity profile, which
was corrected for the angle of observation (70°). The shear rate
distribution followed from the radial derivative of the velocity
profile at each time instant. The maximum value of the first derivative
of the flow profile relative to the radius toward the posterior wall of
the vessel was considered the estimate of the instantaneous
longitudinal WSRp. From the velocity and shear
rate distributions, such parameters as
PWSRP (per second) and the anterior wall
PWSRa (per second) at peak systole,
MWSRP (per second), and
MWSRa (per second), which is the time-averaged
shear rate over 1 cardiac cycle, were determined. The average
intrasubject intrasession variability per measurement is
15% for
PWSR and
12% for MWSR.12 The intrasubject
intersession variability for averaged values (n=15) is
3% for PWSR
and
2% for MWSR.12 Therefore, in the
present study, we averaged >15 measurements.
End-diastolic diameter (d, in
micrometers) and distension (
d, in
micrometers) were assessed simultaneously with
shear rate.
Statistics
To study the differences between a parameter
determined at both locations, a paired t test was used.
Linear regression was used to study the influence of age. A significant
correlation with age was present if the derivative of the linear
regression line was significantly different from zero. Significance was
reached at P
0.05.
| Results |
|---|
|
|
|---|
IMT was larger at location B than at location A (Figure 3A
and 3B
). In all individuals, the
relative spatial difference in IMT was not significantly affected by
age (Figure 3B
). At both locations, IMT, determined at the posterior
wall, increased significantly with age (Figure 3A
).
|
PWSRP (per second) and
MWSRP (per second) were significantly lower at
location B than at location A (Figure 4A
through 4D). The relative spatial differences in
PWSRP and MWSRP were not
affected by age. PWSRP and
MWSRP decreased significantly with age at both
locations (Figure 4A
and 4C
). PWSRP and
MWSRP were not correlated with heart
rate.
|
IMT was negatively correlated with PWSRP
(Figure 5A
) and
MWSRP (Figure 5B
) at location B but not at
location A. At both locations, IMT was positively correlated with age
(Figure 3A
), systolic blood pressure, and body mass index,
which was defined as body weight divided by squared height (data not
shown). At both locations, IMT was not correlated with heart rate.
|
The WSRs determined at locations B and A at the anterior and posterior
walls, as well as the differences between them, are shown in the
Table
. MWSR and PWSR at both the
posterior and anterior walls were significantly higher at location A
than at location B. The differences between the MWSRs and PWSRs between
both locations were a factor of 2 to 3 larger at the posterior wall
than at the anterior wall. At location A, MWSR was significantly higher
at the posterior wall than at the anterior wall, but at location B, the
MWSR values were not significantly different between the walls. At both
locations, PWSR was significantly higher at the posterior than at the
anterior wall. The absolute difference was a factor of 3 larger at
location A than at location B.
|
| Discussion |
|---|
|
|
|---|
Wall shear stress is the product of viscosity and WSR.
Whole-blood viscosity is influenced by wall shear
rate.16 Because of the larger WSR at location A
than at location B, we calculated the viscosity to be, on average,
4.1% (95% confidence interval, =2.8% to 5.4%) smaller at location A
than at location B. Therefore, the difference in near-wall shear stress
between locations A and B, measured to be 23% of that at location A,
is overestimated by 3%. Furthermore, it should be noted that the
maximum first derivative of the flow profile, not the actual WSR, was
considered.6 Shear rate at the wall cannot be
determined because of the low scanning resolution relative to the
arterial diameter.6 Because the
system is activated with a short pulse and because the spatial
resolution, determined by the length of the data window, matches the
system resolution, the eventual resolution along the ultrasound beam is
0.3 mm. The spatial resolution is hardly affected by the
ultrasound beam width (
1 mm) because of the steep observation
angle (70o), but this makes the measurement more
sensitive to minor deviations of the observation angle. To overcome the
problem of low scanning resolution, the maximum shear rate as assessed
was considered as WSR. At the posterior wall, we may assume the maximum
shear rate to be
0.3 mm from the blood-intima boundary. This
implies that the measured shear rate may underestimate WSR when the
actual velocity gradient is steeper distally than at the site of
assessment of maximum shear rate. Therefore, the shear rate as
presented may be considered the lowest possible estimate,
because theoretically it is unlikely that the velocity gradient will be
less steep closer to the wall.
The lower PWSR (Figure 4A
and 4B
) and MWSR (Figure 4C
and 4D
) at
the posterior wall near the bifurcation than at 20 to 30 mm
farther proximally cannot be explained by a difference in diameter or
distension of the artery because these parameters were not
significantly different between both locations. We therefore conclude
that the relatively lower PWSR at the posterior wall near the
bifurcation is related to a difference in shape of the time-dependent
velocity profile. Typical examples of flow profiles at locations A and
B are shown in Figure 2B
. The profiles presented in this figure
show that the velocity profile is more blunted at the more proximal
location in the common carotid artery (location A) than near the
bifurcation (location B). These differences in velocity profiles as
recorded near the bifurcation and more upstream were observed in
all vessels studied. These differences in velocity profile may be
explained by a different influence of reflections originating from the
external carotid artery (Figure 6
). This
explanation, however, is speculative because we were not able to
quantify the shape of the velocity profiles because they change with
time and are asymmetric (Figure 2B
). At location B, only a relatively
small part of the velocity profile is affected, and MWSR is not
significantly different between the posterior and anterior walls (the
Table
). At location A, the reflections are "diluted" and more
spread out over the whole cross section of the vessel, leading to a
more blunted velocity profile and higher PWSRs and MWSRs at both walls
than at location A. The increase is more pronounced at the posterior
wall than at the anterior wall (the Table
), because reflections are not
completely axially oriented but oriented more toward the posterior than
the anterior wall.
|
By relating stationary flow behavior in scale models of human carotid artery bifurcations to intimal thickening in a corresponding series of autopsy specimens, lower WSR, and thus near-wall shear stress, was also found to be correlated inversely with IMT.1,2 In the current study, however, we were able to demonstrate that there is a direct relation between the level of WSR and the degree of intima-media thickening because of the comparable circumstances at both sites of measurement. It is of interest to note that average differences in PWSR and MWSR of only 23% result in a difference in IMT of 29%. The correlation between IMT on one hand and PWSR and MWSR on the other was significant only near the bifurcation, because at this location the increase in IMT with age was more pronounced than at the more proximal location. The observation that arterial distension was not significantly different near the bifurcation and more upstream in the common carotid artery, despite a significant difference in IMT, may be explained by the finding that mainly glycosaminoglycans, a relatively elastic protein, accumulates in the vascular extracellular matrix with increasing age.17 Riley et al18 also observed that thicker vessel walls are not necessarily stiffer than thinner ones.
Two explanations for the negative correlation between WSR or stress and wall thickness may be considered. First, the increased residence time of particles on the inner side of the vessel wall might lead to diffusion of particles into the wall, resulting in a thicker wall. Platelets and macrophages, key elements of atherosclerotic lesions, are more likely to adhere to the arterial wall in regions of increased residence time.19 Second, wall shear stress causes the endothelial cells to orient in the prevailing direction of flow.20 21 Regions of enhanced shear stress are characterized by more elongated endothelial cells, whereas regions of relatively low shear stress are associated with more rounded endothelial cells.22 In the latter regions, cell turnover rates are higher than in the former ones.23 During cell turnover, intercellular junctions become leaky, allowing for an enhanced influx of lipids and other macromolecules. Therefore, the entrance into the wall of macromolecules, as lipids, may be enhanced in low-shear regions.24 Whether accumulation of compounds actually occurs depends on the resistance of the underlying wall, especially the media, to the transport of such molecules across the wall to the adventitia, where they will be removed.25
In summary, in the common carotid artery of humans, WSR is lower near the bifurcation than at more proximal locations. IMT is greater at the site of lower WSR. The relative differences between both locations within a subject are independent of age. We propose that the differences in WSR between both locations are caused by a different effect of reflections originating from the external carotid artery.
Received November 5, 1997; accepted May 19, 1998.
| References |
|---|
|
|
|---|
2.
Ku DN, Giddens DP, Zarins CK, Glagov S.
Pulsatile flow and atherosclerosis in the human carotid
bifurcation. Arteriosclerosis. 1985;5:293302.
3. Friedman MH, Bargeron CB, Deters OJ, Hutchins GM, Mark FF. Correlation between wall shear and intimal thickness at a coronary branch. Atherosclerosis. 1987;68:2733.[Medline] [Order article via Infotrieve]
4. Friedman MH, Hutchins MG, Bargeron CB, Deters OJ, Mark FF. Correlation between intimal thickness and fluid shear in human arteries. Atherosclerosis. 1981;39:425436.[Medline] [Order article via Infotrieve]
5.
Gnasso A, Carallo C, Irace C, Spagnuolo V, De Novara
G, Mattioli PL, Pujia A. Association between intima-media thickness and
wall shear stress in common carotid arteries in healthy male subjects.
Circulation. 1996;94:32573262.
6. Brands PJ, Hoeks APG, Hofstra L, Reneman RS. A noninvasive method to estimate wall shear rate using ultrasound. Ultrasound Med Biol. 1995;21:171185.[Medline] [Order article via Infotrieve]
7. Bots ML, Witteman JCM, Grobbee DE. Carotid intima-media wall thickness in elderly women with and without atherosclerosis of the abdominal aorta. Atherosclerosis. 1993;102:99105.[Medline] [Order article via Infotrieve]
8.
Bonithon-Kopp C, Touboul P, Berr C, Leroux C,
Mainard F, Courbon D, Ducimetiere P. Relation of intima-media thickness
to atherosclerotic plaques in carotid arteries. Arterioscler
Thromb Vasc Biol. 1996;16:310316.
9. Salonen R, Salonen JT. Determinants of carotid intima-media thickness: a population-based ultrasonography study in eastern Finnish men. J Intern Med. 1991;229:225231.[Medline] [Order article via Infotrieve]
10.
Riley WA, Barnes RW, Evans GW, Burke GL. Ultrasonic
measurement of the elastic modulus of the common carotid artery: the
Atherosclerosis Risk in Communities (ARIC) study.
Stroke. 1992;23:952956.
11. Hodges TC, Detmer PR, Dawson DL, Bergelin RO, Beach KW, Hatsukami TS, Zierler BK, Isaacson JA, Strandness DE. Ultrasound determination of total arterial wall thickness. J Vasc Surg. 1994;19:745753.[Medline] [Order article via Infotrieve]
12. Samijo SK, Willigers JM, Brands PJ, Barkhuysen R, Reneman RS, Kitslaar PJEHM, Hoeks APG. Reproducibility of shear rate and shear stress assessment by means of ultrasound in the common carotid artery of young human males and females. Ultrasound Med Biol. 1997;23:583590.[Medline] [Order article via Infotrieve]
13.
Hoeks APG, Samijo SK, Brands PJ, Reneman RS.
Noninvasive determination of shear rate distribution across the
arterial lumen. Hypertension. 1995;26:2633.
14. Hoeks APG, Willekes C, Boutouyrie P, Brands PJ, Willigers JM, Reneman RS. Automated detection of artery wall thickness based on M-line signal processing. Ultrasound Med Biol. 1997;23:10171023.[Medline] [Order article via Infotrieve]
15. Hoeks APG, Samijo SK, Brands PJ, Reneman RS. Assessment of wall shear rate in humans: an ultrasound study. J Vasc Invest. 1995;1:108117.
16. Weaver JPA, Evans A, Walder DN. The effect of increased fibrinogen content of the viscosity of blood. Clin Sci. 1969;36:110.[Medline] [Order article via Infotrieve]
17. Bilato C, Crow MT. Atherosclerosis and the vascular biology of aging. Aging (Milano). 1996;8:221234.[Medline] [Order article via Infotrieve]
18. Riley WAE, Evans GW, Sharrett AR, Burke GL, Barnes RW. Variation of common carotid artery elasticity with intimal-media thickness: the ARIC study. Ultrasound Med Biol. 1997;23:157164.[Medline] [Order article via Infotrieve]
19. Parmentier EM, Morton WA, Petschek HE. Platelet aggregate formation in a region of separated blood flow. J Biochem Eng. 1981;20:20122021.
20. Nerem RM, Girard PR. Hemodynamic influences on vascular endothelial biology. Toxicol Pathol. 1990;18:572582.[Medline] [Order article via Infotrieve]
21.
Flaherty JT, Pierce JE, Ferrans VJ, Patel DJ, Tucker
WK, Fry DL. Endothelial nuclear patterns in the canine
arterial tree with particular reference to
hemodynamic events. Circ Res. 1972;30:2333.
22. Levesque MJ, Nerem RM. The elongation and orientation of cultured endothelial cells in response to shear stress. J Biomech Eng. 1985;107:341347.[Medline] [Order article via Infotrieve]
23. Caplan BA, Schwartz CJ. Increased endothelial cell turnover in areas of in vivo Evans blue uptake in the pig aorta. Atherosclerosis. 1973;17:401417.[Medline] [Order article via Infotrieve]
24. Helmlinger G, Geiger RV, Schreck S, Nerem RM. Effects of pulsatile flow on cultured vascular endothelial cell morphology. J Biomech Eng. 1991;113:123131.[Medline] [Order article via Infotrieve]
25. Nerem RM. Vascular fluid mechanics, the arterial wall, and atherosclerosis. J Biomech Eng. 1992;114:274282.[Medline] [Order article via Infotrieve]
This article has been cited by other articles:
![]() |
M. J. Post and J. Waltenberger Small Is Beautiful: A Miniature Stent Model Arterioscler Thromb Vasc Biol, April 1, 2007; 27(4): 701 - 702. [Full Text] [PDF] |
||||
![]() |
C. A. Ruiz-Feria, Y. Yang, and H. Nishimura Do incremental increases in blood pressure elicit neointimal plaques through endothelial injury? Am J Physiol Regulatory Integrative Comp Physiol, December 1, 2004; 287(6): R1486 - R1493. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. L. Schott, R. P. Wildman, S. Brockwell, L. R. Simkin-Silverman, L. H. Kuller, and K. Sutton-Tyrrell Segment-Specific Effects of Cardiovascular Risk Factors on Carotid Artery Intima-Medial Thickness in Women at Midlife Arterioscler Thromb Vasc Biol, October 1, 2004; 24(10): 1951 - 1956. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Tanaka, F. A. Dinenno, K. D. Monahan, C. A. DeSouza, and D. R. Seals Carotid Artery Wall Hypertrophy With Age Is Related to Local Systolic Blood Pressure in Healthy Men Arterioscler Thromb Vasc Biol, January 1, 2001; 21(1): 82 - 87. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Zureik, P. Ducimetiere, P.-J. Touboul, D. Courbon, C. Bonithon-Kopp, C. Berr, and C. Magne Common Carotid Intima-Media Thickness Predicts Occurrence of Carotid Atherosclerotic Plaques : Longitudinal Results From the Aging Vascular Study (EVA) Study Arterioscler Thromb Vasc Biol, June 1, 2000; 20(6): 1622 - 1629. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. S Reneman and A. P.G Hoeks Noninvasive vascular ultrasound: An asset in vascular medicine Cardiovasc Res, January 1, 2000; 45(1): 27 - 35. [Full Text] [PDF] |
||||
![]() |
L. Kornet, A. P. G. Hoeks, J. Lambregts, and R. S. Reneman In the Femoral Artery Bifurcation, Differences in Mean Wall Shear Stress Within Subjects Are Associated With Different Intima-Media Thicknesses Arterioscler Thromb Vasc Biol, December 1, 1999; 19(12): 2933 - 2939. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
ATVB Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 1998 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |