Analysis of early embryonic great-vessel microcirculation in zebrafish using high-speed confocal μPIV

In the developing cardiovascular system, hemodynamic vascular loading is critical for angiogenesis and cardiovascular adaptation. Normal zebrafish embryos with transgenically-labeled endothelial and red blood cells provide an excellent in vivo model for studying the fluid-flow induced vascular loading. To characterize the developmental hemodynamics of early embryonic great-vessel microcirculation in the zebrafish embryo, two complementary studies (experimental and numerical) are presented. Quantitative comparison of the wall shear stress (WSS) at the first aortic arch (AA1) of wild-type zebrafish embryos during two consecutive developmental stages is presented, using time-resolved confocal micro-particle image velocimetry (μPIV). Analysis showed that there was significant WSS difference between 32 and 48 h post-fertilization (hpf) wild-type embryos, which correlates with normal arch morphogenesis. The vascular distensibility of the arch wall at systole and the acceleration/deceleration rates of time-lapse phase-averaged streamwise blood flow curves were also analyzed. To estimate the influence of a novel intermittent red-blood cell (RBC) loading on the endothelium, a numerical two-phase, volume of fluid (VOF) flow model was further developed with realistic in vivo conditions. These studies showed that near-wall effects and cell clustering increased WSS augmentation at a minimum of 15% when the distance of RBC from arch vessel wall was less than 3 μm or when RBC cell-to-cell distance was less than 3 μm. When compared to a smooth wall, the WSS augmentation increased by a factor of ~1.4 due to the roughness of the wall created by the endothelial cell profile. These results quantitatively highlight the contribution of individual RBC flow patterns on endothelial WSS in great-vessel microcirculation and will benefit the quantitative understanding of mechanotransduction in embryonic great vessel biology, including arteriovenous malformations (AVM).     

Numerical simulation of pulsatile flow in a compliant curved tube model of a coronary artery

The endothelial cells (ECs) lining a blood vessel wall are exposed to both the wall shear stress (WSS) of blood flow and the circumferential strain (CS) of pulsing artery wall motion. These two forces and their interaction are believed to play a role in determining remodeling of the vessel wall and development of arterial disease (atherosclerosis). This study focused on the WSS and CS dynamic behavior in a compliant model of a coronary artery taking into account the curvature of the bending artery and physiological radial wall motion. A three-dimensional finite element model with transient flow and moving boundaries was set up to simulate pulsatile flow with physiological pressure and flow wave forms characteristic of the coronary arteries. The characteristic coronary artery curvature and flow conditions applied to the simulation were: aspect ratio (lambda) = 10, diameter variation (DV) = 6 percent, mean Reynolds number (Re) = 150, and unsteadiness parameter (alpha) = 3. The results show that mean WSS is about 50 percent lower on the inside wall than the outside wall while WSS oscillation is stronger on the inside wall. The stress phase angle (SPA) between CS and WSS, which characterizes the dynamics of the mechanical force pattern applied to the endothelial cell layer, shows that CS and WSS are more out of phase in the coronaries than in any other region of the circulation (-220 deg on the outside wall, -250 deg on the inside wall). This suggests that in addition to WSS, SPA may play a role in localization of coronary atherosclerosis.     

An introduction to biofluid mechanics--basic models and applications

Cardiovascular disease is the primary cause of morbidity and mortality in the western world. Complex hemodynamics play a critical role in the development of atherosclerosis and the processes of aging, as well as many other disease processes. Biofluid mechanics play a major role in the cardiovascular system and it is important to understand the forces and movement of blood cells and whole blood as well as the interaction between blood cells and the vessel wall. Fundamental fluid mechanical, which are important for the understanding of the blood flow in the cardiovascular circulatory system of the human body aspects are presented. Measurement techniques for model studies such as LDA, ultrasound, and MRI studies will be discussed. Viscosity and flow behavior changes specifically the creation of vortices and flow disturbances can be used to show how medication can influence flow behavior. Experiments have shown that hemodynamics may have a strong influence on the creation of aneurysms and varicose veins. Other factors such as vessel wall structure are also important. In preliminary studies, it has been demonstrated that geometry and elasticity of vessel walls help determine flow behavior. High velocity fluctuations indicate flow disturbances that should be avoided. Health care practitioners must understand fluid dynamic factors such as flow rate ratio, pressure and velocity gradients, and flow behavior, velocity distribution, shear stress on the wall and on blood cells. These mechanical factors are largely responsible for the deposit of blood cells and lipids, a leading cause of atherosclerosis. The interaction between blood cells and of the cells with the vessel, leads to the formation of plaques and agglomerations. These deposits are found predominantly at arterial bends and bifurcations where blood flow is disturbed, where a secondary flow is created, and where flow separation regions are found. Experiments on hemodynamic effects in elastic silicon rubber models of the cardiovascular system with flow wire, stents, or patches for vessel surgery will be discussed. These studies can be important in improving diagnostics and therapeutic applications.     

