|
|
||||||||
1From the Research & Technology Directorate and National Center for Space Exploration Research, NASA Glenn Research Center, Cleveland, Ohio; and the 2Cole Eye Institute, Cleveland Clinic Foundation, Cleveland, Ohio.
| Abstract |
|---|
|
|
|---|
METHODS. Increasing concentrations of TA (0–16 ng/mL) were applied topically on embryonic day (E) 7 to the chorioallantoic membrane (CAM) of quail embryos cultured in petri dishes and incubated for an additional 24 or 48 hours until fixation. Binary (black/white) microscopic images of arterial end points were quantified by generational analysis of vessel branching (VESGEN) software to obtain major vascular parameters that include vessel diameter (Dv), fractal dimension (Df), tortuosity (Tv), and densities of vessel area, length, number, and branch point (Av, Lv, Nv, and Brv). For assessment of specific changes in vascular morphology induced by TA, the VESGEN software automatically segmented the vascular tree into branching generations (G1... G10) according to changes in vessel diameter and branching.
RESULTS. Vessel density decreased significantly up to 34% as the function of increasing concentration of TA according to Av, Lv, Brv, Nv, and Df. TA selectively inhibited the growth of new, small vessels because Lv decreased from 13.14 ± 0.61 cm/cm2 for controls to 8.012 ± 0.82 cm/cm2 at 16 ng TA/mL in smaller branching generations (G7–G10) and for Nv from 473.83 ± 29.85 cm–2 to 302.32 ± 33.09 cm–2. In contrast, vessel diameter (Dv) decreased throughout the vascular tree (G1–G10).
CONCLUSIONS. By VESGEN analysis, TA selectively inhibited the angiogenesis of smaller blood vessels, but decreased the vessel diameter of all vessels within the vascular tree.
Although the clinical effect of TA in these diseases is established, the effect of TA on vascular morphology is not well understood. Improved understanding of how TA affects angiogenesis and vascular morphology would be helpful for therapeutic optimization. For this study, site-specific changes within the vascular tree induced by TA were quantified using an in vivo model of angiogenesis.
For evaluation of the effects of TA on the angiogenic vascular tree, the quail chorioallantoic membrane (CAM) is a highly useful model during mid-development, when the rate of CAM angiogenesis is at a maximum.13 14 15 16 17 As described previously, the complex spatial patterns of the branching vascular tree and the associated capillary network can be easily visualized by light and fluorescence/confocal microscopy.13 14 15 16 17 Using fractal analysis, this vascular pattern can be precisely analyzed. Fractal analysis is a recent non-Euclidean mathematical innovation18 that quantifies the space-filling patterns of complex objects. Fractal geometry is common in nature and includes botanical and vascular trees, snowflakes, coastline topography, and even the spatiotemporal scaling of vascular-based physiological metabolism.19 20 As a fractional, nonintegral number that increases according to the increasing density of a space-filling pattern, the fractal dimension (Df) is statistically sensitive to small, early-stage changes in the vascular tree.13 14 15 17 21 A fractal object typically reaches its greatest space-filling capacity using self-similarity, the geometric property by which a pattern such as vascular bifurcational branching is repeated iteratively at continuously decreasing length scales.
The computer program VESGEN (abbreviated from generational analysis of vessel branching) was developed by the National Aeronautics and Space Administration (NASA) as a fully automated, user-interactive program that quantifies major vascular branching parameters using a single, user-provided image of two-dimensional (2D) vascular pattern. The fractal-based VESGEN analysis segments vessels within vascular trees into branching generations (G1, G2, ... Gx) according to changes in vessel diameter and branching. Site-specific changes within the vascular tree induced by angiogenic cytokines or other molecular regulators, such as TA, can then quantified.14 15 16 Thus, the purpose of this study was to evaluate the effects of the TA on branching morphology within the microvascular tree of the quail CAM model using fractal-based VESGEN analysis.
| Materials and Methods |
|---|
|
|
|---|
Culture, Assay, and Mounting
Fertilized eggs of Japanese quail (Coturnix coturnix japonica, Boyds Bird Co., Pullman, WA) were incubated at 37.6° ± 0.2°C under ambient atmosphere, cracked at embryonic day 3 (E3; after incubation of eggs for 56 hours), and cultured further in six-well petri dishes (cross-sectional area, 10 cm2). Quail egg culture and experimental protocols were in accordance with National Institutes of Health guidelines and approved by the Chief Veterinarian Officer of NASA (Ames Research Center). At E7 (after incubation for an additional 96 hours), prewarmed PBS solution containing TA (0.5 mL at 0–16 ng/mL) was applied dropwise to the surface of each CAM. Sterile-filtered stock solutions of TA (T-6501; Sigma, St. Louis, MO) were prepared at 10 mg/mL in 100% ethanol. Additional concentrations of TA were tested to determine the range of dose effectiveness. The total amount of the molecular regulator, rather than its concentration, is the governing parameter because solutions are quickly absorbed into CAM tissue. Quantities of TA are therefore reported as 0 to 8 ng/CAM. After treatment with TA and further incubation for 24 or 48 hours, the embryos (with CAMs) were fixed in 4% paraformaldehyde/2% glutaraldehyde/PBS for several days before dissection and mounting for microscopic analysis.13
Imaging
Aldehyde fixation of the CAM results in high contrast of the arterial tree because of retention of erythrocytes (red blood cells [RBCs]) within arteries but low contrast of the venous tree resulting from evacuation of RBC from veins during dissection (Fig. 1) .13 Digital images (1392 x 1040 pixels) of (terminal) arterial end point vessels from the middle region of the CAM were acquired in grayscale (0–255 intensity) at total 12.5x magnification and resolution of 7.32 µm/pixel (DM4000B microscope [Leica, Wetzlar, Germany] attached to a Retiga EXi CCD camera [Qimaging; Image Pro Plus software]). Previous studies showed that the CAM arterial end point regions analyzed by us are representative of changes induced by topical application of angiogenesis regulators throughout the CAM arterial trees,13 and that there is no significant increase in vascular density at a total magnification of 20x in comparison with 12.5x. Grayscale images were converted to binary (black/white) images of vascular morphology (Fig. 1) by semiautomatic computer processing using Adobe Photoshop 7.0 and NIH ImageJ software (http://rsb.info.nih.gov/ij/). The accuracy of vascular image binarization was confirmed by a second independent, experienced operator.
|
Analysis by VESGEN.
The NASA Glenn computer code VESGEN (Fig. 2) was used to measure parameters of vascular morphology that included vessel length density (Lv), vessel area density (Av), vessel branch point density (Brv), vessel number density (Nv), vessel tortuosity (Tv) and vessel diameter (Dv) for each branching generation G1 through G10. For example, Dv1–2 denotes Dv with respect to branching generations G1–G2. By VESGEN analysis, only one image was found to contain 11 branching generations (i.e., a few vessels of G11), which were therefore merged into image results for G7–10 (G
7). Lv, Av, Nv, and Brv were expressed as density functions by normalization to the area of the image containing the major arterial tree extracted as region of interest (ROI) or the entire image (Fig. 2) . Vessel diameter was calculated as Dv = Av/Lv. A trimmed skeleton was used to obtain accurate measurements for Lv in specific branching generations such as Lv1.16 Tortuosity (Tv) was estimated by the ratio of the length of a trimmed vessel (Lv) determined by VESGEN to the shortest distance between the vessel end points.
|
2) of the parent vessel diameter; therefore, the decrease of vessel diameter to 71% was used as the primary determinant of a new branching generation. As can be seen in biological branching trees (Fig. 2) , however, the branching of relatively symmetric offspring vessels is not perfectly symmetric, and the diameters of few offspring vessels are of the ideal 71% value. In addition, vessels tend to taper. To accommodate a range of vessel diameters within a branching generation, VESGEN contains a 15% default tolerance factor that is user adjustable. For the TA study, the tolerance factor was left at ±15%. For a small number of vessel segments (only five vessel segments within 39 total images), the automatic segmentation by VESGEN was incorrect according to the established criteria of generational classification based on vessel diameter and branching. The user-interactive features of VESGEN were therefore used to override and correct these few inaccuracies. A further consideration is that the most frequent branching event in a vascular tree is the asymmetric offshoot branching of a much smaller vessel from a larger vessel, presumably because of space-filling requirements of vascular branching (Fig. 2) . It is important to note that although the branching pattern of each vascular tree in a CAM or a human retina is unique, the space-filling properties of the vascular trees are remarkably uniform.13 21 VESGEN now functions as a fully automated Java-based software operating as a plug-in to NIH ImageJ. A binary vascular image is the single input required by VESGEN to quantify vascular trees, networks, or tree-network composites of highly 2D tissues such as the avian CAM, human retina, and rodent retina. VESGEN will be publicly available in the near future, and its full capabilities are being described elsewhere.
Fractal Analysis.
To support fractal analysis by the box-counting algorithm,13 each binary image was rescaled slightly to 1370 x 1024. A left- and right-most square image of 1024 x 1024 was extracted from each binary image and skeletonized (i.e., linearized; Fig. 1 ) using the NIH ImageJ skeletonizing algorithm. The fractal dimension (Df) was calculated for binary and skeletonized images by implementation of box-counting at a power of 2 using ImageJ.13 Values of Df for the left and right 1024 x 1024 images were averaged to obtain an overall Df for each original image. The fractal box-counting algorithm has since been incorporated into VESGEN to confirm the TA fractal results and for use in future studies. We consider VESGEN to be a fractal-based analysis of vascular trees because of the complex, non-Euclidean space-filling geometry of the vascular branching structures and VESGEN assignment of vascular parameters to specific, self-similar generations of vascular branching.
| Results |
|---|
|
|
|---|
|
|
Decreases in vascular density as measured by Df and Av in binary images can result from several morphologic changes that include decreased vessel length and number of vessels and/or decreased vessel diameter (Figs. 1B 1E) . Quantitative analysis by VESGEN confirmed visual observations that the morphologic mechanisms of decreased vascular density induced by TA included the decreased number densities restricted to smaller vessels and the overall thinning of vessel diameter. Vessel tortuosity as measured by Tv was unaffected, varying typically between 1.04 and 1.17 within a specimen. By Nv and Lv (Fig. 5) , vessel density decreased significantly for the smallest vessels of G7-G10 but not for large and medium-sized vessels of G1–G2, G3–G4, and G5–G6. Nv7–10 decreased significantly from 474 ± 30 cm–2 in controls to 302 ± 33 cm–2 at 8 ng TA/CAM (P = 0.0017), whereas Nv1–2, Nv3–4, and Nv5–6 remained relatively constant (Fig. 5 ; P = 0.69, 0.53, and 0.61, respectively). For example, Nv1–2 was 5.04 ± 0.32 cm–2 in controls compared with 5.15 ± 0.05 cm–2 at 8 ng TA/CAM. Small but consistent decreases in vessel diameter (Dv) were induced by TA throughout the vascular tree (Fig. 6) . Dv7–10 decreased from 28.1 ± 1.0 µm in controls to 25.5 ± 0.4 µm at 8 ng TA/CAM (P = 0.03), Dv5–6 from 61.1 ± 3.1 µm to 53.0 ± 3.1 µm (P = 0.08), Dv3–4 from 118.1 ± 6.0 µm to 101.6 ± 5.8 µm (P = 0.07), and Dv1–2 from 228.4 ± 12.1 µm to 199.0 ± 11.2 µm (P = 0.09).
|
|
0.10) and confidence intervals (P
0.05). Nonetheless, some binning of generational results, as performed for this study, improves and smooths the data because of increased statistical sampling14 17 and provides greater clarity in the presentation of results. | Discussion |
|---|
|
|
|---|
In this study, TA was applied topically to the quail CAM during mid-embryonic development, when angiogenesis is occurring at its maximum rate and angiogenic cytokines and regulators can be applied easily and uniformly in solution.13 14 15 16 17 The transparent CAM membrane develops rapidly, is highly vascularized, and is essentially 2D. Complex spatial patterns of the branching vascular tree and associated capillary network are easily visualized by light and fluorescence/confocal microscopy and quantified by fractal-based VESGEN analysis. As a relatively convenient experimental model, the angiogenic CAM exhibits some useful morphologic and functional similarities to angiogenic diseases of the quasi-2D retinal vascular tree. For example, region-based fractal methods developed in the CAM were successfully extended to the quantification of progression in diabetic vascular disease using clinical images of the human retina.21
We are using the computer software VESGEN to quantify major vessel parameters of blood and lymphatic vascular remodeling. Output parameters of VESGEN include vessel diameter, tortuosity, fractal dimension, and densities of vessel number, branch point, length, and area. Most parameters can be obtained by the user for the overall vascular image, the major vascular tree, and individual or merged branching generations. Vascular trees are decomposed into branching generations so that cytokine- or therapeutic-regulated modifications can be quantified according to site-specific vessel location. For this study on TA, we chose to analyze generational branching within the major vascular tree. Focusing on branching relationships within a single tree can support precise conclusions about regulator-induced changes in vessel branching relationships. This is the first technical report of results generated by the fully mature, newly automated VESGEN software (version 1.0) that now analyzes 2D vascular trees, networks, and tree-network composites for a number of experimental and clinical tissue applications in angiogenesis and lymphangiogenesis, including the avian CAM and yolk sac, the human retina, rodent retina, and developing coronary vessels in the embryonic heart.
As demonstrated by VESGEN analysis, it now appears that perturbation of rapid, ongoing angiogenesis during the middle stages of CAM development by TA occurs primarily by the selective inhibition or stimulation of new, small blood vessels. As a generalized result for the CAM model, this conclusion appears important, though not particularly surprising, because the molecular and cellular characteristics of angiogenic vascular tissues differ from those of more mature, stable vascular tissues.25 26 Angiogenic perturbants quantified to date in this fractal-based CAM model include the inhibitors TA, transforming growth factor β-1 (TGF-β1)14 and angiostatin,13 and stimulators basic fibroblast growth factor (bFGF)15 and vascular endothelial growth factor-165 (VEGF165).17 Additional regulator-specific effects on vascular morphology, such as overall vessel thinning by TA or thickening of larger blood vessels by VEGF165, have also been quantified. As for TA, morphologic response of the vascular tree to VEGF165 measured by VESGEN was multimodal. Increased vessel density and increased vessel diameter reached maximal frequencies at lower and higher VEGF concentrations, respectively. The study of TGF-β1 in particular considered tissue growth (rescaling) of the entire CAM vascular tree. Each molecular perturbant of angiogenesis has therefore elicited a unique "fingerprint" response that is spatiotemporally distinct and quantifiable, despite the apparent generality of inhibition and stimulation of angiogenesis in the CAM at the level of new, small vessels within the growing vascular tree.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication October 16, 2007; accepted January 21, 2008.
Disclosure: T.L. McKay, None; D.J. Gedeon, None; M.B. Vickerman, P; A.G. Hylton, None; D. Ribita, None; H.H. Olar, None; P.K. Kaiser, P; P. Parsons-Wingerter, P
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Patricia Parsons-Wingerter, Research & Technology Directorate and National Center for Space Exploration Research, NASA Glenn Research Center, MS 110-3, Cleveland, OH 44135; patricia.parsons{at}grc.nasa.gov.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
N. Cheung, K. C. Donaghue, G. Liew, S. L. Rogers, J. J. Wang, S.-W. Lim, A. J. Jenkins, W. Hsu, M. L. Lee, and T. Y. Wong Quantitative Assessment of Early Diabetic Retinopathy Using Fractal Analysis Diabetes Care, January 1, 2009; 32(1): 106 - 110. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |