High wall shear stress and spatial gradients in vascular pathology: a review

Uncertainty quantification of wall shear stress in intracranial aneurysms using a data-driven statistical model of systemic blood flow variability

Adverse wall shear stress (WSS) patterns are known to play a key role in the localisation, formation, and progression of intracranial aneurysms (IAs). Complex region-specific and time-varying aneurysmal WSS patterns depend both on vascular morphology as well as on variable systemic flow conditions. Computational fluid dynamics (CFD) has been proposed for characterising WSS patterns in IAs; however, CFD simulations often rely on deterministic boundary conditions that are not representative of the actual variations in blood flow. We develop a data-driven statistical model of internal carotid artery (ICA) flow, which is used to generate a virtual population of waveforms used as inlet boundary conditions in CFD simulations. This allows the statistics of the resulting aneurysmal WSS distributions to be computed. It is observed that ICA waveform variations have limited influence on the time-averaged WSS (TAWSS) on the IA surface. In contrast, in regions where the flow is locally highly multidirectional, WSS directionality and harmonic content are strongly affected by the ICA flow waveform. As a consequence, we argue that the effect of blood flow variability should be explicitly considered in CFD-based IA rupture assessment to prevent confounding the conclusions.

C-Type Natriuretic Peptides in Coronary Disease

AIMS

C-type natriuretic peptide (CNP) is a paracrine growth factor expressed in the vascular endothelium. Although upregulated in atheromatous arteries, the predictive value of plasma CNP products for outcome in coronary disease is unknown. This study aimed to compare the prognostic value of plasma CNP products with those of other natriuretic peptides in individuals with coronary artery disease, and investigate their associations with cardiac and renal function.

METHODS AND RESULTS

Plasma concentrations of CNP and amino-terminal proCNP (NT-proCNP) were measured at baseline in 2129 individuals after an index acute coronary syndrome admission and related to cardiac and renal function, other natriuretic peptides [atrial NP (ANP) and B-type NP (BNP)] and prognosis (primary end point, mortality; secondary end point, cardiac readmission). Median follow-up was 4 years. At baseline, and in contrast to CNP, ANP, and BNP, plasma NT-proCNP was higher in males and weakly related to cardiac function but strongly correlated to plasma creatinine. All NPs were univariately associated with mortality. Resampling at 4 and 12 months in survivors showed stable concentrations of NT-proCNP whereas all other peptides declined. When studied by diagnosis (myocardial infarction, unstable angina) at index admission using a multivariate model, NT-proBNP predicted mortality and readmission in myocardial infarction. In unstable angina, only NT-proCNP predicted both mortality and cardiac readmission.

CONCLUSIONS

In contrast to the close association of NT-proBNP with cardiac function, and predictive value for outcome after myocardial infarction, plasma NT-proCNP is highly correlated with renal function and is an independent predictor of mortality and cardiac readmission in individuals with unstable angina.

Characterizations and Correlations of Wall Shear Stress in Aneurysmal Flow

Wall shear stress (WSS) is one of the most studied hemodynamic parameters, used in correlating blood flow to various diseases. The pulsatile nature of blood flow, along with the complex geometries of diseased arteries, produces complicated temporal and spatial WSS patterns. Moreover, WSS is a vector, which further complicates its quantification and interpretation. The goal of this study is to investigate WSS magnitude, angle, and vector changes in space and time in complex blood flow. Abdominal aortic aneurysm (AAA) was chosen as a setting to explore WSS quantification. Patient-specific computational fluid dynamics (CFD) simulations were performed in six AAAs. New WSS parameters are introduced, and the pointwise correlation among these, and more traditional WSS parameters, was explored. WSS magnitude had positive correlation with spatial/temporal gradients of WSS magnitude. This motivated the definition of relative WSS gradients. WSS vectorial gradients were highly correlated with magnitude gradients. A mix WSS spatial gradient and a mix WSS temporal gradient are proposed to equally account for variations in the WSS angle and magnitude in single measures. The important role that WSS plays in regulating near wall transport, and the high correlation among some of the WSS parameters motivates further attention in revisiting the traditional approaches used in WSS characterizations.

Effect of Different Rotational Directions of BJUT-II VAD on Aortic Swirling Flow Characteristics: A Primary Computational Fluid Dynamics Study

Background

The BJUT-II VAD is a novel left ventricular assist device (LVAD), which is thought to have significant effects on the characteristics of aortic swirling flow. However, the precise mechanism of the rotational direction of BJTU-II VAD in the aortic swirling flow is unclear.

Material/Methods

A patient-specific aortic geometric model was reconstructed based on the CT data. Three pump’s output flow profiles with varied rotational direction, termed “counterclockwise”, “flat profile”, and “clockwise”, were used as the boundary conditions. The helicity density, area-weighted average of helicity density (Ha), localized normalized helicity (LNH), wall shear stress (WSS), and WSS spatial gradient (WSSG) were calculated to evaluate the swirling flow characteristics in the aorta.

Results

The results demonstrated that the swirling flow characteristics in the aorta and 3 branches are directly affected by the output blood flow of BJUT-II VAD. In the aortic arch, the helicity density, supported by the clockwise case, achieved the highest value. In the 3 branches, the flat profile case achieved the highest helicity density, whereas the maximum WSS and WSSG generated by clockwise case were lower than in other cases.

Conclusions

The outflow of the BJUT-II VAD has significant effects on the aortic hemodynamics and swirling flow characteristics. The helical blood profiles can enhance the strength of aortic swirling flow, and reduce the areas of low WSS and WSSG regions. The clockwise case may have a benefit for preventing development of atherosclerosis in the aorta.

Correlation Between Aneurysm Size and Hemodynamics in One Individual with Multiple Small Intracranial Aneurysms

Objective

A large number of cases are needed in the patient-specific modeling of intracranial aneurysms to establish the statistical significance due to individual variation of risk factors that are difficult to account for. However, these risk factors are critical in hemorrhage risk as demonstrated in large clinical studies. Rupture risks for aneurysms in an individual are easier to compare because these aneurysms are under the same physiological environment, and their only differences are the local hemodynamic factors associated with their anatomic locations.

Methods

Eight small aneurysms (< 7 mm) from one individual were analyzed using patient-specific hemodynamic modeling. Four scenarios with different perfusion assumptions were performed to account for the flow rate at two smaller communicating arteries. Wall shear stresses (WSS) at these aneurysms were compared to determine their relationship with the aneurysm size.

