Article Text

Download PDFPDF

Original research
Hemodynamic conditions that favor bleb formation in cerebral aneurysms
  1. Seyedeh Fatemeh Salimi Ashkezari1,
  2. Fernando Mut1,
  3. Bong Jae Chung2,
  4. Anne M Robertson3,4,
  5. Juan R Cebral1,5
  1. 1 Department of Bioengineering, George Mason University, Fairfax, Virginia, USA
  2. 2 Department of Mathematical Sciences, Montclair State University, Montclair, New Jersey, USA
  3. 3 Department of Mechanical Engineering and Material Science, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
  4. 4 Department of Bioengineering, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
  5. 5 Department of Mechanical Engineering, George Mason University, Fairfax, Virginia, USA
  1. Correspondence to Ms Seyedeh Fatemeh Salimi Ashkezari, Department of Bioengineering, George Mason University, Fairfax, VA 22030, USA; ssalimia{at}gmu.edu

Abstract

Background Although it is generally believed that blebs represent weaker spots in the walls of intracranial aneurysms (IAs), it is largely unknown which aneurysm characteristics favor their development.

Objective To investigate possible associations between aneurysm hemodynamic and geometric characteristics and the development of blebs in intracranial aneurysms.

Methods A total of 270 IAs in 199 patients selected for surgical clipping were studied. Blebs were visually identified and interactively marked on patient-specific vascular models constructed from presurgical images. Blebs were then deleted from the vascular reconstruction to approximate the aneurysm before bleb formation. Computational fluid dynamics studies were performed in these models and in cases without blebs. Hemodynamic and geometric characteristics of aneurysms with and without blebs were compared.

Results A total of 173 aneurysms had no blebs, while 97 aneurysms had a total of 122 blebs. Aneurysms favoring bleb formation had stronger (p<0.0001) and more concentrated inflow jets (p<0.0001), higher flow velocity (p=0.0061), more complex (p<0.0001) and unstable (p=0.0157) flow patterns, larger maximum wall shear stress (WSS; p<0.0001), more concentrated (p=0.0005) and oscillatory (p=0.0004) WSS distribution, and a more heterogeneous WSS field (p<0.0001), than aneurysms without blebs. They were also larger (p<0.0001), more elongated (p<0.0001), had wider necks (p=0.0002), and more distorted and irregular shapes (p<0.0001).

Conclusions Strong and concentrated inflow jets, high-speed, complex, and unstable flow patterns, and concentrated, oscillatory, and heterogeneous WSS patterns favor the formation of blebs in IAs. Blebs are more likely to form in large, elongated, and irregularly shaped aneurysms. These adverse characteristics could be considered signs of aneurysm instability when evaluating aneurysms for conservative observation or treatment.

  • aneurysm
  • hemorrhage

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Introduction

The presence of blebs or secondary focalized bulging of the wall has been included as a risk factor in scoring scales used to assess the risk of rupture of intracranial aneurysms (IAs),1 2 and blebs have been shown to be risk indicators independent of aneurysm size and location.3 Although, it is generally believed that blebs represent weaker spots in the aneurysm wall that are vulnerable to failure, it is largely unknown which patient and aneurysm characteristics favor the development of blebs. Identifying and understanding these characteristics, which could be considered biomarkers for aneurysm instability, is extremely valuable for improving aneurysm evaluation and treatment recommendations.

The progressive weakening of the aneurysm wall is thought to be strongly related to flow-induced degenerative remodeling of the wall through inflammatory processes,4 5 among other multifactorial processes.6 Many studies have investigated associations between abnormal flow conditions and aneurysm rupture and instability.7 Similar associations with aneurysm size and shape characteristics have also been analyzed,8 since geometric characteristics are thought to influence the hemodynamic environment and also the biomechanical stresses within the aneurysm wall, and many of them are relatively simple to measure or calculate from medical images. Therefore, it is logical to investigate possible associations between aneurysm hemodynamic and geometric characteristics and the development of blebs in IAs, which is the objective of this study.

Methods

Aneurysm data and bleb identification

A total of 270 IAs in 199 patients selected for surgical clipping and previously used for studying associations of patient and clinical characteristics with the presence of blebs were further analyzed.9 Deidentified aneurysm vascular geometries were obtained from our database.

The overall methodology for identifying the bleb region in the aneurysm sac is illustrated in figure 1A–D and online supplementary figure 1. Using previously described methods,10 patient-specific vascular models were constructed from three-dimensional rotational angiography or computed tomographic angiography images acquired prior to surgery. In these aneurysms, blebs were identified by visual inspection of corresponding volume rendered 3D images (figure 1A), reconstructed vascular models (figure 1B), and aneurysm Gaussian curvature maps (figure 1C), and interactively marked on the vascular models (figure 1D), using a tool called ChePen3D, as detailed by Salimi Ashkezari et al.9 An illustration of the use of the curvature map for guidance in the definition and marking of blebs is presented in online supplementary figure 1.

