ORIGINAL RESEARCH article

Front. Med. Technol., 03 August 2021 | https://doi.org/10.3389/fmedt.2021.719909

Impact of the Internal Carotid Artery Morphology on in silico Stent-Retriever Thrombectomy Outcome

  • 1Laboratory of Biological Structure Mechanics, Department of Chemistry, Materials and Chemical Engineering, Politecnico di Milano, Milan, Italy
  • 2School of Engineering, National University of Ireland Galway, Galway, Ireland
  • 3Department of Biomedical Engineering and Physics, Amsterdam University Medical Center, Location Academic Medical Center, Amsterdam, Netherlands
  • 4Department of Radiology and Nuclear Medicine, Amsterdam University Medical Center, Amsterdam, Netherlands

The aim of this work is to propose a methodology for identifying relationships between morphological features of the cerebral vasculature and the outcome of in silico simulations of thrombectomy, the mechanical treatment for acute ischemic stroke. Fourteen patient-specific cerebral vasculature segmentations were collected and used for geometric characterization of the intracranial arteries mostly affected by large vessel occlusions, i.e., internal carotid artery (ICA), middle cerebral artery (MCA) and anterior cerebral artery (ACA). First, a set of global parameters was created, including the geometrical information commonly provided in the clinical context, namely the total length, the average diameter and the tortuosity (length over head-tail distance) of the intracranial ICA. Then, a more exhaustive geometrical analysis was performed to collect a set of local parameters. A total of 27 parameters was measured from each patient-specific vascular configuration. Fourteen virtual thrombectomy simulations were performed with a blood clot with the same length and composition placed in the middle of the MCA. The model of TREVO ProVue stent-retriever was used for all the simulations. Results from simulations produced five unsuccessful outcomes, i.e., the clot was not removed from the vessels. The geometric parameters of the successful and unsuccessful simulations were compared to find relations between the vascular geometry and the outcome. None of the global parameters alone or combined proved able to discriminate between positive and negative outcome, while a combination of local parameters allowed to correctly identify the successful from the unsuccessful simulations. Although these results are limited by the number of patients considered, this study indicates a promising methodology to relate patient-specific geometry to virtual thrombectomy outcome, which might eventually guide decision making in the treatment of acute ischemic stroke.

Introduction

The intra-arterial thrombectomy is a minimally invasive procedure based on stent-retriever technology for acute ischemic stroke (AIS) patients. An AIS arises when a large vessel occlusion caused by a thromboembolus (clot) prevents blood supply to the brain tissues. Despite the promising and encouraging results, complications may occur, such as thrombus fragmentation, vascular damage, or microembolization (1, 2).

The most likely location of the AIS is the intracranial internal carotid artery (ICA) and its bifurcated downstream segments. The main individual variability in ICA geometry can be found in the so-called carotid siphon (3). The carotid siphon lies between the carotid canal and the T-junction, where the ICA bifurcates into the middle cerebral artery (MCA) and the anterior cerebral artery (ACA) (Figure 1A). Four different bends can be recognized in the ICA (4, 5): (from the carotid canal) inferior, posterior, anterior, and superior (Figure 1B).

FIGURE 1
www.frontiersin.org

Figure 1. (A) Cerebral arteries mostly affected by AIS: the intracranial internal carotid artery (ICA) which bifurcates at the T-junction into the middle cerebral artery (MCA) and the anterior cerebral artery (ACA). (B) The carotid siphon with the superior, anterior, posterior, and inferior bends. (C) Classes of siphon shape.

Although the way the ICA geometry affects the hemodynamics and different pathologies onset (aneurysms and stenosis) has been already investigated (610), the correlation between the ICA geometry and thrombectomy procedure outcome remains under-studied. Mokin et al. analyzed the correlation between the ICA geometric parameters and the outcome of the thrombectomy procedure, finding that the intracranial tortuosity did not influence the needed number of passes with thrombectomy devices (11). Srivatsa et al. investigated the influence of vessel morphology on the fist-pass success of mechanical thrombectomy in patients with MCA occlusions (12). Despite the limitation of performing manual measurements on 2D angiograms, they were able to conclude that a larger MCA diameter associated with a lower rate of tapering along the ICA was associated with high probability of first-pass recanalization. Schwaiger et al. observed an influence of MCA curvature on the success of recanalization performed with stent-retrievers (13). The curvature was expressed as angle between the ICA and the proximal part of the MCA and angle between proximal and distal MCA segments. Results showed that patients with unsuccessful recanalization had significantly higher angles. However, this study did not differentiate the patients based on other features like the thrombus length and composition or the type of stent-retriever used in the procedure, hence the causality between MCA curvature and successful recanalization could not be clearly assessed. Moreover, the tortuous geometry of the ICA was not considered in the study.

