Shear stress and pressure modulate saphenous vein remodeling ex vivo

Vein graft failure remains an important clinical challenge, but factors contributing to vein graft failure have not clearly been defined. We investigated the role of the mechanical environment in vein remodeling in an ex vivo perfusion system. Porcine saphenous veins were subjected to five different ex vivo hemodynamic environments, including one mimicking an arterial bypass graft, for one week in order to independently assess the effects of shear stress and pressure on vein remodeling. The extent of intimal hyperplasia decreased with culture under increasing shear stress, with veins cultured under the lowest levels of shear stress exhibiting the greatest ratio of intimal/medial area, 0.15+/-0.03, which was greater than that of fresh veins (0.06+/-0.01, p<0.05). All perfused veins displayed characteristics of both medial hypertrophy and eutrophic remodeling, with those veins cultured under elevated pressures showing greater increases in mass and area than those cultured under venous pressures. Medial area correlated with the average pressure under which veins were cultured (R2=0.95, p<0.001), with veins cultured under bypass graft conditions, which were exposed to the greatest pressure during the one week culture, exhibiting the largest medial area (1.69+/-0.15 mm2), which was significantly greater than that of fresh veins (1.08+/-0.05 mm2, p<0.05). However, pulsatility was not a necessary stimulus for medial growth, as increases in medial area were observed in culture conditions in which steady flow and pressure were present. Our results suggest that pressure and shear stress act independently to regulate vein remodeling, influencing changes in vessel size as well as the nature of the remodeling.     

Non-Newtonian effects of blood flow on hemodynamics in distal vascular graft anastomoses

Non-Newtonian fluid flow in a stenosed coronary bypass is investigated numerically using the Carreau-Yasuda model for the shear thinning behavior of the blood. End-to-side coronary bypass anastomosis is considered in a simplified model geometry where the host coronary artery has a 75% severity stenosis. Different locations of the bypass graft to the stenosis and different flow rates in the graft and in the host artery are studied. Particular attention is given to the non-Newtonian effect of the blood on the primary and secondary flow patterns in the host coronary artery and the wall shear stress (WSS) distribution there. Interaction between the jet flow from the stenosed artery and the flow from the graft is simulated by solving the three-dimensional Navier-Stokes equation coupled with the non-Newtonian constitutive model. Results for the non-Newtonian flow, the Newtonian flow and the rescaled Newtonian flow are presented. Significant differences in axial velocity profiles, secondary flow streamlines and WSS between the non-Newtonian and Newtonian fluid flows are revealed. However, reasonable agreement between the non-Newtonian and the rescaled Newtonian flows is found. Results from this study support the view that the residual flow in a partially occluded coronary artery interacts with flow in the bypass graft and may have significant hemodynamic effects in the host vessel downstream of the graft. Non-Newtonian property of the blood alters the flow pattern and WSS distribution and is an important factor to be considered in simulating hemodynamic effects of blood flow in arterial bypass grafts.     

Laminar-to-turbulent transition in pulsatile flow through a stenosis

Laminar-to-turbulent transition in pulsatile flow through a stenosis is studied by means of three-dimensional numerical simulations. The flow transition is associated with the occurrence of a flow instability initiating in the stenosis region. The instability is manifested by a three-dimensional symmetry-breaking and leads to asymmetric separation and intense swirling motion downstream of the stenosis. The above have profound effects on the wall shear stress (WSS). The simulations reveal that the asymmetric separation is extended several radii downstream of the stenosis with substantial WSS fluctuations, in both space and time, occurring in the poststenotic region.     

A novel vein graft model: adaptation to differential flow environments

