Extended Eversion Carotid Endarterectomy: Computation of Hemodynamics

Extended Eversion Carotid Endarterectomy: Computation of Hemodynamics

Accepted Manuscript Extended Eversion Carotid Endarterectomy: Computation of Hemodynamics Tomislav Istvanić, Zvonimir Vrselja, Hrvoje Brkić, Radivoje ...

861KB Sizes 9 Downloads 285 Views

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