In the clinical application, the ICA is usually characterized by its average diameter, length, tortuosity (artery length over head-tail distance), or by other mean values (14). A classification based on the siphon shape is also used to evaluate the artery accessibility (15). In particular, the combined shape of the anterior and posterior bends defines four classes of siphon shape (U, C, V and S) (Figure 1C). Differently, two studies proposed an objective and extensive geometric characterization of carotid siphons, parametrized by a set of local geometric parameters (4, 7). Zhang et al. proposed to characterize the carotid siphon by identifying the two main bends and measuring their curvature, the lumen diameters at the end-points of each bend and the angle and distance between them. Bogunović et al. further enriched the carotid siphon analysis by dividing the ICA in four bends and measuring for each of them the curvature, the vessel diameter and the angles between adjacent bends.

From a technological perspective, the regulatory process for biomedical devices recently started receiving and accepting in silico evidence from modeling and numerical simulations (16). Recently our group developed and validated a finite-element model of the stent-retriever thrombectomy (17) for in silico clinical trials in AIS (18); this methodology was subsequently applied to make a replica of the procedure in a patient-specific case (19).

In this work, we hypothesize that local geometric features of the ICA are related to the virtual stent-retriever thrombectomy outcome. In particular, (i) a geometric characterization of 14 patient-specific cerebral vasculatures is carried out with 4 global and 27 local parameters for each patient; (ii) the thrombectomy procedure is simulated for each patient with the same settings (thrombus characteristics and stent-retriever); (iii) a preliminary post-processing correlation between the global and local parameters and the virtual thrombectomy outcome is proposed as methodology for future investigations.

Methods

Geometric Characterization

We selected patients from the MR CLEAN (Multicenter Randomized Clinical Trial of Endovascular Treatment for Acute Ischemic Stroke in the Netherlands) Registry, which is an ongoing, prospective, observational, multicenter study, including data of stroke patients treated in 16 intervention hospitals in the Netherlands (20). The study population consisted of 14 patients who presented with an AIS due to a large vessel occlusion in the anterior circulation. The patients underwent a non-contrast computed tomography (NCCT) and a computed tomography angiography (CTA) scan before treatment. The segmentations of intracranial vessels were performed with StrokeViewer (NICO.LAB, Amsterdam, The Netherlands). iCAFE (14) (© 2016-2018 University of Washington. Used with permission.), a semi-automated software, was used to extract centerlines, local radius, and to label arterial segments. For a more detailed description of the vascular segmentation and centerline computation, refer to (21). The centerlines of the arteries of interest, i.e. ICA, MCA, and ACA, were extracted (Figure 2) and analyzed to obtain the values of the geometric features describing each patient-specific vessel geometry.

FIGURE 2
www.frontiersin.org

Figure 2. The reconstructed vessel surfaces and centerlines of the 14 studied patients: the vascular segments of interest are ICA (green), MCA (red) and ACA (blue).

First, the geometric parameters routinely considered in the clinical context were measured for each patient geometry. These parameters are the total length of the intracranial ICA just distal form the carotid canal to the T-junction (LICA), the average diameter of the ICA (DICA) and the ICA tortuosity (torICA), calculated as (22):

torICA=LICAdICA-1

where dICA is the Euclidean distance of the ICA endpoints. Additionally, each ICA siphon was visually classified into U, C, V, or S shapes. In this work, this set of parameters is labeled as Global parameters since they are referred to the entire ICA geometry.

Successively, the patient-specific vessel geometries were analyzed to extract the set of Local parameters. This process was made fully automatic by means of a MATLAB (R2020a, The MathWorks, Natick, MA, USA) script. First, the angles created at the T-junction were measured (αICAMCA, αICAACA, αACAMCA). Each angle was measured by fitting the two segments of the vessels converging into the T-junction with a linear segment and by measuring the angle formed in the plane containing both the linear segments (Figure 3). The ICA shape was then analyzed. Following the method proposed in (4), each ICA centerline was divided into four segments: superior, anterior, posterior, and inferior bends, progressively located from the T-junction to the end closest to the heart (Figure 4A). The landmarks dividing adjacent bends were identified with a combined use of the Frenet-Serret triad and a parallel transport frame following the ICA curve (4). The four bends and couples of adjacent bends were analyzed to obtain the following geometric parameters fully describing the ICA geometry:

• bends lengths: Lsup, Lant, Lpos, Linf;

• bends diameters: Dsup, Dant, Dpos, Dinf, measured at the point of maximum bend curvature;

• curvature radii: rsup, rant, rpos, rinf, calculated as the radius of the circle best fitting the bend (Figure 4B);

• bends tortuosity: torsup, torant, torpos, torinf, calculated with the formula:

torbend=Lbenddbend-1

• where Lbend is the length of the bend and dbend is the Euclidean distance of the bend endpoints;

• angles between adjacent bends: α supant, α antpos, α posinf, measured between the planes containing the fitting circle of each bend (Figure 4C);

• distance between the bifurcation point and the starting point of the superior bend (considered as the more distal intersection point between the superior bend and its fitting circle): dTsup (Figure 4D);

• distances between adjacent bends: dsupant, dantpos, dposinf, calculated as the distance between the ending point of the more distal bend and the starting point of the more proximal one, where the starting and ending points of each bend are considered as the distal and proximal intersection point of the bend and the relative fitting circle, respectively (Figure 4D).

FIGURE 3
www.frontiersin.org

Figure 3. Measurement of the angles at the T-junction: from left to right, αICA−MCA, αICA−ACA, αACA−MCA.

FIGURE 4
www.frontiersin.org

Figure 4. ICA analysis: (A) division of the ICA in the 4 bends (superior, anterior, posterior, inferior) separated by the landmarks; (B) calculation of curvature of each bend as the radius of the fitting circle; (C) angles between bends; (D) distance between bends calculated at adjacent intersection points between the bend and the relative fitting circle; for the superior bend, the distance is calculated between the distal intersection point and the ICA endpoint.

Finally, the average diameter of the MCA (DMCA) was calculated. The Local parameters set thus consists of 27 parameters for each patient-specific vessel geometry.

Virtual Thrombectomy Procedure

To study the impact of the geometric features of the intracranial ICA on the virtual thrombectomy procedures, the same thrombus, stent-retriever, and thrombectomy procedure were implemented in the 14 patient-specific vessel geometries.

The patient-specific domains were obtained from the computed centerlines with the aid of level-set images (21). In particular, triangular shell elements with a mean element size of 0.2 mm were used to discretized the vessel walls using the open-source MATLAB toolbox GIBBON (23). Rigid wall material was assumed.

Thrombi with the same length and composition (in terms of red blood cells (RBC) and fibrin contents) were placed in the central part of the MCA (before its bifurcation) for each patient (Figure 5), as this is the most common AIS location (24). The thrombus characteristics were chosen as mean values obtained from the literature. In particular with a length of 14 mm (25) and 35% RBC - 65% fibrin composition (26). The thrombus diameter was set to occlude 90% of the patient-specific MCA diameter where it is located. Thrombi were discretized using ANSA Pre Processor (BETA CAE System, Switzerland) with linear tetrahedral elements with a mean element size of 0.2 mm on the surface. The thrombus non-linear material was modeled with a quasi-hyperelastic foam formulation by fitting stress-strain curves from in vitro tests on ex vivo clots (19).

FIGURE 5
www.frontiersin.org

Figure 5. A clot of 14mm of length placed in the central MCA segment of the 14 patients.

The geometry of the TREVO ProVue 4–20mm stent (Stryker, US) was obtained using a Python code and discretized with linear beam elements with a mean element length of 0.2 mm (Figure 6A). Hughes-Liu formulation with cross-section integration was used. A shape memory alloy material model was used to simulate the NiTi superelastic behavior. Details about the NiTi material parameter definitions can be found in (17).

FIGURE 6
www.frontiersin.org

Figure 6. (A) Finite-element model of the TREVO ProVue 4–20mm stent-retriever. (B) The main steps of the thrombectomy simulation.

The virtual thrombectomy procedure was simulated in the 14 patients with the same settings (Figure 6B): the stent-retriever is crimped in a 0.5mm diameter catheter and tracked to the clot location where it is deployed; finally, it is retrieved to extract the trapped clot. For a more detailed description of the contact algorithms, time step, damping system, and other technical details refer to (17), where a validation of the thrombectomy simulation with in vitro experiments was carried out. In Luraghi et al. (19) the same simulation settings were used to successfully reproduce a clinical case. The simulations were performed on a system featuring 40 CPUs (Intel Xeon64) and 256 GB of RAM, using the commercial finite-element solver LS-DYNA (ANSYS, USA) and lasted about 24 h.