Results

Each of the three largest aneurysms is either the most proximal or distal aneurysm in a given artery so that blood pressure does not have a direct influence on aneurysm size. No wall shear stress-derived hemodynamic variables are found to be related to aneurysm size.

Discussion

A study of multiple aneurysms from one individual offers a unique opportunity to examine various hemodynamic factors without selection biases. Aneurysms greater than 4 mm (Group 1) have a higher product of maximum WSS and area of low WSS; aneurysms smaller than 4 mm (Group 2) have a lower product of maximum WSS and area of low WSS. In addition, aneurysm size is linearly correlated with the flow rate at the parent artery in each group.

Sustained expression of MCP-1 by low wall shear stress loading concomitant with turbulent flow on endothelial cells of intracranial aneurysm

INTRODUCTION

Enlargement of a pre-existing intracranial aneurysm is a well-established risk factor of rupture. Excessive low wall shear stress concomitant with turbulent flow in the dome of an aneurysm may contribute to progression and rupture. However, how stress conditions regulate enlargement of a pre-existing aneurysm remains to be elucidated.

RESULTS

Wall shear stress was calculated with 3D-computational fluid dynamics simulation using three cases of unruptured intracranial aneurysm. The resulting value, 0.017 Pa at the dome, was much lower than that in the parent artery. We loaded wall shear stress corresponding to the value and also turbulent flow to the primary culture of endothelial cells. We then obtained gene expression profiles by RNA sequence analysis. RNA sequence analysis detected hundreds of differentially expressed genes among groups. Gene ontology and pathway analysis identified signaling related with cell division/proliferation as overrepresented in the low wall shear stress-loaded group, which was further augmented by the addition of turbulent flow. Moreover, expression of some chemoattractants for inflammatory cells, including MCP-1, was upregulated under low wall shear stress with concomitant turbulent flow. We further examined the temporal sequence of expressions of factors identified in an in vitro study using a rat model. No proliferative cells were detected, but MCP-1 expression was induced and sustained in the endothelial cell layer.

CONCLUSIONS

Low wall shear stress concomitant with turbulent flow contributes to sustained expression of MCP-1 in endothelial cells and presumably plays a role in facilitating macrophage infiltration and exacerbating inflammation, which leads to enlargement or rupture.

Impact of bifurcation angle and other anatomical characteristics on blood flow - A computational study of non-stented and stented coronary arteries

The hemodynamic influence of vessel shape such as bifurcation angle is not fully understood with clinical and quantitative observations being equivocal. The aim of this study is to use computational modeling to study the hemodynamic effect of shape characteristics, in particular bifurcation angle (BA), for non-stented and stented coronary arteries. Nine bifurcations with angles of 40°, 60° and 80°, representative of ±1 SD of 101 asymptomatic computed tomography angiogram cases (average age 54±8 years; 57 females), were generated for (1) a non-stented idealized, (2) stented idealized, and (3) non-stented patient-specific geometry. Only the bifurcation angle was changed while the geometries were constant to eliminate flow effects induced by other vessel shape characteristics. The commercially available Biomatrix stent was used as a template and virtually inserted into each branch, simulating the T-stenting technique. Three patient-specific geometries with additional shape variation and ±2 SD BA variation (33°, 42° and 117°) were also computed. Computational fluid dynamics (CFD) analysis was performed for all 12 geometries to simulate physiological conditions, enabling the quantification of the hemodynamic stress distributions, including a threshold analysis of adversely low and high wall shear stress (WSS), low time-averaged WSS (TAWSS), high spatial WSS gradient (WSSG) and high Oscillatory Shear Index (OSI) area. The bifurcation angle had a minor impact on the areas of adverse hemodynamics in the idealized non-stented geometries, which fully disappeared once stented and was not apparent for patient geometries. High WSS regions were located close to the carina around peak-flow, and WSSG increased significantly after stenting for the idealized bifurcations. Additional shape variations affected the hemodynamic profiles, suggesting that BA alone has little effect on a patient׳s hemodynamic profile. Incoming flow angle, diameter and tortuosity appear to have stronger effects. This suggests that other bifurcation shape characteristics and stent placement/strategy may be more important than bifurcation angle in atherosclerotic disease development, progression, and stent outcome.

An integrated temperature-compensated flexible shear-stress sensor microarray with concentrated leading-wire

Flexible shear stress sensor is quite important for characterizing curved surface flows. In this work, a novel integrated shear stress sensor microarray is designed with twenty parallel channels, which share the concentrated leading-wire to transmit the ground signal. Electrical pads in rows are easily connected to the circuits with two separate Wheatstone bridges and constant-temperature-difference mode operation is provided for the hot-wires. Temperature crosstalk between adjacent hot-wires is prevented well and the effectiveness of the temperature compensated circuits is verified. Relatively large output response is obtained as the shear stress varies and the sensitivity of the sensors is measured about 0.086 V(2)/Pa(1/3) with nonlinearity lower than 1%, revealing high performance characteristic of the sensors.

Comprehensive Assessment of Coronary Plaque Progression With Advanced Intravascular Imaging, Physiological Measures, and Wall Shear Stress: A Pilot Double-Blinded Randomized Controlled Clinical Trial of Nebivolol Versus Atenolol in Nonobstructive Coronary Artery Disease

Background

We hypothesized that nebivolol, a β‐blocker with nitric oxide–mediated activity, compared with atenolol, a β‐blocker without such activity, would decrease oxidative stress and improve the effects of endothelial dysfunction and wall shear stress (WSS), thereby reducing atherosclerosis progression and vulnerability in patients with nonobstructive coronary artery disease.

Methods and Results

In this pilot double‐blinded randomized controlled trial, 24 patients treated for 1 year with nebivolol 10 mg versus atenolol 100 mg plus standard medical therapy underwent baseline and follow‐up coronary angiography with assessments of inflammatory and oxidative stress biomarkers, microvascular function, endothelial function, and virtual histology intravascular ultrasound. WSS was calculated from computational fluid dynamics. Virtual histology intravascular ultrasound segments were assessed for vessel volumetrics and remodeling. There was a trend toward more low‐WSS segments in the nebivolol cohort (P=0.06). Low‐WSS regions were associated with greater plaque progression (P<0.0001) and constrictive remodeling (P=0.04); conversely, high‐WSS segments demonstrated plaque regression and excessive expansive remodeling. Nebivolol patients had decreased lumen and vessel areas along with increased plaque area, resulting in more constrictive remodeling (P=0.002). There were no significant differences in biomarker levels, microvascular function, endothelial function, or number of thin‐capped fibroatheromas per vessel. Importantly, after adjusting for β‐blocker, low‐WSS segments remained significantly associated with lumen loss and plaque progression.

