AJP - Heart Fuel your research with LabChart
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Heart Circ Physiol 287: H2493-H2500, 2004; doi:10.1152/ajpheart.00371.2004
0363-6135/04 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kalsho, G.
Right arrow Articles by Kassab, G. S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kalsho, G.
Right arrow Articles by Kassab, G. S.

Bifurcation asymmetry of the porcine coronary vasculature and its implications on coronary flow heterogeneity

Ghada Kalsho and Ghassan S. Kassab

Department of Biomedical Engineering, University of California, Irvine, California 92697-2715

Submitted 19 April 2004 ; accepted in final form 30 July 2004


    ABSTRACT
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
The branching pattern of the coronary arteries and veins is asymmetric, i.e., many small vessels branch off of a large trunk such that the two daughter vessels at a bifurcation are of unequal diameters and lengths. One important implication of the geometric vascular asymmetry is the dispersion of blood flow at a bifurcation, which leads to large spatial heterogeneity of myocardial blood flow. To document the asymmetric branching pattern of the coronary vessels, we computed an asymmetry ratio for the diameters and lengths of all vessels, defined as the ratio of the daughter diameters and lengths, respectively. Previous data from silicone elastomer cast of the entire coronary vasculature including arteries, arterioles, venules, and veins were analyzed. Data on smaller vessels were obtained from histological specimens by optical sectioning, whereas data on larger vessels were obtained from vascular casts. Asymmetry ratios for vascular areas, volumes, resistances, and flows of the various daughter vessels were computed from the asymmetry ratios of diameters and lengths for every order of mother vessel. The results show that the largest orders of arterial and venous vessels are most asymmetric and the degree of asymmetry decreases toward the smaller vessels. Furthermore, the diameter asymmetry at a bifurcation is significantly larger for the coronary veins (1.7–6.8 for sinus veins) than the corresponding arteries (1.5–5.8 for left anterior descending coronary artery) for orders 2–10, respectively. The reported diameter asymmetry at a bifurcation leads to significant heterogeneity of blood flow at a bifurcation. Hence, the present data quantify the dispersion of blood flow at a bifurcation and are essential for understanding flow heterogeneity in the coronary circulation.

asymmetry ratio; order number; flow heterogeneity; Poiseuille's resistance; coronary arteries


MYOCARDIAL BLOOD FLOW DISTRIBUTION has been investigated by many authors (see reviews in Refs. 2 and 5). It has been well established that myocardial perfusion is spatially heterogeneous showing a 6- to 10-fold range in normal hearts (2). Bassingthwaighte and colleagues (2) were the first to recognize the fractal nature of coronary flow heterogeneity. They proposed fractal networks to predict the dependence of relative dispersion (SD/mean) on spatial resolution of the perfusion measurements (23). Their model entailed asymmetry of flow distribution in each bifurcation of a fractal tree. It was shown that the simplest fractal model that consists of constant branching asymmetry at each bifurcation predicts the flow data fairly well. If the degree of asymmetry was made to increase with successive branching in the fractal model, however, a better fit to data was obtained. The obvious question is does the coronary arterial tree possess the asymmetry postulated by the fractal models, and, if so, does the asymmetry vary along the bifurcations? The present study was designed to address these issues.

Our hypothesis is that the asymmetry of vascular geometry at a bifurcation is an important determinant of coronary blood flow heterogeneity. We have previously described the branching pattern (connectivity and longitudinal position matrices) and vascular geometry (diameters and lengths) of the coronary vasculature using the diameter-defined Strahler system (7, 10, 22). We found the branching pattern and vascular geometry to be quite asymmetric, which necessitated the use of the diameter-defined Strahler system. The asymmetry was most pronounced for the epicardial vessels that give rise to transmural vessels in "large trunk-small twig" patterns. The asymmetry of the branching pattern was reflected by the pattern of connectivity and longitudinal position matrices and the segments-per-element ratios.

In the present study, we describe the bifurcation asymmetry in terms of local parameters at a bifurcation called asymmetry ratios. For a given mother vessel, we defined the diameter and length asymmetry ratio as the ratio of the daughter diameters and lengths, respectively. The diameter and length ratios were computed for each mother vessel throughout the coronary vasculature (from arteries to arterioles and from venules to veins). The resulting diameter and length asymmetry ratios were used to compute the asymmetry of area, volume, resistance, and flow at various bifurcations throughout the coronary vasculature. These data allow us to predict the flow heterogeneity at a bifurcation along with other important hemodynamic parameters. Furthermore, these data are necessary for the reconstruction of coronary vascular circuits for future hemodynamic analysis.


    METHODS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
