The Role of Hemodynamics through the Circle of Willis in the Development of Intracranial Aneurysm: A Systematic Review of Numerical Models

Intra-Aneurysmal High-Resolution 4D MR Flow Imaging for Hemodynamic Imaging Markers in Intracranial Aneurysm Instability

BACKGROUND AND PURPOSE

Prediction of aneurysm instability is crucial to guide treatment decisions and to select appropriate patients with unruptured intracranial aneurysms (IAs) for preventive treatment. High-resolution 4D MR flow imaging and 3D quantification of aneurysm morphology could offer insights and new imaging markers for aneurysm instability. In this cross-sectional study, we aim to identify 4D MR flow imaging markers for aneurysm instability by relating hemodynamics in the aneurysm sac to 3D morphologic proxy parameters for aneurysm instability.

MATERIALS AND METHODS

In 35 patients with 37 unruptured IAs, a 3T MRA and a 7T 4D MRI flow scan were performed. Five hemodynamic parameters-peak-systolic wall shear stress (WSSMAX) and time-averaged wall shear stress (WSSMEAN), oscillatory shear index (OSI), mean velocity, and velocity pulsatility index-were correlated to 6 3D morphology proxy parameters of aneurysm instability-major axis length, volume, surface area (all 3 size parameters), flatness, shape index, and curvedness-by Pearson correlation with 95% CI. Scatterplots of hemodynamic parameters that correlated with IA size (major axis length) were created.

RESULTS