Conclusion

Nebivolol, compared with atenolol, was associated with greater plaque progression and constrictive remodeling, likely driven by more low‐WSS segments in the nebivolol arm. Both β‐blockers had similar effects on oxidative stress, microvascular function, and endothelial function.

Clinical Trial Registration

URL: https://clinicaltrials.gov/. Unique identifier: NCT01230892.

Flow Dynamics of Aneurysm Growth and Rupture: Challenges for the Development of Computational Flow Dynamics as a Diagnostic Tool to Detect Rupture-Prone Aneurysms

Saccular intracranial aneurysm (sIA) is a relatively common disease that can potentially cause a devastating, life-threatening intracranial hemorrhage. Many sIAs never rupture and thus do not necessitate interventions, making the detection of rupture-prone sIAs a very relevant clinical problem. Moreover, because currently available methods to prevent sIA rupture have significant risks of morbidity and mortality, diagnostic tools that can predict imminent rupture and help plan proper timing of prophylactic interventions, can improve patient care. Hemorrhage from an sIA occurs when hemodynamic stress exceeds sIA wall strength. Computational fluid dynamics (CFD) is a tool with which the hemodynamic stress to which the sIA wall is exposed can be determined non-invasively. Studies using CFD in sIAs have demonstrated associations of wall shear stress (WSS) with aneurysm growth, fragile sIA wall, and sIA rupture; these studies show the potential of CFD as a diagnostic tool. This review discusses the limitations of CFD and of the studies performed, and what needs to be done in order to develop CFD into a useful diagnostic tool to determine aneurysm-specific rupture risk.

Cellular Biomechanics in Drug Screening and Evaluation: Mechanopharmacology

The study of mechanobiology is now widespread. The impact of cell and tissue mechanics on cellular responses is well appreciated. However, knowledge of the impact of cell and tissue mechanics on pharmacological responsiveness, and its application to drug screening and mechanistic investigations, have been very limited in scope. We emphasize the need for a heightened awareness of the important bidirectional influence of drugs and biomechanics in all living systems. We propose that the term 'mechanopharmacology' be applied to approaches that employ in vitro systems, biomechanically appropriate to the relevant (patho)physiology, to identify new drugs and drug targets. This article describes the models and techniques that are being developed to transform drug screening and evaluation, ranging from a 2D environment to the dynamic 3D environment of the target expressed in the disease of interest.

Specific Shrinkage of Carotid Forks in Moyamoya Disease: A Novel Key Finding for Diagnosis

This study was aimed to analyze the outer diameter of the involved arteries in moyamoya disease, using three-dimensional (3D) constructive interference in steady state (CISS) and direct surgical inspection. Radiological evaluation was performed in 64 patients with moyamoya disease. As the controls, six patients with severe middle cerebral artery (MCA) stenosis and 17 healthy subjects were also recruited. On 3D-CISS, the outer diameter was quantified in the supraclinoid portion of internal carotid artery (C1), the horizontal portions of MCA (M1) and anterior cerebral artery (A1), and basilar artery. The involved carotid fork was directly observed during surgery in another series of three adult patients with moyamoya disease. In 53 adult patients with moyamoya disease, the outer diameters of C1, M1, and A1 segments were 2.3 ± 0.7 mm, 1.3 ± 0.5 mm, and 1.0 ± 0.4 mm in the involved side (n = 91), being significantly smaller than the control (n = 17), severe M1 stenosis (n = 6), and non-involved side in moyamoya disease (n = 15, P < 0.01). There were significant correlations between Suzuki’s angiographical stage and the outer diameters of C1, M1, and A1 (P < 0.001). The laterality ratio of C1 and M1 was significantly smaller in unilateral moyamoya disease (n = 20) than the controls and severe MCA stenosis (P < 0.01). Direct observations revealed a marked decrease in the outer diameter of the carotid fork (n = 3). These findings strongly suggest specific shrinkage of the involved arteries in moyamoya disease, which may provide essential information to distinguish moyamoya disease from other intracranial arterial stenosis and shed light on the etiology and novel diagnosis cue of moyamoya disease.

Pressure Changes Within the Sac of Human Cerebral Aneurysms in Response to Artificially Induced Transient Increases in Systemic Blood Pressure

Formation and rupture of cerebral aneurysms have been associated with chronic hypertension. The effect of transient increase in blood pressure and its effect on intra-aneurysmal hemodynamics have not been studied. We examined the effects of controlled increases in blood pressure on different pressure parameters inside the sac of human cerebral aneurysms and corresponding parent arteries using invasive technology. Twelve patients (10 female, 2 male, age 54±15 years) with unruptured cerebral aneurysms undergoing endovascular coiling were recruited. Dual-sensor microwires with the capacity to simultaneously measure flow velocity and pressure were used to measure systolic, diastolic, and mean pressure inside the aneurysm sac and to measure both pressures and flow velocities in the feeder vessel just outside the aneurysm. These pressures were recorded simultaneously with pressures from a radial arterial catheter. Measurements were taken at baseline and then during a gradual increase in systemic systolic blood pressure to a target value of ≈25 mm Hg above baseline, using a phenylephrine infusion. The dose needed to achieve the required increase in radial arterial systolic blood pressure was 0.8±0.2 μg/kg/min. There was a clear linear relationship between changes in radial and aneurysmal pressures with substantial patient-by-patient variation in the slopes of those relationships. The overall increases in systolic and mean pressures in both radial artery and in the aneurysms were similar. Pressures in the aneurysm and in the parent vessels were similar. Peak and mean flow velocities in the parent arteries did not change significantly with phenylephrine infusion, nor did vessel diameters as measured angiographically.

"Do-it-yourself in vitro vasculature that recapitulates in vivo geometries for investigating endothelial-blood cell interactions"

Investigating biophysical cellular interactions in the circulation currently requires choosing between in vivo models, which are difficult to interpret due in part to the hemodynamic and geometric complexities of the vasculature; or in vitro systems, which suffer from non-physiologic assumptions and/or require specialized microfabrication facilities and expertise. To bridge that gap, we developed an in vitro "do-it-yourself" perfusable vasculature model that recapitulates in vivo geometries, such as aneurysms, stenoses, and bifurcations, and supports endothelial cell culture. These inexpensive, disposable devices can be created rapidly (<2 hours) with high precision and repeatability, using standard off-the-shelf laboratory supplies. Using these "endothelialized" systems, we demonstrate that spatial variation in vascular cell adhesion molecule (VCAM-1) expression correlates with the wall shear stress patterns of vascular geometries. We further observe that the presence of endothelial cells in stenoses reduces platelet adhesion but increases sickle cell disease (SCD) red blood cell (RBC) adhesion in bifurcations. Overall, our method enables researchers from all disciplines to study cellular interactions in physiologically relevant, yet simple-to-make, in vitro vasculature models.

Hemodynamics in Idealized Stented Coronary Arteries: Important Stent Design Considerations

Stent induced hemodynamic changes in the coronary arteries are associated with higher risk of adverse clinical outcome. The purpose of this study was to evaluate the impact of stent design on wall shear stress (WSS), time average WSS, and WSS gradient (WSSG), in idealized stent geometries using computational fluid dynamics. Strut spacing, thickness, luminal protrusion, and malapposition were systematically investigated and a comparison made between two commercially available stents (Omega and Biomatrix). Narrower strut spacing led to larger areas of adverse low WSS and high WSSG but these effects were mitigated when strut size was reduced, particularly for WSSG. Local hemodynamics worsened with luminal protrusion of the stent and with stent malapposition, adverse high WSS and WSSG were identified around peak flow and throughout the cardiac cycle respectively. For the Biomatrix stent, the adverse effect of thicker struts was mitigated by greater strut spacing, radial cell offset and flow-aligned struts. In conclusion, adverse hemodynamic effects of specific design features (such as strut size and narrow spacing) can be mitigated when combined with other hemodynamically beneficial design features but increased luminal protrusion can worsen the stent’s hemodynamic profile significantly.

Patient-specific computational modeling of blood flow in the pulmonary arterial circulation

Computational fluid dynamics (CFD) modeling of the pulmonary vasculature has the potential to reveal continuum metrics associated with the hemodynamic stress acting on the vascular endothelium. It is widely accepted that the endothelium responds to flow-induced stress by releasing vasoactive substances that can dilate and constrict blood vessels locally. The objectives of this study are to examine the extent of patient specificity required to obtain a significant association of CFD output metrics and clinical measures in models of the pulmonary arterial circulation, and to evaluate the potential correlation of wall shear stress (WSS) with established metrics indicative of right ventricular (RV) afterload in pulmonary hypertension (PH). Right Heart Catheterization (RHC) hemodynamic data and contrast-enhanced computed tomography (CT) imaging were retrospectively acquired for 10 PH patients and processed to simulate blood flow in the pulmonary arteries. While conducting CFD modeling of the reconstructed patient-specific vasculatures, we experimented with three different outflow boundary conditions to investigate the potential for using computationally derived spatially averaged wall shear stress (SAWSS) as a metric of RV afterload. SAWSS was correlated with both pulmonary vascular resistance (PVR) (R(2)=0.77, P<0.05) and arterial compliance (C) (R(2)=0.63, P<0.05), but the extent of the correlation was affected by the degree of patient specificity incorporated in the fluid flow boundary conditions. We found that decreasing the distal PVR alters the flow distribution and changes the local velocity profile in the distal vessels, thereby increasing the local WSS. Nevertheless, implementing generic outflow boundary conditions still resulted in statistically significant SAWSS correlations with respect to both metrics of RV afterload, suggesting that the CFD model could be executed without the need for complex outflow boundary conditions that require invasively obtained patient-specific data. A preliminary study investigating the relationship between outlet diameter and flow distribution in the pulmonary tree offers a potential computationally inexpensive alternative to pressure based outflow boundary conditions.

The role of wall shear stress in the assessment of right ventricle hydraulic workload

Pulmonary hypertension (PH) is a devastating disease affecting approximately 15-50 people per million, with a higher incidence in women. PH mortality is mostly attributed to right ventricle (RV) failure, which results from RV hypotrophy due to an overburdened hydraulic workload. The objective of this study is to correlate wall shear stress (WSS) with hemodynamic metrics that are generally accepted as clinical indicators of RV workload and are well correlated with disease outcome. Retrospective right heart catheterization data for 20 PH patients were analyzed to derive pulmonary vascular resistance (PVR), arterial compliance (C), and an index of wave reflections (Γ). Patient-specific contrast-enhanced computed tomography chest images were used to reconstruct the individual pulmonary arterial trees up to the seventh generation. Computational fluid dynamics analyses simulating blood flow at peak systole were conducted for each vascular model to calculate WSS distributions on the endothelial surface of the pulmonary arteries. WSS was found to be decreased proportionally with elevated PVR and reduced C. Spatially averaged WSS (SAWSS) was positively correlated with PVR (R (2) = 0.66), C (R (2) = 0.73), and Γ (R (2) = 0.5) and also showed promising preliminary correlations with RV geometric characteristics. Evaluating WSS at random cross sections in the proximal vasculature (main, right, and left pulmonary arteries), the type of data that can be acquired from phase-contrast magnetic resonance imaging, did not reveal the same correlations. In conclusion, we found that WSS has the potential to be a viable and clinically useful noninvasive metric of PH disease progression and RV health. Future work should be focused on evaluating whether SAWSS has prognostic value in the management of PH and whether it can be used as a rapid reactivity assessment tool, which would aid in selection of appropriate therapies.

Morphology and related hemodynamics of the internal carotid arteries of moyamoya patients

BACKGROUND

Morphological studies investigating the intracranial-extradural internal carotid artery with moyamoya disease have not been reported. We designed this case-control study to investigate the morphological differences of the internal carotid artery with moyamoya disease, and to clarify the contributions of these differences to the resultant fluid dynamics.

METHODS

Patients with moyamoya disease and normal controls were assigned to each group. The vascular tortuosity of internal carotid artery was measured with three-dimensional rendering using magnetic resonance angiography. By computational fluid dynamics, hemodynamic characteristics were simulated and compared between two groups.

RESULTS