Post-processing Analysis

The virtual thrombectomy procedure was defined to have a positive outcome if the clot reaches the end of the ICA segment with the stent, and a negative outcome if the clot remains inside the vessel. The outcomes of the 14 simulations were investigated and correlated with the corresponding ICA geometric parameters and their combinations, with both the global and local sets of parameters.

Results

ICA Geometric Features

The measured mean and local sets of geometric parameters are listed in Table 1. A few patients missed some of the bends to perform the local geometric analysis because the most proximal part of the ICA was not included or not discernible in the CTA images: the models of patients 3, 5, 9, and 10 do not include the inferior bend, patient 7 does not have the posterior and inferior bends, patients 6 and 13 do not have the anterior, posterior and inferior bends. In Figure 7, the box plots of the local and global geometric parameters are shown.

TABLE 1
www.frontiersin.org

Table 1. Global and Local parameters of the 14 patient-specific geometries.

FIGURE 7
www.frontiersin.org

Figure 7. Box plots representing the distribution of the local and global geometric parameters.

Virtual Thrombectomy Outcome

The thrombectomy simulations presented with following results: nine had a positive outcome and five a negative outcome. In the successful thrombectomy simulations (patients 1, 2, 3, 6, 7, 9, 11, 13, and 14) the clots remained captured in the stent-retriever during the retrieval phase and reached the end of the ICA segment (Figure 8-top). On the contrary, in the failed thrombectomies (patients 4, 5, 8, 10, and 12), the clots lost contact with the stent-retriever and remained in the anterior bend of the ICA (Figure 8-bottom).

FIGURE 8
www.frontiersin.org

Figure 8. Examples of virtual thrombectomies with positive (top) and negative (bottom) outcome.

Correlation Between ICA Geometry and Thrombectomy Outcome

The outcomes of the thrombectomy simulations were analyzed with the associated global and local geometric parameters to find indicators able to determine the positive or negative outcome. Associations with the thrombectomy outcomes were analyzed for each parameter singularly and for combinations of parameters within each group (global or local). Only the significant associations are here presented.

Using the global geometric parameters, the index with the best performance in discriminating positive and negative cases is the ratio between the clot diameter (90% of MCA diameter) and the average ICA diameter, Dclot/DICA (Index 1 in Table 2), with a threshold of 0.67 (calculated as the mean value between the two patients closest to the boundary, i.e., patients 9 and 12) able to correctly classify 13 cases out of 14 (Figure 9-left). However, besides the misclassified case (patient 4), many correctly classified samples are very close to the separating boundary, indicating that the threshold is not robust in discriminating the cases.

TABLE 2
www.frontiersin.org

Table 2. Indices used for correlating the thrombectomy simulation outcome with global and local geometric parameters of the patient (+ positive outcome, –negative outcome).

FIGURE 9
www.frontiersin.org

Figure 9. Best classification indices using the global and local geometric parameters (on the right, only 12 samples are shown because the segmentation for patients 6 and 13 did not include the anterior bend).

Considering the local geometric parameters, a different index looked more promising to our aim. The combination of the ratio between the clot diameter and the diameter of the anterior bend in the point of maximum curvature, Dclot/Dant (Index 2 in Table 2), and the curvature radius of the anterior bend, rant, allowed for the identification of a separation boundary which correctly classifies the cases, with a wide margin between the two classes (Figure 9-right). For the patients considered in this study, thresholds around 0.55 for Dclot/Dant and 4.5 mm for rant were identified. Namely, when the clot diameter is large with respect to the diameter of the anterior bend in the point of its maximum curvature (Dclot/Dant > 0.55) and the anterior bend has a pronounced curvature (rant < 4.5mm), the stent is not able to remove the clot from the anterior bend, which is the failure mode observed in all the 5 cases with a negative outcome.

Discussion

The proposed geometrical analysis of the ICA is aimed at suggesting an association between geometrical features and in silico thrombectomy outcome.