Accelerated intimal hyperplasia in response to altered flow environment is critical to the process of vein bypass graft failure. Lack of a reproducible animal model for dissecting the mechanisms of vein graft (VG) remodeling has limited progress toward solving this clinically significant problem. Combining a cuffed anastomotic technique with other surgical manipulations, we developed a well-defined, more robust method for studying hemodynamic factors in VG arterialization. VG with fistula placement, complete occlusion, or partial distal branch ligation (DBL) was performed in the carotid artery of 56 rabbits. Extensive hemodynamic and physiological analyses were performed to define the hemodynamic forces and histological adaptations of the wall at 1-28 days. Anastomotic time averaged 12 min, with 100% patency of bilateral grafts and unilateral grafts plus no adjunct or delayed fistula. Bilateral VG-DBL resulted in an immediate disparity in wall shear (0.8 +/- 0.1 vs. 12.4 +/- 1.1 dyn/cm2, ligated vs. contralateral graft). Grafts exposed to low shear stress responded primarily through enhanced intimal thickening (231 +/- 35 vs. 36 +/- 18 microm, low vs. high shear). High-shear-stress grafts adapted through enhanced outward remodeling, with a 24% increase in lumen diameter at 28 days (3.0 +/- 0.1 vs. 3.7 +/- 0.2 mm, low vs. high shear). We have taken advantage of the cuffed anastomotic technique and combined it with a bilateral VG-DBL model to dissect the impact of hemodynamic forces on VG arterialization. This novel model offers a robust, clinically relevant, statistically powerful small animal model for evaluation of high- and low-shear-regulated VG remodeling.     

Mapping vascular response to in vivo Hemodynamics: application to increased flow at the basilar terminus

Hemodynamic forces play critical roles in vascular pathologies such as atherosclerosis, aneurysms, and stenosis. However, detailed relationships between the specific in vivo hemodynamic microenvironment and vascular responses leading to the triggering or exacerbation of pathological remodeling of the vessel remain elusive. We have developed a hemodynamics-biology co-mapping technique that enables in situ correlation between the in vivo blood flow field and vascular changes secondary to hemodynamic insult. The hemodynamics profile is obtained from computational fluid dynamics simulation within the vascular geometry reconstructed from three-dimensional in vivo images, whereas the vascular response is obtained from histology or immunohistochemistry on harvested vascular tissue. The hemodynamics field is virtually sectioned in the histological slicing planes and digitally co-mapped with the histological images, thereby enabling correlation of the specific local vascular responses with the inciting hemodynamic stresses. We demonstrate application of this technique to rabbit basilar terminus subjected to elevated flow. Morphological changes at the basilar terminus 5 days after the flow increase were co-mapped with the initial wall shear stress and wall shear stress gradient distributions, from which localization of destructive remodeling in a specific hemodynamic zone was noticed. This method paves the way for further investigations to determine the connection between in vivo mechanical stimuli and biological responses, such as initiation of aneurysmal remodeling.     

Numerical and experimental models of post-operative realistic flows in stenosed coronary bypasses.

By means of both experimental and finite element methods, we simulated three-dimensional unsteady flows through coronary bypass anastomosis. The host artery includes a stenosis shape located at two different distances of grafting. The inflow rates are issued from in vivo measurements in patients who had undergone coronary bypass surgery a few days before. We provide a comparison between experimental and numerical velocity profiles coupled with the numerical analysis of spatial and temporal wall shear stress evolution. The interaction between the graft and coronary flows has been demonstrated. The phase inflow difference can partly be responsible for specific flow phenomena: jet deflection towards a preferential wall or feedback phenomenon that causes the flapping of the post-stenotic jet during the cardiac cycle. In conclusion, we showed the sensitivity of these typical flows to distance of grafting, inflows waveforms but also to their phase difference.     

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.     

Rule-Based Model of Vein Graft Remodeling

When vein segments are implanted into the arterial system for use in arterial bypass grafting, adaptation to the higher pressure and flow of the arterial system is accomplished thorough wall thickening and expansion. These early remodeling events have been found to be closely coupled to the local hemodynamic forces, such as shear stress and wall tension, and are believed to be the foundation for later vein graft failure. To further our mechanistic understanding of the cellular and extracellular interactions that lead to global changes in tissue architecture, a rule-based modeling method is developed through the application of basic rules of behaviors for these molecular and cellular activities. In the current method, smooth muscle cell (SMC), extracellular matrix (ECM), and monocytes are selected as the three components that occupy the elements of a grid system that comprise the developing vein graft intima. The probabilities of the cellular behaviors are developed based on data extracted from in vivo experiments. At each time step, the various probabilities are computed and applied to the SMC and ECM elements to determine their next physical state and behavior. One- and two-dimensional models are developed to test and validate the computational approach. The importance of monocyte infiltration, and the associated effect in augmenting extracellular matrix deposition, was evaluated and found to be an important component in model development. Final model validation is performed using an independent set of experiments, where model predictions of intimal growth are evaluated against experimental data obtained from the complex geometry and shear stress patterns offered by a mid-graft focal stenosis, where simulation results show good agreements with the experimental data.     

