Hemodynamic findings associated with intraoperative appearances of intracranial aneurysms

Fast simulation of hemodynamics in intracranial aneurysms for clinical use

BACKGROUND

A widely accepted tool to assess hemodynamics, one of the most important factors in aneurysm pathophysiology, is Computational Fluid Dynamics (CFD). As current workflows are still time consuming and difficult to operate, CFD is not yet a standard tool in the clinical setting. There it could provide valuable information on aneurysm treatment, especially regarding local risks of rupture, which might help to optimize the individualized strategy of neurosurgical dissection during microsurgical aneurysm clipping.

METHOD

We established and validated a semi-automated workflow using 3D rotational angiographies of 24 intracranial aneurysms from patients having received aneurysm treatment at our centre. Reconstruction of vessel geometry and generation of volume meshes was performed using AMIRA 6.2.0 and ICEM 17.1. For solving ANSYS CFX was used. For validational checks, tests regarding the volumetric impact of smoothing operations, the impact of mesh sizes on the results (grid convergence), geometric mesh quality and time tests for the time needed to perform the workflow were conducted in subgroups.

RESULTS

Most of the steps of the workflow were performed directly on the 3D images requiring no programming experience. The workflow led to final CFD results in a mean time of 22 min 51.4 s (95%-CI 20 min 51.562 s-24 min 51.238 s, n = 5). Volume of the geometries after pre-processing was in mean 4.46% higher than before in the analysed subgroup (95%-CI 3.43-5.50%). Regarding mesh sizes, mean relative aberrations of 2.30% (95%-CI 1.51-3.09%) were found for surface meshes and between 1.40% (95%-CI 1.07-1.72%) and 2.61% (95%-CI 1.93-3.29%) for volume meshes. Acceptable geometric mesh quality of volume meshes was found.

CONCLUSIONS

We developed a semi-automated workflow for aneurysm CFD to benefit from hemodynamic data in the clinical setting. The ease of handling opens the workflow to clinicians untrained in programming. As previous studies have found that the distribution of hemodynamic parameters correlates with thin-walled aneurysm areas susceptible to rupture, these data might be beneficial for the operating neurosurgeon during aneurysm surgery, even in acute cases.

Intracranial aneurysm CTA images and 3D models dataset with clinical morphological and hemodynamic data

Intracranial aneurysm is a cerebrovascular disease associated with a high rupture risk, often resulting in death or severe disability. Recent advances in AI enable the prediction of intracranial aneurysm initiation, progression, and rupture through medical image analysis. Despite growing research interest, there is a shortage of publicly available datasets for training and validating AI models. This paper presents a comprehensive dataset comprising high-resolution CTA images of 99 patients with 105 MCA aneurysms and 44 normal healthy controls, along with their respective clinical data and 3D models of aneurysms and the parent arteries derived from the CTA images. Furthermore, recognizing the significance of blood hemodynamics on aneurysm development, this dataset also included the morphological and hemodynamic parameters obtained by computational fluid dynamics (CFD) for each patient and healthy control, which can be utilized by researchers without prior CFD experience. This dataset will facilitate hypothesis-driven or data-driven research on intracranial aneurysms, and has the potential to deepen our understanding of this disease.

Transition of intracranial aneurysmal wall enhancement from high to low wall shear stress mediation with size increase: A hemodynamic study based on 7T magnetic resonance imaging

Background

Wall shear stress (WSS) has been proved to be related to the formation, development and rupture of intracranial aneurysms. Aneurysm wall enhancement (AWE) on magnetic resonance imaging (MRI) can be caused by inflammation and have confirmed its relationship with low WSS. High WSS can also result in inflammation but the research of its correlation with AWE is lack because of the focus on large aneurysms limited by 3T MRI in most previous studies.This study aimed to assess the potential association between high or low WSS and AWE in different aneuryms. Especially the relationship between high WSS and AWE in small aneurysm.

Methods

Forty-three unruptured intracranial aneurysms in 42 patients were prospectively included for analysis. 7.0 T MRI was used for imaging. Aneurysm size was measured on three-dimensional time-of-flight (TOF) images. Aneurysm-to-pituitary stalk contrast ratio (CRstalk) was calculated on post-contrast black-blood T1-weighted fast spin echo sequence images. Hemodynamics were assessed by four-dimensional flow MRI.