Figure 1

Vascular modeling: (A) volume rendering of three-dimensional rotational angiography image, (B) patient-specific vascular model, (C) curvature map highlighting bleb (red=positive Gaussian curvature, blue=negative Gaussian curvature), (D) marking of bleb using ChePen3D tool, (E) bleb removal, (F) closure of bleb orifice, (G) meshing. Hemodynamics visualization at peak systole after bleb removal (surrogate for aneurysm prior to bleb formation): (H) inflow jet, (I) wall shear stress magnitude, (J) oscillatory shear index, (K) flow streamlines, (L) vortex corelines, (M) swirling flow around vortex corelines.

Hemodynamic modeling

In order to approximate the hemodynamic conditions and geometry prevalent prior to bleb formation, new vascular models were automatically created by virtually deleting the blebs in aneurysms with identified blebs as shown in figure 1E–G. Briefly, the mesh elements labeled as belonging to the blebs were deleted from the vascular model triangulation (figure 1E), and the hole corresponding to the bleb neck was automatically closed with a surface triangulation algorithm that preserves surface curvature at the edges11 (figure 1F). Finally, a volumetric grid of tetrahedral elements was generated inside this new vascular model for subsequent computational fluid dynamics (CFD) analysis (figure 1G). In aneurysms without blebs, the CFD meshes were generated directly from the original vascular model. Mesh resolution was set to a minimum of 0.02 cm, resulting in grids ranging from approximately 2 to 7 million elements.

Hemodynamic simulations were then performed on both sets of aneurysms: aneurysms with their blebs removed, and aneurysms without blebs. In these simulations, blood was approximated as a Newtonian incompressible fluid, and the unsteady Navier-Stokes equations were numerically solved with a finite element code.12 Vascular walls were approximated as rigid. Pulsatile inflow conditions were imposed by scaling a representative flow waveform with an empirical law relating flow rate and cross-sectional area in internal carotid and vertebral arteries.13 Outflow conditions were imposed by splitting flows consistent with Murray’s law. Simulations were run for two cardiac cycles, and data from the second cycle were used to characterize the hemodynamic conditions in the aneurysm. Example flow visualizations are shown in figure 1H–M.

Data analysis

The aneurysm hemodynamics was characterized by computing flow variables that measure different properties of the intra-aneurysmal flow,14 including: strength (Q) and concentration (ICI) of the inflow jet, speed (VE), viscous dissipation (VD), complexity (corelen) and stability (podent) of the flow pattern, and strength (WSSmean, WSSmax), concentration (SCI) and spread (LSA), and oscillation (OSImean, OSImax) of the wall shear stress (WSS). The mean and maximum WSS normalized with the mean WSS value over the near (<1 cm away from the neck) parent artery were also computed as measures of the strength of the WSS relative to the WSS in the parent artery (see online supplementary table 1).

In addition, critical points in the WSS vector field were detected, counted, and averaged over the cardiac cycle. A critical point is mathematically defined as a location on the aneurysm sac where the shear stress vector is zero (stationary point). Different kinds of critical points, such as saddle points, node sources or sinks, focus sources or sinks, are distinguished by the local derivatives around the stationary point.15 Critical points are valuable for characterizing the local WSS field, which can presumably affect the endothelial function. For example, the point on the aneurysm surface associated with the end of a vortex coreline is a critical point. The mean number of critical points (nCrPoints) is thus a variable that measures the total number of regions with 'insult' to the endothelium by the abnormal WSS field during the cardiac cycle.

Similarly, the geometry of the aneurysm was characterized by computing variables that capture different aspects of the aneurysm geometry,16 including aneurysm sac size (Asize, SR, GAA), aneurysm neck size (Nsize), aneurysm elongation (AR, VOR, BF), and aneurysm shape distortion relative to a spherical shape (NSI, CR) and surface irregularity (UI). More details are presented in online supplementary table 1.

The hemodynamic and geometric characteristics of aneurysms that developed blebs (with their blebs removed to approximate their characteristics before bleb formation) were compared with those of aneurysms without blebs. Specifically, the mean values of variables were statistically compared between these two groups using the two-sample Mann-Whitney U test coded into R scripts. The p values were then adjusted for multiple testing using the Benjamini and Hochberg method available in R. Differences were considered statistically significant for p values <0.05 after adjustment.

Results

A total of 173 aneurysms in 121 patients were classified as having no blebs. Blebs were identified in the remaining 97 aneurysms in 78 patients, for a total of 122 blebs. Comparisons of the patient and clinical characteristics between aneurysms with and without blebs are reported in a separate paper.9 Hemodynamic and geometric differences between aneurysms with blebs (removed) and aneurysms without blebs are presented in table 1. As shown in this table, most of the hemodynamic and geometric variables were significantly different between the two groups of aneurysms, even after adjustment for multiple testing.