Computational simulation of flow-induced arterial remodeling of the pancreaticoduodenal arcade associated with celiac artery stenosis

Arterial remodeling of the pancreaticoduodenal arcade, which enables collateral flow to the liver, spleen, and stomach, is a well-recognized clinical sign of celiac artery (CA) stenosis. However, the hemodynamic changes due to remodeling are poorly understood, despite their importance in surgical procedures such as pancreaticoduodenectomy. In this study, a framework to simulate remodeling of the arterial network following pathological flow alterations was developed and applied to investigate the hemodynamic characteristics of patients with CA stenosis. A one-dimensional–zero-dimensional cardiovascular model was used for blood flow simulation. After introducing CA stenosis into the normal network, arterial remodeling was simulated by iteratively changing the diameter of each artery until time-averaged wall shear stress reached its value under normal conditions. A representative case was simulated to validate the present framework, followed by simulation cases to investigate the impact of stenosis severity on remodeling outcome. A markedly dilated arcade was observed whose diameter agreed well with the corresponding values measured in subjects with CA stenosis, confirming the ability of the framework to predict arterial remodeling. A series of simulations clarified how the geometry and hemodynamics after remodeling change with stenosis severity. In particular, the arterial remodeling and resulting blood flow redistribution were found to maintain adequate organ blood supply regardless of stenosis severity. Furthermore, it was suggested that flow conditions in patients with CA stenosis could be estimated from geometric factors, namely, stenosis severity and arcade diameter, which can be preoperatively and non-invasively measured using diagnostic medical images.     

Numerical analysis of hemodynamics in spastic middle cerebral arteries

Cerebral vasospasm (CVS) is the most common serious complication of subarachnoid hemorrhage. Among the many factors that are associated with the pathogenesis of CVS, hemodynamics plays an important role in the initiation and development of CVS. Numerical simulation was carried out to obtain the flow patterns and wall shear stress (WSS) distribution in spastic middle cerebral arteries. The blood was assumed to be incompressible, laminar, homogenous, Newtonian, and steady. Our simulations reveal that flow velocity and WSS level increase at the stenosis segment of the spastic vessels, but further downstream of stenosis, the WSS significantly decreases along the inner wall, and flow circulation and stagnation are observed. The hydrodynamic resistance increases with the increase of vessel spasm. Moreover, the change of flow field and hydrodynamic forces are not linearly proportional to the spasm level, and the rapid change of hemodynamic parameters is observed as the spasm is more than 50%. Accordingly, in the view of hemodynamic physiology, vessels with less than 30% stenosis are capable of self-restoration towards normal conditions. However, vessels with more than 50% stenosis may eventually lose their capacity to adapt to differing physiologic conditions due to the extreme non-physilogic hemodynamic environment, and the immediate expansion of the vessel lumen might be needed to minimize serious and non-reversible effects.     

Three-dimensional numerical simulations of flow through a stenosed coronary bypass

This work analyzes the flow patterns at the anastomosis of a stenosed coronary bypass. Three-dimensional numerical simulations are performed using a finite elements method. We consider a geometrical model of the host coronary artery with and without a 75% severity stenosis for three different locations from the anastomosis. The flow features - velocity profiles, secondary motions and wall shear stresses - are compared for different configurations of the flow rate and of the distance of the anastomosis from the site of occlusion (called distance of grafting). The combination of the junction flow effects - counter rotating vortices - with the stenosis effects - confined jet flow - is particularly important when the distance of grafting is short. Given that the residual flow issued from the pathologic stenosis being non-negligible after two weeks grafting, models without stenosis cannot predict the evolution of the wall shear stress in the vicinity of the anastomosis.     

Relative contribution of wall shear stress and injury in experimental intimal thickening at PTFE end-to-side arterial anastomoses

