Blood flow shapes intravascular pillar geometry in the chick chorioallantoic membrane
© Lee et al; licensee BioMed Central Ltd. 2010
Received: 30 March 2010
Accepted: 7 July 2010
Published: 7 July 2010
Skip to main content
© Lee et al; licensee BioMed Central Ltd. 2010
Received: 30 March 2010
Accepted: 7 July 2010
Published: 7 July 2010
The relative contribution of blood flow to vessel structure remains a fundamental question in biology. To define the influence of intravascular flow fields, we studied tissue islands--here defined as intravascular pillars--in the chick chorioallantoic membrane. Pillars comprised 0.02 to 0.5% of the vascular system in 2-dimensional projection and were predominantly observed at vessel bifurcations. The bifurcation angle was generally inversely related to the length of the pillar (R = -0.47, P < .001). The pillar orientation closely mirrored the axis of the dominant vessel with an average variance of 5.62 ± 6.96 degrees (p = .02). In contrast, the variance of pillar orientation relative to nondominant vessels was 36.78 ± 21.33 degrees (p > .05). 3-dimensional computational flow simulations indicated that the intravascular pillars were located in regions of low shear stress. Both wide-angle and acute-angle models mapped the pillars to regions with shear less than 1 dyn/cm2. Further, flow modeling indicated that the pillars were spatially constrained by regions of higher wall shear stress. Finally, the shear maps indicated that the development of new pillars was limited to regions of low shear stress. We conclude that mechanical forces produced by blood flow have both a limiting and permissive influence on pillar development in the chick chorioallantoic membrane.
The mechanical influence of blood flow on vessel structure is a fundamental question in developmental  and adaptive  biology. In the chick chorioallantoic membrane, a common model of microvascular network development, the extra-embryonic area undergoes limited development in the absence of blood flow [3, 4]. In later embryogenesis, the onset of a heartbeat and active blood flow is associated with dramatic changes in both embryonic and extra-embryonic vessels [5–7]. In humans, physiological conditions such as growth and exercise lead to adjustments in the structural properties of the vascular network . In pathologic conditions such as inflammation [9, 10] and ischemia [11, 12], structural adaptations appear to be essential for tissue repair and regeneration. Despite these convincing network-level observations, there is little in vivo data on the local interaction between blood flow and vessel structure.
Attempts to clarify the mechanical influence of flow on blood vessels have focused on endothelial cell responses to varying flow patterns in vitro. Flow chamber studies have demonstrated that mechanical forces, such as wall shear stress, have a profound effect on gene transcriptional activity [13–15] and endothelial phenotype [16–18]. In vitro umbilical vein endothelial cells exposed to laminar shear stress reorganize and elongate their cytoskeletal axes in the direction of flow . The response of cultured endothelial cells to in vitro shear stress can also include lamellipodial protrusion and mechanotaxis in the direction of flow [20, 21]. The translation of these in vitro endothelial cell observations to in vivo structural change is less clear.
To define the local influence of intravascular flow fields on vessel structure, we have studied 2nd and 3rd order extra-embryonic microvessels in the chick chorioallantoic membrane (CAM). The 2nd and 3rd order CAM vessels are part of an experimentally accessible planar network that, in contrast to the complex gas exchange and nutrient function of the CAM capillaries, has a simple transport function. More importantly, the 2nd and 3rd order CAM vessels have a unique morphologic feature; namely, intravascular tissue islands or "pillars" . Discrete structures within the blood stream, pillars have several potential advantages in evaluating the effect of blood flow on local vessel development: 1) pillars are lined with normal-appearing endothelium  suggesting a normal responsiveness to intraluminal flow fields, 2) pillars are discrete structures indicating that local changes in pillar geometry have a minimal effect on global blood flow, and 3) pillars can be identified by time-series intravital 2D imaging providing a simultaneous assessment of pillar geometry and surrounding blood flow.
In this report, we used geometry and blood flow measurements derived from intravital microscopy imaging to map the mechanical forces within the CAM vessels--including wall shear stress and blood pressure--using 3D computational flow simulations. Pillar geometry suggested the spatial constraint of high wall shear stress. Further, the development of new pillars was limited to regions with low shear stress. The result suggests both a limiting and permissive influence of wall shear stress on pillar development in the CAM.
Specific pathogen-free, fertilized White Leghorn chicken eggs (G. gallus domesticus) were obtained from Charles River Laboratories (Franklin, CT). The care of the animals was consistent with guidelines of the American Association for Accreditation of Laboratory Animal Care (Bethesda, MD).
For all experiments, a modified, ex ovo (shell-less) culture method was used . Briefly, the eggs were kept in an R-COM 20 digital incubator (GimHae, Korea) at 37.5°C and 70% humidity with automatic turning for 3 days. On embryonic development day (EDD) 3, the eggs were sprayed with 70% ethanol, air-dried in a laminar flow hood and explanted into a 20 × 100 mm Petri dish (Falcon, BD Biosciences, San Jose, CA). The ex ovo cultures were maintained in a humidified 2% CO2 incubator at 37.5°C. To optimize the selective examination of the 2nd and 3rd order conducting, as well as facilitate intravital microscopy identification of the intravascular pillars, intravital microscopy was performed on EDD 13-16.
The CAM was imaged using a Nikon Eclipse TE2000 inverted epifluorescence microscope using Nikon Plan Apo 10x and Plan Fluor 20x objectives. The microscope was custom-fitted with an insulated 37°C convective warming unit with moderate relative humidity . An X-Cite (EXFO, Vanier, Canada) 120 watt metal halide light source and a liquid light guide were used to illuminate the CAM. Excitation and emission filters (Chroma, Rockingham, VT) in separate LEP motorized filter wheels were controlled by a MAC5000 controller (Ludl) and MetaMorph software 7.5 (Molecular Devices, Downingtown, PA). The 14-bit fluorescent images were digitally recorded with an electron multiplier CCD (EMCCD) camera (C9100-02, Hamamatsu, Japan). Images were routinely obtained at frame rates exceeding 50 fps with 2 × 2 binning. The images were recorded in image stacks comprising 100 to 500 frames of video sequences on a Dell Precision workstation (3.06 Ghz dual Xeon processors, 15,000 rpm ultra-SCSI hard drives, 4 gb RAM and an Nvidia Quadro 3450 graphics card with 512 mb memory). The CAMs at EDD13 thru EDD16 were imaged with intermittent time-lapse videos over 6 to 24 hour time period. Selection of vessels was based on an initial visual survey; identified intravascular pillars were then studied in detail. There was no attempt at uniform or random sampling.
The fluorescent plasma marker used for intravital imaging was a 5% fluorescein isothiocyanate (FITC)-dextran (2,000,000 MW; Sigma-Aldrich, St. Louis MO) solution prepared in normal saline immediately prior to injection. In some intravital microscopy experiments, green fluorescent (ex 430 nm; em 510), neutrally-charged, polystyrene spheres (108 beads/ml) were injected with the plasma marker . The 0.5 um microspheres were labeled with derivatives of the BODIPY fluorochrome (ex 488 nm, em 510 nm) using organic solvents (Invitrogen, Eugene, OR). The plasma marker and intravascular tracer solution were injected into the CAM circulation using a micro-fine 0.3 ml insulin syringe with a 30G needle (BD, Franklin Lakes, NJ).
Analysis of video images was performed with MetaMorph (Molecular Devices). Image stacks were created from the 100 to 500 frame sequences. The image stacks were processed with standard MetaMorph filters. After routine thresholding, the image sequences were measured using MetaMorph's integrated morphometry applications. Morphometric measurements such as area, length, orientation, perimeter, hole area and prolate volume were routinely obtained.
The stream-acquired images were stacked to create a time series of 100 or 500 consecutive frames. The stacks were systematically analyzed to ensure the absence of motion artifact. The stack "maximum" operation selected the highest intensity value for each pixel location throughout the time series. Conversely, the stack "minimum" operation selected the lowest intensity value for each pixel location. Other filters such as the "median" operation were similarly applied. The resultant image produced a time series reconstruction of the vessel during the time interval of the image stack.
Video mapping was performed by the serial acquisition of 3 × 3 image stacks. The image stacks were acquired using MetaMorph to control of an LEP XY motorized stage (Ludl); stage accuracy and computer control minimized the time interval between image stacks (<2 sec). After the image stacks at the grid positions were acquired, the images were stitched together for a panoramic video, using the MetaMorph stack montage application.
Similar to other branch angle studies emphasizing localized geometry , the method of branch angle morphometry was designed to be sensitive to variation at the apex of the bifurcation. In layered images, maximally-sized spheres were inscribed in each vessel at the bifurcation. Sequential spheres within each vessel were inscribed so that the surface of the sphere intersected the centerpoint of the preceding sphere. The centerline track of the first two spheres was used to define the vessel coordinates. After routine calibration and thresholding, the pillar axis was determined by the MetaMorph integrated morphometry application. The bifurcation angle was measured as the angle between the pillar and vessel axis.
To develop a finite element mesh of the vessel bifurcation, models were constructed with and without the pillar. The construction of 3D FEM for non-symmetrical, irregularly shaped 3D bifurcating vessel was based on geometry derived from intravital microscopy. Further, the transition at the pillar-wall interface required curved distortion of the mesh to match the morphology of the images from intravital microscopy . Otherwise, geometric measurements provided by intravital microscopy and digital image analysis were represented as faithfully as possible in the FEM model. As automation of the FEM generation was not possible at this stage, a customized approach was required for each model. Segmentation was performed on the original 2D intravital microscopy image to obtain 2D polylines of the countours. Subsequently, 2D splines were constructed from the polylines. Finally, a 3D FEM model made of 6 blocks for smooth mesh continuity was created using 3D NURBS (Non-Uniform Rational B-Spline) . In addition, two blocks of finer mesh were used around the pillar region to detect subtle changes in wall shear stress.
where was the blood velocity in three different directions vx,vy, and vz, ρ was the fluid density (1.05 g/cm3), and p was pressure, μ was the dynamic viscosity (0.03675 g/cm/sec). Equation (1) represented the balance of linear momentum, while equation (2) expressed the incompressibility condition. The code was validated using the analytical solution for shear stress and velocities through a straight expanding tube .
where was the mean velocity of the prescribed inlet parabolic velocity profile, and the vessel diameter D was measured by intravital microscopy. Resulting Re was always less than 1.0, suggesting that the blood flow modeled here was in a viscous flow regime .
The wall shear stress is calculated as , where denotes the tangential velocity immediately to the walls, and is the normal direction at the vessel wall. We first calculate the tangential velocity at the integration points near the wall surface, and then numerically evaluate the velocity gradient ; finally, we obtained 3 components (τ x , τ y , τ z ,) of the wall shear stress vector by multiple the velocity gradient by the viscosity coefficient μ. The effective value of the wall shear stress τ eff at the finite element mesh nodes on the surface was calculated as .
The numerical model was developed using custom code using the C++ object oriented programming language and OpenGL graphic library [32, 34–36]. The brick finite elements with 8 nodes used for our numerical simulation was more efficient than 4 node linear tetrahedral elements . The code was validated using the analytical solution for shear stress and velocities through a straight tube [32, 35]. The system of equations (1) was nonlinear due to the convective term and an unsymmetric Gaussian solver was implemented. An eight node finite element was employed with eight unknown velocities and constant pressure over the element which was recovered in the postprocessing calculation. Mesh independence was reached at 70,000 to 80,000 finite elements with error less then 0.1% for shear stress distribution. A smooth boundary surface was created to maintain a high resolution finite element mesh. For each calculation`, a parallel version of the solver  required 2 hours on 10 parallel processors with 2 GB RAM.
Significance estimates were based on multiple comparisons of paired data by Student-Newman-Keuls or Mann-Whitney test for non-parametric analysis of variance. The values for vessel and pillar orientation for each bifurcation were exported from MetaMorph and plotted in Excel 2007 (Microsoft, Redmond WA). Pearson correlations to the unamplified control were determined using Systat 12 statistical software (Chicago, IL). The significance level for the sample distribution was defined as P < .05.
In this report, we studied the geometry of intraluminal tissue islands--here referred to as intravascular pillars--and their surrounding blood flow in the CAM. Intravital microscopy of vessel bifurcations demonstrated marked pillar anisotropy with pillars orienting along the streamlines of the dominant vessel. Mechanical forces potentially influencing pillar geometry were mapped using 3D computational simulations. Shear maps demonstrated that pillars were predictably constrained by regions of high wall shear stress. Further, the development of new pillars was limited to regions with low shear stress. We conclude that mechanical forces have both a limiting and permissive influence on pillar development in the CAM.
Defining the relationship between mechanical forces and blood flow is experimentally challenging. Focal structural changes can be difficult to recognize in 2D intravital imaging . Definitive structure can be revealed by techniques such as corrosion casting and 3D scanning electron microscopy , but static imaging cannot provide a simultaneous assessment of blood flow. Global morphologic alterations, more readily identified by intravital microscopy, are confounded by the inter-relationship between blood flow and structure. In contrast, intravascular pillars provide a unique opportunity to assess both structure and blood flow. Pillars are discrete structures within the blood stream, relatively isolated from extravascular soluble tissue factors, that are not only identifiable by intravital 2D imaging, but amenable to detailed morphometric analysis. The pillars are lined with normal-appearing endothelium  suggesting a normal responsiveness to intraluminal flow fields.
Computational flow modeling provides important insights into flow-associated mechanical forces. The mechanical forces associated with blood flow include wall shear stress (the frictional force tangential to the vessel wall) and circumferential strain (a blood pressure-related force perpendicular to the direction of flow); wall shear stress being the dominant mechanical force in the smooth and continuous flow of the peripheral CAM microcirculation. Although wall shear stress cannot be directly measured, numerical simulations enable not only the calculation of wall shear stress, but the mapping of these mechanical forces to the vessel wall. The spatial relationships revealed by these maps provide important insights into the role of mechanical forces in shaping vessel structure and the contribution of mechanical forces in localized disease processes such as atherosclerosis.
Despite the value of a numerical analysis, there were also limitations of our computational approach. First, our computational simulation treated the blood as a Newtonian fluid and discounted the mechanical effects of blood cells. We have used a continuum approach in this study, because the chick chorioallantoic membrane has a blood cell concentration that is significantly lower, and more variable, than adults . This observation suggests that the structural modification is relatively insensitive to blood cell concentration. An alternative to continuum modeling is a discretized approach such as discrete particle dynamics (DPD) [41–43]. We have successfully applied DPD to the analysis of blood components in a parallel-plate flow chamber ; however, the computational demands of this promising technique currently preclude its application to the complex geometry of the CAM. Despite the theoretical limitations of a continuum approach, we suspect that our computational model closely approximates the distribution of forces in vivo.
A second limitation of our computational approach is that we do not consider extravascular tissue resistance. Since tissue viscoelasticity can neither be derived from first principles nor measured in vivo, our simplifying assumption is that the vessels function as rigid tubes. Our rationale is that we are modeling smooth and continuous flow in the peripheral extra-embryonic vascular network; not the pulsatile flow in the peri-embryonic vessels. Furthermore, we are not attempting to define the overall energetics of the extra-embryonic tissues, but rather describe the spatial distribution of the mechanical forces within the microvessels. The value of defining force distribution is that it can provide not only a mechanical explanation for structural changes, but a predictive map of endothelial changes. We anticipate that shear maps will facilitate metabolic and transcriptional studies both in vitro (e.g. flow chambers) and in vivo (e.g. laser capture microdissection).
The spatial relationship of pillar anisotropy and wall shear stress suggests that blood flow shapes pillar geometry. The pillar axis was almost uniformly oriented along the streamlines of the dominant vessel in a bifurcation. Consistent with these findings, Thoma observed more than 100 years ago the importance of the "stromrichtung" in determining vessel structure . The relevance of wall shear stress in influencing pillar structure is underscored by the absence of any comparable spatial relationship: the pillars were oriented randomly with respect to 1) different levels of the microcirculation, 2) the longitudinal axis of the embryo, 3) the vitelline (omphalomesenteric) vessels of the yolk sac and 4) the flow streamlines of neighboring nondominant vessels. The strong spatial coincidence of the streamlines and wall shear stress suggests that mechanical forces shape intravascular pillars in the CAM. We speculate that this observation may reflect a more general influence of mechanical forces on all microvessel endothelium.
The influence of blood flow on the lining of CAM microvessels is consistent with many observations of endothelial cells in vitro and intravascular pillars in vivo. The orientation and elongation of the intravascular pillars observed in our study is consistent with the sensitivity of endothelial cells to the direction of flow. Our observations are compatible with the cytoskeletal rearrangement, lamelipodial protrusion and mechanotaxis observed in culture [19–21]. The influence of intravascular pillars is also consistent with the branch angle remodeling and vascular pruning observed in vivo . In each case, the elongation of the pillar likely reflects the dominant streamlines in the flow field. A notable exception is the observation of pillars in divergent flow streams . In our study, new pillars were commonly observed at convergent flow streams or in locations "downstream" from pre-existing pillars. One explanation is simply sampling and/or technical limitations of our study. Another explanation is that pillars are typically formed at convergent bifurcations but are later observed at divergent flow streams after spontaneous flow reversal--an infrequent, but real, observation in the CAM microcirculation.
Because new pillars were not predictable from our shear maps, we suspect that mechanical forces have a permissive role in pillar development; that is, regions of low shear stress permit pillar development stimulated by other growth or developmental signals . A diffusible endothelial activation signal provides one explanation for the vessel irregularities observed during pillar development. A competing hypothesis is that mechanical forces stimulate new pillar formation and may even initiate the related process of intussusceptive angiogenesis . These are both intriguing and plausible possibilities that deserve further study.
embryonic development day
finite element model, FITC: fluorescein isothiocyanate
The authors wish to thank Mr. Zarko Milosevic for his expert technical assistance. Supported in part by NIH Grant HL47078, HL75426 and HL94567
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.