The first aspect this study has highlighted is the importance of analyzing the local parameters rather than the global ones to discriminate the different morphologies of the ICA. For example, while the global parameters (siphon shape and ICA tortuosity) for patients 8 and 9 are the same, the morphology of the ICA bends is clearly different, as shown in Figure 10. Focusing on the anterior bend, the curvature radius appears to be a better index of the bend morphology with respect to the tortuosity index. Again, considering patients 8 and 9 as examples, the difference in the curvature radius (2.7 vs. 5.1 mm) is much more important than the difference in the tortuosity index (1.54 vs. 1.10). It is to be noted that patients 8 and 9 present a different virtual thrombectomy outcome.

FIGURE 10
www.frontiersin.org

Figure 10. Example of patients, 8 and 9, presenting very similar global parameters (green boxes), but different local parameters (red boxes), with a particular focus on the anterior bend (red tract of the ICA centerline), which resulted in different virtual thrombectomy outcomes.

As regards the in silico thrombectomies, the set of local geometric parameters proved to be more efficient in identifying indices able to discriminate positive and negative simulation outcomes with respect to the global parameters. In particular, the anterior bend was identified as the segment of the ICA determining the possible failure of the procedure. A combination of a high ratio between the diameters of the blood clot and of the anterior bend and a small curvature radius of the anterior bend appears to be a decisive feature for negative in silico thrombectomy outcomes. Comparing again patients 8 and 9, from Table 2 it can be observed that since these cases have the same value for the index Dclot/Dant, the outcome cannot be correctly distinguished from this index alone. However, when the radius of the anterior bend curvature in patient 8 is also considered, failure to remove the blood clot appears to be predictable (Figures 9, 10). On the other hand, patients 8 and 11 present with the same value of anterior bend curvature, however, here the difference in the outcome (Figure 9) appears to be due to the relatively large size of the clot with respect to the anterior bend diameter in patient 8. In this study, the indicated threshold values represent only an indication and they should be refined considering more patients with geometric features close to the separating boundaries. Furthermore, the findings presented may be limited to the particular type of thrombus and stent-retriever modeled. Indeed, this study was conducted using a single type of blood clot (65% fibrin content and 14 mm length) and one single stent-retriever design (TREVO ProVue).

A limitation of this work is the small number of patients: future work will gather more patient-specific cerebral vasculature geometries to refine this study. Nonetheless, the 14 patients had substantially different siphon shapes and geometric features, which therefore represented a broad spectrum of vascular geometries. Furthermore, the calculated global and local parameters shown in Figure 7 are in line with what was reported in the literature (4, 7, 24), indicating that the 14 patients of this study are a good representation of the AIS population. The lack of some bends in a subset of patients, due to truncated segmentations of the ICA, is another limitation. In addition, it has not been possible to compare the outcomes of the simulations with the actual outcomes of the clinical procedures. This is because in this study, to focus solely on the impact of the ICA geometry, all simulations were conducted using the same type of retriever device and the same clot position and composition. In reality, however, each patient may have clinically presented with a different type of clot in a different location and may have been treated with a different retriever device. Therefore, the correlations reported here can only be associated with in silico thrombectomy. However, the results presented do shed light on the potential role that the local geometric parameters of the ICA, together with the geometric and mechanical characteristics of the clot and the retriever, may play on predicting the outcome of thrombectomy in clinic. Additionally, this study focused on the ability of the considered stent-retriever to remove the thrombus from the cerebral vessels. However, the recanalization may also be incomplete due to thrombus fragmentation, which can cause embolization. The fragmentation of the clot was not modeled in the simulations. Finally, the vessel was modeled with a rigid material, hence vascular damage due to the procedure could not be assessed. The integration in the simulation of mechanisms of thrombus fragmentation, deformable vessel walls and the introduction of the blood flow to simulate procedures without the use of a balloon guide catheter will be integrated in the future to improve the impact of similar studies.

Future work should investigate the influence of the clot features and the effect of using a different stent-retriever or a different thrombectomy technique on the conducted analyses. Once this analysis is made, it might be used in the clinic in the pre-planning phase to select the best stent-retriever, catheters and procedure for a selected patient. Indeed, once the vessel segmentation process is performed, the implemented automatic script can provide in real-time the patient's set of local parameters needed to assess the most suitable thrombectomy procedure. In this process, the main user-time consuming step would be the image segmentation. The results of this work pave the way to find an association between the virtual thrombectomy outcome and the local set of parameters, rather than the global parameters. This study is not intended to be conclusive, but it can be considered as a first step toward a classification strategy to be used for decision making in the treatment of AIS patients.

Data Availability Statement

The datasets presented in this article are not readily available because this dataset is property of the INSIST consortium. Requests to access the datasets should be directed to https://www.insist-h2020.eu/index.php/contact-us.

Author Contributions

SB and GL designed the study and drafted the manuscript. JR and GD supported the implementation of the study and participated the analysis and the discussion of results. GG and GM carried out the thrombectomy simulations. JK implemented the automatic script to analyze the geometry. KM processed the segmentations and provided the discretized patient-specific geometries with the supervision of PM. PK and NA performed the segmentation process with the supervision of HM. HM, EB, and CM helped in supervising the research project and gave insights for the clinical aspects and revise the manuscript. FM was in charge of overall direction and planning. All authors contributed to the article and approved the submitted version.

Funding

This study has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement No 777072. GL, JM, and FM also received founding from the MIUR FISR-FISR2019_03221 CECOMES. NA also received funding from the AMC medical Research BV, Amsterdam UMC, location AMC, under project No 21937.

Conflict of Interest

HM reports co-founder and shareholder of Nico.lab, a company that focuses on the use of artificial intelligence for medical image analysis. CM received funds from the European Commission (related to this project, paid to institution); and from CVON/Dutch Heart Foundation, Stryker, TWIN Foundation, Health Evaluation Program Netherlands (unrelated; all paid to institution). CM is shareholder of Nico.lab, a company that focuses on the use of artificial intelligence for medical imaging analysis.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Munich SA, Vakharia K, Levy EI. Overview of mechanical thrombectomy techniques. Clin Neurosurg. (2019) 85:S60–S7. doi: 10.1093/neuros/nyz071

CrossRef Full Text | Google Scholar

2. Kühn AL, Vardar Z, Kraitem A, King RM, Anagnostakou V, Puri AS, et al. Biomechanics and hemodynamics of stent-retrievers. J Cereb Blood Flow Metab. (2020) 40:2350–65. doi: 10.1177/0271678X20916002

CrossRef Full Text | Google Scholar

3. Thomas JB, Antiga L, Che SL, Milner JS, Steinman DAH, Spence JD, et al. Variation in the carotid bifurcation geometry of young versus older adults: implications for geometric risk of atherosclerosis. Stroke. (2005) 36:2450–6. doi: 10.1161/01.STR.0000185679.62634.0a

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Bogunović H, Pozo JM, Cárdenes R, Villa-Uriol MC, Blanc R, Piotin M, et al. Automated landmarking and geometric characterization of the carotid siphon. Med Image Anal. (2012) 16:889–903. doi: 10.1016/j.media.2012.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Sanders-Taylor C, Kurbanov A, Cebula H, Leach JL, Zuccarello M, Keller JT. The carotid siphon: a historic radiographic sign, not an anatomic classification. World Neurosurg. (2014) 82:423–7. doi: 10.1016/j.wneu.2013.09.029

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Phan TG, Beare RJ, Jolley D, Das G, Ren M, Wong K, et al. Carotid artery anatomy and geometry as risk factors for carotid atherosclerotic disease. Stroke. (2012) 43:1596–601. doi: 10.1161/STROKEAHA.111.645499

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Zhang C, Pu F, Li S, Xie S, Fan Y, Li D. Geometric classification of the carotid siphon: association between geometry and stenoses. Surg Radiol Anat. (2013) 35:385–94. doi: 10.1007/s00276-012-1042-8

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Spanos K, Petrocheilou G, Karathanos C, Labropoulos N, Mikhailidis D, Giannoukas A. Carotid bifurcation geometry and atherosclerosis. Angiology. (2017) 68:757–64. doi: 10.1177/0003319716678741

CrossRef Full Text | Google Scholar