BACKGROUND: Intimal hyperplastic thickening (IHT) is a frequent cause of prosthetic bypass graft failure. Induction and progression of IHT is thought to involve a number of mechanisms related to variation in the flow field, injury and the prosthetic nature of the conduit. This study was designed to examine the relative contribution of wall shear stress and injury to the induction of IHT at defined regions of experimental end-to-side prosthetic anastomoses. METHODS AND RESULTS: The distribution of IHT was determined at the distal end-to-side anastomosis of seven canine Iliofemoral PTFE grafts after 12 weeks of implantation. An upscaled transparent model was constructed using the in vivo anastomotic geometry, and wall shear stress was determined at 24 axial locations from laser Doppler anemometry measurements of the near wall velocity under conditions of pulsatile flow similar to that present in vivo. The distribution of IHT at the end-to-side PTFE graft was determined using computer assisted morphometry. IHT involving the native artery ranged from 0.0+/-0.1 mm to 0.05+/-0.03 mm. A greater amount of IHT was found on the graft hood (PTFE) and ranged from 0.09+/-0.06 to 0.24+/-0.06 mm. Nonlinear multivariable logistic analysis was used to model IHT as a function of the reciprocal of wall shear stress, distance from the suture line, and vascular conduit type (i.e. PTFE versus host artery). Vascular conduit type and distance from the suture line independently contributed to IHT. An inverse correlation between wall shear stress and IHT was found only for those regions located on the juxta-anastomotic PTFE graft. CONCLUSIONS: The data are consistent with a model of intimal thickening in which the intimal hyperplastic pannus migrating from the suture line was enhanced by reduced levels of wall shear stress at the PTFE graft/host artery interface. Such hemodynamic modulation of injury induced IHT was absent at the neighboring artery wall.     

The Role of Shear Stress in Arteriovenous Fistula Maturation and Failure: A Systematic Review

Introduction

Non-maturation and post-maturation venous stenosis are the primary causes of failure within arteriovenous fistulae (AVFs). Although the exact mechanisms triggering failure remain unclear, abnormal hemodynamic profiles are thought to mediate vascular remodelling and can adversely impact on fistula patency.

Aim

The review aims to clarify the role of shear stress on outward remodelling during maturation and evaluate the evidence supporting theories related to the localisation and development of intimal hyperplasia within AVFs.

Methods

A systematic review of studies comparing remodelling data with hemodynamic data obtained from computational fluid dynamics of AVFs during and after maturation was conducted.

Results

Outward remodelling occurred to reduce or normalise the level of shear stress over time in fistulae with a large radius of curvature (curved) whereas shear stress was found to augment over time in fistulae with a small radius of curvature (straight) coinciding with minimal to no increases in lumen area. Although this review highlighted that there is a growing body of evidence suggesting low and oscillating shear stress may stimulate the initiation and development of intimal medial thickening within AVFs. Further lines of evidence are needed to support the disturbed flow theory and outward remodelling findings before surgical configurations and treatment strategies are optimised to conform to them. This review highlighted that variation between the time of analysis, classification of IH, resolution of simulations, data processing techniques and omission of various shear stress metrics prevented forming pooling of data amongst studies.

Conclusion

Standardised measurements and data processing techniques are needed to comprehensively evaluate the relationship between shear stress and intimal medial thickening. Advances in image acquisition and flow quantifications coupled with the increasing prevalence of longitudinal studies commencing from fistula creation offer viable techniques and strategies to robustly evaluate the relationship between shear stress and remodelling during maturation and thereafter.     

Mechanical properties of native and ex vivo remodeled porcine saphenous veins

When grafted into an arterial environment in vivo, veins remodel in response to the new mechanical environment, thereby changing their mechanical properties and potentially impacting their patency as bypass grafts. Porcine saphenous veins were subjected for one week to four different ex vivo hemodynamic environments in which pressure and shear stress were varied independently, as well as an environment that mimicked that of an arterial bypass graft. After one week of ex vivo culture, the mechanical properties of intact saphenous veins were evaluated to relate specific aspects of the mechanical environment to vein remodeling and corresponding changes in mechanics. The compliance of all cultured veins tended to be less than that of fresh veins; however, this trend was more due to changes in medial and luminal areas than changes in the intrinsic properties of the vein wall. A combination of medial hypertrophy and eutrophic remodeling leads to significantly smaller (p<0.05) wall stresses measured in all cultured veins except those subjected to bypass graft conditions relative to stresses measured in fresh veins at corresponding pressures. Our results suggest that the mechanical environment effects changes in vessel size, as well as the nature of the remodeling, which contribute to altering vein mechanical properties.     

