Skip to main content

Wall shear stress gradient is independently associated with middle cerebral artery aneurysm development: a case-control CFD patient-specific study based on 77 patients

Abstract

Background

Previously published computational fluid dynamics (CFD) studies regarding intracranial aneurysm (IA) formation present conflicting results. Our study analysed the involvement of the combination of high wall shear stress (WSS) and a positive WSS gradient (WSSG) in IA formation.

Methods

We designed a case-control study with a selection of 38 patients with an unruptured middle cerebral artery (MCA) aneurysm and 39 non-aneurysmal controls to determine the involvement of WSS, oscillatory shear index (OSI), the WSSG and its absolute value (absWSSG) in aneurysm formation based on patient-specific CFD simulations using velocity profiles obtained from transcranial colour-coded sonography.

Results

Among the analysed parameters, only the WSSG had significantly higher values compared to the controls (11.05 vs − 14.76 [Pa/mm], P = 0.020). The WSS, absWSSG and OSI values were not significantly different between the analysed groups. Logistic regression analysis identified WSS and WSSG as significant co-predictors for MCA aneurysm formation, but only the WSSG turned out to be a significant independent prognosticator (OR: 1.009; 95% CI: 1.001–1.017; P = 0.025). Significantly more patients (23/38) in the case group had haemodynamic regions of high WSS combined with a positive WSSG near the bifurcation apex, while in the control group, high WSS was usually accompanied by a negative WSSG (14/39). From the analysis of the ROC curve for WSSG, the area under the curve (AUC) was 0.654, with the optimal cut-off value −0.37 Pa/mm. The largest AUC was recognised for combined WSS and WSSG (AUC = 0.671). Our data confirmed that aneurysms tend to form near the bifurcation apices in regions of high WSS values accompanied by positive WSSG.

Conclusions

The development of IAs is determined by an independent effect of haemodynamic factors. High WSS impacts MCA aneurysm formation, while a positive WSSG mainly promotes this process.

Peer Review reports

Background

Intracranial aneurysms (IAs) are defined as pathological dilations of intracranial arteries caused by a weakness within the arterial wall. They are a relatively common disease with a worldwide prevalence of approximately 3% [1]. While often asymptomatic, cerebral aneurysms may rupture, causing a life-threatening condition – subarachnoid haemorrhage (SAH). The vast majority of aneurysms are located in the anterior part of the circle of Willis (the anterior communicating artery, ACoA 18–24%; the internal carotid artery, ICA 37–42%) and in the middle cerebral artery (MCA, 27–35%) [1,2,3]. As IAs tend to develop near certain bifurcations and along vessel curvatures it is hypothesised that the local haemodynamic environment may play a key role in IA formation [4, 5].

Although the morphometric and demographic risk factors responsible for aneurysm growth and rupture were thoroughly tested, researched and established [6,7,8], the exact pathophysiology of aneurysm formation initiation, including haemodynamic factors, remains unclear. Previously published data have shown that the vascular segments mentioned above are affected by complex haemodynamic forces, such as: wall shear stress (WSS), the wall shear stress gradient (WSSG), and temporal fluctuations in WSS [4, 5, 9,10,11,12,13,14,15,16].

Histologic findings revealed that both high and low WSS values may be responsible for initiating the formation of IAs. High WSS induces wall remodelling by the migration of smooth muscle cells, the secretion of inflammatory mediators and endothelial injury [17,18,19]. Low WSS may lead to localised degeneration of the arterial wall by favouring disorganisation and apoptosis of the endothelium, as well as oxidative stress and additional inflammation [20,21,22,23]. Additionally, a positive WSSG downregulates expression of several genes inhibiting proliferation, promoting inflammation and increasing secretion of ADAMTS1 protease [24, 25]. The above changes may result in increased susceptibility of the arterial wall to mechanical damage caused by local blood flow.

A significant number of image-based computational fluid dynamics (CFD) studies connecting specific haemodynamic parameters with intracranial aneurysm growth and rupture have been published in recent years but they present conflicting results, correlating both high and low WSS values with aneurysm growth and rupture or denying the importance of previously analysed parameters, introducing new ones in their place [12, 13]. The above complications may arise from conducting research on relatively small sample groups of patients or using the same averaged velocity profile in every CFD simulation [14]. Therefore, against this background, we designed a case-control study with a selection of patients with an MCA aneurysm and non-aneurysmal controls to determine the involvement of WSS, the WSSG and the oscillatory shear index (OSI) in aneurysm formation based on patient-specific CFD modelling. To the best of our knowledge, this is one of the first CFD studies in which numerical simulations were carried out based on patient-specific velocity profiles obtained from a transcranial Doppler ultrasound. This study is part of a wide range of research studies aimed at confirming the results of clinical observations, experimental research on animals and numerical analyses explaining haemodynamic involvement, mainly a combination of high WSS values and positive WSSG, in aneurysm formation.

Methods

Data source

The inclusion criteria for the case group were as follows: male or female adults, aged 18–75 years, and diagnosed with an unruptured MCA aneurysm. Between June 2013 and June 2017, 38 patients with unruptured MCA aneurysms diagnosed using three-dimensional computed tomography angiography (3D CTA) were enrolled in the case group. The case group consisted of seven males and 31 females ranging from 37 to 75 years old, with a mean age of 56 years. All patients had saccular and unruptured aneurysm located on the MCA bifurcation. The control group consisted of 39 cases of patients with no cerebrovascular diseases, 19 males and 20 females with a mean age of 50 years (ranging from 20 to 72). Every patient in the control group had a CTA done in order to exclude a suspected intracranial aneurysm or to establish the aetiology of minor symptoms, such as headache or vertigo.