Computational study of pulsatile blood flow in prototype vessel geometries of coronary segments

The spatial and temporal distributions of wall shear stress (WSS) in prototype vessel geometries of coronary segments are investigated via numerical simulation, and the potential association with vascular disease and specifically atherosclerosis and plaque rupture is discussed. In particular, simulation results of WSS spatio-temporal distributions are presented for pulsatile, non-Newtonian blood flow conditions for: (a) curved pipes with different curvatures, and (b) bifurcating pipes with different branching angles and flow division. The effects of non-Newtonian flow on WSS (compared to Newtonian flow) are found to be small at Reynolds numbers representative of blood flow in coronary arteries. Specific preferential sites of average low WSS (and likely atherogenesis) were found at the outer regions of the bifurcating branches just after the bifurcation, and at the outer-entry and inner-exit flow regions of the curved vessel segment. The drop in WSS was more dramatic at the bifurcating vessel sites (less than 5% of the pre-bifurcation value). These sites were also near rapid gradients of WSS changes in space and time – a fact that increases the risk of rupture of plaque likely to develop at these sites. The time variation of the WSS spatial distributions was very rapid around the start and end of the systolic phase of the cardiac cycle, when strong fluctuations of intravascular pressure were also observed. These rapid and strong changes of WSS and pressure coincide temporally with the greatest flexion and mechanical stresses induced in the vessel wall by myocardial motion (ventricular contraction). The combination of these factors may increase the risk of plaque rupture and thrombus formation at these sites.     

A Eulerian method to analyze wall shear stress fixed points and manifolds in cardiovascular flows

Biomechanics and Modeling in Mechanobiology - Based upon dynamical systems theory, a fixed point of a vector field such as the wall shear stress (WSS) at the luminal surface of a vessel is a point...     

Numerical investigation of the non-Newtonian pulsatile blood flow in a bifurcation model with a non-planar branch

The pulsatile flow of non-Newtonian fluid in a bifurcation model with a non-planar daughter branch is investigated numerically by using the Carreau-Yasuda model to take into account the shear thinning behavior of the analog blood fluid. The objective of this study is to deal with the influence of the non-Newtonian property of fluid and of out-of-plane curvature in the non-planar daughter vessel on wall shear stress (WSS), oscillatory shear index (OSI), and flow phenomena during the pulse cycle. The non-Newtonian property in the daughter vessels induces a flattened axial velocity profile due to its shear thinning behavior. The non-planarity deflects flow from the inner wall of the vessel to the outer wall and changes the distribution of WSS along the vessel, in particular in systole phase. Downstream of the bifurcation, the velocity profiles are shifted toward the flow divider, and low WSS and high shear stress temporal oscillations characterized by OSI occur on the outer wall region of the daughter vessels close to the bifurcation. Secondary motions become stronger with the addition of the out-of-plane curvature induced by the bending of the vessel, and the secondary flow patterns swirl along the non-planar daughter vessel. A significant difference between the non-Newtonian and the Newtonian pulsatile flow is revealed during the pulse cycle; however, reasonable agreement between the non-Newtonian and the rescaled Newtonian flow is found. Calculated results for the pulsatile flow support the view that the non-planarity of blood vessels and the non-Newtonian properties of blood are an important factor in hemodynamics and may play a significant role in vascular biology and pathophysiology.     

Numerical investigation of the non-Newtonian blood flow in a bifurcation model with a non-planar branch