Procedure for Human Saphenous Veins Ex Vivo Perfusion and External Reinforcement

The mainstay of contemporary therapies for extensive occlusive arterial disease is venous bypass graft. However, its durability is threatened by intimal hyperplasia (IH) that eventually leads to vessel occlusion and graft failure. Mechanical forces, particularly low shear stress and high wall tension, are thought to initiate and to sustain these cellular and molecular changes, but their exact contribution remains to be unraveled. To selectively evaluate the role of pressure and shear stress on the biology of IH, an ex vivo perfusion system (EVPS) was created to perfuse segments of human saphenous veins under arterial regimen (high shear stress and high pressure). Further technical innovations allowed the simultaneous perfusion of two segments from the same vein, one reinforced with an external mesh. Veins were harvested using a no-touch technique and immediately transferred to the laboratory for assembly in the EVPS. One segment of the freshly isolated vein was not perfused (control, day 0). The two others segments were perfused for up to 7 days, one being completely sheltered with a 4 mm (diameter) external mesh. The pressure, flow velocity, and pulse rate were continuously monitored and adjusted to mimic the hemodynamic conditions prevailing in the femoral artery. Upon completion of the perfusion, veins were dismounted and used for histological and molecular analysis. Under ex vivo conditions, high pressure perfusion (arterial, mean = 100 mm Hg) is sufficient to generate IH and remodeling of human veins. These alterations are reduced in the presence of an external polyester mesh.     

Pulsatile flow in an end-to-side vascular graft model: comparison of computations with experimental data

Various hemodynamic factors have been implicated in vascular graft intimal hyperplasia, the major mechanism contributing to chronic failure of small-diameter grafts. However, a thorough knowledge of the graft flow field is needed in order to determine the role of hemodynamics and how these factors affect the underlying biological processes. Computational fluid dynamics offers much more versatility and resolution than in vitro or in vivo methods, yet computations must be validated by careful comparison with experimental data. Whereas numerous numerical and in vitro simulations of arterial geometries have been reported, direct point-by-point comparisons of the two techniques are rare in the literature. We have conducted finite element computational analyses for a model of an end-to-side vascular graft and compared the results with experimental data obtained using laser-Doppler velocimetry. Agreement for velocity profiles is found to be good, with some clear differences near the recirculation zones during the deceleration and reverse-flow segments of the flow waveform. Wall shear stresses are determined from velocity gradients, whether by computational or experimental methods, and hence the agreement for this quantity, while still good, is less consistent than for velocity itself from the wall shear stress numerical results, we computed four variables that have been cited in the development of intiimal hyperplasia-the time-averaged wall shear stress, an oscillating shear index, and spatial and temporal wall shear stress gradients in order to illustrate the versatility of numerical methods. We conclude that the computational approach is a valid alternative to the experimental approach for quantitative hemodynamic studies. Where differences in velocity were found by the two methods, it was generally attributed to the inability of the numerical method to model the fluid dynamics when flow conditions are destabilizing. Differences in wall shear, in the absence of destabilizing phenomena, were more likely to be caused by difficulties in calculating wall shear from relatively low resolution in vitro data.     

Investigation of hemodynamics in the development of dissecting aneurysm within patient-specific dissecting aneurismal aortas using computational fluid dynamics (CFD) simulations