Results

The small aneurysms group had more positive WSS-CRstalk correlation coefficient distribution (dome: 78.6 %, p = 0.009; body: 50.0 %, p = 0.025), and large group had more negative coefficient distribution (dome: 44.8 %, p = 0.001; body: 69.0 %, p = 0.002). Aneurysm size was positively correlated with the significant OSI-CRstalk correlation coefficient at the dome (p = 0.012) and body (p = 0.010) but negatively correlated with the significant WSS-CRstalk correlation coefficient at the dome (p < 0.001) and body (p = 0.017).

Conclusion

AWE can be mediated by both high and low WSS, and translate from high WSS- to low WSS-mediated pathways as size increase. Additionally, AWE may serve as an indicator of the stage of aneurysm development via different correlations with hemodynamic factors.

Enhancing cerebral vasculature analysis with pathlength-corrected 2D angiographic parametric imaging: A feasibility study

BACKGROUND

2D angiographic parametric imaging (API) quantitatively extracts imaging biomarkers related to contrast flow and is conventionally applied to 2D digitally subtracted angiograms (DSA's). In the interventional suite, API is typically performed using 1-2 projection views and is limited by vessel overlap, foreshortening, and depth-integration of contrast motion.

PURPOSE

This work explores the use of a pathlength-correction metric to overcome the limitations of 2D-API: the primary objective was to study the effect of converting 3D contrast flow to projected contrast flow using a simulated angiographic framework created with computational fluid dynamics (CFD) simulations, thereby removing acquisition variability.

METHODS

The pathlength-correction framework was applied to in-silico angiograms, generating a reference (i.e., ground-truth) volumetric contrast distribution in four patient-specific intracranial aneurysm geometries. Biplane projections of contrast flow were created from the reference volumetric contrast distributions, assuming a cone-beam geometry. A Parker-weighted reconstruction was performed to obtain a binary representation of the vessel structure in 3D. Standard ray tracing techniques were then used to track the intersection of a ray from the focal spot with each voxel of the reconstructed vessel wall to a pixel in the detector plane. The lengths of each ray through the 3D vessel lumen were then projected along each ray-path to create a pathlength-correction map, where the pixel intensity in the detector plane corresponds to the vessel width along each source-detector ray. By dividing the projection sequences with this correction map, 2D pathlength-corrected in-silico angiograms were obtained. We then performed voxel-wise (3D) API on the ground-truth contrast distribution and compared it to pixel-wise (2D) API, both with and without pathlength correction for each biplane view. The percentage difference (PD) between the resultant API biomarkers in each dataset were calculated within the aneurysm region of interest (ROI).

RESULTS

Intensity-based API parameters, such as the area under the curve (AUC) and peak height (PH), exhibited notable changes in magnitude and spatial distribution following pathlength correction: these now accurately represent conservation of mass of injected contrast media within each arterial geometry and accurately reflect regions of stagnation and recirculation in each aneurysm ROI. Improved agreement was observed between these biomarkers in the pathlength-corrected biplane maps: the maximum PD within the aneurysm ROI is 3.3% with pathlength correction and 47.7% without pathlength correction. As expected, improved agreement with ROI-averaged ground-truth 3D counterparts was observed for all aneurysm geometries, particularly large aneurysms: the maximum PD for both AUC and PH was 5.8%. Temporal parameters (mean transit time, MTT, time-to-peak, TTP, time-to-arrival, TTA) remained unaffected after pathlength correction.

CONCLUSIONS

This study indicates that the values of intensity-based API parameters obtained with conventional 2D-API, without pathlength correction, are highly dependent on the projection orientation, and uncorrected API should be avoided for hemodynamic analysis. The proposed metric can standardize 2D API-derived biomarkers independent of projection orientation, potentially improving the diagnostic value of all acquired 2D-DSA's. Integration of a pathlength correction map into the imaging process can allow for improved interpretation of biomarkers in 2D space, which may lead to improved diagnostic accuracy during procedures involving the cerebral vasculature.

Hemodynamic Analysis Shows High Wall Shear Stress Is Associated with Intraoperatively Observed Thin Wall Regions of Intracranial Aneurysms

BACKGROUND

Studying the relationship between hemodynamics and local intracranial aneurysm (IA) pathobiology can help us understand the natural history of IA. We characterized the relationship between the IA wall appearance, using intraoperative imaging, and the hemodynamics from CFD simulations.

METHODS

Three-dimensional geometries of 15 IAs were constructed and used for CFD. Two-dimensional intraoperative images were subjected to wall classification using a machine learning approach, after which the wall type was mapped onto the 3D surface. IA wall regions included thick (white), normal (purple-crimson), and thin/translucent (red) regions. IA-wide and local statistical analyses were performed to assess the relationship between hemodynamics and wall type.

RESULTS

Thin regions of the IA sac had significantly higher WSS, Normalized WSS, WSS Divergence and Transverse WSS, compared to both normal and thick regions. Thicker regions tended to co-locate with significantly higher RRT than thin regions. These trends were observed on a local scale as well. Regression analysis showed a significant positive correlation between WSS and thin regions and a significant negative correlation between WSSD and thick regions.

CONCLUSION

Hemodynamic simulation results were associated with the intraoperatively observed IA wall type. We consistently found that elevated WSS and WSS_Norm were associated with thin regions of the IA wall rather than thick and normal regions.

Leveraging Patient-Specific Simulated Angiograms to Characterize Cerebral Aneurysm Hemodynamics using Computational Fluid Dynamics

Cerebral aneurysms (CA) affect nearly 6% of the US population and its rupture is one of the major causes of hemorrhagic stroke. Neurointerventionalists performing endovascular therapy (ET) to treat CA rely on qualitative image sequences obtained under fluoroscopy guidance alone, and do not have access to crucial quantitative information regarding blood flow before, during and after treatment - partially contributing to a failure rate of up to 30%. Computational fluid dynamics (CFD) is a powerful tool that can provide a wealth of quantitative data; however, CFD has found limited utility in the clinic due to the challenges in obtaining hemodynamic boundary conditions for each patient. In this work, we present a novel CFD-based simulated angiogram approach (SAA) that resolves the blood flow physics and interaction between blood and injected contrast agent to extract quantitative hemodynamic parameters which can be used to design real-time parametric imaging analysis. The SAA enables correlating contrast agent transport to the underlying hemodynamic conditions via time-density curves (TDC) obtained at several points in the region of interest. The ability of the TDC and the SAA to provide critical hemodynamic parameters in and around CA anatomies, such as washout and local flow changes is explored and presented. This provides invaluable quantitative data to the clinician at the time of intervention, since it incorporates the physics of blood flow and correlates the contrast transport to hemodynamic parameters quantitatively - thereby enabling the clinician to take informed decisions that improve treatment outcomes.

Rebleeding of Ruptured Intracranial Aneurysm After Admission: A Multidimensional Nomogram Model to Risk Assessment

Objective

Rebleeding is recognized as the main cause of mortality after intracranial aneurysm rupture. Though timely intervention can prevent poor prognosis, there is no agreement on the surgical priority and choosing medical treatment for a short period after rupture. The aim of this study was to investigate the risk factors related to the rebleeding after admission and establish predicting models for better clinical decision-making.

Methods

The patients with ruptured intracranial aneurysms (RIAs) between January 2018 and September 2020 were reviewed. All patients fell to the primary and the validation cohort by January 2020. The hemodynamic parameters were determined through the computational fluid dynamics simulation. Cox regression analysis was conducted to identify the risk factors of rebleeding. Based on the independent risk factors, nomogram models were built, and their predicting accuracy was assessed by using the area under the curves (AUCs).

Result

A total of 577 patients with RIAs were enrolled in this present study, 86 patients of them were identified as undergoing rebleeding after admission. Thirteen parameters were identified as significantly different between stable and rebleeding aneurysms in the primary cohort. Cox regression analysis demonstrated that six parameters, including hypertension [hazard ratio (HR), 2.54; P = 0.044], bifurcation site (HR, 1.95; P = 0.013), irregular shape (HR, 4.22; P = 0.002), aspect ratio (HR, 12.91; P < 0.001), normalized wall shear stress average (HR, 0.16; P = 0.002), and oscillatory stress index (HR, 1.14; P < 0.001) were independent risk factors related to the rebleeding after admission. Two nomograms were established, the nomogram including clinical, morphological, and hemodynamic features (CMH nomogram) had the highest predicting accuracy (AUC, 0.92), followed by the nomogram including clinical and morphological features (CM nomogram; AUC, 0.83), ELAPSS score (AUC, 0.61), and PHASES score (AUC, 0.54). The calibration curve for the probability of rebleeding showed good agreement between prediction by nomograms and actual observation. In the validation cohort, the discrimination of the CMH nomogram was superior to the other models (AUC, 0.93 vs. 0.86, 0.71 and 0.48).

Conclusion

We presented two nomogram models, named CMH nomogram and CM nomogram, which could assist in identifying the RIAs with high risk of rebleeding.

Risk Factors for Recanalization after Coil Embolization

The aim of our study was to identify risk factors for recanalization 6 months after coil embolization using clinical data followed by computational fluid dynamics (CFD) analysis. Methods: Firstly, clinical data of 184 patients treated with coil embolization were analyzed retrospectively. Secondly, aneurysm models for high/low recanalization risk were generated based on ROC curves and their cut-off points. Afterward, CFD was utilized to validate the results. Results: In multivariable analysis, aneurysm filling during the first embolization was an independent risk factor whilst packing density was a protective factor of recanalization after 6 months in patients with aSAH. For patients with unruptured aneurysms, packing density was found to be a protective factor whilst the aneurysm neck size was an independent risk factor. Complex flow pattern and multiple vortices were associated with aneurysm shape and were characteristic of the high recanalization risk group. Conclusions: Statistical analysis suggested that there are various factors influencing recanalization risk. Once certain values of morphometric parameters are exceeded, a complex flow with numerous vortices occurs. This phenomenon was revealed due to CFD investigations that validated our statistical research. Thus, the complex flow pattern itself can be treated as a relevant recanalization predictor.

The association between hemodynamics and wall characteristics in human intracranial aneurysms: a review

Hemodynamics plays a key role in the natural history of intracranial aneurysms (IAs). However, studies exploring the association between aneurysmal hemodynamics and the biological and mechanical characteristics of the IA wall in humans are sparse. In this review, we survey the current body of literature, summarize the studies' methodologies and findings, and assess the degree of consensus among them. We used PubMed to perform a systematic review of studies that explored the association between hemodynamics and human IA wall features using different sources. We identified 28 publications characterizing aneurysmal flow and the IA wall: 4 using resected tissues, 17 using intraoperative images, and 7 using vessel wall magnetic resonance imaging (MRI). Based on correlation to IA tissue, higher flow conditions, such as high wall shear stress (WSS) with complex pattern and elevated pressure, were associated with degenerated walls and collagens with unphysiological orientation and faster synthesis. MRI studies strongly supported that low flow, characterized by low WSS and high blood residence time, was associated with thicker walls and post-contrast enhancement. While significant discrepancies were found among those utilized intraoperative images, they generally supported that thicker walls coexist at regions with prolonged residence time and that thinner regions are mainly exposed to higher pressure with complex WSS patterns. The current body of literature supports a theory of two general hemodynamic-biologic mechanisms for IA development. One, where low flow conditions are associated with thickening and atherosclerotic-like remodeling, and the other where high and impinging flow conditions are related to wall degeneration, thinning, and collagen remodeling.

Elevated Lipid Infiltration Is Associated With Cerebral Aneurysm Rupture

Background: Intracranial aneurysm wall degradation can be associated with lipid infiltration. However, the relationship between lipid infiltration and aneurysm rupture has not been explored quantitatively. To investigate the correlation between lipid infiltration and aneurysm rupture, we utilized patient-specific simulation of low-density lipoprotein (LDL) transport to analyze lipid infiltration in the cerebral aneurysm wall.

Methods: Sixty-two aneurysms were analyzed. Patient blood pressure, plasma LDL concentration, and three-dimensional angiographic images were obtained to simulate LDL transport in aneurysms. Morphological, hemodynamic, and lipid accumulation parameters were compared between ruptures and unruptured groups. Multivariate logistic regression was also performed to determine parameters that are independently associated with rupture.

Results: Size ratio, wall shear stress, low shear area, relative residence time, area-averaged LDL infiltration rate, and maximum LDL infiltration rate were significant parameters in univariate analysis (P < 0.05). Multivariate analysis revealed that only average LDL infiltration remained as a significant variable (P < 0.05). The prediction model derived showed good performance for rupture prediction (AUC, 0.885; 95% CI, 0.794–0.976).

Conclusions: Ruptured aneurysms showed significantly higher LDL infiltration compared to unruptured ones. Our results suggested that lipid infiltration may promote aneurysm rupture. Lipid infiltration characteristics should be considered when assessing aneurysm rupture risk.

4D-CT angiography versus 3D-rotational angiography as the imaging modality for computational fluid dynamics of cerebral aneurysms

Background and purpose

Computational fluid dynamics (CFD) can provide valuable information regarding intracranial hemodynamics. Patient-specific models can be segmented from various imaging modalities, which may influence the geometric output and thus hemodynamic results. This study aims to compare CFD results from aneurysm models segmented from three-dimensional rotational angiography (3D-RA) versus novel four-dimensional CT angiography (4D-CTA).

Methods

Fourteen patients with 16 cerebral aneurysms underwent novel 4D-CTA followed by 3D-RA. Endoluminal geometries were segmented from each modality using an identical workflow, blinded to the other modality, to produce 28 'original' models. Each was then minimally edited a second time to match length of branches, producing 28 additional 'matched' models. CFD simulations were performed using estimated flow rates for 'original' models (representing real-world experience) and patient-specific flow rates from 4D-CTA for 'matched' models (to control for influence of modality alone).

Results

Overall, geometric and hemodynamic results were consistent between models segmented from 3D-RA and 4D-CTA, with correlations improving after matching to control for operator-introduced variability. Despite smaller 4D-CTA parent artery diameters (3.49±0.97 mm vs 3.78±0.92 mm for 3D-RA; p=0.005) and sac volumes (157 (37–750 mm3) vs 173 (53–770 mm3) for 3D-RA; p=0.0002), sac averages of time-averaged wall shear stress (TAWSS), oscillatory shear (OSI), and high frequency fluctuations (measured by spectral power index, SPI) were well correlated between 3D-RA and 4D-CTA 'matched' control models (TAWSS, R2=0.91; OSI, R2=0.79; SPI, R2=0.90).

Conclusions

Our study shows that CFD performed using 4D-CTA models produces reliable geometric and hemodynamic information in the intracranial circulation. 4D-CTA may be considered as a follow-up imaging tool for hemodynamic assessment of cerebral aneurysms.

Hemodynamic characteristics associated with thinner regions of intracranial aneurysm wall

Aneurysm wall thickness is an important determinant of aneurysm progression and intra-procedural rupture. Several previous studies have evaluated the association between hemodynamic stress and aneurysm wall thickness, but conflicting results were obtained and no consensus has been achieved. According to the intraoperative findings, twenty-eight unruptured middle cerebral artery (MCA) aneurysms presented with thin-walled regions were enrolled in our study. Patient-specific 3D aneurysm models were constructed from preoperative computed tomography angiography (CTA) data and computational fluid dynamics (CFD) analyses were performed under pulsatile-flow conditions. Thin-walled regions of aneurysm dome were recognized by two experienced reviewers based on the intraoperative microscopy findings. Hemodynamic parameters derived from CFD analysis, including normalized wall shear stress (NWSS), normalized pressure (NP), the oscillatory shear index (OSI) and relative residence time (RRT), were compared between thin-walled regions and surrounding normal-thickness areas. Of the included aneurysms, twenty-eight pairs of thin-walled and normal surrounding regions were determined. Compared with surrounding tissues, thin-walled regions of aneurysm wall tended to present with higher pressure (1.232 vs 1.043, p < 0.05) and lower wall shear stress (0.693 vs 0.868, p < 0.05). Multivariate analysis revealed that elevated NP was significantly associated with thinning of the local aneurysm wall. Higher pressure and lower WSS were characteristic hemodynamic features associated with thinner regions of the aneurysm wall, elevated NP was an independent risk factor for local aneurysm wall thinning. CFD seems to be a useful method to estimate the location of thin-walled region, which will be helpful in reducing the risk of intraoperative rupture.