Stagnation and complex flow in ruptured cerebral aneurysms: a possible association with hemostatic pattern

Computational fluid dynamics-based virtual angiograms for the detection of flow stagnation in intracranial aneurysms

The goal of this study was to test if CFD-based virtual angiograms could be used to automatically discriminate between intracranial aneurysms (IAs) with and without flow stagnation. Time density curves (TDC) were extracted from patient digital subtraction angiography (DSA) image sequences by computing the average gray level intensity inside the aneurysm region and used to define injection profiles for each subject. Subject-specific 3D models were reconstructed from 3D rotational angiography (3DRA) and computational fluid dynamics (CFD) simulations were performed to simulate the blood flow inside IAs. Transport equations were solved numerically to simulate the dynamics of contrast injection into the parent arteries and IAs and then the contrast retention time (RET) was calculated. The importance of gravitational pooling of contrast agent within the aneurysm was evaluated by modeling contrast agent and blood as a mixture of two fluids with different densities and viscosities. Virtual angiograms can reproduce DSA sequences if the correct injection profile is used. RET can identify aneurysms with significant flow stagnation even when the injection profile is not known. Using a small sample of 14 IAs of which seven were previously classified as having flow stagnation, it was found that a threshold RET value of 0.46 s can successfully identify flow stagnation. CFD-based prediction of stagnation was in more than 90% agreement with independent visual DSA assessment of stagnation in a second sample of 34 IAs. While gravitational pooling prolonged contrast retention time it did not affect the predictive capabilities of RET. CFD-based virtual angiograms can detect flow stagnation in IAs and can be used to automatically identify aneurysms with flow stagnation even without including gravitational effects on contrast agents.

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

Background: The role of regional hemodynamics in the intracranial aneurysmal formation, growth, and rupture has been widely discussed based on numerical models over the past decades. Variation of the circle of Willis (CoW), which results in hemodynamic changes, is associated with the aneurysmal formation and rupture. However, such correlation has not been further clarified yet. The aim of this systematic review is to investigate whether simulated hemodynamic indices of the CoW are relevant to the formation, growth, or rupture of intracranial aneurysm. Methods: We conducted a review of MEDLINE, Web of Science, and EMBASE for studies on the correlation between hemodynamics indices of the CoW derived from numerical models and intracranial aneurysm up to December 2020 in compliance with PRISMA guidelines. Results: Three case reports out of 1046 publications met our inclusion and exclusion criteria, reporting 13 aneurysms in six patients. Eleven aneurysms were unruptured, and the state of the other two aneurysms was unknown. Wall shear stress, oscillatory shear index, von-Mises tension, flow velocity, and flow rate were reported as hemodynamic indices. Due to limited cases and significant heterogeneity between study settings, meta-analysis could not be performed. Conclusion: Numerical models can provide comprehensive information on the cerebral blood flow as well as local flow characteristics in the intracranial aneurysm. Based on only three case reports, no firm conclusion can be drawn regarding the correlation between hemodynamic parameters in the CoW derived from numerical models and aneurysmal formation or rupture. Due to the inherent nature of numerical models, more sensitive analysis and rigorous validations are required to determine its measurement error and thus extend their application into clinical practice for personalized management. Prospero registration number: CRD42021125169.

Complex Intracranial Aneurysms

Complex intracranial aneurysms remain challenging to treat using standard microsurgical or endovascular techniques. These aneurysms often require a combination of deconstructive and reconstructive procedures, such as parent artery occlusion, flow alteration, and blind-alley formation with or without bypass surgery, for effective and enduring therapeutic effects. It is important to determine the type of bypass based on the site of occlusion of the patent artery, anatomical features of the distal vessels, and expected adequate blood flow. In this chapter, we describe the "Standards," "Advances," and "Controversies" in the context of a microsurgical treatment strategy for complex intracranial aneurysms. "Standards" include a combination of frequent and commonly used procedures that have been gathering a certain consensus on their effectiveness. "Advances" include infrequent, demanding, and/or uncertain surgical procedures that are currently under debate. Finally, "Controversies" discuss a number of unsolved issues.