Aortic dissecting aneurysm is one of the most catastrophic cardiovascular emergencies that carries high mortality. It was pointed out from clinical observations that the aneurysm development is likely to be related to the hemodynamics condition of the dissected aorta. In order to gain more insight on the formation and progression of dissecting aneurysm, hemodynamic parameters including flow pattern, velocity distribution, aortic wall pressure and shear stress, which are difficult to measure in vivo, are evaluated using numerical simulations. Pulsatile blood flow in patient-specific dissecting aneurismal aortas before and after the formation of lumenal aneurysm (pre-aneurysm and post-aneurysm) is investigated by computational fluid dynamics (CFD) simulations. Realistic time-dependent boundary conditions are prescribed at various arteries of the complete aorta models. This study suggests the helical development of false lumen around true lumen may be related to the helical nature of hemodynamic flow in aorta. Narrowing of the aorta is responsible for the massive recirculation in the poststenosis region in the lumenal aneurysm development. High pressure difference of 0.21 kPa between true and false lumens in the pre-aneurismal aorta infers the possible lumenal aneurysm site in the descending aorta. It is also found that relatively high time-averaged wall shear stress (in the range of 4-8 kPa) may be associated with tear initiation and propagation. CFD modeling assists in medical planning by providing blood flow patterns, wall pressure and wall shear stress. This helps to understand various phenomena in the development of dissecting aneurysm.     

Velocity measurements in steady flow through axisymmetric stenoses at moderate Reynolds numbers

The velocity field in the neighborhood of axisymmetric constrictions in rigid tubes was investigated using laser Doppler anemometry and flow visualization. Upstream flow conditions were steady; and Reynolds numbers were in the range 500-2000, values which are representative of the larger arteries in humans. Stenoses of 25, 50 and 75% area reduction were studied. Velocity profiles are presented in sufficient detail to allow comparison with computational biofluid dynamics models. Wall shear stresses were estimated from the near wall velocity gradient, and the nature of observed poststenotic flow disturbances is discussed. Results indicate that flow disturbances of discrete oscillation frequency may be more valuable than turbulence as an indicator of early stages of stenosis development. Additionally, despite the fact that poststenotic turbulence exists for the higher degrees of stenosis and Reynolds numbers, the resulting wall shear stresses are only three to four times greater than the Poiseuille value and are considerably less than the wall shear stress within the stenosis itself.     

In vivo demonstration of flow recirculation and turbulence downstream of graded stenoses in canine arteries

The velocity field around arterial stenoses was investigated using a pulsed doppler velocimeter in vivo. Asymmetric zones of recirculation were identified by systolic flow reversal in the post-stenotic field in carotid and iliac arteries of anesthetised dogs. There was a close correlation between shear intensity and turbulence as estimated by the velocity difference between the jet and the recirculation zone and by maximum spectral width respectively. Under the conditions of these experiments, stenosis grade (% diameter reduction) dominated hemodynamic variables such as Reynolds number, oscillation and pulsatility in determining the intensity of turbulence. The method used does not appear to have sufficient resolution to distinguish between turbulence and discrete oscillating velocities (vorticity) nor to allow determination of wall shear stress though the pattern of change of the latter is similar to that found downstream of axisymmetric stenosis in models using steady flow.     

A new imaging technique to study 3-D plaque and shear stress distribution in human coronary artery bifurcations in vivo

OBJECTIVE: Bifurcations of coronary arteries are predilection sites for atherosclerosis and expansive remodeling, the latter being associated with plaque vulnerability. Both are related to blood flow-induced shear stress (SS). We present a new approach to generate 3-D reconstructions of coronary artery bifurcations in vivo and investigate the relationship between SS, wall thickness (WT) and remodeling. METHODS: The patient specific 3-D reconstruction of the main branch of the bifurcation was obtained by combining intravascular ultrasound and biplane angiography, and the 3-D lumen of the side branch was based on biplane angiography only. The two data sets were fused and computational methods were applied to determine the SS distribution, using patient derived flow and viscosity data. The intravascular ultrasound data allowed us to measure local WT and remodeling in the main branch. RESULTS: The lumen reconstruction procedure was successful and it was shown that the impact of the side branch on SS distribution in the main branch diminished within 3mm. Distal to the bifurcation, two continuous regions in the main branch were identified. In the proximal region, we observed lumen preservation, and expansive remodeling. Although a plaque was observed in the low SS region at the non-divider wall, no relationship between SS and WT was found. In the distal region, we observed lumen narrowing and a significant positive relationship between SS and WT. CONCLUSIONS: A new imaging technique was applied to generate a 3-D reconstruction of a human coronary artery bifurcation in vivo. The observed relationship between SS, WT and remodeling in this specific patient illustrates the spatial heterogeneity of the atherosclerosis in the vicinity of arterial bifurcations.