The non-Newtonian fluid flow in a bifurcation model with a non-planar daughter branch is investigated by using finite element method to solve the three-dimensional Navier–Stokes equations coupled with a non-Newtonian constitutive model, in which the shear thinning behavior of the blood fluid is incorporated by the Carreau–Yasuda model. The objective of this study is to investigate the influence of the non-Newtonian property of fluid as well as of curvature and out-of-plane geometry in the non-planar daughter vessel on wall shear stress (WSS) and flow phenomena. In the non-planar daughter vessel, the flows are typified by the skewing of the velocity profile towards the outer wall, creating a relatively low WSS at the inner wall. In the downstream of the bifurcation, the velocity profiles are shifted towards the flow divider. The low WSS is found at the inner walls of the curvature and the lateral walls of the bifurcation. Secondary flow patterns that swirl fluid from the inner wall of curvature to the outer wall in the middle of the vessel are also well documented for the curved and bifurcating vessels. The numerical results for the non-Newtonian fluid and the Newtonian fluid with original Reynolds number and the corresponding rescaled Reynolds number are presented. Significant difference between the non-Newtonian flow and the Newtonian flow is revealed; however, reasonable agreement between the non-Newtonian flow and the rescaled Newtonian flow is found. Results of this study support the view that the non-planarity of blood vessels and the non-Newtonian properties of blood are an important factor in hemodynamics and may play a significant role in vascular biology and pathophysiology.     

Effects of polycythemia on systemic endothelial function in chronic hypoxic lung disease

Chronic obstructive pulmonary disease (COPD) is a major risk factor for cardiovascular disease. Polycythemia, a common complication of hypoxic COPD, may affect systemic vascular function by altering blood viscosity, vessel wall shear stress (WSS), and endothelium-derived nitric oxide (NO) release. Here, we evaluated the effects of hypoxia-related polycythemia on systemic endothelial function in patients with COPD. We investigated blood viscosity, WSS, and endothelial function in 15 polycythemic and 13 normocythemic patients with COPD of equal severity, by recording brachial artery diameter variations in response to hyperemia and by using venous occlusion plethysmography (VOP) to measure forearm blood flow (FBF) responses to a brachial artery infusion of acetylcholine (ACh), bradykinin (BK), sodium nitroprusside (SNP), substance P (SP), isoptin, and N-monomethyl-L-arginine (L-NMMA). At baseline, polycythemic patients had higher blood viscosity and larger brachial artery diameter than normocythemic patients but similar calculated WSS. Flow-mediated brachial artery vasodilation was increased in the polycythemic patients, in proportion to the hemoglobin levels. ACh-induced vasodilation was markedly impaired in the polycythemic patients and negatively correlated with hemoglobin levels. FBF responses to endothelium- (BK, SP) and non-endothelium-dependent (SNP, isoptin) vasodilators were not significantly different between the two groups. L-NMMA infusion induced a similar vasoconstrictor response in both groups, in accordance with their similar baseline WSS. In conclusion, systemic arteries in polycythemic patients adjust appropriately to chronic or acute WSS elevations by appropriate basal and stimulated NO release. Overall, our results suggest that moderate polycythemia has no adverse effect on vascular function in COPD.     

Influence of non-Newtonian properties of blood on the wall shear stress in human atherosclerotic right coronary arteries

The objective of this work is to investigate the effect of non-Newtonian properties of blood on the wall shear stress (WSS) in atherosclerotic coronary arteries using both Newtonian and non-Newtonian models. Numerical simulations were performed to examine how the spatial and temporal WSS distributions are influenced by the stenosis size, blood viscosity, and flow rate. The computational results demonstrated that blood viscosity properties had considerable effect on the magnitude of the WSS, especially where disturbed flow was observed. The WSS distribution is highly non-uniform both temporally and spatially, especially in the stenotic region. The maximum WSS occurred at the proximal side of the stenosis, near the outer wall in the curved artery with no stenosis. The lumen area near the inner wall distal to the stenosis region experienced a lower WSS during the entire cardiac cycle. Among the factors of stenosis size, blood viscosity, and flow rate, the size of the stenosis has the most significant effect on the spatial and temporal WSS distributions qualitatively and quantitatively.     

Flow shear stress regulates endothelial barrier function and expression of angiogenic factors in a 3D microfluidic tumor vascular model