Quantitative assessment of microstructural evolution of intracranial aneurysm wall by vessel wall imaging

PURPOSE

This study aimed to evaluate new quantitative parameters of aneurysm wall enhancement (AWE) on magnetic resonance vessel wall imaging (VWI) in differentiating between the stable and evolving unruptured intracranial aneurysms (UIAs).

METHODS

Thirty-eight consecutive patients with UIAs (27 stable and 11 evolving) underwent VWI with contrast-enhanced 3D T1 volume isotropic turbo spin echo acquisition. The voxel-based enhancement maps were created using pre- and post-contrast images. The aneurysmal lumen with signal suppression by black-blood method was segmented. Then, one voxel outer and inner layers of the lumen contour were automatically segmented. The shape features of the aneurysms and AWE of the two layers were compared between stable and evolving groups.

RESULTS

The shape features, including aneurysm volume, surface, and compacity were significantly different between the stable and evolving groups (P = 0.024, 0.028, and 0.033, respectively). Stable and evolving groups also differed significantly in the AWE at the union of outer and inner layers of the aneurysm wall (P = 0.0082) but not in that of the outer or inner layer alone. Multivariate logistic regression analysis revealed significant differences in aneurysm volume, surface, and AWE at the union of outer and inner layers between the two groups (P = 0.0029, 0.0092, and 0.0033, respectively). Receiver operating characteristics curve analysis revealed that the area under the curve of the logistic regression model was 0.89.

CONCLUSION

Quantitative combined analysis of aneurysm shape features and AWE of the union of outer and inner layers were effective for differentiating between stable and evolving UIAs.

Disturbed flow's impact on cellular changes indicative of vascular aneurysm initiation, expansion, and rupture: A pathological and methodological review

Aneurysms are malformations within the arterial vasculature brought on by the structural breakdown of the microarchitecture of the vessel wall, with aneurysms posing serious health risks in the event of their rupture. Blood flow within vessels is generally laminar with high, unidirectional wall shear stressors that modulate vascular endothelial cell functionality and regulate vascular smooth muscle cells. However, altered vascular geometry induced by bifurcations, significant curvature, stenosis, or clinical interventions can alter the flow, generating low stressor disturbed flow patterns. Disturbed flow is associated with altered cellular morphology, upregulated expression of proteins modulating inflammation, decreased regulation of vascular permeability, degraded extracellular matrix, and heightened cellular apoptosis. The understanding of the effects disturbed flow has on the cellular cascades which initiate aneurysms and promote their subsequent growth can further elucidate the nature of this complex pathology. This review summarizes the current knowledge about the disturbed flow and its relation to aneurysm pathology, the methods used to investigate these relations, as well as how such knowledge has impacted clinical treatment methodologies. This information can contribute to the understanding of the development, growth, and rupture of aneurysms and help develop novel research and aneurysmal treatment techniques.

Computational Fluid Dynamics for Cerebral Aneurysms in Clinical Settings

Hemodynamics is thought to play an important role in the pathogenesis of cerebral aneurysms and recent development of computer technology makes it possible to simulate blood flow using high-resolution 3D images within several hours. A lot of studies of computational fluid dynamics (CFD) for cerebral aneurysms were reported; therefore, application of CFD for cerebral aneurysms in clinical settings is reviewed in this article.CFD for cerebral aneurysms using a patient-specific geometry model was first reported in 2003 and it has been revealing that hemodynamics brings a certain contribution to understanding aneurysm pathology, including initiation, growth and rupture. Based on the knowledge of the state-of-the-art techniques, this review treats the decision-making process for using CFD in several clinical settings. We introduce our CFD procedure using digital imaging and communication in medicine (DICOM) datasets of 3D CT angiography or 3D rotational angiography. In addition, we review rupture status, hyperplastic remodeling of aneurysm wall, and recurrence of coiled aneurysms using the hemodynamic parameters such as wall shear stress (WSS), oscillatory shear index (OSI), aneurysmal inflow rate coefficient (AIRC), and residual flow volume (RFV).

Machine Learning Classification of Cerebral Aneurysm Rupture Status with Morphologic Variables and Hemodynamic Parameters

Purpose

To construct a classification model of rupture status and to clarify the importance of morphologic variables and hemodynamic parameters on rupture status by applying a machine learning (ML) algorithm to morphologic and hemodynamic data of cerebral aneurysms.

Materials and Methods

Between 2011 and 2019, 226 (112 ruptured and 114 unruptured) cerebral aneurysms in 188 consecutive patients were retrospectively analyzed with computational fluid dynamics (CFD). A random forest ML algorithm was applied to the results to create three classification models consisting of only morphologic variables (model 1), only hemodynamic parameters (model 2), and both morphologic variables and hemodynamic parameters (model 3). The accuracy of rupture status classification and the importance of each variable or parameter in the models were computed.

Results

The accuracy was 77.0% in model 1, 71.2% in model 2, and 78.3% in model 3. The three most important features were projection ratio, size ratio, and aspect ratio in model 1; low shear area ratio, oscillatory shear index, and oscillatory velocity index in model 2; and projection ratio, irregular shape, and size ratio in model 3.

Conclusion

Classification models of rupture status of cerebral aneurysms were constructed by applying an ML algorithm to morphologic variables and hemodynamic parameters. The model worked with relatively high accuracy, in which projection ratio, irregular shape, and size ratio were important for the discrimination of ruptured aneurysms.Supplemental material is available for this article.© RSNA, 2020.

Prediction of the extent of thrombus formation in the parent artery after endovascular occlusion of a distal anterior cerebral artery aneurysm using computational fluid dynamics

Endovascular coiling of a cerebral aneurysm and coil occlusion of the parent artery have been occasionally performed to treat cerebral aneurysms; however, it is difficult to predict the accurate extent of thrombus formation in the parent artery proximal to the coiled aneurysm and the coil-occluded parent artery preoperatively, and unexpected occlusion of the arterial branches can occur by thrombus extension into or in the parent artery. The authors describe a case of a distal anterior cerebral artery (ACA) aneurysm treated by endovascular parent artery occlusion (PAO) with preoperative computational fluid dynamics (CFD) prediction of the extent of thrombus formation. A 73-year-old woman presented with subarachnoid hemorrhage and an aneurysm that was located on the right pericallosal artery distal to the paracentral artery bifurcation. Endovascular coiling of the aneurysm and the pericallosal artery was planned. In advance of the treatment, CFD was performed to predict the extent of thrombus formation with specific wall shear stress and shear rate thresholds. The hemodynamic results indicated that coiling of the aneurysm resulted in thrombus formation in the pericallosal artery up to just distal to the paracentral artery ostium; therefore, the treatment was implemented according to the CFD prediction. Postoperative digital subtraction angiography revealed that the extent of thrombus formation was consistent with the preoperative CFD prediction. This technique may prevent unexpected occlusion of arterial branches.

Quantification of hemodynamic irregularity using oscillatory velocity index in the associations with the rupture status of cerebral aneurysms

BACKGROUND

Complex and unstable flow patterns are reported to be associated with the rupture status of cerebral aneurysms, while their evaluation depends on qualitative analysis of streamlines of bloodflow. Oscillatory velocity index (OVI) is a hemodynamic parameter to quantify flow patterns. The aim of this study is to elucidate the associations between OVI and the rupture status of cerebral aneurysms.

METHODS

One hundred and twenty-nine ruptured and unruptured cerebral aneurysms were analyzed with computational fluid dynamics under pulsatile flow conditions. With the use of median value of OVI, all aneurysms were divided into high and low OVI groups. Statistical analysis was performed to compare rupture status, and morphological and hemodynamic parameters between the two groups.

RESULTS

The median value of OVI was 0.006. High OVI was more likely observed in ruptured aneurysms (P=0.028) and associated with irregular shape, complex flow patterns, and unstable flow patterns (P<0.001, respectively). In morphological parameters, maximum size, aspect, projection, size, and volume-to-ostium area ratios were significantly higher in the high OVI group (P<0.001, respectively). In hemodynamic parameters, wall shear stress and wall shear stress gradient were significantly lower, and oscillatory shear index and gradient oscillatory number were significantly higher in the high OVI group (P<0.001, respectively).

CONCLUSION

High OVI was associated with rupture status, and morphological and hemodynamic characteristics of ruptured aneurysms. These results indicate that OVI may serve as a valuable hemodynamic parameter for diagnosing rupture status and risks of aneurysms.

One-stage Stent-assisted Coil Embolization for Rupture-side-unknown Bilateral Vertebral Artery Dissecting Aneurysms in an Acute Stage: A Case Report

Bilateral vertebral artery dissecting aneurysms (VADAs) with subarachnoid hemorrhage (SAH) are rare and their management is still challenging. In this report, we successfully performed one-stage stent-assisted coil embolization (SAC) for bilateral VADAs with SAH in an acute stage, because the ruptured side could not be diagnosed. A 47-year-old woman presented with a sudden onset of headache without laterality, and left-side dominant SAH with bilateral VADAs was noted on computed tomography (CT) scans. The size of aneurysmal dome and neck was similar between the two VADAs, and a bleb was observed only on the right VADA. In computational fluid dynamics (CFD) simulations, findings of wall shear stress (WSS), normalized WSS, and WSS gradient suggested that the left VADA was ruptured, while the oscillatory shear index and aneurysm formation indicator suggested the opposite-side one to be ruptured. Thus, we could not determine which VADA was ruptured by clinical data and CFD analyses. Therefore, we performed simultaneous treatment for the bilateral VADAs by using SAC technique 8 h after the onset under dual antiplatelet and anticoagulation therapies. There was no evidence of rebleeding and stent thrombosis. Stent thrombosis was monitored by duplex color-coded ultrasonography after the intervention. She was discharged without neurological deficits, and 6-month follow-up cerebral angiography demonstrated no recanalization of VADAs. This is the first report showing bilateral VADAs with SAH treated by one-stage SAC within 24 h of SAH, and the potential risks are discussed.

Hemodynamic characteristics of hyperplastic remodeling lesions in cerebral aneurysms

Background & purpose

Hyperplastic remodeling (HR) lesions are sometimes found on cerebral aneurysm walls. Atherosclerosis is the results of HR, which may cause an adverse effect on surgical treatment for cerebral aneurysms. Previous studies have demonstrated that atherosclerotic changes had a correlation with certain hemodynamic characteristics. Therefore, we investigated local hemodynamic characteristics of HR lesions of cerebral aneurysms using computational fluid dynamics (CFD).

Methods

Twenty-four cerebral aneurysms were investigated using CFD and intraoperative video recordings. HR lesions and red walls were confirmed on the intraoperative images, and the qualification points were determined on the center of the HR lesions and the red walls. The qualification points were set on the virtual operative images for evaluation of wall shear stress (WSS), normalized WSS (NWSS), oscillatory shear index (OSI), relative residence time (RRT), and aneurysm formation indicator (AFI). These hemodynamic parameters at the qualification points were compared between HR lesions and red walls.

Results

HR lesions had lower NWSS, lower AFI, higher OSI and prolonged RRT compared with red walls. From analysis of the receiver-operating characteristic curve for hemodynamic parameters, OSI was the most optimal hemodynamic parameter to predict HR lesions (area under the curve, 0.745; 95% confidence interval, 0.603–0.887; cutoff value, 0.00917; sensitivity, 0.643; specificity, 0.893; P<0.01). With multivariate logistic regression analyses using stepwise method, NWSS was significantly associated with the HR lesions.

Conclusions

Although low NWSS was independently associated with HR lesions, OSI is the most valuable hemodynamic parameter to distinguish HR lesions from red walls.