Patients suffering from other vascular malformations that could potentially affect cerebral blood flow, haemodynamically significant carotid artery stenosis or connective tissue disorders (e.g. Marfan syndrome, Loeys-Dietz syndrome and Ehlers-Danlos syndrome types II and IV) were excluded from this study. Additional exclusion criteria from the study were: age under 18 years or over 75 years, presence of multiple cerebral aneurysms (other than mirror MCA aneurysms) or other than aneurysm pathologies in the central nervous system that could have a potential effect on cerebral blood flow, severe systemic disorders, severe heart or multi-organ failures, pregnancy and family history of cerebral aneurysm or genetically determined conditions associated with an increased risk of intracranial aneurysm formation (e.g. autosomal dominant polycystic kidney disease, neurofibromatosis type I, multiple endocrine neoplasia type I, pseudoxanthoma elasticum, hereditary haemorrhagic telangiectasia).

The study’s protocol was approved by the Institutional Review Board at the Medical University of Silesia in Katowice, Poland and all procedures were carried out in accordance with the relevant guidelines and regulations. Each patient was informed about the purpose and course of the research and asked to give their informed consent to participate in the project.

CTA and transcranial colour-coded sonography (TCCS) protocols

All CT studies were performed using a 64-row 128-slice CT system (GE Optima CT 660, GE Healthcare, USA) using the parameters previously described [26]. On average, 50 ml of ionic contrast material (Iomeron 350, Bracco Imaging Deutschland GmbH, Konstanz, Germany) was injected into the basilic vein at 4–4.5 ml/sec via a selectively positioned catheter with the aid of an automatic syringe (OptiVantage DH, Mallincrodt, St Louis, MO, USA).

All TCCS examinations were performed by the same researcher (WK) using a Vivid 3 Pro (GE Healthcare, Chicago, Illinois, USA) equipped with a multi-frequency transcranial probe (1.5–3.6 MHz) through the temporal acoustic window, in accordance with the previously described standards [27]. Anterior cerebral circulation was imaged through the temporal acoustic window with the subject in a supine position. The angle-corrected mean blood-flow velocity (Vm), peak systolic velocity (Vps) and end-diastolic velocity (Ved) were measured for both MCAs (Fig. 1). The spectral Doppler mode was used to graphically display the spectrum of flow velocities that was exported to a graphic file (i.e. BMP format file) and then used in the numerical modelling steps. The groups were analysed in terms of their applied velocity profiles and no statistical differences were found between the case and control groups (Vps: 97.4 vs 99.1 [cm/s], P = 0.665; Ved: 44.0 vs 43.9 [cm/s], P = 0.586; Vm: 66.4 vs 64.8 [cm/s], P = 0.554; respectively).

Fig. 1
figure1

View of a colour-coded image of the middle cerebral artery (MCA) with a corresponding Doppler spectral analysis performed using a transtemporal insonation. The sample volume was placed at a depth of 54 mm in the distal portion of the MCA. ICA – the internal carotid artery; PCA – the posterior cerebral artery; A1 ACA – A1 segment of the anterior cerebral artery

Model construction

The patient-specific models were generated from CTA images using Mimics v16.0 and 3-matic v8.0 (Materialise, Leuven, Belgium). Any touching vessels and lesser branches were removed so the models only consisted of a trunk a bifurcation and the main branches arising from it. Aneurysms were manually removed, and the pre-aneurysmal geometries were recreated. Afterwards, the models were smoothed. Inlet branches were clipped at a manually selected location at the start of the artery perpendicular to its centreline to allow development of the flow. The outlets were clipped as close to the ends of the model as possible. By doing so, we eliminated model inaccuracies that might have affected the final results. Cross-sectional cuts before and after aneurysm removal were inspected to verify if the reconstructed artery was as intended.

The 3D models were pre-discretized and optimized in 3-matic software and subsequently imported as *.igs files to ANSYS Workbench (ANSYS Inc., Canonsburg, PA, USA) and computational meshes were generated using the ANSYS meshing module. The mesh was composed of tetrahedral elements with a fixed side length of 0.3 mm. To ensure an accurate definition of the velocity gradient, three inflation layers were included near the wall. The inflation growth rate was set to 1.2 and the maximum thickness was 0.15 mm. On average, computational meshes in the case group consisted of 810,000 tetrahedral prismatic elements and 740,000 in the control group.

Numerical modelling