Distances were measured from the carotid canal to the siphon. A shorter actual distance was observed in the moyamoya group (p = 0.0170). Vascular tortuosity was significantly low in moyamoya patients showing lower curvature angles in the petrous and intra-cavernous segments (p = 0.0012). Less blood flowed (p < 0.0001) through the narrower internal carotid artery (p < 0.0001) in the moyamoya group at the carotid canal level. The blood flow velocities were not significantly different (p = 0.2332). Faster blood flow and higher wall shear stress in the internal carotid artery bifurcation were verified with computational fluid dynamics.

CONCLUSIONS

Significant morphological differences were confirmed to exist in the intracranial-extradural internal carotid artery of moyamoya patients. These differences might influence the hemodynamics around the bifurcation of the internal carotid artery.

Narrowing the Expertise Gap for Predicting Intracranial Aneurysm Hemodynamics: Impact of Solver Numerics versus Mesh and Time-Step Resolution

BACKGROUND AND PURPOSE

Recent high-resolution computational fluid dynamics studies have uncovered the presence of laminar flow instabilities and possible transitional or turbulent flow in some intracranial aneurysms. The purpose of this study was to elucidate requirements for computational fluid dynamics to detect these complex flows, and, in particular, to discriminate the impact of solver numerics versus mesh and time-step resolution.

MATERIALS AND METHODS

We focused on 3 MCA aneurysms, exemplifying highly unstable, mildly unstable, or stable flow phenotypes, respectively. For each, the number of mesh elements was varied by 320× and the number of time-steps by 25×. Computational fluid dynamics simulations were performed by using an optimized second-order, minimally dissipative solver, and a more typical first-order, stabilized solver.

RESULTS

With the optimized solver and settings, qualitative differences in flow and wall shear stress patterns were negligible for models down to ∼800,000 tetrahedra and ∼5000 time-steps per cardiac cycle and could be solved within clinically acceptable timeframes. At the same model resolutions, however, the stabilized solver had poorer accuracy and completely suppressed flow instabilities for the 2 unstable flow cases. These findings were verified by using the popular commercial computational fluid dynamics solver, Fluent.

CONCLUSIONS

Solver numerics must be considered at least as important as mesh and time-step resolution in determining the quality of aneurysm computational fluid dynamics simulations. Proper computational fluid dynamics verification studies, and not just superficial grid refinements, are therefore required to avoid overlooking potentially clinically and biologically relevant flow features.

The effect of in-plane arterial curvature on blood flow and oxygen transport in arterio-venous fistulae

Arterio-Venous Fistulae (AVF) are the preferred method of vascular access for patients with end stage renal disease who need hemodialysis. In this study, simulations of blood flow and oxygen transport were undertaken in various idealized AVF configurations. The objective of the study was to understand how arterial curvature affects blood flow and oxygen transport patterns within AVF, with a focus on how curvature alters metrics known to correlate with vascular pathology such as Intimal Hyperplasia (IH). If one subscribes to the hypothesis that unsteady flow causes IH within AVF, then the results suggest that in order to avoid IH, AVF should be formed via a vein graft onto the outer-curvature of a curved artery. However, if one subscribes to the hypothesis that low wall shear stress and/or low lumen-to-wall oxygen flux (leading to wall hypoxia) cause IH within AVF, then the results suggest that in order to avoid IH, AVF should be formed via a vein graft onto a straight artery, or the inner-curvature of a curved artery. We note that the recommendations are incompatible-highlighting the importance of ascertaining the exact mechanisms underlying development of IH in AVF. Nonetheless, the results clearly illustrate the important role played by arterial curvature in determining AVF hemodynamics, which to our knowledge has been overlooked in all previous studies.

Human brain microvascular endothelial cells resist elongation due to shear stress

Endothelial cells in straight sections of vessels are known to elongate and align in the direction of flow. This phenotype has been replicated in confluent monolayers of bovine aortic endothelial cells and human umbilical vein endothelial cells (HUVECs) in cell culture under physiological shear stress. Here we report on the morphological response of human brain microvascular endothelial cells (HBMECs) in confluent monolayers in response to shear stress. Using a microfluidic platform we image confluent monolayers of HBMECs and HUVECs under shear stresses up to 16 dyne cm(-2). From live-cell imaging we quantitatively analyze the cell morphology and cell speed as a function of time. We show that HBMECs do not undergo a classical transition from cobblestone to spindle-like morphology in response to shear stress. We further show that under shear stress, actin fibers are randomly oriented in the cells indicating that there is no cytoskeletal remodeling. These results suggest that HBMECs are programmed to resist elongation and alignment under shear stress, a phenotype that may be associated with the unique properties of the blood-brain barrier.

Hydrogel-based methods for engineering cellular microenvironment with spatiotemporal gradients

Natural cellular microenvironment consists of spatiotemporal gradients of multiple physical (e.g. extracellular matrix stiffness, porosity and stress/strain) and chemical cues (e.g. morphogens), which play important roles in regulating cell behaviors including spreading, proliferation, migration, differentiation and apoptosis, especially for pathological processes such as tumor formation and progression. Therefore, it is essential to engineer cellular gradient microenvironment incorporating various gradients for the fabrication of normal and pathological tissue models in vitro. In this article, we firstly review the development of engineering cellular physical and chemical gradients with cytocompatible hydrogels in both two-dimension and three-dimension formats. We then present current advances in the application of engineered gradient microenvironments for the fabrication of disease models in vitro. Finally, concluding remarks and future perspectives for engineering cellular gradients are given.

A comparative review of the hemodynamics and pathogenesis of cerebral and abdominal aortic aneurysms: lessons to learn from each other

OBJECTIVE

Cerebral aneurysms (CAs) and abdominal aortic aneurysms (AAAs) are degenerative vascular pathologies that manifest as abnormal dilations of the arterial wall. They arise with different morphologies in different types of blood vessels under different hemodynamic conditions. Although treated as different pathologies, we examine common pathways in their hemodynamic pathogenesis in order to elucidate mechanisms of formation.

MATERIALS AND METHODS

A systematic review of the literature was performed. Current concepts on pathogenesis and hemodynamics were collected and compared.

RESULTS

CAs arise as saccular dilations on the cerebral arteries of the circle of Willis under high blood flow, high wall shear stress (WSS), and high wall shear stress gradient (WSSG) conditions. AAAs arise as fusiform dilations on the infrarenal aorta under low blood flow, low, oscillating WSS, and high WSSG conditions. While at opposite ends of the WSS spectrum, they share high WSSG, a critical factor in arterial remodeling. This alone may not be enough to initiate aneurysm formation, but may ignite a cascade of downstream events that leads to aneurysm development. Despite differences in morphology and the structure, CAs and AAAs share many histopathological and biomechanical characteristics. Endothelial cell damage, loss of elastin, and smooth muscle cell loss are universal findings in CAs and AAAs. Increased matrix metalloproteinases and other proteinases, reactive oxygen species, and inflammation also contribute to the pathogenesis of both aneurysms.

