- Commentary
- Open Access
- Published:

# Fractal parameters and vascular networks: facts & artifacts

*Theoretical Biology and Medical Modelling*
**volume 5**, Article number: 12 (2008)

## Abstract

### Background

Several fractal and non-fractal parameters have been considered for the quantitative assessment of the vascular architecture, using a variety of test specimens and of computational tools. The fractal parameters have the advantage of being scale invariant, i.e. to be independent of the magnification and resolution of the images to be investigated, making easier the comparison among different setups and experiments.

### Results

The success of several commercial and/or free codes in computing the fractal parameters has been tested on well known exact models. Based on such a preliminary study, we selected the code Frac-lac in order to analyze images obtained by visualizing the angiogenetic process occurring in chick Chorio Allontoic Membranes (CAM), assumed to be paradigmatic of a realistic 2D vascular network. Among the parameters investigated, the fractal dimension D_{f} proved to be the most robust estimator for CAM vascular networks. Moreover, only D_{f} was able to discriminate between effective and elusive increases in vascularization after drug-induced angiogenic stimulations on CAMs.

### Conclusion

The fractal dimension D_{f} is likely to be the most promising tool for monitoring the effectiveness of anti-angiogenic therapies in various clinical contexts.

## Introduction

The concept of fractal dimension was first introduced by Hausdorff [1] as a generalization of the geometrical dimension, and subsequently developed by Kolmogorov & Tihomirov [2]. By introducing a scale ε, according to which the original length of a segment is partitioned, and counting the number N of self-similar parts resulting from the partitioning, the fractal dimension D is defined as:

When D is integer, it reduces to the current geometrical dimension (i.e, by partitioning in 3 equal parts the sides of a square, we obtain 2 = log_{3} (9), doing the same with a cube we obtain 3 = log_{3} (27), etc). However, also non-integer values of D are possible, corresponding to different 'recipes', e.g. by dividing a segment in 3 parts (ε = 3) and considering a new figure in which N = 4 segments form a cusp (such as in the Koch curve). Another example is the branching structure with N = 5 equal branches (as in the 'bush1' configuration, see Figure 1). Such plots produce a 'covering' of the plane which is intrinsically self-similar or 'scaling' invariant.

Although, in mathematical terms, a structure can be analyzed at an arbitrary high resolution, and D is actually defined as a 'limit', when the 'conceptual' structure is plotted as a 'real' image, the original scale-independence is lost. Considering the image of a self-similar structure, two intrinsic limits, i.e. the image resolution and its dimension, define the minimal (ε_{m}) and maximal (ε_{M}) values for the scale parameter ε. When the image of a fractal structure is considered, in order to apply Eq. (1), a practical approach for estimating D would be that of selecting a series of values for ε (ε_{m} < ε < ε_{M}), perform for each of them a tessellation of the image in boxes and compute the number N_{b} of boxes containing the image structures (box-counting method, BCM), introduced by Mandelbrot [3]. For a detailed description of the method see Bunde & Havlin [4].

Results are clearly dependent on the procedure of tessellation and of evaluation of the linear regression. The quantification of such a procedure for different computational algorithms implemented in the most popular codes is the goal of the first part of our study.

In the real world, self-similarity becomes an even less well-defined concept. Accordingly, structures are generally defined as 'fractal-like', or by means of truncated fractals. First, real structures, although similar to segments, possess a 'thickness' (which can be disregarded provided some 'skeletonization' procedure is performed), which adds some 'noise' to the ideal fractal structure to which it resembles. Moreover, the minimal and maximal values within which the scaling behaviour is restricted (i.e. self-similarity is satisfied) can be further reduced by the structure itself. This point is of great importance for the biological systems, which have been recently investigated using scaling relationships by Brown and colleagues [5].

The rationale for an approach based on the assessment of fractal properties resides on previous studies on self-similar architectures observed in many biological structures, such as the bronchial tree [6] and the placental villous tree [7] and on the occurrence of scaling relationships (i.e., when two variables X and Y relate according to a given power law:

As an example, Y can be the basal metabolic rate, X the body mass of living organisms and *p* is a non-integer, fixed value (*p* = 3/4 according to Kleiber [8]and West et al [9]). In another example Y is the basal metabolic rate, X the tumor mass and *p* a non-integer value changing according to the developmental phase of the tumor itself [10].

Such scaling properties are actually assumed to be originated by the microvascular structure, which is responsible for the delivery of nutrients to body cells. In principle, a generic branching system, whose purpose is to exchange nutrients at its endpoints, may satisfactorily develop without restrictive rules, such as self-similarity. However, optimalization (i.e. maximal delivery of nutrients at the endpoints) is reached imposing some constraints on the vessels length and diameters, in order to minimize energy expenditure for flow and metabolic purposes, which can be expressed as a power law relation between parent and daughter vessels' dimensions (see Zamir, [11]). As West pointed out [9], power exponents *p* are expected to be different in the pulsatile regimen (arteries, for which at each branching point, the area is preserved to minimize wave reflection) and in the microcirculation.

Finally, almost no geometrical constraints are expected for vascular networks induced by angiogenesis and elicited by various 'growth factors' produced by proliferating tissues, e.g. tumors. The characteristics of such vascular patterns are mainly due to the interplay with the host, and vascular growth is driven by local information about pressure, blood velocity, and by the presence of some randomness and noise. This point was emphasized by Sandau & Kurz [12, 13], who suggested an extension of the concept of fractal dimension, called 'complexity', which best describes the case of non self-similar networks. Other parameters, related to the vascular network positional, topological and orientational orders were introduced by Guidolin et al [14].

Actually, the former considerations explain why Baish, Jain & Gazit [15] found that the in-vivo estimation of the fractal dimension of planar vascular networks in normal tissues and in four different tumor lines, implanted in the dorsal skinfold chamber of immuno-deficient mice ranges between the value of 2 for normal capillaries, 1.7 for arteries and veins and 1.88 for tumor vessels.

A similar estimation done on casts, derived from the application of other models proposed in the literature, gives values of 2 for the 'space-filling' growth model, 1.71 for the 'diffusion limited aggregation' model (simulating the arterio-venous system) and 1.90 for the 'invasion percolation' model (simulating the tumor neovasculature).

Due to the possible relevance of the fractal parameters for characterizing the neovascular tumor structures, which may be of interest for both diagnostic [16] and therapeutic [17] purposes, we wish here to investigate the vascular fractal dimension in the simplified model of the Chorio-Allontoic membranes (CAM) of chick eggs. A comparison with other parameters currently evaluated is performed. Moreover, since different expression of pro- and anti-angiogenic factors should elicit differences in the microvascular network development, we challenged the CAM with such compounds and evaluated their effects on the investigated parameters.

## Methodology

### Definition of the parameters

As previously stated, many geometrical forms in nature and, in particular, many vascular networks are thought to be 'fractal-like', i.e. they can be subdivided, up to a given scale, into 'self-similar' parts. The most popular (and potentially useful) parameters defined to describe fractals (from here on named 'fractal parameters') are the fractal dimension D_{f}, which is the 'experimental' counterpart of D defined in Eq. (1), and the lacunarity L. Qualitatively, D_{f} specifies how completely a fractal-like structure fills the space for decreasing scales, while the lacunarity assesses its texture, i.e. the distribution and size of the empty domains. The operative definitions given in our paper are the following. Images of the vascular network (of linear dimension Λ), obtained from any technique (microscopy, RMN, etc) are normally pre-processed in order to reduce their grey levels to a dichotomic (black/white) binary figure, or sometimes even skeletonised, i.e. the difference in diameter of the branches is neglected. Then a matrix of squares of side l = εΛ, with ε spanning in a given range (ε_{m} < ε < ε_{M}), is superposed on the binary image and the corresponding number N_{b} of boxes containing at least one black 'pixel' of the image, is computed.

This procedure, using Eq. (2), with Y = N_{b}, Y_{0} = 1 and X = ε, allows to define D_{f} = p. In other words, D_{f} is computed as the slope of the straight line Log N_{b} = D_{f} log ε in a log-log representation. If the variable under investigation is the variance over the square of the mean number of black points in the box of side l = εΛ, being σ the standard deviation and μ the mean, i.e. Y = (σ/μ)^{2} and Y_{0} = 1 then, from Eq. (2) we can argue that the lacunarity of the image is L = *p*. By virtue of their definitions, both D_{f} and L are expected to be scale invariant (**in a given range)** according to Eq.(2).

Among non-fractal parameters, the most diffused is the *Vascular density* V_{d}, which is generally evaluated by estimating the fraction of the image area covered by vessels. Also the distribution of diameters and lengths of the vessels, as well as the number of generations in the tree-like networks, are often evaluated [18]. In order to monitor angiogenesis, the number of neo-formed vessels has been quantified by counting the number of network nodes per unit section, called *angiogenic index* or *forks density F*_{
d
}[19].

### Validation of the parameters estimation methodology by applying various computational codes on 'exact models'

A few examples of 'exact models', i.e. '**fractal-like'** trees generated by dedicated SWs (e.g. the shareware 'treegenerator', http://www.treegenerator.com/) were considered, some of which simulating the arterial tree (as the 'bush1' and 'bush 3' models, see Figure 1) and the capillary network (the so-called 'Hilbert' model, see Figure 1). Their D_{f} is intrinsically defined by the constitutive constructive algorithm, i.e. the 'bush1' structure is obtained by dividing the initial segment into 5 sub-segments of length equal to 1/3 of the initial one, and replicating the above process. Therefore its fractal dimension is D_{f} = ln5/ln3 = 1.465. Analogously, 'bush 3' was proved to have D_{f} = 1.5 and the Hilbert network D_{f} = 2. No 'exact' estimations for the lacunarity L of the same models are available. We will therefore simply check its scaling properties.

It is straightforward to observe that, for what concerns the vascular density V_{d} and the number of forks F_{d}, no scale-invariance is expected. Images of trees extended to different generations were generated, in order to check the ability of various programs, mainly working on the basis of the BCM, to estimate the correct D_{f} and L values given by the exact calculations independently from the network extension. The code 'Winrhizo' was applied to simple fractal models to estimate F_{d}.

The commercial codes compared in this paper are 'Fractalyse', (ThèMA, F), Whinrhizo' (Regent Instruments Inc.) and Image Pro Plus (Media Cybernetics), while FDSURFFT is a MATLAB^{®} routine and FracLac is a freely downloadable Plugin of ImageJ http://rsb.info.nih.gov/ij/.

Table 1 summarizes the results for D_{f} obtained by means of several codes. Data are missing for the cases in which the program was unable to give an estimation.

The lacunarity L was also evaluated for the same exact models using FracLac. After a direct verification that L changes very slightly with the number of generations (according to the expected scaling properties), we found (mean ± Standard Deviation): L = 0.374 ± 0.008 for 'bush 1', L = 0.67 ± 0.04 for 'bush 3' and L = 0.17 ± 0.02 for Hilbert.

To summarize, according to Tab 1, only FracLac could evaluate D_{f} at best for all the exact models. Some of the programs, optimized for branching structures, such as FDSURFFT and Winrhizo, performed well for the 'bush' models but failed with 'Hilbert', while programs optimized for capillary structures, such as Image Pro Plus and Fractalyse, failed in estimating D_{f} for branched structures. Moreover, FracLac can estimate lacunarity. As far as limitations are concerned, we conclude that some computer codes addressed to the evaluation of the fractal dimensions perform well only in limited contexts (for instance for tree-like and/or capillary-like structures), and it is therefore necessary to check the code performance before application.

## CAMs as models for 2d vascular trees

### Estimation of the parameters on the CAMs: scaling properties and robustness

The Chorio-Allontoic Membrane (CAM) in the fertilized hen egg starts to develop on the 5th day of incubation by the fusion of the chorion and the allantois [20]. It is completely formed on the 11th day and lies attached under the inner eggshell membrane (See Figure 2). Throughout CAM development, blood vessels and capillaries expand intensively [20, 21]. The angiogenesis has been monitored using several parameters. For the chick embryo, Kirkner et al [22] showed that D_{f} exhibits an almost linear dependence on the incubation time, from around 1.3 at day 6 up to 1.8 at days 13–14, followed by a decrease in the later stages of maturation. A similar pattern was shown by Parsons-Wingerter for the quail embryo [23]. Their results (obtained by using the MATLAB-supported VESGEN code) have been validated using the FracLac code on the same images.

Six eggs were incubated at 37°C in a humidified environment. After 72 hours eggs were oriented and windowed, 6–10 ml removed from the egg and embryo were visually checked for heart beat (vitality). At day 10 Watman paper discs (6 mm diameter) were exposed to UV, soaked in hydrocortisone (3 mg/ml in ethanol 100%), dried and placed on CAMs. Contaminated CAMs were excluded from the study. After treatments, membranes are fixed with 3.7% PAF in PBS and after 10 min CAMs were removed and placed into a Petri dish containing PAF. Discs were excised and imaging performed a JVC TK-C1380E color video camera (ImageProPlus 4.0 imaging software) connected to a stereomicroscope (model SZX9; Olympus). Pictures, taken at different magnifications, were processed with the dedicated software.

At first the scaling properties were tested among three different magnifications (6×, 10× and 16×) for the parameters V_{d}, L and D_{f}, showing that both V_{d} and D_{f} vary significantly depending on the magnification only when it is larger than 10 (p < 0.05 ANOVA), while no test for L was reliable due to the huge data dispersion (see fig 3).

In conclusion, D_{f} satisfies the scaling properties (as well as V_{d}), provided images are taken at low and intermediate magnification (< 16 x).

In order to test the parameters robustness, i.e. suitability with the (sometimes unpredictable) variations in its operating environment, we applied the code FracLac on the six images (magnification 10×) to estimate D_{f}, L and V_{d}. Each evaluation of D_{f} and L was performed averaging results from four computational procedures, which assumed as starting point for the box counting algorithm the 4 different corners of the image.

Results are given in figure 4. It is apparent that D_{f} is very consistently defined (no statistically difference, p < 0.05), and the final value is 1,733 ± 0,006. On the contrary, the result for L is affected by a large deviation (0,36 ± .0,11), producing a percentage error of about 50% and vanifying any further statistical analysis Also the non-fractal parameter V_{d} shows a much larger, statistically significant variability, i.e. V_{d} = 0,37 ± 0,14 (p < 0.05).

In conclusion, D_{f} proved to be a very robust parameter, showing only very slight changes among the estimated values on different CAMs. Moreover, the value of D_{f} estimated on CAMs confirms its characteristics of being an efficient distributive system.

### Estimation of the parameters on the CAMs: testing sensibility and specificity

In order to test whether the parameters are sensitive (i.e. can correctly identify a condition of enhanced vascularisation without false positives) and specific (i.e. can correctly identify a condition of apparent increase of vascularisation, without false negatives), we performed a drug response analysis. The main application of the studies reported in the literature was the investigation of CAMs response to various proangiogenic drugs (mainly VEGF and FGF) in terms of vascular architecture variation. For instance, when administered at concentrations between 1,25 and 2,5 μg, VEGF 165 induces increases in arterial density and in arterial diameter. Correspondingly, D_{f} was shown to increase from 1,65 ± 0,01 to 1,69 ± 0,01 at its maximal dose [24, 25].

Another commonly used proangiogenic compound is the Fibroblast Growth Factor 2 (FGF). Guidolin's group showed a significant increase in D_{f} from 1,51 ± 0,01 of the control to 1,62 ± 0,04 after treatment at day 12 at the dose of 500 ng. On the contrary, L decreases from 0,288 ± 0,02 to 0,198 ± 0,02 [14]. The same group showed that, provided the angiogenetic process is antagonized with angiostatic factors, the fractal dimension of the vascular network reflects the observed decrease in branching patterns (about 10% from the starting value of 1.20) [14]. Other drugs proved to be far less effective in promoting angiogenesis. For instance Angiopoietin 1 (Ang), a ligand for the Tie2 receptor of endothelial cells, acts as a regulator of angiogenesis in several experimental models, although it is not clear whether it has a pro- or anti-angiogenic effect [26–31].

On six CAM specimens, prepared according to the procedures already described above, the fractal dimension D_{f} (using FracLac), the vascular density V_{d} and the fork density F_{d} (using Winrhizo) were computed after administration of FGF and Ang. Statistical analysis was performed using the Kruskal-Wallis nonparametric test. Following the administration of FGF and Ang, we obtained the following results (Figures 5a and 5b): D_{f} for the control was (1.733 ± 0,006), while after administration of FGF increased to (1,826 ± 0,042, p < 0,05) and after administration of Ang was not statistically different versus the control (1,709 ± 0,061).

As far as the vascular density is concerned, in the control V_{d} was (0,37 ± 0,14), after FGF was (0,33 ± 0,10) and after Ang was (0,26 ± 0,02), with no statistical difference among groups. Finally we examined the third index, i. e. the number of forks per mm^{2}, F_{d} finding that in the control it was (2,3 ± 0,4), after FGF administration became (8,9 ± 1,2) and after Ang (8,2 ± 1,3). The fork density increased significantly for both FGF and Ang treatment (see Figure 5c).

## Conclusion

Fractal parameters, such as the fractal dimension and lacunarity, have been widely used to investigate vascular systems, particularly those formed by neoangiogenesis. In this contribution we have at first validated several computational codes and discussed some of their possible weaknesses in the estimation of the fractal parameters. We then applied the code FracLac to CAM images. We found D_{f} = 1,733 ± 0,006 at day 10, which is consistent with literature reports. We noted that the parameter D_{f} is very robust, reproducible and reliable, in contrast with the other fractal (L) and non-fractal parameters (V_{d} and F_{d}) considered.

The most remarkable result of our analysis is that completely different conclusions can be drawn from the same set of data following drug treatments depending on the chosen parameters.

The density of forks F_{d} showed a marked increase after administration of both Ang and FGF, suggesting that both drugs are equally effective in promoting branching. On the contrary, none of the drugs affected the vascular density V_{d}, with a surprisingly discrepancy with respect to expectations from the literature [14, 32].

On the same set of data D_{f} resolved a finer discrimination between the pro-angiogenic effect of FGF compared to the more controversial regulatory effect of Ang. It proved therefore to be both sensitive for the angiogenic affect of FGF and specific for the ineffectiveness of Ang.

According to Stoeltzing[31]Ang possibly plays a role in a later phase of developmental angiogenesis, such as remodelling and maturation of vessel network. Such effect, although of great importance on a biological ground, is not expected to modify the fractal-like structure of the CAM microcirculation.

Our analysis refers only to 2D structures and needs to be extended to the more realistic case of 3D microvascular networks. It is important, however, to stress that important applications to 2D microcirculatory systems have already been investigated, e.g. by De Felice and colleagues [33]. In particular, they were able to show that D_{f} could discriminate the oral vascular networks (gingival and vestibular oral mucosa) from controls and carriers of hereditary non-polyposis colorectal cancer (Lynch Cancer Family Syndrome II), i.e. D_{f} is a marker for LCFS2.

Our general conclusions is that, if caution is paid in the selection of the images, in their handling, and in the selection of the code, the fractal dimension D_{f} can represent a fast, reliable and robust parameter for evaluating angiogenetic processes. After further validation and refining, such parameters may potentially be useful in vascular and vascular-related diagnostics, for instance in monitoring the effectiveness of tumor anti-angiogenic therapies, anti-proliferative treatments in autoimmune syndromes and retinopathies and for stadiation of tumor developmental phases. Such applications would be an interesting extension to previous results, which have been considered to be relevant for diagnosis and therapy, i.e. that the fractal dimension is a prognostic factor for laryngeal carcinoma [34], for endometrial carcinoma [35] and for ovarian cancer [36]. Even more important, allometric scaling invariance has been observed in tumor growth processes involving lymph nodes [16], and could have a clinical impact in the future.

## References

- 1.
Hausdorff E: Dimension und auβeres Maβ. Mathematische Annalen. 1919, 79: 157-179. 10.1007/BF01457179.

- 2.
Kolmogorov A, Tihomirov V: ε-entropy and ε-capacity of sets of functional spaces. American Mathematical Society Translations. 1961, 17: 277-364.

- 3.
Mandelbrot B: How long is the coast of britain? statistical self-similarity and fractional dimension. Science. 1967, 156: 636-638. 10.1126/science.156.3775.636.

- 4.
Bunde A, Havlin S: A Brief Introduction to Fractal Geometry. 1994, Berlin: Springer

- 5.
Brown JH, Gupta VK, Li BL, Milne BT, Restrepo C, West GB: The fractal nature of nature: power laws, ecological complexity and biodiversity. Philos Trans R Soc Lond B Biol Sci. 2002, 357 (1421): 619-626. 10.1098/rstb.2001.0993.

- 6.
Kitaoka H, Suki B: Branching design of the bronchial tree based on a diameter-flow relationship. J Appl Physiol. 1997, 82: 968-976.

- 7.
Bergman DL, Ullberg U: Scaling properties of the placenta's arterial tree. J Theor Biol. 1998, 193: 731-738. 10.1006/jtbi.1998.0741.

- 8.
Kleiber M: The fire of life: an introduction to animal energetics. 1975, Huntington, NY: Robert E Krieger

- 9.
West GB, Brown JH, Enquist BJ: A general model for ontogenetic growth.[see comment]. Nature. 2001, 413: 628-631. 10.1038/35098076.

- 10.
Guiot C, Delsanto PP, Carpinteri A, Pugno N, Mansury Y, Deisboeck TS: The dynamic evolution of the power exponent in a universal growth model of tumors. J Theor Biol. 2006, 240: 459-463. 10.1016/j.jtbi.2005.10.006.

- 11.
Zamir M: On fractal properties of arterial trees. J Theor Biol. 1999, 197: 517-526. 10.1006/jtbi.1998.0892.

- 12.
Sandau K, Kurz H: Modelling of vascular growth processes: a stochastic biophysical approach to embryonic angiogenesis. J Microsc. 1994, 175: 205-213.

- 13.
Sandau K, Kurz H: Measuring fractal dimension and complexity – an alternative approach with an application. J Microsc. 1997, 186: 164-176. 10.1046/j.1365-2818.1997.1270685.x.

- 14.
Guidolin D, Vacca A, Nussdorfer GG, Ribatti D: A new image analysis method based on topological and fractal parameters to evaluate the angiostatic activity of docetaxel by using the Matrigel assay in vitro. Microvasc Res. 2004, 67: 117-124. 10.1016/j.mvr.2003.11.002.

- 15.
Baish JW, Gazit Y, Berk DA, Nozue M, Baxter LT, Jain RK: Role of tumor vascular architecture in nutrient and drug delivery: an invasion percolation-based network model. Microvasc Res. 1996, 51: 327-346. 10.1006/mvre.1996.0031.

- 16.
Demicheli R, Biganzoli E, Boracchi P, Greco M, Hrushesky WJ, Retsky MW: Allometric scaling law questions the traditional mechanical model for axillary lymph node involvement in breast cancer. J Clin Oncol. 2006, 24: 4391-4396. 10.1200/JCO.2006.05.5988.

- 17.
Jain RK: Normalization of tumor vasculature: an emerging concept in antiangiogenic therapy. Science. 2005, 307: 58-62. 10.1126/science.1104819.

- 18.
Zetter BR: Angiogenesis and tumor metastasis. Annu Rev Med. 1998, 49: 407-424. 10.1146/annurev.med.49.1.407.

- 19.
Cascone I, Napione L, Maniero F, Serini G, Bussolino F: Stable interaction between alpha5beta1 integrin and Tie2 tyrosine kinase receptor regulates endothelial cell response to Ang-1. J Cell Biol. 2005, 170: 993-1004. 10.1083/jcb.200507082.

- 20.
Vargas A, Zeisser-Labouebe M, Lange N, Gurny R, Delie F: The chick embryo and its chorioallantoic membrane (CAM) for the in vivo evaluation of drug delivery systems. Advanced Drug Delivery Reviews. 2007, 59: 1162-1176. 10.1016/j.addr.2007.04.019.

- 21.
Staton CA, Stribbling SM, Tazzyman S, Hughes R, Brown NJ, Lewis CE: Current methods for assaying angiogenesis in vitro and in vivo. Int J Exp Pathol. 2004, 85: 233-248. 10.1111/j.0959-9673.2004.00396.x.

- 22.
Kirkner FJ, Wisham WW, Giedt FH: A report on the validity of the rorschach prognostic rating scale. Journal of Projective Techniques. 1953, 17: 465-470.

- 23.
Parsons-Wingerter P, Lwai B, Yang MC, Elliott KE, Milaninia A, Redlitz A, Clark JI, Sage EH: A novel assay of angiogenesis in the quail chorioallantoic membrane: stimulation by bFGF and inhibition by angiostatin according to fractal dimension and grid intersection. Microvasc Res. 1998, 55: 201-214. 10.1006/mvre.1998.2073.

- 24.
Avakian A, Kalina RE, Sage EH, Rambhia AH, Elliott KE, Chuang EL, Clark JI, Hwang JN, Parsons-Wingerter P: Fractal analysis of region-based vascular change in the normal and non-proliferative diabetic retina. Curr Eye Res. 2002, 24: 274-280. 10.1076/ceyr.24.4.274.8411.

- 25.
Parsons-Wingerter P, Chandrasekharan UM, McKay TL, Radhakrishnan K, DiCorleto PE, Albarran B, Farr AG: A VEGF165-induced phenotypic switch from increased vessel density to increased vessel diameter and increased endothelial NOS activity. Microvasc Res. 2006, 72: 91-100. 10.1016/j.mvr.2006.05.008.

- 26.
Chae JK, Kim I, Lim ST, Chung MJ, Kim WH, Kim HG, Ko JK, Koh GY: Coadministration of angiopoietin-1 and vascular endothelial growth factor enhances collateral vascularization. Arterioscler Thromb Vasc Biol. 2000, 20 (12): 2573-2578.

- 27.
Hangai M, Moon YS, Kitaya N, Chan CK, Wu DY, Peters KG, Ryan SJ, Hinton DR: Systemically expressed soluble Tie2 inhibits intraocular neovascularization. Hum Gene Ther. 2001, 12: 1311-1321. 10.1089/104303401750270968.

- 28.
Hawighorst T, Skobe M, Streit M, Hong YK, Velasco P, Brown LF, Riccardi L, Lange-Asschenfeldt B, Detmar M: Activation of the tie2 receptor by angiopoietin-1 enhances tumor vessel maturation and impairs squamous cell carcinoma growth. Am J Pathol. 2002, 160: 1381-1392.

- 29.
Shim WS, Teh M, Bapna A, Kim I, Koh GY, Mack PO, Ge R: Angiopoietin 1 promotes tumor angiogenesis and tumor vessel plasticity of human cervical cancer in mice. Exp Cell Res. 2002, 279: 299-309. 10.1006/excr.2002.5597.

- 30.
Uemura A, Ogawa M, Hirashima M, Fujiwara T, Koyama S, Takagi H, Honda Y, Wiegand SJ, Yancopoulos GD, Nishikawa S: Recombinant angiopoietin-1 restores higher-order architecture of growing blood vessels in mice in the absence of mural cells.[see comment]. J Clin Invest. 2002, 110: 1619-1628.

- 31.
Stoeltzing O, Ahmad SA, Liu W, McCarty MF, Wey JS, Parikh AA, Fan F, Reinmuth N, Kawaguchi M, Bucana CD, Ellis LM: Angiopoietin-1 inhibits vascular permeability, angiogenesis, and growth of hepatic colon cancer tumors. Cancer Res. 2003, 63: 3370-3377.

- 32.
Parsons-Wingerter P, Chandrasekharan UM, McKay TL, Radhakrishnan K, DiCorleto PE, Albarran B, Farr AG: A VEGF165-induced phenotypic switch from increased vessel density to increased vessel diameter and increased endothelial NOS activity. Microvasc Res. 2006, 72: 91-100. 10.1016/j.mvr.2006.05.008.

- 33.
De Felice C, Latini G, Bianciardi G, Parrini S, Fadda GM, Marini M, Laurini RN, Kopotic RJ: Abnormal vascular network complexity: a new phenotypic marker in hereditary non-polyposis colorectal cancer syndrome.[see comment]. Gut. 2003, 52: 1764-1767. 10.1136/gut.52.12.1764.

- 34.
Delides A, Panayiotides I, Alegakis A, Kyroudi A, Banis C, Pavlaki A, Helidonis E, Kittas C: Fractal dimension as a prognostic factor for laryngeal carcinoma. Anticancer Res. 2005, 25: 2141-2144.

- 35.
Dey P, Rajesh L: Fractal dimension in endometrial carcinoma. Anal Quant Cytol Histol. 2004, 26: 113-116.

- 36.
Kikuchi A, Kozuma S, Sakamaki K, Saito M, Marumo G, Yasugi T, Taketani Y: Fractal tumor growth of ovarian cancer: sonographic evaluation.[see comment]. Gynecol Oncol. 2002, 87: 295-302. 10.1006/gyno.2002.6842.

## Acknowledgements

GF Varetto is recipient of a Progettolagrange CRT grant. We wish to thank PP Delsanto for useful suggestions and comments.

## Author information

### Affiliations

### Corresponding author

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

## About this article

### Cite this article

Mancardi, D., Varetto, G., Bucci, E. *et al.* Fractal parameters and vascular networks: facts & artifacts.
*Theor Biol Med Model* **5, **12 (2008). https://doi.org/10.1186/1742-4682-5-12

Received:

Accepted:

Published:

### Keywords

- Fractal Dimension
- Fibroblast Growth Factor
- Basal Metabolic Rate
- Vascular Network
- Fractal Parameter