All experiments were performed in accordance with national and local ethical guidelines, including the Institute of Laboratory Animal Resources, National Institutes of Health Guide for the Care and Use of Laboratory Animals, and University of California-Irvine policies regarding the use of animals in research.

Available data. We have previously described in detail the methods of preparation, measurement, and morphometric analysis of the right coronary artery (RCA), left anterior descending (LAD) coronary artery, left circumflex (LCx) artery, coronary sinus, and thebesian veins (10, 11). The morphometric data on the coronary vessels of diameters <50 µm were obtained from histological specimens, whereas the morphometric data on vessels of diameters >50 µm were obtained from cast studies as described in Kassab et al. (10). The morphometric measurements (diameters, lengths, and node-to-node connectivity) of the microvasculature were made by changing the focal plane through the thickness of the histological sections. For cast studies, the trunks of major coronary arteries and veins were sketched, and their segments measured down to ~50 µm in diameter. The raw coronary morphometric data showing node-to-node connectivity, diameters, and lengths are available on http://cvbiomech.eng.uci.edu.

Once the morphometric measurements were completed, we organized the segments and their diameters and lengths in terms of the diameter-defined Strahler's ordering scheme as described in Refs. 7 and 10. Briefly, the capillary blood vessels are defined as vessels of order 0. The smallest arterioles supplying blood to the capillaries are assigned an order number of 1. The smallest venules draining the capillaries are assigned an order number of –1. When two arterioles of order 1 meet, the confluent vessel is given an order of number 2 if its diameter exceeds the diameters of the order 1 vessels by an amount specified by a set of formulas or diameter criterion (10) or remain as order 1 if the diameter of the confluent is not larger than the amount specified by the formulas. When an order 2 artery meets another order 1 artery, the order number of the confluent is 3 if its diameter is larger by an amount specified by the diameter criterion or remains at 2 if its diameter does not increase sufficiently. This process was continued until all arterial segments were arranged in increasing diameter and assigned the order numbers 1, 2, 3, ..., n, ... Similarly, the veins are assigned the orders –1, –2, –3..., –n, ... The criterion is a value that is midway between the mean diameter of order n + 1SD and the mean diameter of the order n – 1SD for all n from 1 upward (10).

Asymmetry at a bifurcation. With the existence of a full set of ordered data on the coronary vasculature, we proceeded to define and compute the asymmetry ratios by considering each bifurcation node. The largest and the two smaller vessels meeting at a node are called the mother and daughter vessels, respectively, if the order numbers of these vessels are so arranged that n ≥ m ≥ k with diameters Dn ≥ Dm ≥ Dk. Unlike the coronary arteries, which have cylindrical cross sections and can be quantified by a diameter, the veins have approximately elliptical cross sections. We shall denote the major and minor axes diameters of the coronary veins as Dmaj and Dmin, respectively. For each mother vessel, n, the asymmetry ratio for the major and minor diameters of the daughter vessels, SDmaj and SDmin, are the ratios defined by

(1a)
and

(1b)
where SDmaj≥ 1, by definition. Analogously, the asymmetry ratio for the length of the daughter vessels, Sl, is the ratio defined by

(2)
The asymmetry ratios of the vessel cross-sectional area, A, and volume, v, are similarly defined by

(3)
and

(4)
because the cross-sectional area of an elliptical venous vessel is A = ({pi}/4)DmajDmin, and the volume of a branch is the product of its cross-sectional area and length (11).

The flow resistance, R, for a vein with an elliptical cross section is computed as follows: assuming that the daughter vessels are long and slender, that the flow is laminar and steady, and that the flow disturbances at the entry and exit sections are negligible, then the resistance to flow in the daughter vessels is given by (Ref. 11)

(5)
where µ is the apparent viscosity. Hence, we can define an asymmetry ratio of the flow resistance by the equation

(6a)
where

(6b)
Sµ is defined as the apparent viscosity ratio of the two daughter branches (µmk). The data on the variation of viscosity with vessel diameter and hematocrit are given by Pries et al. (19), which was used in our calculations. They proposed a modified viscosity relationship based on a compilation of literature data on relative blood viscosity in tube flow in vitro and in vivo, which reflects the Fahraeus-Lindqvist effect as given by:

(7)
where D is the vessel diameter. This relationship was used throughout the coronary vasculature for a given vessel diameter D at a hematocrit of 0.45.

For the coronary arteries, the major and minor axes diameters are equal, and hence Eq. 1 reduces to SD = (Dm/Dk). Similarly, Eq. 3 becomes SA = SD2. Finally, Eq. 5 becomes Poiseuille's resistance (i.e., R = 128 µl/{pi}D), and Eq. 6 reduces to SR = SµSl/SD4.

Horsfield et al. (6) have previously defined the difference in order number between two daughter branches at a bifurcation ({Delta} = mk) to describe the asymmetry at a bifurcation of an airway. Additionally, we shall introduce a parameter {varepsilon} ({varepsilon} = nm), which is the difference between the order of the mother vessel and that of the larger daughter.

Data analysis. All asymmetry parameters were curve fitted by a relation of the form ln(Sx) = a + bn + cn2, where ln is the natural logarithm and x corresponds to diameter (D), length (l), cross-sectional area (A), volume (v), or resistance (R). The empirical constants a, b, and c were determined by a least-squares fit. ANOVA was used to detect differences among different orders for the various parameters (i.e., SD, Sl, SA, Sv, SR, {Delta}, and {varepsilon}). Although our data are skewed, it is generally accepted that ANOVA is robust to violations of the normality assumption (4). Statistical significance was accepted at the 0.05 level.


    RESULTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
We obtained data from the entire range of orders throughout the coronary arterial and venous trees. A total of 5,082, 4,014, and 1,946 arterial bifurcations of RCA, LAD, and LCx arteries and 13,778 and 4,505 venous bifurcations of coronary sinus and thebesian veins and their branches, respectively, was analyzed. Probability frequency distribution functions were obtained for SD and Sl for each order of the entire coronary vasculature and were found to be right skewed. Tables 1GoGoGo5 show a summary of the means ± SD of SD and Sl in RCA, LAD, and LCx arteries and sinusal and thebesian veins, respectively. Using ANOVA, we found a significant dependence of SD and Sl on the order number of the mother vessel. We also found that, at a given order number, the SD for the venous trees was significantly larger than those of the corresponding arteries (orders 2–10). The Sl was larger for veins of orders 3–6 but smaller for orders 2 and 7–10 compared with arteries. The data on SD and Sl were represented by a quadratic exponential curve fit and are summarized in Table 6.


View this table:
[in this window]
[in a new window]
 
Table 1. Asymmetry ratios of diameters and lengths of daughter vessels in each order of mother vessels in the RCA and its branches of the pig

 

View this table:
[in this window]
[in a new window]
 
Table 2. Asymmetry ratios of diameters and lengths of daughter vessels in each order of mother vessels in the LAD artery and its branches of the pig

 

View this table:
[in this window]
[in a new window]
 
Table 3. Asymmetry ratios of diameters and lengths of daughter vessels in each order of mother vessels in the LCx artery and its branches of the pig

 

View this table:
[in this window]
[in a new window]
 
Table 4. Asymmetry ratios of major diameters and lengths of daughter vessels in each order of mother vessels in the coronary sinus vein and its branches of the pig

 

View this table:
[in this window]
[in a new window]
 
Table 5. Asymmetry ratios of major diameters and lengths of daughter vessels in each order of mother vessels in the Thebesian veins and their branches of the pig

 

View this table:
[in this window]
[in a new window]
 
Table 6. Regression coefficients for the various asymmetry parameters

 
Although SD ≥ 1, by definition, Sl may be less than or greater than 1. We quantified the fraction of larger daughter vessels that also have longer lengths, i.e., Flm > lk = [N(lm > lk)/N(lm > lk) + N(lm lk)], where N(lm > lk) and N(lm > lk) represent the number of bifurcations where the larger daughter diameter vessel has a longer and shorter length, respectively. Figure 1 shows the relationship between Flm > lk and the order number of the mother vessel throughout the coronary vasculature.



View larger version (16K):
[in this window]
[in a new window]
 
Fig. 1. Relationship between the fraction of larger daughter vessels that also have longer lengths than the smaller daughter vessels and the order number of the mother vessel of the coronary vasculature [Flm > lk = N(lm > lk)/N(lm > lk) + N(lm < lk)], where N(lm > lk) and N(lm > lk) represent the number of bifurcations where the larger daughter diameter vessel has a longer and shorter length, respectively. RCA, right coronary artery; LAD, left anterior descending coronary artery; LCx, left circumflex artery.

 
The SA and Sv of daughter branches of coronary blood vessels of order n were computed from the SD and Sl as given by Eqs. 3 and 4, respectively. Figures 2 and 3 show the relationship between the mean SA and Sv, respectively, for the entire range of orders. It can be seen that SA and Sv increase rapidly with the order number for the coronary arteries and veins, becoming very large at highest order numbers. The SA is significantly larger for the veins compared with the arteries for orders 2–10, which is similar to SD. The Sv for the veins was only significantly larger for orders 2–6. The curve fits of SA and Sv are summarized in Table 6.



View larger version (17K):
[in this window]
[in a new window]
 
Fig. 2. Relationship between the cross-sectional area ratio of daughter vessels and the order number of the mother vessel of the coronary vasculature.

 


View larger version (17K):
[in this window]
[in a new window]
 
Fig. 3. Relationship between the volume ratio of daughter vessels and the order number of the mother vessel of the coronary vasculature.

 
The ratio of the Poiseuille's resistance of the daughter branches, SR, was computed from the diameter (major and minor), length, and viscosity ratios as given by Eq. 6. The relationship between SR and the order number of the mother vessel is shown in Fig. 4. We found a significant effect of the mother vessel, such that as the order of the mother vessel increases, both the arteries and veins show a significant increase in asymmetry. The SR was significantly larger for the veins of orders 4–8 compared with arteries. Furthermore, the SR was significantly smaller (more asymmetric) for the thebesian compared with sinusal veins for orders –6 to –10. Table 6 summarizes the empirical constants of SR for a particular curve fit.



View larger version (19K):
[in this window]
[in a new window]
 
Fig. 4. Relationship between the Poiseuille's resistance ratio of daughter vessels and the order number of the mother vessel of the coronary vasculature.

 
The asymmetry can be described with respect to the order numbers at a bifurcation. The order asymmetry ratio between daughter vessels is given by {Delta} = (mk), as shown in Fig. 5. The deviation from symmetry ({Delta} = 0) is largest for the larger vessels. The value of {Delta} was significantly larger for the veins compared with arteries for orders 3 and 5–11. The order asymmetry between mother and larger daughter vessel is given by {varepsilon} = (n m) and is shown in Fig. 6. Similarly, it is clear that the asymmetry (deviation from {varepsilon} = 1) is greatest for the larger vessels.



View larger version (16K):
[in this window]
[in a new window]
 
Fig. 5. Relationship between the difference in the order between the two daughter vessels, {Delta}, and the order number of the mother vessel of the coronary vasculature.

 


View larger version (23K):
[in this window]
[in a new window]
 
Fig. 6. Relationship between the difference in the order between the mother and larger daughter vessel, {varepsilon}, and the order number of the mother vessel of the coronary vasculature.

 
The coefficient of variance (CV = SD/mean) was determined for the various asymmetry quantities for each order number of mother vessel. Figure 7 shows the CV of SD for the arterial and venous trees. The CV for the SD is significantly larger for the venous trees compared with the arterial trees.



View larger version (17K):
[in this window]
[in a new window]
 
Fig. 7. Relationship between the coefficient of variance, CV (SD/mean), and the order number of the mother vessel of the coronary vasculature.

 

    DISCUSSION
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
Asymmetry at a bifurcation. Circuitry asymmetry has many measures. In this study, we described vascular asymmetry at a bifurcation with the use of ratios of diameter and length of daughter vessels (SD and Sl, respectively). Our data show that the largest mean SD is ~10 and 12 for the largest arterial and venous bifurcations (Tables 1 and 5), respectively, which implies that the larger daughter vessel is on the average 10 and 12 times greater than its sister vessel.

As we have previously stated, SD ≥ 1 by definition. Sl, on the other hand, is not bound in this way because a larger diameter daughter vessel may be shorter than the smaller daughter vessel. The fraction of larger diameter daughter vessels that have longer lengths than the smaller vessels is shown in Fig. 1. It is interesting to note that for the smallest five orders, ~50% of the daughter vessels with larger diameters also have longer lengths. For the larger mother vessels, however, a higher percentage of daughter vessels of larger diameter have longer lengths.

Knowledge of SD and Sl allows the determination of SA, Sv, and SR of the two daughter branches of each mother vessel. The SA is simply the square of SD. However, the SA is not simply the square of mean SD because of the skewed shape of the distribution. The Sv is the product of SA and Sl. The SA and Sv are greatest for the larger vessels and decrease toward the smaller vessels. The mean asymmetry ratios of area and volume at the highest orders are ~150 and 560 (RCA) and 180 and 360 (sinus vein), as shown in Figs. 2 and 3, respectively. The enormity of the asymmetry is remarkable. Finally, the SR considers the SD, Sl, and Sµ as expressed by Eq. 6 and shown in Fig. 4. The SR is < 1 because a larger vessel has a lower resistance than a smaller vessel. The SD dominates over the Sl because of the fourth power exponent.

Comparison with other studies. Kassab et al. (9, 13) have previously described the bifurcation asymmetry of the right ventricular branches in normal and hypertensive right ventricles. They found that the SD and Sl increase significantly with order number for both control and hypertensive hearts. As expected, the control data on the right ventricular branches were similar to the present data for the RCA tree.

Horsfield et al. (6) introduced the concept of {Delta} to describe the branching asymmetry at a bifurcation. They showed that a fixed value of {Delta} implies a property of self-similarity, i.e., each succeeding bifurcation having the branching pattern and relative dimensions of the parent vessel but on a smaller scale. In the present study, we found that the value of {Delta} varies with order number. It is interesting to note that the first six orders of arterial trees have a fairly uniform value of {Delta} of ~1 but increase with higher orders (Fig. 5). We have added a new parameter, {varepsilon}, which describes the order asymmetry between mother and larger daughter vessel (Fig. 6). The collective data on {Delta} and {varepsilon} can be used to reconstruct the order connectivity of the vascular tree.

Relation of asymmetry ratios to design rules. We have previously predicted several structure-structure and structure-function relationships for coronary arterial trees (29). In the process, a vessel segment was defined as a stem that perfuses a tree or a crown (25). The crown includes all vessels proximal to the stem including the first segment of capillary vessels. If A and Q represent the mean cross-sectional area and blood flow of a stem, respectively, and V and L represent the cumulative distal arterial volume and length of a crown, respectively, we predicted the following relationships:

(8)

(9)

(10)
where Dmax, Qmax, Vmax, and Lmax correspond to the diameter and flow of the most proximal stem to the arterial tree, the volume of the entire tree of interest, and the cumulative arterial length of the entire tree. The values of {gamma}, {beta}, and {delta} have been previously reported as 0.41, 1.35, and 2.42 (RCA); 0.41, 1.37, and 2.45 (LAD); and 0.40, 1.38, and 2.47 (LCx) for the arterial trees, respectively (29). These relationships have been determined and validated in the following three ways: 1) a hemodynamic analysis of coronary arterial blood flow based on detailed anatomic data yields these relationships over the entire arterial network (12); 2) a generalization of Murray's cost function and conservation of energy predict the same result over the entire tree (29); and 3) in vivo data on flow, length and volume using digital subtraction angiography (28) verified the relationships for vessels proximal to 0.5 mm in diameter as observed in an angiogram.

To relate the asymmetry ratios of stem quantities (diameter, cross-sectional area, and flow) to crown quantities (crown length and volume), we combined Eq. 8 with Eq. 1 for an artery to yield

(11a)
or

(11b)
since {gamma}{delta} = 1, as shown previously (29). SL is the crown asymmetry length compared with the segment asymmetry length Sl. The exponent of Eq. 11 has a value of ~2.5. Hence, the asymmetry ratio for the crown length increases with SD raised to the power of 2.5. Similarly, the SV and SQ (Eqs. 9 and 10) can be expressed in terms of SD as:

(12)
and

(13)
respectively. The exponents of Eqs. 12 and 13 are ~3.4 and 2.5, respectively. Murray's law predicts a fixed value of 3 for {delta} irrespective of the organ of interest based on the minimum energy hypothesis (18).

We have also previously verified a form for the crown resistance (Rc), i.e., equivalent resistance of vessels in a crown, as

(14)
where kR and {varepsilon} are constants determined from experimental data (10). The reported value of {varepsilon} was approximately –0.55 (29). The asymmetry ratio of the crown resistance can be expressed as

(15)
Equations 11 and 12 can be combined with Eq. 15 to yield

(16)
The exponent is numerically equal to approximately –1.3, i.e., SRc is <1, similar to SR. We can relate the asymmetry of flow to that of crown resistance if we combine Eqs. 13 and 16 and use {gamma}{delta} = 1 as

(17a)
Because {delta} is ~2.5, the exponent of Eq. 17 is –0.52. Hence, for the coronary arterial tree, Eq. 17a can be reduced to

(17b)
which is a remarkably simple relationship between flow asymmetry at a bifurcation and the crown asymmetry.

A relation between SRc and SR can be determined if we combine Eq. 16 with Eq. 6 for a cylindrical vessel, i.e.,

(18a)
The exponent of Eq. 17 is obviously –0.32. Again, this can be simplified as

(18b)
In the microcirculation, at the site of vascular resistance, Sl and Sµ are ~1 (Tables 15) and hence SRc ~ SR(1/3). This remarkably simple relation implies that the crown resistance ratio at a bifurcation can be determined from the resistance ratio of the daughter vessels for coronary microvessels. In summary, the conclusion of the above analysis is that the crown ratios (Eqs. 11, 12, and 16) are significantly more asymmetric than the diameter asymmetry ratio.