Navier-Stokes equations describing blood flow were solved by performing CFD simulations using an ANSYS Fluent, a finite volume-based CFD solver. Blood was considered to be an incompressible Newtonian fluid with a density of 1056 kg/m3 and a constant viscosity of 0.0035 Pa*s. Arterial walls were assumed to be rigid and a no-slip boundary condition was applied at them. Boundary conditions describing patient-specific inlet velocity profiles were based on the MCA blood-flow waveforms measured with TCCS and analysed by a custom-written script (see link: https://github.com/kudlacz964/tccs2ansys). The programme transcribed the BMP graphic files acquired from the TCCS into numerical files readable by the ANSYS computing environment for each patient. The Python Imaging Library used in the script allowed for highly accurate readings. An initial value of 0 was used for the pressure at both outlets.

The length of the time step was 0.001 s and the duration of the simulation was dependent on the duration of the patient’s cardiac cycle. To establish convergence in every time step, two cardiac cycles were calculated and only the data obtained from the second cycle was analysed.

Haemodynamic parameters

The following parameters were analysed within a defined circular patch with a radius of 5 mm and a centre in the bifurcation apex, which was consistent among the subjects and encapsulated the aneurysm initiation sites in the dataset.

WSS is a tangential, frictional stress exerted on the vessel wall by the flow of the blood. During the research, the time-averaged WSS was analysed, which was calculated by averaging the WSS vector over the cardiac cycle. In this study, a high WSS was defined as time-averaged WSS values in the 98th percentile for each case.

$$ WSS=\frac{1}{T}{\int}_0^T\left|{wss}_i\right| dt $$

where wssi is the instantaneous WSS vector and T is the duration of the cardiac cycle.

The WSSG was calculated by taking the spatial derivative of the time-averaged WSS. Apart from the WSSG, additional WSSG-derived parameters were also distinguished: absolute WSSG (absWSSG) and WSSG direction (dirWSSG). The WSSG was defined as the WSSG averaged over the high WSS area. The WSSG measurements were taken from one branch – in the case group, from the branch developing aneurysm or from the branch with the higher WSS values when aneurysm was big enough to cover both branches or the bifurcation apex. In the control group, the WSSG measurements were taken from branch with the higher WSS values.

The dirWSSG was the categorical variable: positive (dirWSSG+) and negative dirWSSG (dirWSSG-). This variable was positive when the WSSG vector direction was the same as the WSS vector originating from the same point, while for dirWSSG, it was the opposite.

$$ WSSG=\sqrt{{\left(\frac{\partial {\tau}_{w,p}}{\partial p}\right)}^2+{\left(\frac{\partial {\tau}_{w,q}}{\partial q}\right)}^2} $$

where τw is the WSS vector, the p-direction corresponds to the time-averaged direction of the WSS, and the q-direction is perpendicular to p.

OSI was introduced to measure changes in the direction of WSS during the cardiac cycle. Its values range from 0.0 in the case of unidirectional flow to 0.5 in regions of high and complete departure of the WSS vector from its original axis. The OSI was defined as the OSI averaged over the high WSS area.

$$ OSI=\frac{1}{2}\left(1-\frac{\left|{\int}_0^T{wss}_i dt\right|}{\int_0^T\left|{wss}_i\right| dt}\right) $$

Statistical analysis

Results were visualised using ParaView (Kitware Inc., Clifton Park, NY, USA) and Matplotlib and all data components were analysed using a Statistica v.13.3 package (StatSoft, Tulsa, OK, USA). The Shapiro-Wilk test was performed to assess the normal distribution of the investigated groups. Depending on the distribution, the parameters between groups were examined using the Student’s t-test or Mann–Whitney U test. To test if the proportion of cases characterised by a combination of high WSS and positive WSSG was significantly different between analysed groups, the Fisher’s exact test was used. P-values < 0.05 were considered statistically significant. For the sake of statistical analysis, dirWSS was encoded as follows: positive dirWSSG – 1, negative dirWSSG – 0. All haemodynamic parameters were subjected to logistic regression analysis with a stepwise addition mode. The potential risk factors for MCA aneurysm formation were identified on the univariate analysis. Intercorrelations between analysed parameters were examined using the Spearman’s rank correlation test. All variables were entered into a logistic regression model. Uncorrelated variables with p-values < 0.1 on univariate analysis were included in the multivariate logistic regression model to identify the independent predictors of MCA aneurysm formation. The results were presented as odds ratios (ORs) and their 95% confidence intervals (CIs). The independent predictors of an MCA aneurysm were subjected to receiver-operating characteristic (ROC) analysis to identify the area under the curve (AUC) values with optimal sensitivity and specificity.

Results

Haemodynamic parameters

As none of the haemodynamic parameters were normally distributed, they were analysed using the Mann–Whitney U test (Table 1). Patients from the case group had positive WSSG with significantly higher values (11.05 vs − 14.76 [Pa/mm], P = 0.020). The WSS, absWSSG and OSI values were not significantly different between the case and control groups (69.36 vs 76.59 [Pa], P = 0.158; 23.63 vs 35.07 [Pa/mm], P = 0.208; 0.175 vs 0.185, P = 0.759; respectively).

Table 1 Comparison of haemodynamic parameters in the analysed groups

Fig. 2 shows the distribution of WSS and WSSG in exemplificatory MCA bifurcation models from both groups. Areas coloured red represent those with high WSS values while areas with low WSS values are blue. It should be noted that those places with the highest WSS were not located at the bifurcation apex but in close proximity. There were usually two places of relatively high WSS – one located on each branch (Fig. 2a, b, d, e).

Fig. 2
figure2

Distribution of WSS (a and d); bifurcation with points representing positive (white) and negative (black) WSSG (b and e); WSS (red) and WSSG (blue) distribution along the wall (c and f), WSSG and WSS measurements taken along the dotted lines (a and d); representative MCA bifurcation models from case (upper row) and control (lower row) groups. WSS – wall shear stress; WSSG – wall shear stress gradient

Charts (Fig. 2c and f) present the characteristic haemodynamic patterns observed in both groups. The WSSG values near the bifurcation apex were positive and tended to decrease with distance, dropping below zero (Fig. 2c and f). Significantly more patients in the case group had haemodynamic regions of high WSS values combined with high and positive WSSG near the bifurcation apex (23/38), while in control group, high WSS values were usually accompanied by negative WSSG (14/39) (Table 2, P = 0.041).

Table 2 Comparison of regions of high WSS values combined with high and positive WSSG near the bifurcation apex in the analysed groups

In both groups, the regions of high WSS were located slightly more often on the larger branch (case group 55.3%, n = 38; control group 56.4%, n = 39). In the vast majority of cases (32/38), it was not possible to define from which branch the aneurysm had originated. In the remaining patients (6/38), the aneurysm locations corresponded with areas of high WSS and a positive WSSG, and overlapped the aneurysm ostia containing the initiation zone.

Predictors of MCA aneurysm – logistic regression and ROC analysis

The univariate logistic regression analysis identified two parameters as significant predictors for MCA aneurysm formation: WSSG (OR: 1.008; 95% CI: 1.0–1.016; P = 0.028) and dirWSSG (OR: 2.738; 95% CI: 1.088–6.889; P = 0.030) (Table 3). Intercorrelations between the analysed parameters were examined using the Spearman’s rank correlation test. For obvious reasons, the WSSG values were significantly correlated with dirWSSG (correlation coefficient: 0.865, P < 0.0001). Additionally, absWSSG was correlated with WSS (correlation coefficient: 0.488, P < 0.0001) and with OSI (correlation coefficient: -0.309, P = 0.006). Other parameters were not correlated.

Table 3 Risk factors of MCA aneurysm formation: univariate and multivariate logistic regression analysis

For the final multivariate logistic regression model, we included two variables – WSS and WSSG – identified as factors significantly influencing MCA aneurysms formation (OR: 0.986; 95% CI: 0.971–1.0; P = 0.049; OR: 1.009; 95% CI: 1.001–1.017; P = 0.025, respectively) (Table 3). The logistic regression formula with the corresponding coefficients for WSS and WSSG is presented below:

$$ \frac{1}{1+{\mathrm{e}}^{-\left({\upbeta}_0+{\upbeta}_1\mathrm{WSS}+{\upbeta}_2\mathrm{WSSG}\right)}} $$

where β0 = 1.133, β1 = − 0.014, β2 = 0.009, WSS – wall shear stress, WSSG – wall shear stress gradient.

The ROC curve for discovered risk factors for MCA aneurysm formation is presented in Fig. 3. From the analysis of the ROC curve for WSSG, the area under the curve (AUC) was 0.654, with the optimal cut-off value −0.37 Pa/mm (sensitivity 0.658 and specificity 0.641). The largest AUC value was recognised for combined WSS and WSSG (AUC = 0.671).

Fig. 3
figure3

ROC curves for all of the most significant risk factors of MCA aneurysm formation. WSS – wall shear stress; WSSG – wall shear stress gradient

Discussion

The main findings of this study are: (1) cerebral aneurysm development might be an independent effect of haemodynamic factors and (2) aneurysms form near the bifurcation apices in regions of high WSS values accompanied by positive WSSG.

Haemodynamic forces are considered to provide a regulatory function for both the physiologic and pathologic aspects of endothelial activity in the blood vessels. Endothelial cells, through molecular sensors such as integrins and mechanosensitive ion channels located on their surface, transform the WSS-induced mechanical signals into biological signals, activating molecular pathways to maintain vascular homeostasis [21, 22, 28, 29]. Endothelial dysfunction is known to contribute to intracranial aneurysm pathogenesis [11]. Data reported from animal experimentation showed that damage or injury to the endothelium is usually the first incident in intracranial aneurysm formation [30]. Haemodynamic stress activates key molecular pathways that result in the upregulation of adhesion molecules and chemotactic cytokines leading to inflammatory zone formation and compromising the endothelial cell-to-cell junctions [22, 29, 31, 32]. These processes are based on the overexpression of a large group of genes, such as ADAMTS1, PLAU, PLAT, and TIMP3, all of which are involved in extracellular matrix (ECM) degradation, induced by very high WSS [17]. In vivo studies have shown that high flow and WSS induce significant proliferation of endothelial cells (ECs) and smooth muscle cells (SMCs), accompanied by significant increases in matrix metallopeptidases (MMPs) [18, 19] responsible for the degradation and digestion of extracellular matrix components. The high level of MMPs degrades cell basement membranes and internal elastic lamina (IEL), allowing proliferating ECs to migrate into the new enlarged areas and allowing proliferating SMCs to reorientate in the media [18, 19] leading, consequently, to aneurysmal wall remodelling. Many CFD-studies have demonstrated a positive correlation with high WSS in bifurcation aneurysms [33,34,35,36,37], whereas a negative correlation was found in sidewall aneurysms [12, 13, 38, 39]. Our study showed a similar increase in WSS at MCA bifurcation in both studied groups of patients (control and aneurysmal), which suggests that high WSS alone cannot initiate IA formation. Univariate logistic regression analysis performed in our study revealed only two significant predictors for MCA aneurysm formation – the WSSG value and its direction (dirWSSG), but only the WSSG turned out to be a significant independent prognosticator. Therefore, in our opinion, high WSS impacts on MCA aneurysm formation, whereas a positive WSSG mainly promotes this process. Additionally, the results of our study confirmed the lack of a correlation between OSI and bifurcation aneurysm formation, which is consistent with the previously published data [12, 40].

In in vivo studies on dogs in which arterial bifurcations were surgically created from native common carotid arteries with increased blood flow (animals were on a high salt diet, and some had one of their renal arteries ligated), combined with CFD simulations performed on the in vivo images and results mapped onto histological images, have suggested the coexistence of increased WSS and positive high WSSG at the bifurcation as triggers for aneurysm initiation [41, 42]. However, little is known about the effect of high WSS and the WSSG on the formation of IAs in humans. The majority of image-based CFD studies investigating quantitative haemodynamic variables and their relationships demonstrated a positive correlation between WSS and/or WSSG and bifurcation aneurysm formation [34,35,36], while studies performed for sidewall aneurysms denied such a correlation [12, 39, 40]. Of note, all the above-mentioned studies were mainly case series reports conducted on a small group of patients, whereas our clinical-control study was based on a large homogeneous MCA bifurcation group of 77 patients, which increased the statistical power of the analyses performed.

The complex influence of several haemodynamic factors, such as high WSS with corresponding positive WSSG values, may induce degenerative changes in the arterial wall and result in aneurysm formation [41,42,43]. According to the work by Dolan et al. [25], a positive WSSG induces gene expression that could promote the proliferation of ECs, apoptosis, and extracellular processing. Specifically, a positive WSSG downregulated RPRM and BMP4 genes which are involved in cell cycle arrest and cell growth and differentiation, respectively [25]. In addition, a positive WSSG leads to downregulation of the genes that inhibit apoptosis (such as: MCP1, CSF2, BMP4, and THBS1) increasing the propensity of ECs to be apoptotic [25]. However, the key factor in the aneurysm formation process is the high WSS and positive WSSG induction of matrix remodelling. A positive WSSG increases significant ADAMTS1 expression in ECs (an ECM protease, which degrades ECM proteins including aggrecan, versican and thrombospondin), and downregulates TAGLN (a repressor of MMP-9 expression) further increasing ECM degradation [25]. All these processes make the arterial wall more vulnerable to mechanical damage [24, 25]. Interestingly, when taken together, a positive WSSG and WSS may augment the effect of high WSS on gene expression, as opposed to a negative WSSG that antagonises the effect of high WSS on gene expression [25].

Previous CFD studies were unable to identify strong risk factors or predictors for aneurysm formation among the known haemodynamic parameters (WSS, WSSG, OSI). This induced a strong need and motivation to find new haemodynamic parameters describing aneurysm initiation regions in a way previously unknown. Two new variables, derived from the previously described WSS and WSSG, were introduced specifically to study aneurysm formation: the aneurysm formation indicator (AFI) [13] and the gradient oscillatory number (GON) [12]. AFI quantifies the change of the WSS vector direction, by calculating the angle between the instantaneous and time-averaged WSS vector at any time point throughout the pulse cycle. Mantha et al. [13] found significant correlations between AFI and the location of sidewall aneurysm formation and concluded that the introduced parameter may be useful in identifying areas for future aneurysm formation. However, Shimogonya et al. [12] could not find significant correlations between previously analysed haemodynamic parameters (including AFI) and the location of aneurysm formation, so they proposed a new index for the initiation of a cerebral aneurysm – GON. GON quantifies fluctuations of the compression forces integrated over one cardiac cycle. Although, correlations between an elevated GON and the location of aneurysm formation were found [12], later studies [38, 40] showed that elevated GON values may be found in many non-aneurysmal sites, indicating that GON might be sensitive but not specific for aneurysm formation. The most recently proposed variable – transWSS (transverse wall shear stress) – which describes the multidirectional nature of disturbed flow, was strongly correlated with atherosclerotic lesions [44]. Due to its unique properties, Geers et al. [34] decided to assess transWSS influence in the context of aneurysm initiation. For both cases and controls, regions of high transWSS were concentrated, suggesting that flow disturbances remained in the same location throughout the cardiac cycle. Although they reported that cases were found to have significantly higher transWSS values, there was no clear correlation with the aneurysm initiation site.

Similarly to other authors, we found that the places of highest WSS were located not at the bifurcation apex, but in close proximity. Following flow impingement at a bifurcation apex, it splits into daughter branches and experiences rapid acceleration followed by deceleration creating areas of positive and negative WSSG. In the analysed bifurcations from the case group, we identified three haemodynamic regions previously discovered and described by Meng et al. [41, 42]: impingement, acceleration and recovery zones. As MCA aneurysms tend to form near the bifurcation apex, on one of the daughter branches, aneurysm ostia containing initiation zones were overlapping with regions of high WSS and positive WSSG.

The WSS experienced at a bifurcation is dependent on its geometry, including the radii of all the vessels involved and the bifurcation angle of the arterial tree [45, 46]. The increase in the asymmetry of the vascular diameters of bifurcations and the values of the bifurcation angles, and the decrease in the bilateral angles formed between the daughter and parent vessels and the inclination angle (formed between the parent vessel axis and the plane containing the axes of the two daughter vessels) likely contribute to abnormally enhanced haemodynamic stresses at the arterial bifurcations [26, 47,48,49,50,51,52,53]. Thus, arterial morphology and local haemodynamics are interrelated: the geometry directly determines the blood-flow parameters, while the flow drives arterial wall remodelling which determines the future geometry, promoting aneurysm formation.

Conclusions

This study demonstrated that the development of cerebral aneurysms might be an independent effect of haemodynamic factors. Furthermore, we also demonstrated that high WSS impacts MCA aneurysm formation, while positive WSSG mainly promotes this process. There is a considerable need to study the role of the interactions between haemodynamic forces, arterial geometry and the arterial wall mechanobiology in the natural history of intracranial aneurysm formation in experimental models. Patient-specific fluid structure interaction simulations, involving interactions between the elastic arterial wall and flowing blood, must be enhanced by histopathology studies as these approaches are complementary and mutually informative. Analysis of such simulations may be helpful in discovering the exact mechanisms of IA formation and may enable us to develop new and more efficient diagnostic tools.

Availability of data and materials

The custom-written script used to generate patient-specific inlet velocity profiles is available online (see link: https://github.com/kudlacz964/tccs2ansys).

Abbreviations

absWSSG:

Absolute wall shear stress gradient

ACoA:

Anterior communicating artery

AFI:

Aneurysm formation index

AUC:

Area under the curve

CFD:

Computational fluid dynamics

CI:

Confidence interval

CTA:

Computed tomography angiography

dirWSSG:

Direction wall shear stress

EC:

Endothelial cell

ECM:

Extracellular matrix

GON:

Gradient oscillatory number

IA:

Intracranial aneurysm

ICA:

Internal carotid artery

IEL:

Internal elastic lamina

MCA:

Middle cerebral artery

MMP:

Matrix metallopeptidase

OSI:

Oscillatory shear index

OR:

Odds ratio

PCA:

Posterior cerebral artery

SAH:

Subarachnoid haemorrhage

SMC:

Smooth muscle cell

TCCS:

Transcranial colour-coded sonography

Ved:

End-diastolic velocity

Vm:

Mean velocity

Vps:

Peak systolic velocity

WSS:

Wall shear stress

WSSG:

Wall shear stress gradient

References

  1. 1.

    Vlak MH, Algra A, Brandenburg R, Rinkel GJ. Prevalence of unruptured intracranial aneurysms, with emphasis on sex, age, comorbidity, country, and time period: a systematic review and meta-analysis. Lancet Neurol. 2011;10(7):626–36. https://doi.org/10.1016/S1474-4422(11)70109-0.

    Article  PubMed  Google Scholar 

  2. 2.

    Rinkel GJ, Djibuti M, Algra A, van Gijn J. Prevalence and risk of rupture of intracranial aneurysms: a systematic review. Stroke. 1998;29(1):251–6. https://doi.org/10.1161/01.STR.29.1.251.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Weir B. Unruptured intracranial aneurysms: a review. J Neurosurg. 2002;96(1):3–42. https://doi.org/10.3171/jns.2002.96.1.0003.

    Article  PubMed  Google Scholar 

  4. 4.

    Alfano JM, Kolega J, Natarajan SK, Xiang J, Paluch RA, Levy EI, et al. Intracranial aneurysms occur more frequently at bifurcation sites that typically experience higher hemodynamic stresses. Neurosurgery. 2013;73(3):497–505. https://doi.org/10.1227/NEU.0000000000000016.

    Article  PubMed  Google Scholar 

  5. 5.

    Alnaes MS, Isaksen J, Mardal KA, Romner B, Morgan MK, Ingebrigtsen T. Computation of hemodynamics in the circle of Willis. Stroke. 2007;38(9):2500–5. https://doi.org/10.1161/STROKEAHA.107.482471.

    Article  PubMed  Google Scholar 

  6. 6.

    Kleinloog R, de Mul N, Verweij BH, Post JA, Rinkel GJE, Ruigrok YM. Risk factors for intracranial aneurysm rupture: a systematic review. Neurosurgery. 2018;82(4):431–40. https://doi.org/10.1093/neuros/nyx238.

    Article  PubMed  Google Scholar 

  7. 7.

    Backes D, Rinkel GJ, Laban KG, Algra A, Vergouwen MD. Patient- and aneurysm-specific risk factors for intracranial aneurysm growth: a systematic review and meta-analysis. Stroke. 2016;47(4):951–7. https://doi.org/10.1161/STROKEAHA.115.012162.

    Article  PubMed  Google Scholar 

  8. 8.

    Müller TB, Vik A, Romundstad PR, Sandvei MS. Risk factors for Unruptured intracranial aneurysms and subarachnoid hemorrhage in a prospective population-based study. Stroke. 2019;50(10):2952–5. https://doi.org/10.1161/STROKEAHA.119.025951.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Cebral JR, Meng H. Counterpoint: realizing the clinical utility of computational fluid dynamics--closing the gap. AJNR Am J Neuroradiol. 2012;33(3):396–8. https://doi.org/10.3174/ajnr.A2994.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Meng H, Tutino VM, Xiang J, Siddiqui A. High WSS or low WSS? Complex interactions of hemodynamics with intracranial aneurysm initiation, growth, and rupture: toward a unifying hypothesis. AJNR Am J Neuroradiol. 2014;35(7):1254–62. https://doi.org/10.3174/ajnr.A3558.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Dolan JM, Kolega J, Meng H. High wall shear stress and spatial gradients in vascular pathology: a review. Ann Biomed Eng. 2013;41(7):1411–27. https://doi.org/10.1007/s10439-012-0695-0.

    Article  PubMed  Google Scholar 

  12. 12.

    Shimogonya Y, Ishikawa T, Imai Y, Matsuki N, Yamaguchi T. Can temporal fluctuation in spatial wall shear stress gradient initiate a cerebral aneurysm? A proposed novel hemodynamic index, the gradient oscillatory number (GON). J Biomech. 2009;42(4):550–4. https://doi.org/10.1016/j.jbiomech.2008.10.006.

    Article  PubMed  Google Scholar 

  13. 13.

    Mantha A, Karmonik C, Benndorf G, Strother C, Metcalfe R. Hemodynamics in a cerebral artery before and after the formation of an aneurysm. AJNR Am J Neuroradiol. 2006;27(5):1113–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  14. 14.

    Can A, Du R. Association of Hemodynamic Factors with Intracranial Aneurysm Formation and Rupture: systematic review and meta-analysis. Neurosurgery. 2016;78(4):510–20. https://doi.org/10.1227/NEU.0000000000001083.

    Article  PubMed  Google Scholar 

  15. 15.

    Staarmann B, Smith M, Prestigiacomo CJ. Shear stress and aneurysms: a review. Neurosurg Focus. 2019;47(1):E2. https://doi.org/10.3171/2019.4.FOCUS19225.

    Article  PubMed  Google Scholar 

  16. 16.

    Ku DN, Giddens DP, Zarins CK, Glagov S. Pulsatile flow and atherosclerosis in the human carotid bifurcation. Positive correlation between plaque location and low oscillating shear stress. Arteriosclerosis. 1985;5(3):293–302. https://doi.org/10.1161/01.ATV.5.3.293.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Dolan JM, Sim FJ, Meng H, Kolega J. Endothelial cells express a unique transcriptional profile under very high wall shear stress known to induce expansive arterial remodeling. Am J Physiol Cell Physiol. 2012;302(8):C1109–18. https://doi.org/10.1152/ajpcell.00369.2011.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Sho E, Sho M, Singh TM, Nanjo H, Komatsu M, Xu C, et al. Arterial enlargement in response to high flow requires early expression of matrix metalloproteinases to degrade extracellular matrix. Exp Mol Pathol. 2002;73(2):142–53. https://doi.org/10.1006/exmp.2002.2457.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Tronc F, Mallat Z, Lehoux S, Wassef M, Esposito B, Tedgui A. Role of matrix metalloproteinases in blood flow-induced arterial enlargement: interaction with NO. Arterioscler Thromb Vasc Biol. 2000;20(12):E120–6.

    CAS  Article  Google Scholar 

  20. 20.

    Malek AM, Alper SL, Izumo S. Hemodynamic shear stress and its role in atherosclerosis. Jama. 1999;282(21):2035–42. https://doi.org/10.1001/jama.282.21.2035.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Nixon AM, Gunel M, Sumpio BE. The critical role of hemodynamics in the development of cerebral vascular disease. J Neurosurg. 2010;112(6):1240–53. https://doi.org/10.3171/2009.10.JNS09759.

    Article  PubMed  Google Scholar 

  22. 22.

    Resnick N, Yahav H, Shay-Salit A, Shushy M, Schubert S, Zilberman LC, et al. Fluid shear stress and the vascular endothelium: for better and for worse. Prog Biophys Mol Biol. 2003;81(3):177–99. https://doi.org/10.1016/S0079-6107(02)00052-4.

    Article  PubMed  Google Scholar 

  23. 23.

    Traub O, Berk BC. Laminar shear stress: mechanisms by which endothelial cells transduce an atheroprotective force. Arterioscler Thromb Vasc Biol. 1998;18(5):677–85. https://doi.org/10.1161/01.ATV.18.5.677.

    CAS  Article  PubMed  Google Scholar 

  24. 24.

    Dolan JM, Meng H, Singh S, Paluch R, Kolega J. High fluid shear stress and spatial shear stress gradients affect endothelial proliferation, survival, and alignment. Ann Biomed Eng. 2011;39(6):1620–31. https://doi.org/10.1007/s10439-011-0267-8.

    Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Dolan JM, Meng H, Sim FJ, Kolega J. Differential gene expression by endothelial cells under positive and negative streamwise gradients of high wall shear stress. Am J Physiol Cell Physiol. 2013;305(8):C854–66. https://doi.org/10.1152/ajpcell.00315.2012.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Kaspera W, Ćmiel-Smorzyk K, Wolański W, Kawlewska E, Hebda A, Gzik M, et al. Morphological and hemodynamic risk factors for middle cerebral artery aneurysm: a case-control study of 190 patients. Sci Rep. 2020;10(1):2016. https://doi.org/10.1038/s41598-019-56061-2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Kaspera W, Majchrzak H, Ladzinski P, Tomalski W. Color Doppler sonographic evaluation of collateral circulation in patients with cerebral aneurysms and the occlusion of the brachiocephalic vessels. Neurosurgery. 2005;57(6):1117–26 discussion −26.

    Article  Google Scholar 

  28. 28.

    Resnick N, Yahav H, Khachigian LM, Collins T, Anderson KR, Dewey FC, et al. Endothelial gene regulation by laminar shear stress. Adv Exp Med Biol. 1997;430:155–64. https://doi.org/10.1007/978-1-4615-5959-7_13.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Resnick N, Yahav H, Schubert S, Wolfovitz E, Shay A. Signalling pathways in vascular endothelium activated by shear stress: relevance to atherosclerosis. Curr Opin Lipidol. 2000;11(2):167–77. https://doi.org/10.1097/00041433-200004000-00010.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Kim C, Kikuchi H, Hashimoto N, Kojima M, Kang Y, Hazama F. Involvement of internal elastic lamina in development of induced cerebral aneurysms in rats. Stroke. 1988;19(4):507–11. https://doi.org/10.1161/01.STR.19.4.507.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Li MH, Li PG, Huang QL, Ling J. Endothelial injury preceding intracranial aneurysm formation in rabbits. West Indian Med J. 2014;63(2):167–71. https://doi.org/10.7727/wimj.2013.129.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Sheinberg DL, McCarthy DJ, Elwardany O, Bryant JP, Luther E, Chen SH, et al. Endothelial dysfunction in cerebral aneurysms. Neurosurg Focus. 2019;47(1):E3. https://doi.org/10.3171/2019.4.FOCUS19221.

    Article  PubMed  Google Scholar 

  33. 33.

    Baek H, Jayaraman MV, Karniadakis GE. Wall shear stress and pressure distribution on aneurysms and infundibulae in the posterior communicating artery bifurcation. Ann Biomed Eng. 2009;37(12):2469–87. https://doi.org/10.1007/s10439-009-9794-y.

    Article  PubMed  Google Scholar 

  34. 34.

    Geers AJ, Morales HG, Larrabide I, Butakoff C, Bijlenga P, Frangi AF. Wall shear stress at the initiation site of cerebral aneurysms. Biomech Model Mechanobiol. 2017;16(1):97–115. https://doi.org/10.1007/s10237-016-0804-3.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Kulcsár Z, Ugron A, Marosfoi M, Berentei Z, Paál G, Szikora I. Hemodynamics of cerebral aneurysm initiation: the role of wall shear stress and spatial wall shear stress gradient. AJNR Am J Neuroradiol. 2011;32(3):587–94. https://doi.org/10.3174/ajnr.A2339.

    Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Le WJ, Zhu YQ, Li MH, Yan L, Tan HQ, Xiao SM, et al. New method for retrospective study of hemodynamic changes before and after aneurysm formation in patients with ruptured or unruptured aneurysms. BMC Neurol. 2013;13(1):166. https://doi.org/10.1186/1471-2377-13-166.

    Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Singh PK, Marzo A, Howard B, Rufenacht DA, Bijlenga P, Frangi AF, et al. Effects of smoking and hypertension on wall shear stress and oscillatory shear index at the site of intracranial aneurysm formation. Clin Neurol Neurosurg. 2010;112(4):306–13. https://doi.org/10.1016/j.clineuro.2009.12.018.

    Article  PubMed  Google Scholar 

  38. 38.

    Ford MD, Hoi Y, Piccinelli M, Antiga L, Steinman DA. An objective approach to digital removal of saccular aneurysms: technique and applications. Br J Radiol. 2009;82 Spec No 1:S55–61.

    CAS  Article  Google Scholar 

  39. 39.

    Lauric A, Hippelheuser J, Safain MG, Malek AM. Curvature effect on hemodynamic conditions at the inner bend of the carotid siphon and its relation to aneurysm formation. J Biomech. 2014;47(12):3018–27. https://doi.org/10.1016/j.jbiomech.2014.06.042.

    Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Chen H, Selimovic A, Thompson H, Chiarini A, Penrose J, Ventikos Y, et al. Investigating the influence of haemodynamic stimuli on intracranial aneurysm inception. Ann Biomed Eng. 2013;41(7):1492–504. https://doi.org/10.1007/s10439-013-0794-6.

    Article  PubMed  Google Scholar 

  41. 41.

    Meng H, Swartz DD, Wang Z, Hoi Y, Kolega J, Metaxa EM, et al. A model system for mapping vascular responses to complex hemodynamics at arterial bifurcations in vivo. Neurosurgery. 2006;59(5):1094–100; discussion 100-1. https://doi.org/10.1227/01.NEU.0000245599.92322.53.

    Article  PubMed  Google Scholar 

  42. 42.

    Meng H, Wang Z, Hoi Y, Gao L, Metaxa E, Swartz DD, et al. Complex hemodynamics at the apex of an arterial bifurcation induces vascular remodeling resembling cerebral aneurysm initiation. Stroke. 2007;38(6):1924–31. https://doi.org/10.1161/STROKEAHA.106.481234.

    Article  PubMed  PubMed Central  Google Scholar 

  43. 43.

    Metaxa E, Tremmel M, Natarajan SK, Xiang J, Paluch RA, Mandelbaum M, et al. Characterization of critical hemodynamics contributing to aneurysmal remodeling at the basilar terminus in a rabbit model. Stroke. 2010;41(8):1774–82. https://doi.org/10.1161/STROKEAHA.110.585992.

    Article  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Peiffer V, Sherwin SJ, Weinberg PD. Computation in the rabbit aorta of a new metric - the transverse wall shear stress - to quantify the multidirectional character of disturbed blood flow. J Biomech. 2013;46(15):2651–8. https://doi.org/10.1016/j.jbiomech.2013.08.003.

    Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Tütüncü F, Schimansky S, Baharoglu MI, Gao B, Calnan D, Hippelheuser J, et al. Widening of the basilar bifurcation angle: association with presence of intracranial aneurysm, age, and female sex. J Neurosurg. 2014;121(6):1401–10. https://doi.org/10.3171/2014.8.JNS1447.

    Article  PubMed  Google Scholar 

  46. 46.

    Zhang XJ, Li CH, Hao WL, Zhang DH, Ren CF, Gao BL. Enlarged anterior cerebral artery bifurcation angles may induce abnormally enhanced hemodynamic stresses to initiate aneurysms. World Neurosurg. 2018;120:e783–e91. https://doi.org/10.1016/j.wneu.2018.08.167.

    Article  PubMed  Google Scholar 

  47. 47.

    Baharoglu MI, Lauric A, Wu C, Hippelheuser J, Malek AM. Deviation from optimal vascular caliber control at middle cerebral artery bifurcations harboring aneurysms. J Biomech. 2014;47(13):3318–24. https://doi.org/10.1016/j.jbiomech.2014.08.012.

    Article  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Kaspera W, Ładziński P, Larysz P, Hebda A, Ptaszkiewicz K, Kopera M, et al. Morphological, hemodynamic, and clinical independent risk factors for anterior communicating artery aneurysms. Stroke. 2014;45(10):2906–11. https://doi.org/10.1161/STROKEAHA.114.006055.

    Article  PubMed  Google Scholar 

  49. 49.

    Zhang XJ, Gao BL, Hao WL, Wu SS, Zhang DH. Presence of anterior communicating artery aneurysm is associated with age, bifurcation angle, and vessel diameter. Stroke. 2018;49(2):341–7. https://doi.org/10.1161/STROKEAHA.117.019701.

    Article  PubMed  Google Scholar 

  50. 50.

    Zhang XJ, Gao BL, Li TX, Hao WL, Wu SS, Zhang DH. Association of Basilar Bifurcation Aneurysms with age, sex, and bifurcation geometry. Stroke. 2018;49(6):1371–6. https://doi.org/10.1161/STROKEAHA.118.020829.

    Article  PubMed  Google Scholar 

  51. 51.

    Zhang XJ, Hao WL, Zhang DH, Gao BL. Asymmetrical middle cerebral artery bifurcations are more vulnerable to aneurysm formation. Sci Rep. 2019;9(1):15255. https://doi.org/10.1038/s41598-019-51734-4.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. 52.

    Zhang XJ, Hao WL, Zhang DH, Gao BL. Asymmetrical than symmetrical cerebral arterial bifurcations are more vulnerable to aneurysm presence. Sci Rep. 2019;9(1):17144. https://doi.org/10.1038/s41598-019-53715-z.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  53. 53.

    Sadatomo T, Yuki K, Migita K, Imada Y, Kuwabara M, Kurisu K. Differences between middle cerebral artery bifurcations with normal anatomy and those with aneurysms. Neurosurg Rev. 2013;36(3):437–45. https://doi.org/10.1007/s10143-013-0450-5.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by funding from the Silesian University of Technology, Zabrze, Poland (Project No. 07/030/RGP19/2041).

Author information

Affiliations

Authors

Contributions

WK conceptualised and supervised the study. MZ, WK, EK and WW acquired the data. WK, MZ and WW analysed data. AH conducted the statistical analysis. MZ and WK wrote the manuscript. WK, WW and PŁ edited the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Wojciech Kaspera.

Ethics declarations

Ethics approval and consent to participate

The study’s protocol was approved by the Institutional Review Board at the Medical University of Silesia in Katowice and all procedures were carried out in accordance with the relevant guidelines and regulations. Each patient was informed about the purpose and course of the research and asked to give their informed consent to participate in the project.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zimny, M., Kawlewska, E., Hebda, A. et al. Wall shear stress gradient is independently associated with middle cerebral artery aneurysm development: a case-control CFD patient-specific study based on 77 patients. BMC Neurol 21, 281 (2021). https://doi.org/10.1186/s12883-021-02251-3

Download citation

Keywords

  • Middle cerebral artery
  • Intracranial aneurysm formation
  • Wall shear stress
  • Wall shear stress gradient
  • Computational fluid dynamics
  • Haemodynamics