Accepted Manuscript Extended Eversion Carotid Endarterectomy: Computation of Hemodynamics Tomislav Istvanić, Zvonimir Vrselja, Hrvoje Brkić, Radivoje Radić, Igor Lekšan, Goran Curic PII:
S0890-5096(15)00620-2
DOI:
10.1016/j.avsg.2015.05.034
Reference:
AVSG 2504
To appear in:
Annals of Vascular Surgery
Received Date: 13 February 2015 Revised Date:
14 May 2015
Accepted Date: 26 May 2015
Please cite this article as: Istvanić T, Vrselja Z, Brkić H, Radić R, Lekšan I, Curic G, Extended Eversion Carotid Endarterectomy: Computation of Hemodynamics, Annals of Vascular Surgery (2015), doi: 10.1016/j.avsg.2015.05.034. This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
ACCEPTED MANUSCRIPT 1
Extended Eversion Carotid Endarterectomy: Computation of Hemodynamics
2
Tomislav Istvanića,* , Zvonimir Vrseljab,c,*, Hrvoje Brkićd, Radivoje Radićc, Igor Lekšanc,
3
Goran Curice,**¶
5
*These authors equally contributed to the manuscript.
6
a
7
b
8
c
9
Croatia
RI PT
4
Department of Vascular Surgery of Clinic for Surgery, University Hospital Osijek, Croatia Department of Radiology, University Hospital Osijek, Croatia
10
d
11
University of Osijek, Croatia
12
e
13
Medicine, University of Osijek, Croatia
M AN U
SC
Department of Anatomy and Neuroscience, Faculty of Medicine, University of Osijek,
Department of Biophysics, Medical Statistics and Medical Informatics, Faculty of Medicine,
Department of Medical Chemistry, Biochemistry and Clinical Chemistry, Faculty of
TE D
14
**Correspondence to Goran Curic, Faculty of Medicine, University of Osijek, J. Huttlera 4,
16
31000 Osijek, Croatia.
[email protected]. Phone +385911612584. Fax +38531505615.
18
AC C
17
EP
15
1
ACCEPTED MANUSCRIPT ABSTRACT
2
Background: Stroke prevention includes surgery for significant stenosis of internal carotid
3
artery. Consensus on a standard approach lacks and one alternative approach is eversion
4
carotid endarterectomy (eCEA). In order to overcome disadvantages of eCEA, we developed
5
extended eversion carotid endarterectomy (exeCEA). Aiming to investigate hemodynamics
6
after different surgical approaches, we created computational fluid dynamics models of
7
exeCEA and eCEA with included progressing lumen narrowing - representation of artery
8
restenosis at the incision line.
9
Methods: Blood flow velocities, volume flow rates and wall shear stress were established in
SC
RI PT
1
carotid arteries from models of eCEA and exeCEA with included increasing groove (1, 1.5, 2
11
and 2.5 mm) at the ‘incision line’, under input pressure of 120 and 150 mmHg.
12
Results: For the corresponding restenosis grade, models of exeCEA had a larger orifice
13
towards internal carotid artery, lower blood flow velocities and higher volume flow rates in
14
internal carotid artery, with lower volume flow rates in external carotid artery. Wall shear
15
stress values in internal carotid artery of exeCEA models were lower than in eCEA models,
16
later reaching the thrombotic range.
17
Conclusions: Computational fluid dynamics showed better hemodynamic properties in
18
exeCEA models, indicating presented approach might be better at preserving brain perfusion.
19
Keywords: carotid endarterectomy, circulation, blood flow, wall shear stress, computational
20
fluid dynamics
TE D
EP
AC C
21
M AN U
10
2
ACCEPTED MANUSCRIPT 1. Introduction
2
Stroke is the second most common cause of death and stroke survivors have the highest risk
3
for a recurrent stroke [1]. About 30% of all strokes is considered to be associated with
4
stenosis of carotid artery (i.e., due to reduction of blood flow, or thrombi developing at
5
atherosclerotic plaque) [1]. In order to reduce absolute risk for stroke, prevention strategy
6
includes eversion carotid endarterectomy (eCEA): surgery of stenotic internal carotid artery
7
(ICA) [2]. Carotid endarterectomy (CEA) is recommended when the lumen of the ICA is
8
reduced more than 70%, as documented by noninvasive imaging, i.e. Doppler
9
ultrasonography. CEA is performed in both (I) symptomatic (history of ipsilateral minor
RI PT
1
stroke or transient ischemic attack) and (II) asymptomatic patients (evident stenosis at carotid
11
bifurcation).
12
Common approaches are conventional endarterectomy with patch angioplasty and eCEA;
13
conventional CEA is performed more often, although some large studies slightly favor eCEA
14
[2-6]. These procedures have similar proportions of postoperative stroke, myocardial
15
infarction and recurrent stenosis [4, 6, 7]. Restenosis develops primarily because of intimal
16
hyperplasia at the site of vessel injury [8, 9]. (Re)stenosis increases blood velocity that can
17
lead to increased turbulence. Altered hemodynamics corresponds with sites of development of
18
vascular wall pathology (intimal thickening and atherosclerosis) and it is associated with
19
altered properties of blood (hypercoagulability) [10]. CEA alters the endothelial lining and
20
changes flow pattern; blood velocity, volume flow rate and wall shear stress (WSS), although
21
effects on hemodynamics are partially understood [11].
22
The eCEA has several disadvantages; difficulty of intraoperative shunt insertion (if needed),
23
possible baroreceptor damage, lack of association of reduced restenosis rates and stroke, and
24
inadequate endarterectomy when plaque extends well into CCA [4]. Guided by surgical
25
experience, knowledge of anatomy and hemodynamics, we projected that modification of
26
incision might overcome some disadvantages of eCEA. ‘Hockey stick’ shaped incision in
27
extended eCEA (exeCEA) involves transection of CCA, from bifurcation distally, until
28
reaching proximal end of plaque, where incision ends by extending under angle of 45°
29
towards ICA (Fig. 1).
30
We modeled 3D flow simulations in carotid vessels after exeCEA and eCEA procedure, and
31
tested hemodynamics after lumen narrowing (restenosis at site of vessel cut) to assess the
32
effects of local geometry changes on volume flow rates. Here, we present results of
AC C
EP
TE D
M AN U
SC
10
3
ACCEPTED MANUSCRIPT 1
computational fluid dynamics (CFD) models of restenosis at incision line after simulated
2
eCEA and exeCEA.
3
2. Material and methods
5
The 3D carotid arterial tree (physiological, lumen obstruction free) model was created on the
6
basis of MR angiography of a healthy middle age female, after obtaining informed consent.
7
Local ethical committee approved the study. The 3Dslicer software [12] extracted data from
8
MR. Raw model was processed with MeshLab (MeshLab Visual Computing Lab - ISTI -
9
CNR http://meshlab.sourceforge.net/) software and further processed with SolidWorks,
SC
RI PT
4
software for building models and managing flow simulations.
11
The eCEA and exeCEA were mimicked (Fig. 1) on created 3D carotid model. Models
12
included simplified folds that represented intimal hyperplasia or recurrent atherosclerotic
13
disease that narrow the vessel diameter. Increasing folds at ‘incision line’ (1, 1.5, 2 and 2.5
14
mm) – represented restenosis grade I-IV, respectively (accordingly referred in the text).
15
Mathematical approximations of incision cross-sectional areas showed that orifice after
16
exeCEA is approximately 1.6 times larger (Calculation shown in Supplemental Figure 1). In
17
the simulations, the blood was defined as a non-Newtonian liquid (density of 1 g/cm3), wall
18
material was set as rubber (elastic modulus 6.2 N/mm2), temperature was 310.15 K (37°C),
19
surrounding pressure of 101325 Pa, and adiabatic wall conditions were assumed.
20
Volume flow over time at model inlet lid (beginning of CCA) was adjusted to mimic a
21
pulsatile flow (single heart cycle) [13]. To assess the effects of changing local geometry
22
(restenosis) during ‘usual’ daily pressure oscillations in operated patients, outlet boundary
23
condition was set as a constant pressure, with value of 120 or 150 mmHg [14]. Since the
24
volume flow rates between eCEA and exeCEA procedure were of primary interest, flow rate
25
partition at carotid bifurcation was not fixed outlet boundary condition. Although, it is
26
recommended to avoid the coupling of unsteady flow and multiple outflows with constant
27
pressure as outlet boundary condition [15], we chose this setting in order to assess the relative
28
volume flow rate distribution (see Limitations for details). Because the relative values
29
between eCEA and exeCEA were of interest, we find that used simplifications were mitigated
30
- as only differing variable was the shape of incision line. Similarly, differing resistances of
31
ICA and ECA were not included in the flow simulations since only relative hemodynamic
AC C
EP
TE D
M AN U
10
4
ACCEPTED MANUSCRIPT values were of interest (eCEA vs. exeCEA). Data acquisition of hemodynamic parameters
2
was performed in the middle third of CCA, and at middle thirds of modeled ICA and ECA,
3
distally from narrowing. WSS size area was analyzed using open source ImageJ software
4
[16]. According to previous studies, WSS values between 3.5 Pa and 7 Pa were categorized as
5
“high WSS” [17, 18], while WSS values above 7 Pa were categorized as “thrombotic WSS”
6
[10]. Calculation of areas under curve (AUC), using trapezoidal approximation, was used for
7
comparison of hemodynamic data.
8
2.1. Statistical analysis
9
IBM SPSS Statistics v15.0 was used for statistical analyses. Normality of distribution was
SC
RI PT
1
tested with Shapiro – Wilks test. Time dependent blood velocities and volume flow rates for a
11
given locations were compared using repeated measures T-test. A two-sided p<0.05 was
12
considered significant.
M AN U
10
13 14
3. RESULTS
16
3.1. Blood velocities
17
Blood velocities in CCA were similar in both models for corresponding restenosis grades.
18
Overall, blood velocities in ICA of exeCEA models varied less with progressing restenosis
19
and were lower than in corresponding eCEA models (Table 1.). ECA blood velocities
20
increased with progressing restenosis in all models and were higher in eCEA models (Table
21
2.). In models of restenosis grade IV, blood velocities in ICA continued to fall in the same
22
rate, but blood velocities in ECA increased sharply, more prominently in eCEA models
23
(Tables 1 and 2).
24
3.2. Volume flow rates
25
CCA volume flow rates were constant in all models. ICA volume flow rates decreased more
26
with progressing restenosis in eCEA models, so exeCEA models had significantly higher flow
27
rates in ICA, with the exception of restenosis grade III model for 150 mmHg pressure
28
(p=0.373). ICA peak flow rates for restenosis grade II-IV were lower in eCEA models than in
29
physiological and exeCEA models (data not shown).
AC C
EP
TE D
15
5
ACCEPTED MANUSCRIPT ECA volume flow rates increased with increasing systolic pressure and progressing
2
restenosis, and increased more prominently in eCEA restenosis models (Table 2). Similarly to
3
ICA volume flow rates, ECA volume flow rates did not differ only in grade III restenosis
4
models under simulated pressure of 150 mmHg (p=0.131).
5
3.3 Wall shear stress
6
Under the pressure of 120 mmHg and during maximal volume flow, ICA in eCEA models
7
across all restenosis grades, on average, had 43% more vessel’s surface exposed to high WSS,
8
out of which, on average, 11% of the vessel surface was exposed to the thrombotic WSS
9
(Fig. 2 and 3 A). Similar results were obtained under the pressure of 150 mmHg during
SC
RI PT
1
maximal volume flow, where ICA in eCEA models across all restenosis grades, on average,
11
had 50% more vessel’s surface exposed to high WSS, out of which, on average, 40% of
12
vessel surface was exposed to thrombotic WSS (Fig. 3 C).
13
In both models, WSS in ECA increased with increasing grade of restenosis. In ECA of eCEA
14
models, relatively constant proportion of vessel surface was exposed to high WSS, but with
15
progressing restenosis, proportion of vessels surface exposed to thrombotic WSS increased
16
(Fig. 3 B and D). In ECA of exeCEA models, vessel surface exposed to high WSS increased
17
with progressing restenosis (Fig. 3 B and D). Thrombotic WSS (values of 7 Pa or higher)
18
were not detected in ICA nor ECA of exeCEA models (Fig. 2).
19
TE D
M AN U
10
4. DISCUSSION
21
4.1. Surgical techniques and restenosis
22
Leading invasive stroke prevention strategy includes surgical intervention, recommended in
23
asymptomatic and symptomatic patients and in later ‘should be performed sooner rather than
24
later […] within two weeks’ [19]. CEA is a standard surgical approach for carotid
25
revascularization, but alternative to CEA for symptomatic patients is carotid artery stenting,
26
which has advantages for some patients [11, 20]. Conventional CEA includes longitudinal
27
incision from CCA into ICA, while the incision in eCEA is made from carotid bifurcation
28
obliquely, parallel to ECA, downward to CCA lateral border. The later is performed more
29
often and commonly includes larger, oblique transection of CCA. Both techniques have
30
advantages and disadvantages [1-5]. Advantages of eCEA include omitting incision of ICA
AC C
EP
20
6
ACCEPTED MANUSCRIPT (minimizing risk of restenosis of ICA), while disadvantage is restenosis and need for an
2
additional incision of CCA, where atheroma spreads deep in CCA (Fig. 1) [21].
3
Several modifications of eCEA has been described previously; including longitudinal
4
arteriotomy that begins on CCA and skews towards ECA [22], and anterior incision of carotid
5
bulbus [23]. In presented exeCEA approach, the cut is limited to CCA, enables better
6
visualization of proximal end of endarterectomy and length of incision can be adjusted to
7
remove atheroma from proximal parts of CCA. At carotid bifurcation, a typical location of
8
carotid atheroma, the vessel incision line in exeCEA is 1.6 times larger (Supplemental Figure
9
1). As the larger incision means technical ease, vessel injury is also larger - where intimal
RI PT
1
hyperplasia (considered as the most common cause of restenosis) is expected to develop [9].
11
Although incision is larger in exeCEA, the intimal ridge occupies relatively smaller
12
proportion of the orifice of a crucial vessel, the ICA (Fig. 1). Oblique tansection of CCA is
13
often extended during eCEA, in order to create a larger orifice towards ICA, but the orifice is
14
still smaller than in exeCEA (calculation not shown), implying that a larger volume flow
15
towards ICA should be present after exeCEA. Other contributing factors for development of a
16
restenosis is an intimal flap in suture line, continuing atherosclerosis, increased turbulence
17
[24] and change in vessel wall properties [25].
18
4.2. Hemodynamics after carotid endarterectomy
19
As intimal hyperplasia at the site of vessel injury is expected to develop after every CEA, we
20
modeled lumen narrowing at the suture line after original and extended eCEA. Although
21
carotid revascularization is performed on stenotic arteries (arterial lining with irregular
22
atherosclerotic plaques), our simplified model is based on ‘clean’ vessels, since one would
23
expect that after the procedure the arteries would be ‘clean’. As restenosis at suture develops,
24
the lumen of ICA narrows, causing the reduction of volume flow in ICA and shunt of blood
25
into ECA (increase in volume flow in ECA). Beside reduced blood delivery into the brain,
26
altered hemodynamics in narrowed vessel further jeopardizes the brain perfusion, as thrombi
27
that might build-up on the atherosclerotic plaque can cause ipsilateral thromboembolic stroke.
28
Recanalization of ICA leads to the reduction of a peak blood flow velocity and pressure
29
gradients, but effects on WSS and vascular wall properties are less understood [26-28].
30
Intraluminal vascular wall is subjected to physical forces that affect its physiological
31
response, development of atherosclerosis and other wall pathologies [29]. Hemodynamic
32
forces that act on arteries can be perpendicular (i.e., blood pressure) or parallel to the wall
AC C
EP
TE D
M AN U
SC
10
7
ACCEPTED MANUSCRIPT (i.e., WSS) [30]. Both low and high WSS are considered pathological [10, 30, 31]. A low
2
WSS and its high oscillations are linked to the development of local atherosclerosis, through
3
remodeling of the vessel wall [10]. WSS in large arteries ranges from 1 to 10 Pa, reaching
4
peak values during increased cardiac output or hypertension [10, 30]. High WSS leads to
5
development of vortices that entrap non-activated platelets and cause platelet activation and
6
aggregation [31]. Therefore, progressing stenosis, through increased WSS leads to platelet
7
activation and thrombus formation [32]. Lumen narrowing in eCEA models resulted in higher
8
blood velocities causing higher WSS (reaching the thrombotic range; above 7 Pa) (Fig. 2) [10,
9
33]. Accumulated WSS over time is other possible biomechanical parameter of thrombus
RI PT
1
formation [34] and CFD indicate that exeCEA should have less chance of thrombus
11
formation. As stated previously, due to the selected CFD parameters absolute values of WSS
12
might not be similar to experimentally observed values but their relative (eCEA vs. exeCEA)
13
values can be taken into account. Our results show that eCEA models had higher WSS than
14
exeCEA models, indicating that after exeCEA approach thrombotic events might be less
15
frequent, but recurrent atherosclerotic process might be more progressive (as in vessels of
16
exeCEA models part of surface had low WSS). These assumptions remain to be tested in
17
comparative clinical study.
18
The major advantage of exeCEA should be a wider lumen (cross-section area) at the site of
19
restenosis, allowing better flow and pressure pulsation transmission [35]. Due to surgical
20
procedure and subsequent restenosis, altered arterial geometry results with altered flow
21
pulsations [36]. In physiological conditions, 70% of blood from CCA goes into ICA, but as
22
ICA narrows, more blood is shunted towards ECA. ECA blood flow is usually overlooked
23
and considered clinically irrelevant, although increased volume blood flow through ECA
24
means less blood volume in ICA. Volume flow rate in ECA increased with increasing
25
restenosis in both models, but ‘shunt’ was more pronounced in eCEA models (Table 2). WSS
26
in ECA also increased with increasing restenosis, more in eCEA models, where increasingly
27
larger ECA vessel surface was exposed to thrombotic WSS (B and D in Fig. 3). In the light of
28
previous research [37], hemodynamics of ECA indicate that its occlusion after exeCEA might
29
be less frequent. Overall, the hemodynamics of ICA and ECA changed less with progression
30
of restenosis in exeCEA models. Therefore, presented modification of surgical approach
31
might have advantages over the original eCEA.
32
4.3. Study limitations
AC C
EP
TE D
M AN U
SC
10
8
ACCEPTED MANUSCRIPT In order to assess the relative difference in volume flow distribution, velocity and WSS after
2
eCEA and exeCEA, constant pressure was selected as boundary condition (therefore not
3
taking in count peripheral resistance). Such simplification can deviate the results when model
4
includes pulsatile flow and multiple outlets, but in a model of a relatively simple geometry
5
(like the one with single inlet and two outlets), the used boundary condition might be
6
applicable for CFD analysis [15]. We aimed to compare two approaches and, thus, selected
7
CFD parameters were uniformly applied across eCEA and exeCEA models. In this way, the
8
only differing variable was the shape and position of the incision line with its restenosis fold.
9
Therefore, we find that such ‘relative’ experimental setting mitigated our study’s major
10
limitation. CFD studies have become valuable tool in hemodynamic analyses, allowing
11
interpretation of clinical observations and experimental planning based on CFD.
12
4.4. Conclusions
13
Presented CFD models of modification of standard eversive carotid endarterectomy approach
14
showed better hemodynamic parameters with progressing restenosis (lower wall shear stress,
15
higher volume flow rate and lower blood velocities in ICA). These results indicate that, when
16
compared to eCEA, exeCEA might be less thrombotic prone and better at long-term
17
preservation of brain perfusion.
SC
M AN U
TE D EP AC C
18
RI PT
1
9
ACCEPTED MANUSCRIPT 1
Acknowledgements
2
None.
3
Disclosure
5
The authors report no proprietary or commercial interest in any product mentioned or concept
6
discussed in this article.
RI PT
4
AC C
EP
TE D
M AN U
SC
7
10
ACCEPTED MANUSCRIPT
2 3 4 5
REFERENCES 1) Brønnum-Hansen, H., M. Davidsen, and P. Thorvaldsen, Long-term survival and causes of death after stroke. Stroke, 2001. 32(9): p. 2131-2136. 2) Shah, D.M., et al., Carotid endarterectomy by eversion technique: its safety and
RI PT
1
durability. Ann Surg. 1998. 228(4): p. 471.
3) Cao, P., et al., A randomized study on eversion versus standard carotid
7
endarterectomy: study design and preliminary results: the Everest Trial. J Vasc Surg.
8
1998. 27(4): p. 595-605.
10
4) Cao, P., et al., Eversion versus conventional carotid endarterectomy: late results of a
M AN U
9
SC
6
prospective multicenter randomized trial. J Vasc Surg. 2000. 31(1): p. 19-30.
11
5) Cao, P., P. De Rango, and S. Zannetti, Eversion vs conventional carotid
12
endarterectomy: a systematic review. Eur J Vasc Endovasc Surg. 2002. 23(3): p. 195-
13
201.
16
TE D
15
6) Crawford, R.S., et al., Restenosis after eversion vs patch closure carotid endarterectomy. J Vasc Surg. 2007. 46(1): p. 41-48. 7) Mohler, E. and R.M. Fairman, Carotid endarterectomy. [UpToDate - Wolters Kluwer
EP
14
Health].
June
2014.
Available
at:
18
endarterectomy. Accessed June 30, 2014
http://www.uptodate.com/contents/carotid-
AC C
17
19
8) Baan, J., et al., Vessel wall and flow characteristics after carotid endarterectomy:
20
eversion endarterectomy compared with Dacron patch plasty. Eur J Vasc Endovasc
21 22 23 24 25
Surg. 1997. 13(6): p. 583-591.
9) Makihara, N., et al., Characteristic sonographic findings of early restenosis after carotid endarterectomy. Ultrasound Med. 2008. 27(9): p. 1345-1352. 10) Malek, A.M., S.L. Alper, and S. Izumo, Hemodynamic shear stress and its role in atherosclerosis. Jama, 1999. 282(21): p. 2035-2042. 11
ACCEPTED MANUSCRIPT 1 2 3 4
11) Harrison, G.J., et al., Closure technique after carotid endarterectomy influences local hemodynamics. J Vasc Surg. 2014. 60(2): p. 418-427. 12) Fedorov, A., et al., 3D Slicer as an image computing platform for the Quantitative Imaging Network. Magn Reson Imaging. 2012. 30(9): p. 1323-1341. 13) Wake, A.K., et al., Choice of in vivo versus idealized velocity boundary conditions
6
influences physiologically relevant flow patterns in a subject-specific simulation of
7
flow in the human carotid bifurcation. J Biomech Eng. 2009. 131(2): p. 021013.
RI PT
5
14) Morbiducci, U., et al., Outflow conditions for image-based hemodynamic models of
9
the carotid bifurcation: implications for indicators of abnormal flow. J Biomech Eng.
12 13 14
M AN U
11
2010. 132(9): p. 091005.
15) Grinberg, L. and G.E. Karniadakis, Outflow boundary conditions for arterial networks with multiple outlets. Ann Biomed Eng. 2008. 36(9): p. 1496-1514. 16) Schneider, C.A., W.S. Rasband, and K.W. Eliceiri, NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012. 9(7): p. 671-675.
TE D
10
SC
8
17) Augst, A., et al., Analysis of complex flow and the relationship between blood
16
pressure, wall shear stress, and intima-media thickness in the human carotid artery.
17
Am J Physiol Heart Circ Physiol. 2007. 293(2): p. H1031-H1037.
EP
15
18) Kamenskiy, A.V., et al., A mathematical evaluation of hemodynamic parameters after
19
carotid eversion and conventional patch angioplasty. Am J Physiol Heart Circ
20
AC C
18
Physiol. 2013. 305(5): p. H716-H724.
21
19) Howell, S., Carotid endarterectomy. Br J Anaesth. 2007. 99(1): p. 119-131.
22
20) Almekhlafi, M., et al., When Is Carotid Angioplasty and Stenting the Cost-Effective
23
Alternative for Revascularization of Symptomatic Carotid Stenosis? A Canadian
24
Health System Perspective. AJNR Am J Neuroradiol. 2014. 35(2): p. 327-332.
12
ACCEPTED MANUSCRIPT 1 2
21) Findlay, J.M., et al., Carotid endarterectomy: a review. Can J Neurol Sci. 2004. 31(01): p. 22-36. 22) Okur, F.F., et al., Experience of 352 Carotid Endarterectomy Cases Operated on
4
without Any Shunt and a Technical Modification: Thirty-Day Results. Turkiye
5
Klinikleri J Med Sci. 2013. 33(5): p. 1216-1223.
8 9 10 11
27(2): p. 178-185.
24) Fietsam, R., et al., Hemodynamic effects of primary closure versus patch angioplasty
SC
7
23) Kumar, S., et al., Modified eversion carotid endarterectomy. Ann Vasc Surg. 2013.
of the carotid artery. Ann Vasc Surg. 1992. 6(5): p. 443-449.
25) Friedman, M.H., Some atherosclerosis may be a consequence of the normal adaptive
M AN U
6
RI PT
3
vascular response to shear. Atherosclerosis 1990. 82(3): p. 193-196. 26) Harloff, A., et al., Comparison of blood flow velocity quantification by 4D flow MR
13
imaging with ultrasound at the carotid bifurcation. AJNR Am J Neuroradiol. 2013.
14
34(7): p. 1407-1413.
TE D
12
27) Sachar, R., et al., Severe bilateral carotid stenosis: the impact of ipsilateral stenting on
16
Doppler-defined contralateral stenosis. J Am Coll Cardiol. 2004. 43(8): p. 1358-1362.
17
28) Aleksic, M., et al., Changes in internal carotid blood flow after CEA evaluated by
19 20 21 22
transit-time flowmeter. Eur J Vasc Endovasc Surg. 2006. 31(1): p. 14-17. 29) Resnick, N., et al., Fluid shear stress and the vascular endothelium: for better and for
AC C
18
EP
15
worse. Prog Biophys Mol Biol. 2003. 81(3): p. 177-199.
30) Davies, P.F., Flow-mediated endothelial mechanotransduction. Physiol Rev. 1995. 75(3): p. 519-560.
23
31) Biasetti, J., F. Hussain, and T.C. Gasser, Blood flow and coherent vortices in the
24
normal and aneurysmatic aortas: a fluid dynamical approach to intra-luminal
25
thrombus formation. J. R. Soc. Interface. 2011: p. rsif20110041.
13
ACCEPTED MANUSCRIPT 1
32) Holme, P.A., et al., Shear-induced platelet activation and platelet microparticle
2
formation at blood flow conditions as in arteries with a severe stenosis. Arterioscler
3
Thromb Vasc Biol. 1997. 17(4): p. 646-653.
6 7 8 9
disturbed laminar blood flow. Thromb Haemost. 1996. 75(5): p. 827-832.
RI PT
5
33) Barstad, R., et al., Reduced effect of aspirin on thrombus formation at high shear and
34) Hellums, J.D., 1993 Whitaker Lecture: biorheology in thrombosis research. Ann Biomed Eng. 1994. 22(5): p. 445-455.
35) Young, D., Fluid mechanics of arterial stenoses. J Biomech Eng. 1979. 101(3): p.
SC
4
157-175.
36) Farrar, D.J., H.D. Green, and D.W. Peterson, Noninvasively and invasively measured
11
pulsatile haemodynamics with graded arterial stenosis. Cardiovasc Res. 1979. 13(1):
12
p. 45-57.
TE D
15
artery following carotid endarterectomy? BMC Surg. 2008. 8(1): p. 20.
EP
14
37) Abbas, S.M., D. Adams, and P. Vanniasingham, What happens to the external carotid
AC C
13
M AN U
10
14
ACCEPTED MANUSCRIPT Table 1. Internal carotid artery blood velocities and volume flow rates in eCEA and exeCEA models. BLOOD VELOCITY Restenosis
VOLUME FLOW RATE p
eCEA [AUC]
exeCEA [AUC]
p
Grade I
0.332
0.304
<0.001
5.28
5.43
<0.001
Grade II
0.327
0.303
<0.001
5.13
5.37
<0.001
Grade III
0.281
0.282
0.734
4.71
5.05
<0.001
Grade IV
0.220
0.234
0.003
3.61
3.93
0.005
Grade I
0.415
0.385
<0.001
6.59
6.72
<0.001
Grade II
0.404
0.385
<0.001
6.29
6.66
<0.001
Grade III
0.343
0.331
0.034
5.93
6.01
0.373
Grade IV
0.267
0.287
0.001
4.49
4.86
0.003
120 mmHg
M AN U
150 mmHg
RI PT
eCEA [AUC]
exeCEA [AUC]
SC
Input pressure
AUC - area under curve.
Table 2. External carotid artery blood velocities and volume flow rates in eCEA and exeCEA models. BLOOD VELOCITY
Restenosis
TE D
Systolic pressure
VOLUME FLOW RATE
eCEA [AUC]
exeCEA [AUC]
p
eCEA [AUC]
exeCEA [AUC]
p
Grade I
0.318
0.271
<0.001
3.21
3.13
0.002
Grade II
0.328
0.275
<0.001
3.33
3.20
<0.001
Grade III
0.364
0.307
<0.001
3.71
3.59
<0.001
Grade IV
0.455
0.377
<0.001
4.82
4.54
0.014
AC C
EP
120 mmHg
Grade I
0.388
0.336
<0.001
4.02
3.98
0.048
Grade II
0.406
0.340
<0.001
4.19
4.05
<0.001
Grade III
0.448
0.403
<0.001
4.63
4.75
0.131
Grade IV
0.562
0.472
<0.001
6.11
5.79
0.018
150 mmHg
AUC - area under curve.
15
ACCEPTED MANUSCRIPT FIGURE TITLES AND CAPTIONS
Figure 1. Incision line in eversion and modified eversion carotid endarterectomy and
RI PT
their cross–sectional areas. Scheme of carotid vessels with incision line after eversion (eCEA) and modified eversion carotid endarterectomy (exeCEA). exeCEA involves transection of CCA, from bifurcation
SC
distally, until reaching proximal end of plaque, where incision ends by extending under angle
M AN U
of 45° towards ICA.
Figure 2. Wall shear stress in modified and eversion carotid endarterectomy. Results of computation of wall shear stress in carotid arteries after eversion (eCEA) and modified eversion carotid endarterectomy (exeCEA) during single heart cycle under pressure
restenosis grade I-IV.
TE D
of 120 and 150 mmHg. Simulation included increasing groove at ‘incision line’ - representing
EP
Figure 3. Relative WSS in ICA and ECA after eCEA and exeCEA procedure. Since the relative values were of interest, WSS values in each chart are normalized to eCEA
AC C
grade I. Line representing WSS higher than 7 Pa in exeCEA is not present since thrombotic range with this procedure is never reached. Supplemental Figure 1. Representation of incision line in eversion and modified eversion carotid endarterectomy and calculation of their cross–sectional areas.
16
AC C
EP
TE D
M AN U
SC
RI PT
ACCEPTED MANUSCRIPT
AC C
EP
TE D
M AN U
SC
RI PT
ACCEPTED MANUSCRIPT
AC C
EP
TE D
M AN U
SC
RI PT
ACCEPTED MANUSCRIPT
ACCEPTED MANUSCRIPT SUPPLEMENT METHODS Approximation of cross-sectional areas in eCEA and exeCEA Several assumptions were made in order to simplify calculation of cross-sectional areas. A
AC C
EP
TE D
M AN U
SC
shape composed of rectangle and half of the ellipse.
RI PT
cross-sectional incision areas; eCEA has a elliptical shape, while the exeCEA has a complex
dCCA=1 2beCEA=2.05xdCCA 2aeCEA=dCCA
ACCEPTED MANUSCRIPT The cross-sectional area of ellipse in eCEA is: PeCEA=axbxπ PeCEA= (1.02xdCCA)x(0.5xdCCA)xπ
The cross-sectional area of ellipse in exeCEA is:
PexeCEA=Prectangle+Phalf-ellipse
l1=2.02xdCCA
M AN U
l2= dCCA
SC
P=(a x b)+(a x b x π)/2
RI PT
PeCEA=1.61xdCCA2
Prectangle= l1 x l2
Prectangle = 2.02dCCA2
2bexeCEA = dCCA/cos 45°
TE D
2bexeCEA= 1.42x dCCA bexeCEA= 0.71 x dCCA
PexeCEA=Prectangle+Phalf-ellipse
EP
PexeCEA=(a x b)+ (a x b x π)/2
AC C
PexeCEA=2.02dCCA2+(0.5xdCCAx0.71xdCCAx π)/2 PexeCEA=2.57xdCCA2
Ratio of exeCEA to eCEA cross-sectional areas:
PexeCEA:PeCEA=2.57x dCCA2 : 1.61xdCCA2 PexeCEA:PeCEA=1.6