CONCLUSION

Our review revealed similar pathways in seemingly different pathologies. We also highlight the need for cross-disciplinary studies to aid in finding similarities between pathologies.

Piezo1 integration of vascular architecture with physiological force

The mechanisms by which physical forces regulate endothelial cells to determine the complexities of vascular structure and function are enigmatic. Studies of sensory neurons have suggested Piezo proteins as subunits of Ca(2+)-permeable non-selective cationic channels for detection of noxious mechanical impact. Here we show Piezo1 (Fam38a) channels as sensors of frictional force (shear stress) and determinants of vascular structure in both development and adult physiology. Global or endothelial-specific disruption of mouse Piezo1 profoundly disturbed the developing vasculature and was embryonic lethal within days of the heart beating. Haploinsufficiency was not lethal but endothelial abnormality was detected in mature vessels. The importance of Piezo1 channels as sensors of blood flow was shown by Piezo1 dependence of shear-stress-evoked ionic current and calcium influx in endothelial cells and the ability of exogenous Piezo1 to confer sensitivity to shear stress on otherwise resistant cells. Downstream of this calcium influx there was protease activation and spatial reorganization of endothelial cells to the polarity of the applied force. The data suggest that Piezo1 channels function as pivotal integrators in vascular biology.

Exacerbation of intracranial aneurysm and aortic dissection in hypertensive rat treated with the prostaglandin F-receptor antagonist AS604872

Intracranial aneurysm (IA) and aortic dissection are both complications of hypertension and characterized by degeneration of the media. Given the involvement of prostaglandin F2α and its receptor, FP, in extracellular matrix remodeling in a mouse model of pulmonary fibrosis, here we induced hypertension and IA in rats by salt loading and hemi-lateral ligation of renal and carotid arteries and examined effects of a selective FP antagonist, AS604872, on these vascular events. AS604872 significantly accelerated degeneration of the media in both cerebral artery and aorta as evidenced by thinning of the media and disruption of the elastic lamina and promoted IA and aortic dissection. Notably, AS604872 induced expression of pro-inflammatory genes such as E-selectin in lesions and significantly enhanced macrophage infiltration. Suppression of surface expression of E-selectin with cimetidine prevented macrophage infiltration and aortic dissection. Thus, AS604872 exacerbates vascular inflammation in hypertensive rats and facilitates IA and aortic dissection. These results demonstrate that both IA and aortic dissection are caused by chronic inflammation of the arterial wall, which is worsened by AS604872, cautioning that other FP antagonists may share such deleterious actions in vascular homeostasis and suggesting that AS604872 can be used to make models of these vascular diseases with extensive degeneration.

Rho/ROCK signal cascade mediates asymmetric dimethylarginine-induced vascular smooth muscle cells migration and phenotype change

Asymmetric dimethylarginine (ADMA) induces vascular smooth muscle cells (VSMCs) migration. VSMC phenotype change is a prerequisite of migration. RhoA and Rho-kinase (ROCK) mediate migration of VSMCs. We hypothesize that ADMA induces VSMC migration via the activation of Rho/ROCK signal pathway and due to VSMCs phenotype change. ADMA activates Rho/ROCK signal pathway that interpreted by the elevation of RhoA activity and phosphorylation level of a ROCK substrate. Pretreatment with ROCK inhibitor, Y27632 completely reverses the induction of ADMA on ROCK and in turn inhibits ADMA-induced VSMCs migration. When the Rho/ROCK signal pathway has been blocked by pretreatment with Y27632, the induction of ERK signal pathway by ADMA is completely abrogated. Elimination of ADMA via overexpression of dimethylarginine dimethylaminohydrolase 2 (DDAH2) and L-arginine both blocks the effects of ADMA on the activation of Rho/ROCK and extra cellular signal-regulated kinase (ERK) in VSMCs. The expression of differentiated phenotype relative proteins was reduced and the actin cytoskeleton was disassembled by ADMA, which were blocked by Y27632, further interpreting that ADMA inducing VSMCs migration via Rho/ROCK signal pathway is due to its effect on the VSMCs phenotype change. Our present study may help to provide novel insights into the therapy and prevention of atherosclerosis.

In vitro shear stress measurements using particle image velocimetry in a family of carotid artery models: effect of stenosis severity, plaque eccentricity, and ulceration

Atherosclerotic disease, and the subsequent complications of thrombosis and plaque rupture, has been associated with local shear stress. In the diseased carotid artery, local variations in shear stress are induced by various geometrical features of the stenotic plaque. Greater stenosis severity, plaque eccentricity (symmetry) and plaque ulceration have been associated with increased risk of cerebrovascular events based on clinical trial studies. Using particle image velocimetry, the levels and patterns of shear stress (derived from both laminar and turbulent phases) were studied for a family of eight matched-geometry models incorporating independently varied plaque features - i.e. stenosis severity up to 70%, one of two forms of plaque eccentricity, and the presence of plaque ulceration). The level of laminar (ensemble-averaged) shear stress increased with increasing stenosis severity resulting in 2-16 Pa for free shear stress (FSS) and approximately double (4-36 Pa) for wall shear stress (WSS). Independent of stenosis severity, marked differences were found in the distribution and extent of shear stress between the concentric and eccentric plaque formations. The maximum WSS, found at the apex of the stenosis, decayed significantly steeper along the outer wall of an eccentric model compared to the concentric counterpart, with a 70% eccentric stenosis having 249% steeper decay coinciding with the large outer-wall recirculation zone. The presence of ulceration (in a 50% eccentric plaque) resulted in both elevated FSS and WSS levels that were sustained longer (∼20 ms) through the systolic phase compared to the non-ulcerated counterpart model, among other notable differences. Reynolds (turbulent) shear stress, elevated around the point of distal jet detachment, became prominent during the systolic deceleration phase and was widely distributed over the large recirculation zone in the eccentric stenoses.

Endothelial nitric oxide synthase and superoxide mediate hemodynamic initiation of intracranial aneurysms

Background