Implications of SD on flow heterogeneity. The bifurcation asymmetry of resistance may not be directly proportional to the flow asymmetry because the flow distribution depends on the total distal resistance (crown resistance) of the two daughter vessels, as shown in Eq. 17. Hence, it is desirable to have a more direct relationship between the geometric asymmetry ratios and the flow distribution. Such relationship is expressed by Eq. 13, which shows a power-law dependence of flow asymmetry on the diameter asymmetry. The implication of Eq. 13 is that a small increase in diameter asymmetry ratio can cause a large increase in flow asymmetry. Hence, the heterogeneity of flow at a bifurcation is very sensitive to the diameter ratio of daughter vessels.

The CV can be estimated for the asymmetry of flow. If we take the differential of Eq. 13 and divide it by Eq. 13, we obtain

(19a)
or

(19b)
Hence, the CV(SQ) is proportional to the CV(SD) with {delta} (a value of 2.5) as the proportionality constant. We found that CV(SD) of veins was significantly larger than that of arteries. Hence, the CV(SQ) at a venous bifurcation is expected to be greater than that at an arterial bifurcation. This result has important implication for coronary venous retroperfusion where blood supply to the ischemic myocardium may be provided by making an arterial anastomosis with the coronary venous system (15). It suggests that the flow delivery to the capillary vessels will be more heterogeneous than antegrade flow. We have previously shown that the capillary network homogenizes the blood flow supplied by the arterial tree (14).

Critique of methods. Although the coronary arterial tree is dominated by bifurcations (98% bifurcations, 2% trifurcations) (10), the coronary veins consists of 86%, 12.8%, 1%, and 0.2% bifurcations, trifurcations, quadrifications, and quintifurcations, respectively (11). The present data apply only to the bifurcations but can be easily extended to trifurcations, etc. Furthermore, although the arteries have a fairly uniform diameter along their segment (i.e., diameters only change upon branching), the veins have more complex variations within segments. Because the measurement of venous diameter represents an average over the segment, the asymmetry ratio for the vein also represents an average in that sense. Finally, the asymmetry ratios may change with pressure because the coronary vessels are distensible. The present data correspond to a particular distribution of pressure through the diastolic state.

The present asymmetry ratios were expressed in terms of order numbers as determined by the diameter-defined Strahler system. Numerous ordering schemes exist, however, which may result in somewhat different asymmetry ratios. The Strahler's scheme was used because it handles the asymmetry of the system very well (22). The diameter-defined modification was proposed to eliminate the overlap in the diameters between successive orders (7, 10). In a previous study (9), we considered the diameter of the mother vessel as an independent variable regardless of the order number. We found that the diameter asymmetry ratios were similar in both cases.

The asymmetry in diameters also affects the asymmetry of bifurcation angles. This issue has been examined in detail by Zamir and colleagues (27). Zamir (26) showed that in a nonsymmetrical bifurcation, the branching angle that the larger branch makes relative to the parent vessel decreases with an increase in asymmetry ratio. The angle data affect the optimality of vascular junctions as flow conduits and dictate the spatial position of the vessels relative to the myocardium.

Although asymmetry at bifurcation is an important determinant of flow distribution, it is unlikely to be the only determinant of flow heterogeneity. In the beating heart, a heterogeneous distribution of stresses can affect the perfusion pattern (1). The addition of coronary tone plays a major role in modulating regional flow (1). Although the present anatomy was obtained under dilated state, this case is applicable when the heart is stressed physiologically or in the presence of coronary artery disease. Although the coronary anatomy establishes the potential distribution of coronary flow, mechanical, metabolic (3, 21), neural (24), and signal transduction (20) dictate local perfusion of myocardium.

Significance of study. The present study shows that the large asymmetry at a bifurcation can result in significant asymmetry of flow and potentially large spatial heterogeneity of blood flow. We have previously shown that such geometric asymmetry can remodel in right ventricular hypertrophy during flow overload (9, 13). Hence, the asymmetry parameters may be important indices of flow distribution in the heart in health and disease.

The lumen cross-sectional area and length of coronary arteries and their blood volume and flow rate can be determined from an angiogram using digital subtraction angiography (8, 16, 17). Hence, the SD, SL, SV, and SQ can be quantified and their relationships (Eqs. 1113) can be established as signatures of normal patients where coronary disease can be studied in relation to changes in these relationships. Furthermore, knowledge of SD or SQ allows the determination of SRc through Eqs. 16 or 17, respectively. The assessment of ratio of total arterial or crown resistance in a normal region compared with that in an ischemic or infarct region has obvious clinical significance.