Table 1

Hemodynamic and geometric characteristics of aneurysms with blebs (after bleb removal, mimicking conditions before bleb formation) and aneurysms without blebs. See online supplementary table 1 for more details on the hemodynamic and geometric variables

Hemodynamically, aneurysms with blebs (with their blebs deleted) had stronger (Q, p<0.0001) and more concentrated inflow jets (ICI, p<0.0001), higher intra-aneurysmal flow velocity (VE, p=0.0061), more complex (corelen, p<0.0001) and unstable (podent, p=0.0157) flow patterns than aneurysms without blebs. In addition, they also had larger maximum wall shear stress (WSSmax, p<0.0001; MaxWSSnorm, p=0.0005), more concentrated (SCI, p=0.0005) and oscillatory (OSImax, p=0.0004) WSS distribution, and a larger number of critical points of the wall shear stress field (nCrPoints, p<0.0001).

Geometrically, aneurysms that developed blebs were larger (Asize, p<0.0001; SR, p<0.0001; GAA, p<0.0001), more elongated (AR, p<0.0001; VOR, p<0.0001; BF, p<0.0001), had wider necks (Nsize, p=0.0002), and more shape distortion (CR, p=0.0016, NSI, p<0.0001) and more surface irregularity (UI, p=0.0016) than aneurysms without blebs.

Discussion

Aneurysm shape irregularity and, in particular, the presence of blebs has been considered a risk factor for rupture of cerebral aneurysms.1–3 However, little is known about the conditions that predispose aneurysms to develop blebs and place them at immediate risk of rupture. As such, the focus of this work was to investigate what flow conditions and aneurysm geometries could favor bleb formation, which would subsequently predispose the wall for failure and aneurysm rupture. For this purpose, a series of aneurysms with and without blebs were compared. The geometry and hemodynamics of the aneurysms with blebs were analyzed by first virtually removing the bleb from the vascular model in order to approximate the conditions prior to bleb formation, which were then compared with those of aneurysms without blebs. In this manner, the comparison is essentially between those aneurysms that are prone to develop blebs and those that are unlikely to develop blebs. There are two underlying assumptions in this approach. First, is the assumption that the aneurysm shape and hemodynamics before the blebs were formed are reasonably approximated by deleting the blebs from vascular models reconstructed from 3D images obtained after the blebs were formed. This, in turn, presupposes that blebs develop and grow more rapidly than the rest of the aneurysm, which is not unreasonable since longitudinal studies have suggested two modes of growth, one associated with overall enlargement of the aneurysm and another associated with focalized growth and bleb development.17 However, there may be other modes of aneurysm growth and bleb development—for example, early bleb formation followed by aneurysm enlargement, or progressive growth of the aneurysm and simultaneous bleb development. Although it is unknown how frequently these growth modes occur in the aneurysm population, previous studies have shown in a few cases rapid development of blebs without major enlargement of the rest of the aneurysm sac in longitudinal images,18 19 suggesting that our assumption is reasonable. The second underlying assumption is that aneurysms without blebs at the time of imaging have had a lower likelihood of developing blebs than those that have already developed a bleb. This conjecture seems reasonable, but it is not known if some of these aneurysms would develop blebs in the future. In spite of these limitations, the findings of this study are intriguing as further discussed below.

The methodology used for defining blebs in cerebral aneurysms involves the manual step of visually identifying the blebs and interactively marking them on the vascular models, with guidance from volume-rendered 3D angiographic images and aneurysm surface curvature maps. It has previously been shown that this approach is quite consistent between different operators.9 Furthermore, all the subsequent steps (bleb removal, hole triangulation, meshing, flow simulation, and hemodynamic and geometric characterization) are fully automatic, thus making the entire process quite objective and reproducible.

Our results suggest that in general aneurysm flow conditions characterized by high and concentrated inflow jets, high speed, complex and unstable flow patterns with an associated concentrated, complex, and oscillatory wall shear stress distribution seem more favorable for the development of blebs in intracranial aneurysms. Furthermore, blebs seem to be in general more likely to develop in larger, more elongated, and more irregular aneurysms. To illustrate these differences, examples of flow visualizations in two aneurysms at the middle cerebral artery (MCA) bifurcation are presented in figure 2, one without a bleb (left) and another with a bleb (right, in this case flow simulation corresponds to the aneurysm prior to bleb formation—that is, with the bleb virtually deleted). Further examples are presented in online supplementary figure 2–5 for aneurysms at the posterior communicating artery (PCOM), anterior communicating artery (ACOM), M1 segment of the MCA, and MCA bifurcation, respectively. These visualizations show, in general, stronger inflow jets and impingement regions, higher and more heterogeneous wall shear stress, and more complex flow patterns in aneurysms that developed blebs than in those without blebs.