Hemodynamic insults at arterial bifurcations are believed to play a critical role in initiating intracranial aneurysms. Recent studies in a rabbit model indicate that aneurysmal damage initiates under specific wall shear stress conditions when smooth muscle cells (SMCs) become pro-inflammatory and produce matrix metalloproteinases (MMPs). The mechanisms leading to SMC activation and MMP production during hemodynamic aneurysm initiation are unknown. The goal is to determine if nitric oxide and/or superoxide induce SMC changes, MMP production and aneurysmal remodeling following hemodynamic insult.

Methods

Bilateral common carotid artery ligation was performed on rabbits (n = 19, plus 5 sham operations) to induce aneurysmal damage at the basilar terminus. Ligated animals were treated with the nitric oxide synthase (NOS) inhibitor LNAME (n = 7) or the superoxide scavenger TEMPOL (n = 5) and compared to untreated animals (n = 7). Aneurysm development was assessed histologically 5 days after ligation. Changes in NOS isoforms, peroxynitrite, reactive oxygen species (ROS), MMP-2, MMP-9, and smooth muscle α-actin were analyzed by immunohistochemistry.

Results

LNAME attenuated ligation-induced IEL loss, media thinning and bulge formation. In untreated animals, immunofluorescence showed increased endothelial NOS (eNOS) after ligation, but no change in inducible or neuronal NOS. Furthermore, during aneurysm initiation ROS increased in the media, but not the intima, and there was no change in peroxynitrite. In LNAME-treated animals, ROS production did not change. Together, this suggests that eNOS is important for aneurysm initiation but not by producing superoxide. TEMPOL treatment reduced aneurysm development, indicating that the increased medial superoxide is also necessary for aneurysm initiation. LNAME and TEMPOL treatment in ligated animals restored α-actin and decreased MMPs, suggesting that eNOS and superoxide both lead to SMC de-differentiation and MMP production.

Conclusion

Aneurysm-inducing hemodynamics lead to increased eNOS and superoxide, which both affect SMC phenotype, increasing MMP production and aneurysmal damage.

Extent of flow recirculation governs expression of atherosclerotic and thrombotic biomarkers in arterial bifurcations

AIMS

Atherogenesis, evolution of plaque, and outcomes following endovascular intervention depend heavily on the unique vascular architecture of each individual. Patient-specific, multiscale models able to correlate changes in microscopic cellular responses with relevant macroscopic flow, and structural conditions may help understand the progression of occlusive arterial disease, providing insights into how to mitigate adverse responses in specific settings and individuals.

METHODS AND RESULTS

Vascular architectures mimicking coronary and carotid bifurcations were derived from clinical imaging and used to generate conjoint computational meshes for in silico analysis and biocompatible scaffolds for in vitro models. In parallel with three-dimensional flow simulations, geometrically realistic scaffolds were seeded with human smooth muscle cells (SMC) or endothelial cells and exposed to relevant, physiological flows. In vitro surrogates of endothelial health, atherosclerotic progression, and thrombosis were locally quantified and correlated best with an quantified extent of flow recirculation occurring within the bifurcation models. Oxidized low-density lipoprotein uptake, monocyte adhesion, and tissue factor expression locally rose up to three-fold, and phosphorylated endothelial nitric oxide synthase and Krüppel-like factor 2 decreased up to two-fold in recirculation areas. Isolated testing in straight-tube idealized constructs subject to static, oscillatory, and pulsatile conditions, indicative of different recirculant conditions corroborated these flow-mediated dependencies.

CONCLUSIONS

Flow drives variations in vascular reactivity and vascular beds. Endothelial health was preserved by arterial flow but jeopardized in regions of flow recirculation in a quasi-linear manner. Similarly, SMC exposed to flow were more thrombogenic in large recirculating regions. Health, thrombosis, and atherosclerosis biomarkers correlate with the extent of recirculation in vascular cells lining certain vascular geometries.

Regional variation of wall shear stress in ascending thoracic aortic aneurysms

The development of an ascending thoracic aortic aneurysm is likely caused by excessive hemodynamic loads exerted on the aneurysmal wall. Computational fluid-dynamic analyses were performed on patient-specific ascending thoracic aortic aneurysms obtained from patients with either bicuspid aortic valve or tricuspid aortic valve to evaluate hemodynamic and wall shear parameters, imparting aneurysm enlargement. Results showed an accelerated flow along the outer aortic wall with helical flow in the aneurysm center for bicuspid aortic valve ascending thoracic aortic aneurysms. In a different way, tricuspid aortic valve ascending thoracic aortic aneurysms exhibited normal systolic flow without substantial secondary pattern. Analysis of wall shear parameters evinced a high and locally varying wall shear stress on the outer aortic wall and high temporal oscillations in wall shear stress (oscillatory shear index) on either left or right side of aneurysmal aorta. These findings may explain the asymmetric dilatation typically observed in ascending thoracic aortic aneurysms. Simulations of a hypertensive scenario revealed an increase in wall shear stress upon 44% compared to normal systemic pressure models. Computational fluid-dynamics–based analysis may allow identification of wall shear parameters portending aneurysm dilatation and hence guide preventative intervention.

A multiple-channel, multiple-assay platform for characterization of full-range shear stress effects on vascular endothelial cells

Vascular endothelial cells (VECs), which line blood vessels and are key to understanding pathologies and treatments of various diseases, experience highly variable wall shear stress (WSS) in vivo (1-60 dyn cm(-2)), imposing numerous effects on physiological and morphological functions. Previous flow-based systems for studying these effects have been limited in range, and comprehensive information on VEC functions at the full spectrum of WSS has not been available yet. To allow rapid characterization of WSS effects, we developed the first multiple channel microfluidic platform that enables a wide range (~15×) of homogeneous WSS conditions while simultaneously allowing trans-monolayer assays, such as permeability and trans-endothelial electrical resistance (TEER) assays, as well as cell morphometry and protein expression assays. Flow velocity/WSS distributions between channels were predicted with COMSOL simulations and verified by measurement using an integrated microflow sensor array. Biomechanical responses of the brain microvascular endothelial cell line bEnd.3 to the full natural spectrum of WSS were investigated with the platform. Under increasing WSS conditions ranging from 0 to 86 dyn cm(-2), (1) permeabilities of FITC-conjugated dextran and propidium iodide decreased, respectively, at rates of 4.06 × 10(-8) and 6.04 × 10(-8) cm s(-1) per dyn cm(-2); (2) TEER increased at a rate of 0.8 Ω cm(2) per dyn cm(-2); (3) increased alignment of cells along the flow direction under increasing WSS conditions; and finally (4) increased protein expression of both the tight junction component ZO-1 (~5×) and the efflux transporter P-gp (~6×) was observed at 86 dyn cm(-2) compared to static controls via western blot. We conclude that the presented microfluidic platform is a valid approach for comprehensively assaying cell responses to fluidic WSS.