9. Gregg S, Li TY, Hétu MF, Pang SC, Ewart P, Johri AM. Relationship between carotid artery atherosclerosis and bulb geometry. Int J Cardiovasc Imaging. (2018) 34:1081–90. doi: 10.1007/s10554-018-1319-z

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Ngo MT, Kwak HS, Ho CG, Koh EJ. Longitudinal study of carotid artery bifurcation geometry using magnetic resonance angiography. Vascular. (2019) 27:312–7. doi: 10.1177/1708538119828262

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Mokin M, Waqas M, Chin F, Rai H, Senko J, Sparks A, et al. Semi-automated measurement of vascular tortuosity and its implications for mechanical thrombectomy performance. Neuroradiology. (2020) 63:381–9. doi: 10.1007/s00234-020-02525-6

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Srivatsa S, Duan Y, Sheppard JP, Pahwa S, Pace J, Zhou X, et al. Cerebral vessel anatomy as a predictor of first-pass effect in mechanical thrombectomy for emergent large-vessel occlusion. J Neurosurg. (2021) 134:576–84. doi: 10.3171/2019.11.JNS192673

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Schwaiger BJ, Gersing AS, Zimmer C, Prothmann S. The curved MCA: influence of vessel anatomy on recanalization results of mechanical thrombectomy after acute ischemic stroke. Am J Neuroradiol. (2015) 36:971–6. doi: 10.3174/ajnr.A4222

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Chen L, Mossa-Basha M, Balu N, Canton G, Sun J, Pimentel K, et al. Development of a quantitative intracranial vascular features extraction tool on 3D MRA using semiautomated open-curve active contour vessel tracing. Magn Reson Med. (2018) 79:3229–38. doi: 10.1002/mrm.26961

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Krayenbuehl H, Yasargil M, Huber P. Cerebral Angiography. 2nd ed. New York, NY: Thieme (1982).

Google Scholar

16. The American Society of Mechanical Engineers. Assessing Credibility of Computational Modeling Through Verificationa and Validation: Application to Medical Devices. (2018). Available at: https://www.asme.org/codes-standards/find-codes-standards/v-v-40-assessing-credibility-computational-modeling-verification-validation-application-medical-devices

17. Luraghi G, Rodriguez Matas JF, Dubini G, Berti F, Bridio S, Duffy S, et al. Applicability assessment of a stent-retriever thrombectomy finite-element model. Interface Focus. (2021) 11:20190123. doi: 10.1098/rsfs.2019.0123

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Konduri PR, Marquering HA, van Bavel EE, Hoekstra A, Majoie CB, INSIST Investigators. In-silico trials for treatment of acute ischemic stroke. Front Neurol. (2020) 11:1062. doi: 10.3389/fneur.2020.558125

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Luraghi G, Bridio S, Rodriguez Matas JF, Dubini G, Boodt N, Gijsen FJH, et al. The first virtual patient-specific thrombectomy procedure. J Biomech. (2021) (In press). doi: 10.1016/j.jbiomech.2021.110622

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Jansen IG, Mulder MJ, Goldhoorn RJ. Endovascular tre-atment for acute ischaemic stroke in routine clinical practice: prospective, observational cohort study (MR CLEAN Registry). BMJ. (2018) 9;360:k949. doi: 10.1136/bmj.k949

CrossRef Full Text | Google Scholar

21. Moerman KM, Konduri PR, Fereidoonnezhad B, Marquering HA, van der Lugt A, Luraghi G, et al. Development of a patient-specific cerebral vasculature fluid-structure-interaction model. J Biomech. (2021) (In press). doi: 10.31224/osf.io/qaujs

CrossRef Full Text | Google Scholar

22. Piccinelli M, Veneziani A, Steinman DA, Remuzzi A, Antiga L. A framework for geometric analysis of vascular structures: application to cerebral aneurysms. IEEE Trans Med Imaging. (2009) 28:1141–55. doi: 10.1109/TMI.2009.2021652

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Moerman KM. GIBBON: the geometry and image-based bioengineering add-on. J Open Source Softw. (2018) 3:506. doi: 10.21105/joss.00506

CrossRef Full Text | Google Scholar

24. Dutra BG, Tolhuisen ML, Alves HC, Treurniet KM, Kappelhof M, Yoo AJ, et al. Thrombus imaging characteristics and outcomes in acute ischemic stroke patients undergoing endovascular treatment. Stroke. (2019) 50:2057–64. doi: 10.1161/STROKEAHA.118.024247

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Boodt N, Compagne KC, Dutra BG, Samuels N, Tolhuisen ML, Alves HC, et al. Stroke etiology and thrombus computed tomography characteristics in patients with acute ischemic stroke: a MR CLEAN Registry Substudy. Stroke. (2020) 51:1727–35. doi: 10.1161/STROKEAHA.119.027749

CrossRef Full Text | Google Scholar

26. Staessens S, Fitzgerald S, Andersson T, Clarençon F, Denorme F, Gounis MJ, et al. Histological stroke clot analysis after thrombectomy: technical aspects and recommendations. Int J Stroke. (2020) 15:467–76. doi: 10.1177/1747493019884527