WSSMAX and WSSMEAN correlated negatively with all 3 size parameters (strongest for WSSMEAN with volume (r = -0.70, 95% CI -0.83 to -0.49) and OSI positively (strongest with major axis length [r = 0.87, 95% CI 0.76-0.93]). WSSMAX and WSSMEAN correlated positively with shape index (r = 0.61, 95% CI 0.36-0.78 and r = 0.49, 95% CI 0.20-0.70, respectively) and OSI negatively (r = -0.82, 95% CI -0.9 to -0.68). WSSMEAN and mean velocity correlated negatively with flatness (r = -0.35, 95% CI -0.61 to -0.029 and r = -0.33, 95% CI -0.59 to 0.007, respectively) and OSI positively (r = 0.54, 95% CI 0.26-0.74). Velocity pulsatility index did not show any statistically relevant correlation.

CONCLUSIONS

Out of the 5 included hemodynamic parameters, WSSMAX, WSSMEAN, and OSI showed the strongest correlation with morphologic 3D proxy parameters of aneurysm instability. Future studies should assess these promising new imaging marker parameters for predicting aneurysm instability in longitudinal cohorts of patients with IA.

Knowledge domain and emerging trends in the rupture risk of intracranial aneurysms research from 2004 to 2023

BACKGROUND

Intracranial aneurysms (IAs) pose significant health risks, attributable to their potential for sudden rupture, which can result in severe outcomes such as stroke and death. Despite extensive research, the variability of aneurysm behavior, with some remaining stable for years while others rupture unexpectedly, remains poorly understood.

AIM

To employ bibliometric analysis to map the research landscape concerning risk factors associated with IAs rupture.

METHODS

A systematic literature review of publications from 2004 to 2023 was conducted, analyzing 3804 documents from the Web of Science Core Collection database, with a focus on full-text articles and reviews in English. The analysis encompassed citation and co-citation networks, keyword bursts, and temporal trends to delineate the evolution of research themes and collaboration patterns. Advanced software tools, CiteSpace and VOSviewer, were utilized for comprehensive data visualization and trend analysis.

RESULTS

Analysis uncovered a total of 3804 publications on IA rupture risk factors between 2006 and 2023. Research interest surged after 2013, peaking in 2023. The United States led with 28.97% of publications, garnering 37706 citations. Notable United States-China collaborations were observed. Capital Medical University produced 184 publications, while Utrecht University boasted a citation average of 69.62 per publication. "World Neurosurgery" published the most papers, contrasting with "Stroke", the most cited journal. The PHASES score from "Lancet Neurology" emerged as a vital rupture risk prediction tool. Early research favored endovascular therapy, transitioning to magnetic resonance imaging and flow diverters. "Subarachnoid hemorrhage" stood out as a recurrent keyword.

CONCLUSION

This study assesses global IA research trends and highlights crucial gaps, guiding future investigations to improve preventive and therapeutic approaches.

Geometry and Symmetry of Willis' Circle and Middle Cerebral Artery Aneurysms Development

Background: A relationship between the geometry and symmetry of Willis' circle and intracranial aneurysms was reported for anterior communicating and posterior communicating (PCom) aneurysms. A similar association with the middle cerebral artery (MCA) aneurysms instead appeared weaker. Methods: We reviewed 432 patients from six Italian centers with unilateral MCA aneurysms, analyzing the relationship between the caliber and symmetry of Willis' circle and the presence of ruptured and unruptured presentation. CT-angiograms were evaluated to assess Willis' circle geometrical characteristics and the MCA aneurysm side, dimension and rupture status. Results: The hypoplasia of the first segment of the anterior cerebral artery (A1) was in approximately one-quarter of patients and PCom hypoplasia was in almost 40%. About 9% had a fetal PCom ipsilaterally to the aneurysm. By comparing the aneurysmal and healthy sides, only the PCom hypoplasia appeared significantly higher in the affected side. Finally, the caliber of the internal carotid artery (ICA) and the first segment of MCA (M1) caliber were significantly greater in patients with unruptured aneurysms, and PCom hypoplasia appeared related to the incidence of an ipsilateral MCA aneurysm and its risk of rupture. Conclusions: Although according to these findings asymmetries of Willis' circle are shown to be a risk factor for MCA aneurysm formation and rupture, the indifferent association with ipsilateral or contralateral hypoplasia remains a datum of difficult hemodynamic interpretation, thereby raising the concern that this association may be more casual than causal.

Influence of circle of Willis modeling on hemodynamic parameters in anterior communicating artery aneurysms and recommendations for model selection

Computational fluid dynamics (CFD) has been utilized to calculate hemodynamic parameters in anterior communicating artery aneurysm (AComA), which is located at a junction between left and right A1 and A2 segments. However, complete or half circle of Willis (CoW) models are used indiscriminately. This study aims to suggest recommendations for determining suitable CoW model. Five patient-specific CoW models with AComA were used, and each model was divided into complete, left-half, and right-half models. After validating the CFD using a flow experiment, the hemodynamic parameters and flow patterns in five AComAs were compared. In four out of five cases, inflow from one A1 side had a dominant influence on the AComA, while both left and right A1 sides affected the AComA in the remaining case. Also, the average difference in time-averaged wall shear stress between the complete and half models for four cases was 4.6%, but it was 62% in the other case. The differences in the vascular resistances of left and right A1 and A2 segments greatly influenced the flow patterns in the AComA. These results may help to enhance clinicians' understanding of blood flow in the brain, leading to improvements in diagnosis and treatment of cerebral aneurysms.

Image-based hemodynamic simulations for intracranial aneurysms: the impact of complex vasculature

PURPOSE

Hemodynamics play an important role in the assessment of intracranial aneurysm (IA) development and rupture risk. The purpose of this study was to examine the impact of complex vasculatures onto the intra-vessel and intra-aneurysmal blood flow.

METHODS

Complex segmentation of a subject-specific, 60-outlet and 3-inlet circle of Willis model captured with 7T magnetic resonance imaging was performed. This model was trimmed to a 10-outlet model version. Two patient-specific IAs were added onto both models yielding two pathological versions, and image-based blood flow simulations of the four resulting cases were carried out. To capture the differences between complex and trimmed model, time-averaged and centerline velocities were compared. The assessment of intra-saccular blood flow within the IAs involved the evaluation of wall shear stresses (WSS) at the IA wall and neck inflow rates (NIR).

RESULTS

Lower flow values are observed in the majority of the complex model. However, at specific locations (left middle cerebral artery 0.5 m/s, left posterior cerebral artery 0.25 m/s), higher flow rates were visible when compared to the trimmed counterpart. Furthermore, at the centerlines the total velocity values reveal differences up to 0.15 m/s. In the IAs, the reduction in the neck inflow rate and WSS in the complex model was observed for the first IA (IA-A δNIRmean =  - 0.07ml/s, PCA.l δWSSmean =  - 0.05 Pa). The second IA featured an increase in the neck inflow rate and WSS (IA-B δNIRmean = 0.04 ml/s, PCA.l δWSSmean = 0.07 Pa).

CONCLUSION

Both the magnitude and shape of the flow distribution vary depending on the model's complexity. The magnitude is primarily influenced by the global vessel model, while the shape is determined by the local structure. Furthermore, intra-aneurysmal flow strongly depends on the location in the vessel tree, emphasizing the need for complex model geometries for realistic hemodynamic assessment and rupture risk analysis.

Hemodynamic Parameters in the Parent Arteries of Unruptured Intracranial Aneurysms Depend on Aneurysm Size and Are Different Compared to Contralateral Arteries: A 7 Tesla 4D Flow MRI Study

BACKGROUND

Different Circle of Willis (CoW) variants have variable prevalences of aneurysm development, but the hemodynamic variation along the CoW and its relation to presence and size of unruptured intracranial aneurysms (UIAs) are not well known.

PURPOSE

Gain insight into hemodynamic imaging markers of the CoW for UIA development by comparing these outcomes to the corresponding contralateral artery without an UIA using 4D flow magnetic resonance imaging (MRI).

STUDY TYPE

Retrospective, cross-sectional study.

SUBJECTS

Thirty-eight patients with an UIA, whereby 27 were women and a mean age of 62 years old.

FIELD STRENGTH/SEQUENCE

Four-dimensional phase-contrast (PC) MRI with a 3D time-resolved velocity encoded gradient echo sequence at 7 T.

ASSESSMENT

Hemodynamic parameters (blood flow, velocity pulsatility index [vPI], mean velocity, distensibility, and wall shear stress [peak systolic (WSSMAX ), and time-averaged (WSSMEAN )]) in the parent artery of the UIA were compared to the corresponding contralateral artery without an UIA and were related to UIA size.

STATISTICAL TESTS

Paired t-tests and Pearson Correlation tests. The threshold for statistical significance was P < 0.05 (two-tailed).

RESULTS

Blood flow, mean velocity, WSSMAX , and WSSMEAN were significantly higher, while vPI was lower, in the parent artery relative to contralateral artery. The WSSMAX of the parent artery significantly increased linearly while the WSSMEAN decreased linearly with increasing UIA size.

CONCLUSIONS

Hemodynamic parameters and WSS differ between parent vessels of UIAs and corresponding contralateral vessels. WSS correlates with UIA size, supporting a potential hemodynamic role in aneurysm pathology.

LEVEL OF EVIDENCE

2 TECHNICAL EFFICACY: Stage 2.

Personalised simulation of hemodynamics in cerebrovascular disease: lessons learned from a study of diagnostic accuracy

Intracranial atherosclerotic disease (ICAD) poses a significant risk of subsequent stroke but current prevention strategies are limited. Mechanistic simulations of brain hemodynamics offer an alternative precision medicine approach by utilising individual patient characteristics. For clinical use, however, current simulation frameworks have insufficient validation. In this study, we performed the first quantitative validation of a simulation-based precision medicine framework to assess cerebral hemodynamics in patients with ICAD against clinical standard perfusion imaging. In a retrospective analysis, we used a 0-dimensional simulation model to detect brain areas that are hemodynamically vulnerable to subsequent stroke. The main outcome measures were sensitivity, specificity, and area under the receiver operating characteristics curve (ROC AUC) of the simulation to identify brain areas vulnerable to subsequent stroke as defined by quantitative measurements of relative mean transit time (relMTT) from dynamic susceptibility contrast MRI (DSC-MRI). In 68 subjects with unilateral stenosis >70% of the internal carotid artery (ICA) or middle cerebral artery (MCA), the sensitivity and specificity of the simulation were 0.65 and 0.67, respectively. The ROC AUC was 0.68. The low-to-moderate accuracy of the simulation may be attributed to assumptions of Newtonian blood flow, rigid vessel walls, and the use of time-of-flight MRI for geometric representation of subject vasculature. Future simulation approaches should focus on integrating additional patient data, increasing accessibility of precision medicine tools to clinicians, addressing disease burden disparities amongst different populations, and quantifying patient benefit. Our results underscore the need for further improvement of mechanistic simulations of brain hemodynamics to foster the translation of the technology to clinical practice.

Validation of a cerebral hemodynamic model with personalized calibration in patients with aneurysmal subarachnoid hemorrhage

This study aims to validate a numerical model developed for assessing personalized circle of Willis (CoW) hemodynamics under pathological conditions. Based on 66 computed tomography angiography images, investigations were obtained from 43 acute aneurysmal subarachnoid hemorrhage (aSAH) patients from a local neurovascular center. The mean flow velocity of each artery in the CoW measured using transcranial Doppler (TCD) and simulated by the numerical model was obtained for comparison. The intraclass correlation coefficient (ICC) over all cerebral arteries for TCD and the numerical model was 0.88 (N = 561; 95% CI 0.84-0.90). In a subgroup of patients who had developed delayed cerebral ischemia (DCI), the ICC had decreased to 0.72 but remained constant with respect to changes in blood pressure, Fisher grade, and location of ruptured aneurysm. Our numerical model showed good agreement with TCD in assessing the flow velocity in the CoW of patients with aSAH. In conclusion, the proposed model can satisfactorily reproduce the cerebral hemodynamics under aSAH conditions by personalizing the numerical model with TCD measurements. Clinical trial registration: [http://www.trialregister.nl/], identifier [NL8114].