Critical role of TNF-alpha-TNFR1 signaling in intracranial aneurysm formation

Background

Intracranial aneurysm (IA) is a socially important disease due to its high incidence in the general public and the severity of resultant subarachnoid hemorrhage that follows rupture. Despite the social importance of IA as a cause of subarachnoid hemorrhage, there is no medical treatment to prevent rupture, except for surgical procedures, because the mechanisms regulating IA formation are poorly understood. Therefore, these mechanisms should be elucidated to identify a therapeutic target for IA treatment. In human IAs, the presence of inflammatory responses, such as an increase of tumor necrosis factor (TNF)-alpha, have been observed, suggesting a role for inflammation in IA formation. Recent investigations using rodent models of IAs have revealed the crucial role of inflammatory responses in IA formation, supporting the results of human studies. Thus, we identified nuclear factor (NF)-kappaB as a critical mediator of inflammation regulating IA formation, by inducing downstream pro-inflammatory genes such as MCP-1, a chemoattractant for macrophages, and COX-2. In this study, we focused on TNF-alpha signaling as a potential cascade that regulates NF-kappaB-mediated IA formation.

Results

We first confirmed an increase in TNF-alpha content in IA walls during IA formation, as expected based on human studies. Consistently, the activity of TNF-alpha converting enzyme (TACE), an enzyme responsible for TNF-alpha release, was induced in the arterial walls after aneurysm induction in a rat model. Next, we subjected tumor necrosis factor receptor superfamily member 1a (TNFR1)-deficient mice to the IA model to clarify the contribution of TNF-alpha-TNFR1 signaling to pathogenesis, and confirmed significant suppression of IA formation in TNFR1-deficient mice. Furthermore, in the IA walls of TNFR1-deficient mice, inflammatory responses, including NF-kappaB activation, subsequent expression of MCP-1 and COX-2, and infiltration of macrophages into the IA lesion, were greatly suppressed compared with those in wild-type mice.

Conclusions

In this study, using rodent models of IAs, we clarified the crucial role of TNF-alpha-TNFR1 signaling in the pathogenesis of IAs by inducing inflammatory responses, and propose this signaling as a potential therapeutic target for IA treatment.

Aneurysmal remodeling in the circle of Willis after carotid occlusion in an experimental model

Carotid occlusions are associated with de novo intracranial aneurysm formation in clinical case reports, but this phenomenon is not widely studied. We performed bilateral carotid ligation (n=9) in rabbits to simulate carotid occlusion, and sham surgery (n=3) for control. Upon euthanasia (n=3 at 5 days, n=6 at 6 months post ligation, and n=3 at 5 days after sham operation), vascular corrosion casts of the circle of Willis (CoW) were created. Using scanning electron microscopy, we quantified gross morphologic, macroscopic, and microscopic changes on the endocasts and compared findings with histologic data. At 5 days, CoW arteries of ligated animals increased caliber. The posterior communicating artery (PCom) increased length and tortuosity, and the ophthalmic artery (OA) origin presented preaneurysmal bulges. At 6 months, calibers were unchanged from 5 days, PComs further increased tortuosity while presenting segmental dilations, and the OA origin and basilar terminus presented preaneurysmal bulges. This exploratory study provides evidence that flow increase after carotid occlusion produces both compensatory arterial augmentation and pathologic remodeling such as tortuosity and saccular/fusiform aneurysm. Our findings may have considerable clinical implications, as these lesser-known consequences should be considered when managing patients with carotid artery disease or choosing carotid ligation as a therapeutic option.

Propose a wall shear stress divergence to estimate the risks of intracranial aneurysm rupture

Although wall shear stress (WSS) has long been considered a critical indicator of intracranial aneurysm rupture, there is still no definite conclusion as to whether a high or a low WSS results in aneurysm rupture. The reason may be that the effect of WSS direction has not been fully considered. The objectives of this study are to investigate the magnitude of WSS (|WSS|) and its divergence on the aneurysm surface and to test the significance of both in relation to the aneurysm rupture. Patient-specific computational fluid dynamics (CFD) was used to compute WSS and wall shear stress divergence (WSSD) on the aneurysm surface for nineteen patients. Our results revealed that if high |WSS| is stretching aneurysm luminal surface, and the stretching region is concentrated, the aneurysm is under a high risk of rupture. It seems that, by considering both direction and magnitude of WSS, WSSD may be a better indicator for the risk estimation of aneurysm rupture (154).

Differential gene expression by endothelial cells under positive and negative streamwise gradients of high wall shear stress

Flow impingement at arterial bifurcations causes high frictional force [or wall shear stress (WSS)], and flow acceleration and deceleration in the branches create positive and negative streamwise gradients in WSS (WSSG), respectively. Intracranial aneurysms tend to form in regions with high WSS and positive WSSG. However, little is known about the responses of endothelial cells (ECs) to either positive or negative WSSG under high WSS conditions. We used cDNA microarrays to profile gene expression in cultured ECs exposed to positive or negative WSSG for 24 h in a flow chamber where WSS varied between 3.5 and 28.4 Pa. Gene ontology and biological pathway analysis indicated that positive WSSG favored proliferation, apoptosis, and extracellular matrix processing while decreasing expression of proinflammatory genes. To determine if similar responses occur in vivo, we examined EC proliferation and expression of the matrix metalloproteinase ADAMTS1 under high WSS and WSSG created at the basilar terminus of rabbits after bilateral carotid ligation. Precise hemodynamic conditions were determined by computational fluid dynamic simulations from three-dimensional angiography and mapped on immunofluorescence staining for the proliferation marker Ki-67 and ADAMTS1. Both proliferation and ADAMTS1 were significantly higher in ECs under positive WSSG than in adjacent regions of negative WSSG. Our results indicate that WSSG elicits distinct EC gene expression profiles and particular biological pathways including increased cell proliferation and matrix processing. Such EC responses may be important in understanding the mechanisms of intracranial aneurysm initiation at regions of high WSS and positive WSSG.