An additional significance of the asymmetry ratios relates to the reconstruction of the full vascular arterial tree for blood flow analysis. The process of obtaining statistical morphometric data from a tree is unique, given that the same set of rules are followed each time (e.g., the diameter-defined Strahler system). The creation of a circuit from morphometric data, however, is a nonunique process. An infinite number of circuits can be created corresponding to a set of morphometric data but not uniquely specified by them. In addition to the connectivity and longitudinal position matrixes and the diameters, lengths, and number of vessels previously reported, the proposed asymmetry ratios impose additional constraints on the reconstruction of tree circuits.


    GRANTS
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 
This research was supported by the National Heart, Lung, and Blood Institute Grant 2 R01 HL-055554-06. G. S. Kassab is a recipient of an American Heart Association Established Investigator Award.


    ACKNOWLEDGMENTS
 
We thank Winny Tan and Marli Amin for excellent technical assistance.


    FOOTNOTES
 

Address for reprint requests and other correspondence: G. S. Kassab, Dept. of Biomedical Engineering, Univ. of California-Irvine, 204 Rockwell Engineering Center, Irvine, CA 92697-2715 (E-mail: gkassab{at}uci.edu)

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.


    REFERENCES
 TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 GRANTS
 REFERENCES
 

  1. Austin RE, Smedira NG, Squiers TM, and Hoffman JIE. Influence of cardiac contraction and coronary vasomotor tone on regional myocardial blood flow. Am J Physiol Heart Circ Physiol 266: H2542–H2553, 1994.[Abstract/Free Full Text]
  2. Bassingthwaighte J, RB King, and SA Roger. Fractal nature of regional myocardial blood flow heterogeneity. Circ Res 65: 578–590, 1989.[Abstract/Free Full Text]
  3. Caldwell JH, Martin GV, Raymond GM, and Bassingthwaighte JB. Regional myocardial flow and capillary permeability-surface area products are nearly proportional. Am J Physiol Heart Circ Physiol 267: H654–H666, 1994.[Abstract/Free Full Text]
  4. Glass G, Peckham PD, and Sanders JR. Consequences of failure to meet assumptions underlying the analysis of variance and covariance. Rev Educ Psychol 42: 237–288, 1972.
  5. Hoffman JIE. Heterogeneity of myocardial blood flow. Basic Res Cardiol 90: 103–111, 1995.[CrossRef][ISI][Medline]
  6. Horsfield K, Dart G, Olson DE, Filley GF, and Cumming G. Models of the human bronchial tree. J Appl Physiol 31: 201–217, 1971.
  7. Kassab GS. Morphometry of the Coronary Vasculature in the Pig (PhD dissertation). San Diego, CA: Univ. of California, 1990.
  8. Kassab GS and Molloi S. Cross-Sectional Area and volume compliance of the porcine left coronary arteries. Am J Physiol Heart Circ Physiol 281: H623–H628, 2001.[Abstract/Free Full Text]
  9. Kassab GS, Schatz A, Imoto K, and Fung YC. Remodeling of the bifurcation asymmetry of right ventricular branches in hypertrophy. Ann Biomed Eng 28: 424–430, 2000.[CrossRef][ISI][Medline]
  10. Kassab GS, Rider CA, Tang NJ, and Fung YCB. Morphometry of pig coronary arterial trees. Am J Physiol Heart Circ Physiol 265: H350–H365, 1993.[Abstract/Free Full Text]
  11. Kassab GS, Lin D, and Fung YC. Morphometry of the pig coronary venous system. Am J Physiol Heart Circ Physiol 267: H2100–H2113, 1994.[Abstract/Free Full Text]
  12. Kassab GS, Berkley J, and Fung YC. Analysis of porcine coronary arterial blood flow with detailed anatomical data. Ann Biomed Eng 25: 204–217, 1997.[ISI][Medline]
  13. Kassab GS, Imoto K, White FC, Rider CA, Fung YC, and Bloor CM. Coronary arterial tree remodeling in right ventricular hypertrophy. Am J Physiol Heart Circ Physiol 265: H366–H375, 1993.[Abstract/Free Full Text]
  14. Kassab GS, Lee KN, and Fung YC. A hemodynamic analysis of coronary capillary blood flow based on detailed anatomical and distensibility data. Am J Physiol Heart Circ Physiol 277: H2158–H2166, 1999.[Abstract/Free Full Text]
  15. Mohl W. The relevance of coronary sinus interventions in cardiac surgery. J Card Surg 39: 245–250, 1991.
  16. Molloi S, Zhou Y, and Kassab GS. Regional volumetric coronary blood flow measurement by digital angiography: in-vivo validation. Acad Radiol 11: 757–866, 2004.[ISI][Medline]
  17. Molloi S, Kassab GS, and Zhou Y. Quantification of coronary arterial volume by digital angiography: In vivo validation. Circulation 104: 2351–2357, 2001.[Abstract/Free Full Text]
  18. Murray CD. The physiological principle of minimum work. I. The vascular system and the cost of blood volume. Proc Natl Acad Sci USA 12: 207–214, 1926.[Free Full Text]
  19. Pries AR, Neuhaus D, and Gaehtgens P. Blood viscosity in tube flow: dependence on diameter and hematocrit. Am J Physiol Heart Circ Physiol 263: H1770–H1778, 1992.[Abstract/Free Full Text]
  20. Ramirez J, Rodriguez E, and Weiss HR. Effect of isoproterenol on the relationship between regional cyclic AMP and local coronary blood flow in the canine myocardium. Pharmacology 48: 41–48, 1994.[ISI][Medline]
  21. Sonntag M, Deussen A, Schultz J, Loncar R, Hort W, and Schrader J. Spatial heterogeneity of blood flow in the dog heart. I: glucose uptake, free adenosine and oxidative/glycolytic enzyme activity. Pflügers Arch 432: 439–450, 1996.[CrossRef][ISI][Medline]
  22. Strahler AN. Hyposmetric (area altitude) analysis of erosional topology. Bull Geol Soc Am 63: 1117–1142, 1952.[Abstract]
  23. Van Beek JHGM, Roger SA, and Bassingthwaighte JB. Regional myocardial flow heterogeneity explained with fractal networks. Am J Physiol Heart Circ Physiol 257: H1670–H1680, 1989.[Abstract/Free Full Text]
  24. Verhoeven L, Gerlock T, and Weiss HR. Relationship between myocardial tissue norepinephrine and coronary flow heterogeneity. Microvasc Res 41: 289–298, 1991.[CrossRef][ISI][Medline]
  25. Wahle A, Wellnhofer E, Mugaragu I, Sauer HU, Oswald H, and Fleck E. Quantitative Volume Analysis of Coronary Vessel Systems by 3-D Reconstruction From Biplane Angiograms. San Francisco, CA: IEEE Medical Imaging Conference, 1993, p. 1217–1221.
  26. Zamir M. Nonsymmetrical bifurcations in arterial branching. J Gen Physiol 72: 837–845, 1978.[Abstract/Free Full Text]
  27. Zamir M, Phipps S, Langille BL, and Wonnacott TH. Branching characteristics of coronary arteries in rats. Can J Physiol Pharmacol 62: 1453–1459, 1984.[ISI][Medline]
  28. Zhou Y, Kassab GS, and Molloi S. In vivo validation of the design rules of the coronary arteries and their application in the assessment of diffuse disease. Phys Med Biol 47: 977–993, 2002.[ISI][Medline]
  29. Zhou Y, Kassab GS, and Molloi S. On the design of the coronary arterial tree: a generalization of Murray's law. Phys Med Biol 44: 2929–2945, 1999.[CrossRef][ISI][Medline]



This article has been cited by other articles:


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
B. Kaimovitz, Y. Huo, Y. Lanir, and G. S. Kassab
Diameter asymmetry of porcine coronary arterial trees: structural and functional implications
Am J Physiol Heart Circ Physiol, February 1, 2008; 294(2): H714 - H723.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
Y. Huo and G. S. Kassab
A hybrid one-dimensional/Womersley model of pulsatile blood flow in the entire coronary arterial tree
Am J Physiol Heart Circ Physiol, June 1, 2007; 292(6): H2623 - H2633.
[Abstract] [Full Text] [PDF]


Home page
Physiol. Rev.Home page
N. Westerhof, C. Boer, R. R. Lamberts, and P. Sipkema
Cross-talk between cardiac muscle and coronary vasculature.
Physiol Rev, October 1, 2006; 86(4): 1263 - 1308.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
G. S. Kassab
Functional hierarchy of coronary circulation: direct evidence of a structure-function relation
Am J Physiol Heart Circ Physiol, December 1, 2005; 289(6): H2559 - H2565.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (6)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kalsho, G.
Right arrow Articles by Kassab, G. S.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kalsho, G.
Right arrow Articles by Kassab, G. S.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online
Copyright © 2004 by the American Physiological Society.