PubMed Abstract | CrossRef Full Text | Google Scholar

Appendix

INSIST Investigators

Charles Majoie1, Henk Marquering1,2, Ed van Bavel2, Alfons Hoekstra3, Diederik W.J. Dippel4, Hester L. Lingsma5, Aad van der Lugt6, Noor Samuels4,5,6, Nikki Boodt4,5,6, Yvo Roos7, Simon de Meyer8, Senna Staessens8, Praneeta Konduri1,2, Nerea Arrarte Terreros2,Bastien Chopard9, Franck Raynaud9, Remy Petkantchin9, Vanessa Blanc-Guillemaud10, Mikhail Panteleev11,12, Alexey Shibeko11, Francesco Migliavacca13, Gabriele Dubini13, Giulia Luraghi13, Jose Felix Rodriguez Matas13, Sara Bridio13, Patrick Mc Garry14, Kevin Moerman14, Behrooz Fereidoonnezhad14, Michael Gilvarry15, Ray McCarthy15, Sharon Duffy15, Stephen Payne16, Tamas Jozsa16, Wahbi El-Bouri13, Sissy Georgakopoulou2, Victor Azizi3, Raymond Padmos3, Anushree Dwivedi15

1 Department of Radiology and Nuclear Medicine, Amsterdam UMC, location AMC, Amsterdam, Netherlands

2 Biomedical Engineering and Physics, Amsterdam UMC, location AMC, Amsterdam, Netherlands

3 Computational Science Lab, Faculty of Science, Institute for Informatics, University of Amsterdam, Amsterdam, Netherlands

4 Department of Neurology, Erasmus MC University Medical Center, Rotterdam, Netherlands

5 Department of Public Health, Erasmus MC University Medical Center, Rotterdam, Netherlands

6 Department of Radiology & Nuclear Medicine, Erasmus MC University Medical Center, Rotterdam, Netherlands

7 Department of Neurology, Amsterdam UMC, location AMC, Amsterdam, Netherlands

8 Laboratory for Thrombosis Research, KU Leuven Campus Kulak Kortrijk, Kortrijk, Belgium

9 Computer Science Department, University of Geneva, CUI, Carouge, Switzerland

10 Institut de Recherches Internationales Servier, Coubevoie Cedex, France

11 Center for Theoretical Problems of Physicochemical Pharmacology RAS, Moscow, Russia

12 Dmitry Rogachev National Research Center of Pediatric Hematology, Oncology and Immunology, Moscow, Russia; Faculty of Physics, Lomonosov Moscow State University, Moscow, Russia

13 Laboratory of Biological Structure Mechanics (LaBS), Department of Chemistry, Materials and Chemical Engineering “Giulio Natta”, Politecnico di Milano, Milano, Italy

14 College of Engineering and Informatics, National University of Ireland Galway, Ireland; National Centre for Biomedical Engineering Science, National University of Ireland Galway, Galway, Ireland

15 Cerenovus, Galway Neuro Technology Centre, Galway, Ireland

16 Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Oxford, United Kingdom

Keywords: insist, finite element analysis, carotid siphon, acute ischemic stroke, internal carotid artery, digital twin

Citation: Bridio S, Luraghi G, Rodriguez Matas JF, Dubini G, Giassi GG, Maggio G, Kawamoto JN, Moerman KM, McGarry P, Konduri PR, Arrarte Terreros N, Marquering HA, van Bavel E, Majoie CBLM and Migliavacca F (2021) Impact of the Internal Carotid Artery Morphology on in silico Stent-Retriever Thrombectomy Outcome. Front. Med. Technol. 3:719909. doi: 10.3389/fmedt.2021.719909

Received: 03 June 2021; Accepted: 09 July 2021;
Published: 03 August 2021.

Edited by:

Stephane Avril, Ecole des Mines de Saint-Étienne, France

Reviewed by:

Nenad Filipovic, University of Kragujevac, Serbia
Miquel Aguirre, Ecole des Mines de Saint-Étienne, France

Copyright © 2021 Bridio, Luraghi, Rodriguez Matas, Dubini, Giassi, Maggio, Kawamoto, Moerman, McGarry, Konduri, Arrarte Terreros, Marquering, van Bavel, Majoie and Migliavacca. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Sara Bridio, sara.bridio@polimi.it

These authors have contributed equally to this work and share first authorship

The INSIST investigators are given in the Appendix