Endothelial cells lining blood vessels are exposed to various hemodynamic forces associated with blood flow. These include fluid shear, the tangential force derived from the friction of blood flowing across the luminal cell surface, tensile stress due to deformation of the vessel wall by transvascular flow, and normal stress caused by the hydrodynamic pressure differential across the vessel wall. While it is well known that these fluid forces induce changes in endothelial morphology, cytoskeletal remodeling, and altered gene expression, the effect of flow on endothelial organization within the context of the tumor microenvironment is largely unknown. Using a previously established microfluidic tumor vascular model, the objective of this study was to investigate the effect of normal (4 dyn/cm²), low (1 dyn/cm²), and high (10 dyn/cm²) microvascular wall shear stress (WSS) on tumor-endothelial paracrine signaling associated with angiogenesis. It is hypothesized that high WSS will alter the endothelial phenotype such that vascular permeability and tumor-expressed angiogenic factors are reduced. Results demonstrate that endothelial permeability decreases as a function of increasing WSS, while co-culture with tumor cells increases permeability relative to mono-cultures. This response is likely due to shear stress-mediated endothelial cell alignment and tumor-VEGF-induced permeability. In addition, gene expression analysis revealed that high WSS (10 dyn/cm²) significantly down-regulates tumor-expressed MMP9, HIF1, VEGFA, ANG1, and ANG2, all of which are important factors implicated in tumor angiogenesis. This result was not observed in tumor mono-cultures or static conditioned media experiments, suggesting a flow-mediated paracrine signaling mechanism exists with surrounding tumor cells that elicits a change in expression of angiogenic factors. Findings from this work have significant implications regarding low blood velocities commonly seen in the tumor vasculature, suggesting high shear stress-regulation of angiogenic activity is lacking in many vessels, thereby driving tumor angiogenesis.     

Circumferential vascular deformation after stent implantation alters wall shear stress evaluated with time-dependent 3D computational fluid dynamics models.

The success of vascular stents in the restoration of blood flow is limited by restenosis. Recent data generated from computational fluid dynamics (CFD) models suggest that stent geometry may cause local alterations in wall shear stress (WSS) that have been associated with neointimal hyperplasia and subsequent restenosis. However, previous CFD studies have ignored histological evidence of vascular straightening between circumferential stent struts. We tested the hypothesis that consideration of stent-induced vascular deformation may more accurately predict alterations in indexes of WSS that may subsequently account for histological findings after stenting. We further tested the hypothesis that the severity of these alterations in WSS varies with the degree of vascular deformation after implantation. Steady-state and time-dependent simulations of three-dimensional CFD arteries based on canine coronary artery measurements of diameter and blood flow were conducted, and WSS and WSS gradients were calculated. Circumferential straightening introduced areas of high WSS between stent struts that were absent in stented vessels of circular cross section. The area of vessel exposed to low WSS was dependent on the degree of circumferential vascular deformation and axial location within the stent. Stents with four vs. eight struts increased the intrastrut area of low WSS in vessels, regardless of cross-sectional geometry. Elevated WSS gradients were also observed between struts in vessels with polygonal cross sections. The results obtained using three-dimensional CFD models suggest that changes in vascular geometry after stent implantation are important determinants of WSS distributions that may be associated with subsequent neointimal hyperplasia.     

Numerical investigation of physiologically realistic pulsatile flow through arterial stenosis

Numerical simulations of pulsatile blood flow in straight tube stenosis models were performed to investigate the poststenotic flow phenomena. In this study, three axisymmetrical and three asymmetrical stenosis models with area reduction of 25%, 50% and 75% were constructed. A measured human common carotid artery blood flow waveform was used as the upstream flow condition which has a mean Reynold's number of 300. All calculations were performed with high spatial and temporal resolutions. Flow features such as velocity profiles, flow separation zone (FSZ), and wall shear stress (WSS) distributions in the poststenotic region for all models are presented. The results have demonstrated that the formation and development of FSZs in the poststenotic region are very complex, especially in the flow deceleration phase. In axisymmetric stenoses the poststenotic flow is more sensitive to changes in the degree of stenosis than in asymmetric models. For severe stenoses, the stenosis influence length is shorter in asymmetrical models than in axisymmetrical cases. WSS oscillations (between positive and negative values) have been observed at various downstream locations in some models. The amplitude of the oscillation depends strongly on the axial location and the degree of stenosis.     

Numerical modeling of simulated blood flow in idealized composite arterial coronary grafts: steady state simulations

This paper presents a comparative study of simulated blood flow in different configurations of simplified composite arterial coronary grafts (CACGs). Even though the composite arterial grafting is increasingly used in cardiac surgery, it is still questionable whether or not the blood flow in such grafts can adequately meet the demands of the native myocardial circulation. A computational fluid dynamics (CFD) model was developed to conduct computer-based studies of simulated blood flow in four different geometric configurations of CACGs, corresponding to routinely used networks in cardiac surgery coronary grafts (T, Y, Pi and sequential). The flow was assumed three-dimensional, laminar and steady and the fluid as Newtonian, while the vessel walls were considered as inelastic and impermeable. It was concluded that local haemodynamics, practically described by velocity, pressure drop, wall shear stress (WSS) and flow rates, may be strongly influenced by the local geometry, especially at the anastomotic sites. The computations were made at mean flow rates of 37.5, 75 and 150ml/min. The side-branch outflow rates, computed for each bypass graft, showed noticeable differences. The results, which were found both qualitatively and quantitatively consistent with other studies, indicate that the Pi-graft exhibits significantly less uniform distribution of outflow rates than the other geometric configurations. Moreover, prominent variations in WSS and velocity distribution among the assessed CACGs were predicted, showing remarkable flow interactions among the arterial branches. The lowest shear stress regions were found on the lateral walls of bifurcations, which are predominantly susceptible to the occurrence of coronary artery disease (CAD). In contrast, the highest WSS were observed at the turn of the arterial branches.     

Alterations in regional vascular geometry produced by theoretical stent implantation influence distributions of wall shear stress: analysis of a curved coronary artery using 3D computational fluid dynamics modeling

Background  

The success of stent implantation in the restoration of blood flow through areas of vascular narrowing is limited by restenosis. Several recent studies have suggested that the local geometric environment created by a deployed stent may influence regional blood flow characteristics and alter distributions of wall shear stress (WSS) after implantation, thereby rendering specific areas of the vessel wall more susceptible to neointimal hyperplasia and restenosis. Stents are most frequently implanted in curved vessels such as the coronary arteries, but most computational studies examining blood flow patterns through stented vessels conducted to date use linear, cylindrical geometric models. It appears highly probable that restenosis occurring after stent implantation in curved arteries also occurs as a consequence of changes in fluid dynamics that are established immediately after stent implantation.     

Epigenetic response of endothelial cells to different wall shear stress magnitudes: A report of new mechano-miRNAs

Endothelial cells (ECs) respond to flow stress via a variety of mechanisms, leading to various intracellular responses that can modulate the vessel wall and lead to diseases if the flow is disturbed. Mechano-microRNAs (miRNAs) are a subset of miRNAs in the ECs that are flow responsive. Mechano-miRNAs were shown to be related to atherosclerosis pathophysiology, and a number of them were identified as pathologic. Here, we exposed human carotid ECs to different wall shear stresses (WSS), high and low, and evaluated the response of miRNAs by microarray and quantitative polymerase chain reaction analysis. We discovered five new mechano-miRNAs that were not reported in that context previously to the best of our knowledge. Moreover, functional pathway analysis revealed that under low WSS conditions, several pathways regulating apoptosis are affected. In addition, KLF2 and KLF4, known atheroprotective genes, were downregulated under low WSS and upregulated under high WSS. KLF2 and VCAM1, both angiogenic, were upregulated under high WSS. NOS3, which is vascular protective, was also upregulated with higher WSS. On the contrary, ICAM-1 and E-selectin, both atherogenic and proinflammatory, were upregulated with high WSS. Collectively, the epigenetic landscape with the gene expression analysis reveals that low WSS is associated with a proapoptotic state, while high WSS is associated with a proliferative and proinflammatory state.     

Mis-sizing of stent promotes intimal hyperplasia: impact of endothelial shear and intramural stress

Stent can cause flow disturbances on the endothelium and compliance mismatch and increased stress on the vessel wall. These effects can cause low wall shear stress (WSS), high wall shear stress gradient (WSSG), oscillatory shear index (OSI), and circumferential wall stress (CWS), which may promote neointimal hyperplasia (IH). The hypothesis is that stent-induced abnormal fluid and solid mechanics contribute to IH. To vary the range of WSS, WSSG, OSI, and CWS, we intentionally mismatched the size of stents to that of the vessel lumen. Stents were implanted in coronary arteries of 10 swine. Intravascular ultrasound (IVUS) was used to size the coronary arteries and stents. After 4 wk of stent implantation, IVUS was performed again to determine the extent of IH. In conjunction, computational models of actual stents, the artery, and non-Newtonian blood were created in a computer simulation to yield the distribution of WSS, WSSG, OSI, and CWS in the stented vessel wall. An inverse relation (R(2) = 0.59, P < 0.005) between WSS and IH was found based on a linear regression analysis. Linear relations between WSSG, OSI, and IH were observed (R(2) = 0.48 and 0.50, respectively, P < 0.005). A linear relation (R(2) = 0.58, P < 0.005) between CWS and IH was also found. More statistically significant linear relations between the ratio of CWS to WSS (CWS/WSS), the products CWS × WSSG and CWS × OSI, and IH were observed (R(2) = 0.67, 0.54, and 0.56, respectively, P < 0.005), suggesting that both fluid and solid mechanics influence the extent of IH. Stents create endothelial flow disturbances and intramural wall stress concentrations, which correlate with the extent of IH formation, and these effects were exaggerated with mismatch of stent/vessel size. These findings reveal the importance of reliable vessel and stent sizing to improve the mechanics on the vessel wall and minimize IH.     

PEG-albumin supraplasma expansion is due to increased vessel wall shear stress induced by blood viscosity shear thinning

We studied the extreme hemodilution to a hematocrit of 11% induced by three plasma expanders: polyethylene glycol (PEG)-conjugated albumin (PEG-Alb), 6% 70-kDa dextran, and 6% 500-kDa dextran. The experimental component of our study relied on microelectrodes and cardiac output to measure both the rheological properties of plasma-expander blood mixtures and nitric oxide (NO) bioavailability in vessel walls. The modeling component consisted of an analysis of the distribution of wall shear stress (WSS) in the microvessels. Our experiments demonstrated that plasma expansion with PEG-Alb caused a state of supraperfusion with cardiac output 40% above baseline, significantly increased NO vessel wall bioavailability, and lowered peripheral vascular resistance. We attributed this behavior to the shear thinning nature of blood and PEG-Alb mixtures. To substantiate this hypothesis, we developed a mathematical model of non-Newtonian blood flow in a vessel. Our model used the Quemada rheological constitutive relationship to express blood viscosity in terms of both hematocrit and shear rate. The model revealed that the net effect of the hemodilution induced by relatively low-viscosity shear thinning PEG-Alb plasma expanders is to reduce overall blood viscosity and to increase the WSS, thus intensifying endothelial NO production. These changes act synergistically, significantly increasing cardiac output and perfusion due to lowered overall peripheral vascular resistance.     

Increased shear stress with upregulation of VEGF-A and its receptors and MMP-2, MMP-9, and TIMP-1 in venous stenosis of hemodialysis grafts

Venous injury and subsequent venous stenosis formation are responsible for hemodialysis graft failure. Our hypothesis is that these pathological changes are in part related to changes in wall shear stress (WSS) that results in the activation of matrix regulatory proteins causing subsequent venous stenosis formation. In the present study, we examined the serial changes in WSS, blood flow, and luminal vessel area that occur subsequent to the placement of a hemodialysis graft in a porcine model of chronic renal insufficiency. We then determined the corresponding histological, morphometric, and kinetic changes of several matrix regulatory proteins including VEGF-A, its receptors, matrix metalloproteinase (MMP)-2, MMP-9, tissue inhibitor of matrix metalloproteinase (TIMP)-1, and TIMP-2. WSS was estimated by obtaining blood flow and luminal vessel area by performing phase-contrast MRI with magnetic resonance angiography in 21 animals at 1 day after graft placement and prior to death on day 3 (n = 7), day 7 (n = 7), and day 14 (n = 7). At all time points, the mean WSS at the vein-to-graft anastomosis was significantly higher than that at the control vein (P < 0.05). WSS had a bimodal distribution with peaks on days 1 and 7 followed by a significant reduction in WSS by day 14 (P < 0.05 compared with day 7) and a decrease in luminal vessel area compared with control vessels. By day 3, there was a significant increase in VEGF-A and pro-MMP-9 followed by, on day 7, increased pro-MMP-2, active MMP-2, and VEGF receptor (VEGFR)-2 (P < 0.05) and, by day 14, increased VEGFR-1 and TIMP-1 (P < 0.05) at the vein-to-graft anastomosis compared with control vessels. Over time, the neointima thickened and was composed primarily of alpha-smooth muscle actin-positive cells with increased cellular proliferation. Our data suggest that hemodialysis graft placement leads to early increases in WSS, VEGF-A, and pro-MMP-9 followed by subsequent increases in pro-MMP-2, active MMP-2, VEGFR-1, VEGFR-2, and TIMP-1, which may contribute to the development of venous stenosis.