Figure 2

Examples of hemodynamics (at peak systole) for middle cerebral artery (MCA) aneurysms with and without blebs. Left panel – aneurysm with bleb removed (surrogate for IA prior to bleb formation): (A) geometry of MCA aneurysm after removal of one bleb (insert shows marked bleb in red), (B) inflow jet (iso-velocity surface), (C) wall shear stress (WSS) magnitude, (D) flow pattern (streamlines), (E) vortex corelines. Right panel – representative MCA aneurysm that presented without bleb: (F) geometry of MCA aneurysm without bleb, (G) inflow jet, (H) WSS magnitude, (I) flow pattern, (J) vortex corelines.

Since treatment decisions are more difficult for small aneurysms, the analysis was repeated for aneurysms smaller than 7 mm. The results are presented in online supplementary table 2. It can be seen that most differences remained significant for this subset. Interestingly, viscous dissipation became significantly different between small aneurysms with and without blebs, whereas measures of flow instability (podent) and oscillation (OSI) became not significantly different. Similarly, neck size and measures of aneurysm shape distortion (NSI, CR, UI) were no longer significantly different in this subset. This suggests that flow instabilities and shape irregularity are less important factors for bleb formation in small aneurysms than measures of flow strength and WSS heterogeneity.

Similar hemodynamic and geometric characteristics have been previously associated with aneurysm rupture,20 as well as aneurysm growth and instability,21 and have been used in aneurysm rupture probability models.22 However, other studies have suggested connections for low and oscillatory WSS with aneurysm rupture.23 24 In those studies, the aneurysm WSS was normalized with the average WSS in the parent artery. In the current study, the normalized mean WSS tended to be lower in aneurysms that developed blebs than in aneurysms without blebs, but this difference did not reach statistical significance. In contrast, the maximum normalized WSS was significantly larger. Note that these differences cannot in general be explained by the difference in aneurysm size. In addition, high flow conditions and irregular aneurysm shapes have also been related to increased aneurysm wall inflammation25 and decreased aneurysm wall strength.26 Bleb formation seems to be a highly localized phenomenon, and the flow conditions identified here as favorable for bleb formation are characterized by strong and concentrated flow characteristics (inflow jet, WSS distribution, critical points). Therefore, it seems reasonable to propose that these adverse flow and geometric characteristics may be signs of aneurysm instability and could be used to better understand which aneurysms are more likely to develop blebs and progress to a more fragile and rupture-prone status. Moreover, once the blebs are formed, they become high-risk features that could increase local wall stress27 and permeability, leading to inflammation, which has been linked to wall enhancement in MR vessel wall imaging.28

Conclusions

Hemodynamic conditions characterized by strong and concentrated inflow jets, high speed, complex and unstable flow patterns, and concentrated, oscillatory, and heterogeneous wall shear stress patterns favor the formation of blebs in cerebral aneurysms. Similarly, blebs are more likely to form in large, elongated, and irregularly shaped aneurysms. These adverse characteristics could be considered signs of vulnerability to aneurysm instability when evaluating aneurysms for conservative observation or treatment.

References

Footnotes

  • Contributors SFSA, AMR, and JRC designed the study. FM contributed to development of the methodology. FM and JRC designed the software tools. BJC simulated vascular reconstructions. SFSA curated the data. SFSA and JRC performed the data analysis. SFSA, FM, AMR, and JRC contributed to interpretation of the results. AMR and JRC aquired funding, supervised students and coordinated the project.SFSA and JRC drafted the manuscript. All authors contributed to manuscript edition and approved the final manuscript.

  • Funding This work was supported by NIH grant R01NS097457.

  • Competing interests None declared.

  • Patient consent for publication Not required.

  • Ethics approval The protocols for patient consent, handling of patient data and analysis were approved by the institutional review board (IRB) at the University of Pittsburgh (Protocol # STUDY20020015), University of Illinois at Chicago (Protocol # 2015-0322), Allegheny General Hospital (Protocol # RC-5141), and Helsinki University Hospital. The whole study's IRB is overseen by the University of Pittsburgh's IRB.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data availability statement Data are available upon reasonable request. The data that support the findings of this study are available from the corresponding author, upon request.

Linked Articles

  • Hemorrhagic Stroke
    Seyedeh Fatemeh Salimi Ashkezari Felicitas J Detmer Fernando Mut Bong Jae Chung Alexander K Yu Christopher J Stapleton Alfred P See Sepideh Amin-Hanjani Fady T Charbel Behnam Rezai Jahromi Mika Niemelä Juhana Frösen Ji Zhou Spandan Maiti Anne M Robertson Juan R Cebral