Am J Physiol Heart Circ Physiol 292: H3006-H3018, 2007.
First published February 16, 2007; doi:10.1152/ajpheart.01012.2006
0363-6135/07 $8.00
Shear stress paradigm for perinatal fractal arterial network remodeling in lambs with pulmonary hypertension and increased pulmonary blood flow
Zahra Ghorishi,1
Jay M. Milstein,1
Francis R. Poulain,1
Anita Moon-Grady,1,2
Theresa Tacy,2
Stephen H. Bennett,1
Jeffery R. Fineman,2 and
Marlowe W. Eldridge3
1Neonatology Division, Department of Pediatrics, University of California, Davis and 2Department of Pediatrics, University of California, San Francisco, California; and 3Department of Pediatrics, Population Health Sciences and Biomedical Engineering, University of Wisconsin, Madison, Wisconsin
Submitted 14 September 2006
; accepted in final form 12 February 2007
 |
ABSTRACT
|
|---|
Congenital heart disease with increased blood flow commonly leads to the development of increased pulmonary vascular reactivity and pulmonary arterial hypertension by mechanisms that remain unclear. We hypothesized a shear stress paradigm of hemodynamic reactivity and network remodeling via the persistence and/or exacerbation of a fetal diameter bifurcation phenotype [parent diameter d0 and daughters d1
d2 with
< 2 in (d1/d0)
+ (d2/d0)
and area ratio
< 1 in
= (d12+ d22)/ d02] that mechanically acts as a high resistance magnifier/shear stress amplifier to blood flow. Evidence of a hemodynamic influence on network remodeling was assessed with a lamb model of high-flow-induced secondary pulmonary hypertension in which an aortopulmonary graft was surgically placed in one twin in utero (Shunt twin) but not in the other (Control twin). Eight weeks after birth arterial casts were made of the left pulmonary arterial circulation. Bifurcation diameter measurements down to 0.010 mm in the Shunt and Control twins were then compared with those of an unoperated fetal cast. Network organization, cumulative resistance, and pressure/shear stress distributions were evaluated via a fractal model whose dimension D0
delineates hemodynamic reactivity. Fetus and Control twin D0 differed: fetus D0 = 1.72, a high-resistance/shear stress amplifying condition; control twin D0 = 2.02, an area-preserving transport configuration. The Shunt twin (D0 = 1.72) maintained a fetal design but paradoxically remodeled diameter geometry to decrease cumulative resistance relative to the Control twin. Our results indicate that fetal/neonatal pulmonary hemodynamic reactivity remodels in response to shear stress, but the response to elevated blood flow and pulmonary hypertension involves the persistence and exacerbation of a fetal diameter bifurcation phenotype that facilitates endothelial dysfunction/injury.
pulmonary arterial morphometry; branching complexity; fractals
THE PULMONARY CIRCULATION undergoes a dramatic hemodynamic change at birth, from a state of high pressure/low flow to one of low pressure/high flow (50). This change occurs in conjunction with a decrease in pulmonary arterial resistance (8) and a transformation in the complexity of arterial diameter network organization (3). However, in cases of congenital heart disease with a systemic-to-pulmonary communication, an abnormally increased flow is superimposed on the usual decrease in resistance and commonly leads to development of pulmonary hypertension with functional and structural complications (25). While the mechanisms are not completely understood, the present consensus is that endothelial injury and dysfunction increase vascular tone and reactivity, leading to an increase in pulmonary artery pressure that is initially considered reversible (20). Pulmonary hypertension becomes problematic when vascular wall remodeling leads to the appearance of smooth muscle and alterations of medial wall thickness that encroach upon the lumen diameter (20). Since the forces of pressure and shear stress are known to be transduced and translated into molecular mechanisms that are instrumental in regulating diameter, vascular tone, and the remodeling of the vascular wall (11), further understanding of the factors aggravating these forces within the pulmonary arterial circulation should aid in understanding patterns of reactivity and the demarcation of reversibility in the disease process.
In general, mammalian vascular networks are thought to possess a universal fractal pattern of diameter network organization that appears independent of organ system, body size, and hemodynamic state (32, 33, 5557). However, the fetal pulmonary circulation is extraordinary. Fetal arterial bifurcations appear to possess a hemodynamic phenotype that is distinguished by an area ratio <1 (3), not >1 as typically found in the adult pulmonary circulation (28). This curious vessel diameter topology has a design that mechanically aggravates hemodynamic reactivity per unit blood flow, leading to an elevated driving pressure, a magnified resistance, and the concomitant amplification of shear stress in the course of branching (3). Although intended to be rapidly adapting to an adult pattern via perturbations in blood flow and shear stress in the fetal-neonatal transition after birth, the shear stress-amplifying characteristics of the fetal network design may be precarious to the pulmonary circulation in the presence of an elevated blood flow. Shear stress is a force known to influence the remodeling of vascular diameter and/or vascular wall properties in response to transient and chronic changes in blood flow (35). However, acute elevations beyond critical levels can induce endothelial dysfunction and/or injury (19, 35, 54). Thus the interaction of the fetal phenotype with amplified shear stress secondary to elevated blood flow, or any other source of fetal induced endothelial injury and/or pattern of dysfunction (7), may cause the fetal phenotype to persist or become exacerbated in the presence of elevated blood flow after birth (25).
While the changes in vascular diameter geometry and vascular wall morphology of the pulmonary circulation are well understood, the structural remodeling patterns of arterial network organization in response to the hemodynamic forces of pressure and shear stress are not. We hypothesized that hemodynamic reactivity and network remodeling via the persistence and/or exacerbation of a fetal diameter bifurcation phenotype might mechanically act as a high resistance magnifier/shear stress amplifier to blood flow. Therefore, one purpose of this study was to test the hypothesis that a fetal phenotype of arterial branching persists and/or is exacerbated after birth. We utilized a well-established lamb model of congenital heart disease with increased pulmonary blood flow, in which an aortopulmonary graft is surgically placed in one twin in utero (Shunt twin) but not in the other twin (Control twin) (40, 44). Diameter measurements in lung casts were made to characterize longitudinal branching patterns and design of bifurcations from the left pulmonary artery (LPA) in a Shunt twin and its Control twin, relative to an unoperated fetal lung (3). In addition, as the hemodynamic forces of pressure and shear stress are linked to multiple mechanisms of vascular diameter and wall remodeling, our second objective was to evaluate the fractal geometry, topology, and design of the fetal, Shunt, and Control pulmonary arterial networks as a means of predicting the distribution of pressure and shear stress, and to determine whether or not their respective fractal dimensions conform to a universal design independent of hemodynamic forces.
 |
METHODS
|
|---|
Cast preparation.
The respective studies were approved by animal use protocols issued by the University of California Davis (3) and the University of California San Francisco (40). Lung casts of 8-wk-old twin lambs and a term fetus (140 days gestation, 145 days term) were made from experimental animals originating from previous studies (3, 40). In the aortopulmonary shunt study, mixed-breed pregnant Western ewes with twin lambs (between 137 and 141 days gestation, term = 145 days) were operated on under sterile conditions (44). One twin had an aortopulmonary shunt surgically implanted as previously described (40) (Shunt twin). The other lamb was not operated on (Control twin). As previously described (40), after spontaneous birth the lambs were kept with their mothers and weighed daily, and respiratory rate and heart rate were measured. In Shunt lambs furosemide (1 mg/kg im) was administered daily and elemental iron (50 mg im) was given weekly. Eight weeks after delivery the lambs were euthanized with an intravenous injection of pentobarbital sodium (Euthanasia CII; Central City Medical, Union City, CA) and subjected to bilateral thoracotomy. In one set of twins, lung casts were made of the left pulmonary arterial circulation, while the right pulmonary arterial circulation was harvested for tissue. Lung casts were made as follows: a 4.5-Fr cannula was inserted into the proximal LPA. The airways were expanded with saline under a hydrostatic pressure gradient of 20 cmH2O via a tracheal cannula. The pulmonary arterial circulation was washed free of blood by saline perfusion. To remain consistent with the morphometric diameter measurement methods utilizing gelatin-barium-filled arteries, a high-viscosity methyl methacrylate plastic (Coe Tray Plastic, GC America, Chicago, IL) was injected into the pulmonary arterial system slowly with a syringe over a 2-min period under a maximum pressure of 75 mmHg. The plastic was allowed to polymerize overnight, whereupon the lung tissue was macerated in a 20% KOH bath for 35 days. The remaining tissue was washed away gently with distilled water, and the cast was allowed to dry before bifurcation diameter measurements.
Diameter measurement.
A large number of diameters within bifurcations of the left pulmonary arterial tree were sampled according to previous methods (3). In the Shunt and Control twins, the samples consisted of diameters of bifurcations along the LPA branch, plus an extensive sampling of bifurcations within the upper cranial lobe. Bifurcation diameters were measured with a methodology described previously (3). Briefly, a video microscope system (Infinivar Video Microscope with Zoom, Infinity Photo-optical, Boulder, CO) connected to a Cohu model 2600 solid-state camera (Cohu, San Diego, CA) and to a Macintosh computer (Cupertino, CA) was used to acquire and analyze images via a program (Object Image) developed by Norbert Vischer of the Faculty of Science, University of Amsterdam, Amsterdam, The Netherlands (http://simon.bio.uva.nl). Object Image is an extended version of the program NIH Image (http://rsb.info.nih.gov/nih-image), which allows diameter measurements and derived calculations to be recorded into an extensive database. To measure diameters, the three-dimensional branching aspect of each cast was broken into pieces to facilitate placement of bifurcations onto a plane. Assuming vessels were circular, vessel diameter was calculated as the average of two diameter measurements per vessel segment. The three averaged diameter measurements, consisting of a parent vessel, d0, and two daughter vessels, d1 and d2, where d0 > d1
d2, were stored in a database and used to make functional morphometric calculations within bifurcations.
Functional morphometric relationships.
The measurements d0, d1, and d2 were utilized to calculate functional morphometric parameters that characterize local hemodynamic reactivity, via diameter geometry, topology, and design of bifurcations within the pulmonary arterial tree (APPENDIX; Ref. 3). Diameter branching geometry is defined by the diameter ratio
 | (1) |
along with the asymmetry ratio
= d2/d1
1. The area ratio
 | (2) |
designates the topology between parent and daughter diameters in terms of the relative cross-sectional area available for perfusion. The design parameter
was iteratively solved (28) via the equation
 | (3) |
The values of
and/or
have hemodynamic significance conferring reactivity to bifurcation network organization (Fig. 1; APPENDIX): a high-resistance/shear stress-amplifying fetal design (3) occurs when
< 2 and/or
< 1; a lower resistance transport design (57) is evident when
2 and/or
1.

View larger version (24K):
[in this window]
[in a new window]
|
Fig. 1. Impact of diameter network organization on magnifying cumulative resistance (Rcum) and amplifying shear stress ( ) to blood flow with branching order as seen in a symmetric dichotomous fractal arterial network via different area ratios ( ) and fractal dimensions (D0). A: hemodynamic reactivity ( ), defined as the tendency to magnify and amplify properties with arterial branching, has interrelationships with the scaling of diameter geometry (via parent d0 and daughter d1 = d2 vessels), topology [via the area ratio = (d
12+ d22)/d02= 2d12/d02], and fractal design D0 = 2/(1 ln /ln2) according to = 2R1/R0 = 2 1/ 0 = ( )3D0/2D0(APPENDIX). Hemodynamic reactivity spans a spectrum of network designs and topologies that encompass different theoretical laws of network form and function: West-Brown-Enquist law D0 = 2.0 (57); Kurz-Sandau law D0 = 2.7 (32); Murray's law D0 = 3.0 (33, 51). Reactivity is neutral ( = 1) when the critical area ratio c = 21/3 and D0 = 3 is achieved; increases ( > 1) when < c; and decreases ( < 1) when < c. B: influence of on the cumulative resistance/branching order Rcum(k) in successive orders of branching (k = 1, 2,..., 5, designated by different symbols: , order 1; , 2; , 3; , 4; , 5), normalized by parent resistance R0: Rcum(k)/R0 = 1 + + 2 + ... d. C: influence of on amplifying shear stress in successive vessel orders k = 1, 2,..., 5 via (k)/ 0 = k.
|
|
Diameter ordering/analysis of arterial tree.
We categorized the parent diameters (d0) of bifurcations into orders within the arterial tree via a diameter ranking method that calculates an average diameter per order,
k, as previously described (3). The diameter data set for each bifurcation (d0, d1, d2,
, and
) was first ranked from the largest parent diameter
0 (designated order 0) down to the smallest diameter, by orders k = 0, 1... n categorized by the overall average diameter ratio
d according to
dk
(di/
0) >
d(k+1). Here, di is a parent diameter of a bifurcation in the sorted list, d0 is the diameter of the largest (main or left) pulmonary artery, and
d represents the average diameter ratio, computed from the pool of Rd calculated from bifurcations. Once orders were assigned to each parent diameter, the average values of diameter
k, log
k, and log
k of each order in the arterial tree were calculated. The ranking procedure results in a nonoverlapping range of average parent diameters for each vessel diameter order, where each order is assumed to be statistically independent from another but where the log-transformed values of
k and
k are assumed to possess a statistically self-similar Gaussian distribution with bifurcation level independent of the average value of
or
(3).
Fractal model.
The average values of
k and
k for parent diameters within orders, from the largest diameter artery
0 down to terminal arteries
0.010 mm, were used to synthesize an equivalent homogeneous fractal branching tree to summarize and compare the hemodynamic reactive properties of arterial network organization (APPENDIX). The fractal properties of the arterial tree are summarized by global ratios: a network branching ratio
b, a network diameter ratio
d, and a fractal dimension D0, derived from the distribution of average
k and
k in the following way. To determine the number of vessels with each order of branching, a local branching ratio for each order was determined by
 | (4) |
where Rdk is the local diameter ratio between orders derived from
 | (5) |
and
k is the average value of
for the order k. The number of vessels per order was calculated according to Nd = N0Rb1Rb2... Rbd. The global network branching ratio
b was determined from the numbers of vessels per order via Nk = N0
bk, computed via the slope of the linear regression relationship between vessel number N and order
 | (6) |
where
b is the antilog of the slope and the intercept was comprised of the main parent N0 = 1. The global network diameter ratio
d, determining the decrement in average diameter with order
k =
0
dk and representative of an average of all local diameter ratios Rdk between orders (Eq. 5), was evaluated by linear regression of the logarithm of average diameter
k vs. order k,
 | (7) |
where
d is the negative of the antilog of the regression slope. The fractal dimension of the model, D0, designating the design of the arterial tree, was then evaluated from the power law relationship between vessel number Nk = N0
dkD0 and average diameter
k =
0
dk via the slope of the linear regression of
 | (8) |
where the fractal dimension D0 is equal to the negative of the regression slope. The values of the fractal dimension D0, the network diameter ratio
d, and the network branching ratio
b in the fetal, Control twin, and Shunt twin lung casts were then used to calculate their respective cumulative resistance and shear stress distribution (APPENDIX; Fig. 1). The cumulative resistance (Rcum) for each order k (0k+1) down to vessels
0.020 mm was calculated according to
 | (9) |
The hydraulic resistance of the LPA to blood flow was computed via Poiseuille's law,
 | (10) |
where µ is blood viscosity (assumed to be 3.0 cP), l0 is vessel length, and
0 is the LPA vessel diameter. Equation 9 assumes the presence of a length-diameter scaling relationship
k = a
k
2 where a and
2 are constants (30). We assumed that
2 = 1 (30) and that a = 1.7, where the constant a was evaluated on the basis of human pulmonary data (52) ordered according to the Horsfield method (27). The rectilinear pressure response to an increase in flow (13) was calculated according to
 | (11) |
where q is pulmonary arterial flow in the parent vessel normalized to body weight and pint is the intercept pressure, set to 25 mmHg in the fetal state (46) and 6 mmHg in the postnatal state. The associated pressure distribution in the fractal continuum equivalent tree at a steady-state flow qss was evaluated according to
 | (12) |
The distribution of shear stress with branching order k was calculated according to
 | (13) |
where
 | (14) |
and qss is the absolute steady-state blood flow (ml/min) with
0 the diameter of the LPA. Steady-state values for flow q used in the calculations were taken from previously published studies in which LPA flow was measured: 139-day fetal lambs/145 days term (with average body wt 5.0 kg), where q = 9.5 ml·min1·kg1 or q = 38 ml/min (46); 8-wk-old Control twins (average body wt 22.6 kg), where q = 33 ml·min1·kg1 or 743 ml/min; and 8-wk-old Shunt twins (Ref. 40; average body wt 18.9 kg), where shunt open q =115 ml·min1·kg1 or 2,200 ml/min. The Reynolds number per order k associated with a given flow to the network was calculated according to Re(k) = Re(0)Rd(1D0)k, where Re(0) = q
/µd0, with the density of blood
= 1.05 g/cm and blood viscosity µ = 0.03 g·cm1·s1 (APPENDIX).
Statistical analysis.
The distributions of
and
exhibit log-normal behavior; therefore, a two-way ANOVA was performed on log-transformed values of
and
to evaluate differences between the fetal, Shunt twin, and Control twin data. One factor was the experimental group (Control twin, Shunt twin, fetus), and the other factor was branching order k. Our null hypothesis was that the mean values of
and
, over all orders, are equal to a transport design (
= 2 and
= 1), considered to be a universal theoretical pattern of arterial vascular network organization (57). Results of the analysis of the log-transformed values of
and
are reported as means ± 95% confidence intervals (CI) (58), expressed in terms of their geometric mean adjusted values (37). The regression slopes D0 of Eq. 8 derived from the lung casts are also reported in terms of 95% CI. Significant differences comparing branching order design, topology, and group fractal dimensions were distinguished by nonoverlapping CI, indicative of a P > 0.05 level of significance. Statistical tests were performed with JMP 6.0.
 |
RESULTS
|
|---|
A total of n = 46,673 bifurcation diameter measurements for
and
were evaluated in the three lung casts (Control twin n = 16,443 bifurcations, Shunt twin n = 19,816 bifurcations, and fetus n = 10,414 bifurcations), where 1 bifurcation = 1 parent diameter connected to 2 daughter diameters.
Distributions of arterial design (
) and topology (
) with diameter and daughter asymmetry (
).
ANOVA revealed significant differences in the averages between experimental groups and within orders for (
,
) at a P < 0.0001. In the fetal arterial cast
=1.736 (CI 1.661.815) and
= 0.846 (CI 0.8160.878); in the Control twin cast
=1.94 (CI 1.882.00) and
= 0.934 (CI 0.9110.957); and in the Shunt twin cast
= 1.671 (CI 1.601.744) and
= 0.850 (CI 0.8210.881). These results indicate that, relative to the fetal cast, the control cast altered its bifurcation design and topology toward a West-Brown-Enquist design. However, the Shunt twin demonstrated the persistence of a fetal design and topology.
Figure 2 compares conventional histologically derived indices of vascular-diameter reactivity [percentage of smooth muscle and medial thickness-diameter ratio measured previously in this experimental model (40)] relative to indices of hemodynamic reactivity, average
and
, relative to average diameter and branching order. The fetal-reactive regions (
< 1,
< 2) are colored gray in Fig. 2, while the white area (
1,
2) indicates a hemodynamically low-reactive adult transport criterion of bifurcation branching. Statistically significant differences of
and
between groups within orders according to ANOVA are earmarked within the Shunt twin markers as a white circle within a black circle in Fig. 2. Significant differences between any two groups within an order are characterized by nonoverlapping CI. Overlapping CI were observed for orders 06 in the Control twin and Shunt twin and orders 02 in the fetus, and they are not shown in Fig. 2. Figure 2 illustrates that the fetal pulmonary circulation exhibits a regional mechanical source of hemodynamic reactivity in intra-acinar arteries [0.020-mm to 0.100-mm parent inner diameter (ID)], a region we previously showed to be devoid of smooth muscle(40). The Shunt twin lamb demonstrated a similar fetuslike design profile of
< 2 and
< 1 in preacinar arteries >0.150 mm. In addition, intra-acinar Shunt twin arteries (0.020- to 0.040-mm ID), which we previously showed to demonstrate the appearance of increased smooth muscle and medial wall thickness (40), were associated with significant differences in hemodynamic reactivity (smaller
and
) relative to both the Control twin and the fetus.

View larger version (36K):
[in this window]
[in a new window]
|
Fig. 2. Comparison of classic histological indices of vascular reactivity (A) and hemodynamic reactivity (B and C) in 8-wk-old Control twin and Shunt twin lambs, both indices relative to the distribution of vascular diameter. A: vascular reactivity histologically is defined by appearance and/or elevation in % smooth muscle and/or increased % medial wall thickness relative to diameters ordered via accompanying airways classified as terminal bronchiole (TB), respiratory bronchioles (RB), alveolar ducts (AD), and the alveolar wall (AW) (40). B and C: hemodynamic reactivity expressed by degree that the average values of (B) and (C) [±95% confidence intervals (CI)] deviate from a universal transport design state is demarcated by = 2, = 1. ANOVA indicated that the average and were significantly different between groups (fetus, Shunt twin, and Control twin) and between branching orders, indicating regional variations of hemodynamic reactivity. Significant differences within a branching order are indicated by a white dot superimposed upon the Shunt indicator (black dots). Hemodynamically reactive values ( < 2, < 1) are evident before the histological appearance of vascular reactivity in the fetal condition. Differences in and between fetal and Control twin indicate that the pulmonary arterial circulation adapts to hemodynamic conditions of a high-flow/low-pressure state after birth by remodeling area ratios in a universal direction, toward area ratios equal to 1, which mitigates the hemodynamic reactivity of the fetal phenotype. However, the shunt exacerbates the fetal design in preacinar arteries and in small intra-acinar arteries, whose overall average indicates a persistence of the fetal configuration.
|
|
The fetal cast demonstrated significantly different average values of diameter asymmetry
relative to the Control twin and the Shunt twin: Fetus 0.62 (SD 0.20), Control twin 0.71 (SD 0.18); Shunt twin 0.70 (SD 0.19). Figure 3 shows the interrelationship between area ratio
, bifurcation design exponent
, and diameter asymmetry
as a bivariate frequency distribution. The overall range and spread of
between the casts are similar and substantial. The fetal lung cast bifurcations show a predominant frequency toward acute asymmetry (
0.5), with area ratios favoring 0.8 <
< 1 rather than
> 1. In the Control twin and Shunt twin casts, the greatest frequency of asymmetry occurred in a higher range, 0.6 <
< 0.8. The area ratios of the Control twin demonstrate their greatest frequencies about 0.8 <
< 1.2, while the area ratios in the Shunt twin included relatively larger numbers of bifurcations 0.6 <
< 0.8.
Fractal power law properties.
The power law behavior synthesized via Eqs. 47 for the fetal, Control twin, and Shunt twin lung casts is shown in Fig. 4.
In Fig. 4A, the diameter ratio (
d) of the Shunt twin was significantly different from that of the fetal and Control twin casts: Shunt twin
d = 1.41 ± 0.015 (CI), fetus
d = 1.30 ± 0.003 (CI), and Control twin
d = 1.31 ± 0.001 (CI). While the difference between the fetus and the Control twin was statistically significant, it was small. The relatively small difference between the fetus and Control twin cast
d values, plus the parallel lines in average diameter with order in Fig. 4A, indicate that growth and development scaled diameters upward with branching order. However, in the Shunt twin the average arterial diameters for orders 08 were larger relative to the Control twin or fetus cast, while orders 910 of the Shunt and Control twins were similar. Diameters of vessel orders 11 and larger were smaller in the Shunt twin compared with the Control twin. Alternatively, the fractal regression indicates that the intersection between Control twin and Shunt twin occurs at order 11, at a diameter of 0.5 mm. Thus the differences between the Shunt twin and other casts indicate that the Shunt twin remodeled and rescaled diameter geometry about a diameter pivot region of 0.51.0 mm, leading to smaller resistances in larger arteries and elevated resistances in smaller arteries.
In Fig. 4B, the branching ratios
b and their respective 95% CI were statistically different from one another: fetus
b = 1.572 ± 0.03 (CI), Shunt twin
b = 1.79 ± 0.03 (CI), and Control twin
b = 1.73 ± 0.002 (CI).
In Fig. 4C, the fractal dimension ± CI of the fetus was D0 = 1.72 ± 0.069, Shunt twin D0 =1.70 ± 0.038, and Control twin D0 = 2.02 ± 0.021. Thus the Shunt twin was not statistically different from the fetus, but both were both significantly different from the Control twin, indicating the persistence of a fetal fractal network design in the arterial tree at 8 wk of shunting.
Together these results indicate that the influence of a chronic shunt interrupts an apparent transition in arterial network organization at birth: the fetal-Control transition maintained the scaling of arterial diameter geometry uniformly increasing diameters but remodeled arterial design and topology to an adult transport design. Shunting remodels diameter geometry by rescaling large and small diameters to where large arterial vessels increase diameter to decrease resistance and small arterial vessels decrease diameter to increase resistance, while leaving the topology and fractal design of the pulmonary circulation in a hyperreactive fetuslike configuration.
Cumulative resistance distribution.
Figure 5 illustrates the influence of fractal network organization in the fetal, Control twin, and Shunt twin casts on the cumulative resistance distribution versus diameter, as calculated from the regression values of
d and D0 according to Eqs. 9 and 10. In general, the smaller fractal dimension (D0 = 1.7) of the fetus and the Shunt twin demonstrates greater resistance reactivity with arterial branching, evidenced by a larger slope of cumulative resistance with decreasing diameter. However, the Shunt twin demonstrates an absolute resistance that is less than the fetus and the Control twin: the remodeling of arterial geometry (R0,
b) counteracts the increased hemodynamic reactivity of a smaller fractal dimension by decreasing the absolute resistance per order, leaving a smaller cumulative arterial resistance.

View larger version (13K):
[in this window]
[in a new window]
|
Fig. 5. Cumulative resistance of the fetus, Control twin, and Shunt twin fractal network in terms of diameter (d0) or order calculated from power law relationships of Fig. 4 according to Eq. 9. The Shunt and fetal lambs have similar slopes of Rcum change/order relative to the Control lamb because of smaller fractal dimensions of design. The values of R0 and Rd determine the vertical offsets of the Rcum graph: the fetus has the largest cumulative resistance with branching because of the size scale of the left pulmonary artery (LPA) and increased reactivity due to the smaller fractal dimension with arterial branching. The Shunt twin has a smaller cumulative resistance for all parent diameters of branching, which arises from the geometric scaling leading to a larger LPA diameter d0 = 2.12 cm and larger diameter ratio (Rd = 1.41) increasing diameter in orders 010 (Fig. 3A) relative to Control (LPA d0 = 1.15 cm and Rd = 1.3) and fetus (LPA d0 = 0.53 cm and Rd = 1.31). This graph indicates that by 8 wk hemodynamic forces of hypertension and increased blood flow of the shunt remodel vascular arterial diameter geometry (R0 and Rd) in a manner that paradoxically decreases the cumulative resistance relative to control (44), but maintains the hemodynamic resistance reactivity of the fetal state (slope), via the persistence of the fetal topology and design (fractal dimensions of fetus = Shunt).
|
|
Pressure-flow curves/pressure distribution.
Figure 6 illustrates the influence of Rcum on the predicted pressure-flow curves (Fig. 6A) of the LPA circulation in the fetal, Control twin, and Shunt twin fractal configurations and (Fig. 6B) the corresponding distribution of pressure versus diameter at a given steady states of flow. The fetal state is the most pressure reactive, yielding hypertensive pressures in all arterial regions over a small flow range. Fractal network organization in the Shunt twin configuration is the least pressure reactive to flow of the three states, with hypertension arising secondary to the increased steady-state blood flow. As Fig. 6 indicates, the Control twin state also predicts pulmonary hypertension at a flow comparable to the Shunt twin, but with higher driving pressures and higher values in vessels >0.060 mm. Conversely, arterial network remodeling in the shunt condition predicts that shunt closure will result in nonhypertensive pressures with smaller pulmonary vascular resistances than the corresponding Control twin.

View larger version (19K):
[in this window]
[in a new window]
|
Fig. 6. A: influence of Rcum on pressure-flow relationships and pressure distribution vs. diameter. Rcum influences the slope of pressure-flow curves determining the pressure reactivity to a given increment in flow about its steady-state value [LPA steady-state flows Qlpa: fetus = 10 (46), Control twin = 33 (40), Shunt twin = 100 (40)]. The fetus has the greatest reactivity, followed by the Control twin and then the Shunt twin. Paradoxically, the remodeling of diameter geometry, topology, and design in the Shunt twin decreases the pressure response to all increments in flow. Increasing steady-state flow in the pulmonary arterial circulation to some threshold results in a clinical definition of pulmonary hypertension (Ppa = 25 mmHg). Structural remodeling in the Shunt twin elevates the steady-state flow threshold to 80 ml·min1·kg1, compared with the lower threshold of 50 ml·min1·kg1 in the Control twin. B: comparison of steady-state pressure-diameter profiles computed at low flow (33 ml·min1·kg1) and high flow (100 ml·min1·kg1) in the Control and Shunt twins. The structural remodeling in the Shunt twin lamb diminishes the hypertensive pressure stimulus compared with the Control twin at high flow.
|
|
Shear stress distribution.
Figure 7 illustrates the influence of diameter network organization on the amplification of shear stress at steady-state blood flows in the fractal model. The calculated values for
0 shear stress and Reynolds numbers in the LPA for the different states were as follows: fetus, LPA d0 = 5.4 mm at a flow of q = 38 ml/min yields
0 = 1.25 dyn-cm2, Re(0) = 82; Shunt twin with the shunt-open condition, LPA d0 = 21.2 mm and flow q0 = 2,200 ml/min yields
0 = 1.16 dyn-cm2 at Re(0) = 774; Control twin, LPA d0 = 11.5 mm and q0 =743 ml/min yields
0 = 2.51 dyn-cm2, Re(0) = 486. The common feature in all groups is that extremely large shear stress avalanches are developed with low Reynolds numbers (Re < 10) in arteries of 0.020 mm. For comparison purposes, the gray area in Fig. 7 represents the reported ranges of critical and/or erosion yield shear stresses (1658,000 dyn-cm2) of endothelial dysfunction/damage to an acute increase in shear stress as it occurs in large-diameter conduit vessels over a range of Reynolds numbers (0.212,300) (19, 35, 54). Interestingly, the fetal fractal model (
0.020 = 1,200 dyn-cm2 at Re0.020 = 1.4) and Control twin model (
0.020 = 1,100 dyn-cm2 at Re0.020 = 0.43) demonstrated similar magnitudes of shear stress experienced in arterial vessels of 0.020-mm ID as derived from their respective steady-state blood flows. However, the Shunt twin steady-state values of blood flow yielded a magnitude of shear stress at
0.010 =7,100 dyn-cm2 at Re0.020 = 2.6, which extends into the upper region of thresholds reported to induce endothelial injury.

View larger version (20K):
[in this window]
[in a new window]
|
Fig. 7. Predicted arterial shear stress distribution in fetal, Control twin, and Shunt twin fractal models. Shear stress amplified with branching order k according to k = 0Rd(3D0)k using steady-state blood flows in the parent arteries to solve for shear stress in the parent 0 = 32µq0/ d03. Here, shear stress reactivity, the increase in shear stress with branching order, is determined primarily by the degree that D0 deviates from D0 = 3 (Murray's law). If the distribution with order is = D0 = 3, then shear stress is expected to be constant, with branching and equal to 0. Otherwise, if = D0 < 3, shear stress is amplified in the course of branching, to a degree determined by how small = D0 becomes. A: gray area demarcates reported ranges of critical values of acute shear stress within large-conduit arteries that induce endothelial dysfunction as a permeability change or endothelial mechanical injury (19, 35, 54). At steady-state flow rates, the predicted shear stress in intra-acinar arteries in the fetus and the Control twin at 8 wk are predicted to be similar. The shear stress in the Shunt twin is predicted to be greater for all arterial diameters. B: comparison of Shunt and Control twin shear stress distributions at corresponding steady-state flow rates, high-flow equivalent (100 ml·min1·kg1) vs. low flow (33 ml·min1·kg1). Increasing or decreasing flow increases/decreases amplitude of shear stress in small vessels, but amplification is always greater with smaller fractal dimensions (Shunt twin > Control twin) for any given flow rate.
|
|
 |
DISCUSSION
|
|---|
The present study is the first to demonstrate the existence of an unusual fetuslike pattern of arterial diameter network organization after birth in an experimental model of pulmonary hypertension with increased pulmonary blood flow. Persistence and/or exacerbation of the fetal pattern after birth suggests a mechanical basis of hemodynamic reactivity that may be linked to the origin of the disease and implicates a shear stress paradigm of structural remodeling in the pathology and progression of pulmonary hypertension in high-flow lesions. Although this study is limited to a morphometric analysis based on a sample of three lung casts only and not a population, our results indicate that a fetal phenotype may be a major factor that influences the state of structural remodeling in intra-acinar vessels when viewed in the context of classic grades of vascular remodeling determined by vessel geometry alone and is involved in two divergent patterns of adaptation in arterial network organization in a model of neonatal congenital heart disease with increased pulmonary blood flow.
Previous work regarding pulmonary vascular abnormalities in patients with pulmonary hypertension has focused primarily on histological appearance of the arterial wall. In 1958, Heath and Edwards (24a) described a classification sequence that ranged from reversible medial hypertrophy to terminal changes such as angiomatoid formation and fibrinoid necrosis. Rabinovitch et al. in 1978 (42a) described a modification that featured alterations in normal remodeling and growth in three progressive stages (AC), characterized by progressive medial hypertrophy, abnormal muscle extension, and reduced artery size and concentration, that have been correlated with perioperative hemodynamics. The histological pattern of vascular remodeling in the Shunt twin at 8 wk of age occurring in this experimental model has been previously characterized as grade B alterations, with medial hypertrophy and abnormal appearance of muscle in the walls of intra-acinar arteries (40). A fundamental limitation of histologically derived indices of vascular reactivity is that it is not clear whether their impact on arterial diameter comprises a mechanical cause of pulmonary hypertension or represents an adaptation in response to it (47). Alternatively, our study was concerned with further characterization of the high-flow secondary hypertension experimental model with specific regard to the influence of elevated flow complicating the structural remodeling of a preexisting fetal phenotype of arterial diameter network organization to a postnatal adult pattern in the Shunt twin, but not in the Control twin. In this study we show that a fetal phenotypic source of hemodynamic reactivity (
< 2,
< 1) is evident in arterial diameters of fetal and Shunt twin intra-acinar arteries. We note that fetal intra-acinar arteries are typically devoid of smooth muscle (2, 41), in which case the fetal phenotype represents a preexisting mechanical source and state of aggravated hemodynamic reactivity contributing to the high-resistance, hypertensive condition normally seen in the late-term fetal pulmonary circulation (50). Consequently, as Fig. 1 indicates, we would expect that the structural remodeling of the fetal phenotype to smaller area ratios would be of greater hemodynamic significance to the progression and timing of hypertension than that expected from the structural remodeling and encroachment of vessel diameter alone fixed at a given area ratio (47). Moreover, because the unusual fetal topological pattern between diameters represents a source of hemodynamic reactivity not considered in the classic vascular remodeling perspective, it is likely that changes in the topological form of hemodynamic reactivity in the fetal-neonatal state may precede the timing and appearance of more complex lesions after birth (43). Indeed, preferential patterns of structural remodeling of bifurcation topology and design in relationship to the appearance of smooth muscle and wall thickness are seen in our data in two ways. First, as seen in the 8-wk Control twin lung cast (Fig. 3), the appearance of a normal pattern of smooth muscle is related to the timing of significant increases in the area ratio. In this case, the migration of
and
toward values of
2 and
1, and to the fractal dimension of D0
2, indicates that the degree of hemodynamic reactivity is significantly diminished relative to that of the fetus. Under these circumstances, the structural and functional impact of the residual fetal configuration in the Control twin on overall hemodynamic reactivity at this point may be comparatively minimal relative to the emerging influence of smooth muscle on vascular reactivity and tone. Second, Fig. 3 also illustrates that the abnormally increased muscle mass in Shunt twin intra-acinar vessels (0.020- to 0.040-mm ID) relative to the Control twin is associated with diameters of parent/daughter vessels with significantly smaller
and
, relative to both the fetal and Control twin lungs. Under these circumstances, the increased muscle mass and the aggravated hemodynamic reactive design and topology indicate a preferential action of smooth muscle vascular remodeling that can exacerbate the fetal phenotype in the course of hypertension disease. Thus the regional pulmonary patterns of
and
we observed in this study emphasize that there is a preexisting mechanical source of hemodynamic reactivity and remodeling in the fetal pulmonary circulation that extends beyond the apparently latent remodeling patterns of vascular wall properties at a vessel level of organization.
The conceptual advantage of the fractal model is that it simplifies the complexity of arterial network organization into a simple static parametric characterization of arterial geometry (Rd), topology (area ratio
), and design (fractal dimension D0), whose impact on hemodynamics can be inferred from an equivalent symmetric bifurcation-branching arterial network (24, 30). The fractal dimension and/or area ratio have structural and functional significance in that their values define the global hemodynamic reactivity of the arterial tree, along with the pattern of transformation of pressure and shear stress with arterial branching (Fig. 1; APPENDIX) (4, 30). Our morphometric results demonstrated differences in the fractal dimension between the fetus (D0 = 1.72) and the Control twin (D0 = 2.02), suggesting a transformation in fractal arterial network organization at birth; furthermore, our results indicate that this transformation did not take place in the Shunt twin (D0 = 1.70). However, we emphasize that our study is limited in the sense that it is based on the results of three lung casts only, which are isolated samples of a term lamb fetus and a postnatal state at 8 wk of birth. Unfortunately, the limited number of lung casts we studied does not identify the exact timing of when the adult bifurcation transport design is achieved after birth. Also, the limited number of lung casts we utilized may not necessarily reflect the developmental dynamics of change in network properties representative of the population of human subjects or other species (39) experiencing the same hemodynamic conditions. Here, genotypic variation and heterogeneity in patterns of vascular wall remodeling may exist within a population of individuals and between species that also influences the distribution and reactivity of the fetal phenotype of arterial branching and its responsiveness to hemodynamic forces for arterial network structural remodeling (53). In addition, at a given developmental time period lung casts are static pictures of network organization that neglect the dynamic information about the influence of vasomotor activity on bifurcation constriction or dilation. Here it should be appreciated that the actions of vasomotion and vascular tone are now more complex beyond their influence on vascular diameter (12, 17), because they also modulate bifurcation design and topology, which can also influence the fractal dimension of design. Thus a spectrum of hemodynamic reactivity in network organization may be present at different times under different hemodynamic circumstances in the perinatal period than those observed in this study. Despite limitations, the pervasive nature of the unusual fetal design phenotype and its divergent directions of structural remodeling seen here after birth indicate that persistence and/or aggravation of the fetal phenotype after birth may be important in developing further understanding of the spectrum of abnormal patterns in hemodynamic reactivity and vascular wall remodeling seen in patients with pulmonary hypertension secondary to elevated blood flow (43).
The fractal dimensions observed in the fetus and the Shunt twin are curious because they are at odds with what is predicted to be universal for arterial network organization in mammalian vascular systems in general (32, 33, 5557). While hemodynamic forces, such as cyclic or mean pressure and shear stress secondary to increased flow, are stimuli known to influence endothelial dysfunction and vascular wall remodeling (10), various theories predict that arterial vascular network organization maintains a universal design, or a target configuration, with arterial branching ratios always
1, with the condition of
< 1 improbable. An assertion made by West-Brown-Enquist (57) is that arterial network organization has a universal design based on a widely known metabolic scaling principle in adult species, which uniquely links the fractal dimensions of arterial network organization to the exponent of the 3/4 law metabolic scaling law [Kleiber's relationship (29)]. Their theory maintains that a fractal design of D0 = 2 and a diameter topology of
= 1 exist for all arterial transport and capillary exchange networks, and that this design remains independent of animal size, ontogenic state, and hemodynamic condition (56, 57). Another classic design, called Murray's law, dictates that arterial network organization has a universal configuration of D0 = 3.0, a condition optimized for constant shear stress and resistance with arterial branching (33). Alternatively, Kurz and Sandau (32) predicted that the bifurcation scaling exponents of arterial network organization remodel during ontogenesis and converge to a target configuration of D0 = 2.7, a design indicative of a vascular system of constant wall tension. Whatever the target, the assumption of a fixed, universal design independent of hemodynamic conditions is relevant and crucial to the study of structural remodeling in pulmonary hypertension. Design invariance with respect to hemodynamic conditions is an operational premise that serves to justify the idea that the complexity of network organization is not a precursor to the generation of hemodynamic forces that are instrumental to the remodeling process of pulmonary hypertension disease. In addition, network design invariance provides a necessary theoretical justification that the complexity of the pulmonary arterial system can be exclusively reduced to a vessel level of organization to elucidate mechanisms controlling diameter restriction in pulmonary hypertension. What the theories of arterial network organization overlook is that the fetal pulmonary circulation does not function as a generalized transport/exchange organ system in utero. In addition, hemodynamic conditions in the fetal pulmonary circulation are also much different from those in the target adult pulmonary circulation: the fetal pulmonary bed is hypertensive in relation to the systemic bed and is subject to a relatively small blood flow (50). At birth the fetal configuration switches in response to an increase in blood flow, via the removal of the placenta and the closure of the ductus arteriosis and foramen ovale, in conjunction with a decrease in capillary resistance with airway expansion. Our morphometric results indicate that the hemodynamic transition at birth is associated with an important, necessary transformation in arterial network organization, such that the fetal configuration functions as a rapidly adapting high-impedance transformer, which is capable of dynamically altering its hemodynamic reactivity in response to shear stress perturbations via changes in blood flow.
Although our data are limited in not allowing us to generalize to the fractal state of arterial network organization in other fetal and neonatal states as a population or conditions of hypertension, we can address the general question of whether fractal network organization is independent of hemodynamic conditions in another way, via the influence of fractal network organization on the scaling of cumulative arterial resistance and slopes of pressure-flow curves relative to body weight. As Eq. 12 indicates, Rcum functionally dictates the slope of the pressure-flow curve (13), with increases in Rcum (Eq. 10) accomplished via smaller fractal dimensions and larger diameter ratios Rd, while a compensating decrease in Rcum can be accomplished via a smaller R0. However, arterial resistance Rcum is predicted to scale with body weight independent of species, state of development, and hemodynamic condition via a fixed fractal dimension according to different patterns: the West-Brown-Enquist law predicts resistance to universally scale according to Rcum
M3/4 with arterial D0 = 2 (57), while Murray's law predicts Rcum
M1 with arterial D0 = 3 (51). Consequently, if vascular network organization is independent of hemodynamic conditions and state of pulmonary development, then the slopes of pulmonary pressure-flow curves sampled from different species, ages, and conditions of pulmonary hypertension would be anticipated to be devoid of outliers and invariant about a given power law. Figure 8 summarizes the scaling of Rcum estimated via the slope of pulmonary pressure-flow curves drawn from the literature (1, 15, 16, 22, 31, 36, 38, 42, 45) and indicates that the resistance scaling law for the pulmonary circulation of the normotensive adult does not adhere to a Rcum
M3/4 as predicted by the West-Brown-Enquist theory, but instead migrates to an M1 power law only in the adult case free of pulmonary hypertension. At the same time, under conditions of pulmonary hypertension in the fetal-neonatal state (8, 46) and the adult state with primary pulmonary hypertension (36), the deviations about the M1 power law indicate that resistance does not adhere to any particular law at all. Under these circumstances, the degrees of freedom possible in a fractal pattern of network organization to influence resistance suggest that the elevated cumulative resistance seen in the fetus and adult with pulmonary hypertension may possibly share a common topology (
< 1) and design (D0 < 2) in arterial network organization. At the same time, if arterial network organization is, in general, not independent of hemodynamic conditions, then there must be a common hemodynamic feedback signal to condition arterial design to target the fractal dimensions of arterial network organization from one resistance state of organization to another. As for targets and their feedback signals, LaBarbera (33) proposed constant shear stress for a target of D0 = 3 while Kurz and Sandau (32) theorized constant wall tension for a target of D0 = 2.7. Regardless of the signal and the target design, intermediate states of Rcum and its associated fractal dimension seen in Fig. 8 can be understood exclusively in terms of a network design that is dynamic, with a capacity to change its organization to act as a variable hemodynamic transformer (3, 4) that changes its reactivity under different hemodynamic conditions.

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 8. A: Rcum power laws. Adult pulmonary cumulative resistance in absence of hypertension adheres to Rcum M1 power law (51), not a universal Rcum M3/4 law (57). Pressure-flow data for Rcum derived from Refs. 1, 15, 16, 22, 31, 36, 38, 42, 45. B: Rcum does not adhere to an invariant power law over developmental/physiological timescales. Rcum·M M0 median trend line (0.08 mmHg·kg·ml1·min1) is value of impedance predicted by universal power law and allometric result of a mean pressure of 15 mmHg divided by allometric scaling of cardiac output with body size. Outliers to the law include a vertical allometric comparison between species at same body weight demonstrating significantly different resistances (39), not a singular impedance scaled to M3/4 design. Also, high-resistance outliers to Rcum M1 occur in fetal hypertension states (8, 46) and in the adult with primary pulmonary hypertension (PPHN) (36).
|
|
One major pattern of structural remodeling observed in this study was the transformation in fractal dimension between the fetal and postnatal states (fetus and Control twin), which was not seen in the Shunt twin. There is likely more than one contribution responsible for the structural remodeling of arterial network organization. As the hemodynamic transition at birth is associated with changes in vascular tone (17), we suspect that vascular tone is involved in remodeling of bifurcation design. Indeed, Griffith et al. (23) demonstrated in isolated bifurcations that vasodilation selectively increases the area ratio in bifurcations from a baseline state of
= 2 to
3, a constant shear stress design, via a nitric oxide (NO) mechanism, while vasoconstriction from a baseline state to a fetal design of
1.25 and
< 1 is induced by endothelin (ET)-1. Another factor expected to influence the structural remodeling of arterial design and topology is the influence of shear stress. In general, the endothelium responds to shear stress by altering its functional phenotype and endothelial gene expression via several thousand genes over different timescales (21). The endothelium incorporates machinery for a rapid transduction-translation response in function and geometry coupled to shear stress (11), which includes potassium channels (9), vasodilation via NO and ET-1 (5, 18), and possibly lung Krüppel-like transcription factors to gradually remodel the state of tone (14). Over hours to days shear stress is capable of signaling a spectrum of remodeling changes within the vessel wall (endothelium, intima, smooth muscle, and adventitia) working as a syncytium that add or detract from vessel wall mass to either encroach or increase diameter (21, 34). Thus, insofar as the unusual configuration of the term fetal state is already preconditioned to hypertension, and small increments in flow lead to large transient avalanches of shear stress, causal reasoning implicates that it is the transduction of shear stress via the endothelium and its translation of endothelial/smooth muscle function to reorganize bifurcation topology, design, and hemodynamic reactivity that is a principal hemodynamic stimulus for remodeling to an adult configuration at birth, rather than the vasodilation of vascular diameter alone (17). We also suspect that shear stress is implicated in the abnormal structural remodeling of arterial diameter geometry, topology, and design in the Shunt twin. In postnatal experimental models of congenital heart disease, where the reduction in resistance has occurred in conjunction with an adult pulmonary arterial configuration, we would have anticipated that the chronically elevated pressure would have increased pulmonary vascular resistance (6, 48). However, while the reactive branching pattern persisted and was aggravated in smaller arterial vessels, the Shunt twin demonstrated a paradoxical decrease in absolute cumulative resistance relative to the Control twin, which is consistent with our previous findings (40, 44). In this example, the decrease in resistance can be attributed to geometric remodeling leading to a significant increase in larger diameter vessels. Figure 5 illustrates the regional influence of arterial diameter scaling in the Shunt twin relative to the fetus and the Control twin: about a pivot point of order 11, orders 08 possess larger diameters and orders 911 approximate Control twin diameters, while orders 1222 are smaller than Control twin diameters. The impact of the increased diameters in the larger arteries is to increase arterial volume, relative to control, a common finding in this model and in patients with chronic high-flow hypertension (44). This geometric feature of structural remodeling may have arisen from the increased flow, rather than the pressure alone, via a negative feedback mechanism of diameter structural remodeling dependent on an endothelial shear stress mechanism, mediated by NO (34). However, the persistent pattern of fetal arterial network organization in the Shunt twin in combination with elevated blood flow suggests a potential initiating physical-mechanical cause and site of action of endothelial injury/dysfunction due to amplified shear stress (25). In Fig. 7 the flow in the Shunt twin may be elevated to a sufficient extent to lead to excessive chronic levels of shear stress known to transiently induce permeability changes and mechanical injury (19, 35, 54). Unfortunately, our study is limited in that it does not identify whether the amplitude of shear stress alone is an effective cause of the persistence of the fetal bifurcation design. Alternatively, other factors operating within the fetal state besides shear stress-induced endothelial dysfunction and/or injury, such as the disruption of endothelial/smooth muscle function and communication (10), may be instrumental effective causes interrupting the fetal-neonatal transformation to a low-impedance transport design in which elevated shear stress is a complication. Under these conditions, amplified shear stress itself may represent an additional complication factor that interferes with the structural remodeling of the fetal bifurcation design to an adult configuration. We submit that the nature and direction of structural remodeling in arterial networks with area ratios <1 under a shear stress paradigm may have implications for understanding the origins and course of vascular reactivity in pulmonary hypertension in the perinatal period.
In summary, our results indicate that the unusual hemodynamic reactivity of the perinatal pulmonary arterial circulation in response to elevated blood flow leading to pulmonary hypertension can be potentially explained via a shear stress paradigm of arterial diameter network remodeling. In an elevated pulmonary blood flow lamb model, the increased resistance and reactivity predicted in pulmonary hypertension has a preexisting mechanical basis evidenced by a fetal diameter branching phenotype that acts as a transformer/amplifier of hemodynamic forces that persists and/or becomes exacerbated after birth. This fetal phenotypic pattern appears before the emergence of abnormal vascular histological patterns of remodeling seen in the course of secondary pulmonary hypertension and serves to maintain and/or exacerbate the hemodynamic reactivity of the fetal phenotype. The levels of shear stress predicted in preacinar and intra-acinar arterial vessels in our model indicate that the critical levels of shear stress leading to endothelial dysfunction/injury may arise mechanically from the interaction of blood flow and the fetal phenotype of arterial network organization itself. While the actual mechanisms responsible for arterial network remodeling and adaptation in the perinatal period are not yet identified, our shear stress paradigm suggests possible remodeling pathways of arterial network diameter topology and design that extend beyond present understanding of a stationary pattern of arterial diameter geometry under a transmural pressure paradigm of vascular remodeling diameter and wall properties. Given the pressure-flow properties of the pulmonary circulation in patients with other forms of pulmonary hypertension, and the theoretical dependence of the area ratio on hemodynamic reactivity, we speculate that diverse forms of pulmonary arterial hypertension may share a common origin and formal cause embedded in replication of a high-impedance fetal pattern of branching complexity.
 |
APPENDIX: ROLE OF ARTERIAL FRACTAL NETWORK ORGANIZATION ON HEMODYNAMIC REACTIVITY
|
|---|
Our present concept of vascular reactivity is based primarily on factors that influence vascular diameter impacting vessel hemodynamic resistance via Poiseuille's law via Eq. 10. However, an additional source of hemodynamic reactivity arises from an influence of arterial network topology and design, which affects the distribution and scaling of hemodynamic forces via the manner in which diameters are interconnected. The purpose of this APPENDIX is to summarize the functional morphometric relationships of topology and design used to characterize hemodynamic reactivity in terms of dimensionless constants, such as the area ratio
, the bifurcation scaling exponent
, the fractal dimension D0, and the resistance/shear stress reactivity constant
with arterial branching.
Bifurcation topology.
The area ratio expresses the cross-sectional area of perfusion with arterial branching that can be approximated in terms of network topology and geometry
 | (A1) |
This relationship states that the change in cross-sectional area of perfusion to blood flow is determined by the product of the relative increase in the number of vessels with arterial branching (via the branching ratio
b = Nk+1/Nk) and the relative decrement in average diameter (
d =
k+1/
k), where the subscripts k and k+1 refer to the order of branching, with k the parent to daughters k+1. A bifurcation is the simplest network topology in which the area ratio can be explicitly expressed in terms of diameter geometry. In general, arterial networks branch asymmetrically: branching typically consists of a parent diameter d0 and daughter diameters d1 and d2 such that d0 > d1
d2. Here bifurcation summarizes the bifurcation's topology via a geometric relationship
 | (A2) |
where the local bifurcation diameter ratio is defined as Rd = d1/d0 and the asymmetry ratio is defined as
= d2/d1. Thus, in the context of the definition in Eq. A1, the diameter ratio and daughter asymmetry of the bifurcation are associated with hemodynamically equivalent average branching ratio and diameter ratio. In the case of a symmetric bifurcation
= 1 the area ratio has a simple topological-geometric interpretation, as
 | (A3) |
where
b = 2 and
d = Rd. Alone, the area ratio alone says nothing about its hemodynamic impact in response to flow.
Bifurcation design.
A structure-function scaling relationship between flow and diameter within a vessel leads to the concept of a hemodynamic design for a given diameter network topology that can be used for hemodynamic analysis and prediction. Roy and Woldenberg (49) argued that the steady-state flow (q) in the vasculature has a general scaling relationship with diameter (d) according to
 | (A4) |
where k is a constant. The value of the scaling constant
has hemodynamic significance, imparting a design or plan for the steady components of flow and pressure that are optimal in some sense (3). In the case of a bifurcation, the conservation of flow dictates that q0 = q1 + q2. With the fraction of flow distributed in daughter vessels according to (p1 = q1/q0, p2 = q2/q0 = 1 p1) the flow-diameter scaling relationship in the bifurcation is
 | (A5) |
In the case of an asymmetric bifurcation the area ratio is related to the diameter (49)
 | (A6) |
where the value for a symmetric bifurcation network is noted by
= 1.
The value of
, although dimensionless, has hemodynamic significance that reflects an underlying law or principle associated with the influence of topology of network organization on the distribution of blood flow (3):
= 4,
> 1 confers constant resistance with branching;
= 3,
> 1 constant shear stress (Murray's law);
= 2.7,
> 1 constant wall tension (Kurz-Sandau law); and
= 2,
= 1 constant velocity (West-Brown-Enquist law). The hemodynamic influence of a given value of
also impacts the scaling of hemodynamic properties associated with flow q, such as the Reynolds number (Re = q
/µd), flow velocity (
= q/
d2), and/or shear stress (
= 32q/
d2). In essence, arterial network organization can influence a given hemodynamic property depending on the value of the design parameter
. In a symmetric bifurcation (Rb = 2 and Rd = 21/
):
 | (A7) |
 | (A8) |
 | (A9) |
These indices ignore the development of flow, the formation of fluid dynamic boundary layers, and flow-diameter asymmetry. However, what these scaling relationships have in common is that there is a critical value of
crit that influences whether network topology either attenuates or amplifies a given hemodynamic property.
Network topology and design.
Horsfield and Woldenberg (28) argued that, akin to flow-diameter scaling in a bifurcation, there is a fractal flow-diameter scaling relationship in arterial trees. The number of branches in successive orders down a fractal tree, such as the pulmonary circulation, increases by a constant proportion (
b = Nk+1/Nk), while the mean flow in branches in successive orders will decrease by a proportion
=
k+1/
k that is inversely proportional to the branching ratio,
b =
1. The scaling relationship between the average decrease in diameter with branching,
d =
k+1/
k in conjunction with the vessel number or inverse fraction of flow distributed, leads to the seed of a fractal power law scaling relationship,
 | (A10) |
If Rb and Rd are constant with branching order k, then the arterial network can be formally considered as self-similar, with the scaling constant
approximately equal to a fractal dimension
D0 (30). Under these circumstances, the topology network organization via the area ratio is related to the fractal network scaling property according to
 | (A11) |
Consequently, network organization and design has an exponential scaling influence with branching order k on the attenuation/amplification properties according to
 | (A12) |
 | (A13) |
 | (A14) |
The idea that network organization influences hemodynamic properties in response to blood flow suggests that the diameter scaling properties can be summarized in terms of a reactivity parameter
that demonstrates its dependence on arterial topology and design with arterial branching. Thus, for shear stress reactivity, the relative increment of shear stress from the first order 0 to the first branching level 1 is given by
 | (A15) |
Figure 1 illustrates the relative influence of
and D0 for
relative to one level of branching, as well as the amplification effects over several orders of branching.
Resistance reactivity.
Under the assumptions that Poiseuille's law applies (Eq. 10), the cumulative resistance of a complex network depends on how resistance changes with arterial branching
 | (A16) |
where
k is the average resistance (Eq. 10) of the vessels of order k with average diameter
k and length
k, and Nk is the number of vessels within order k. Expressing Eq. A16 in terms of the presence of self-similar fractal scaling relationships, a constant diameter ratio
d =
k+1/
k and a constant length ratio
l =
k+1/
k, Rcum becomes
 | (A17) |
In asymmetric trees, the relative values of
d and
b depend on the methods of diameter ordering: Horsfield ordering yields
bH
2, and a smaller
dH(26, 27) with more orders for the same hemodynamically equivalent tree than that ordered by a Strahler method with
bS
2 (30). The branching ratios are related by 1 = (RdH)1 + (RdS)1 = (RdH)1 + (RdH)(1+
) via a symmetry factor
(26, 27). If an additional scaling relationship for vessel length and diameter is taken into account, assuming that arterial network organization has a common fractal length-diameter scaling relationship
k = a
k
2(
2 = 1 and a is an arbitrary scaling constant that depends on the method of vessel ordering) and fractal scaling condition
b =
dD0, then Rcum can be expressed more simply as
 | (A18) |
In effect, the resistance of order k is R(k) = R0
d(3
)k, in which the fractal network organization imparts a contributing reactive factor
R =
d3D0 within each order. Thus, under the assumptions of network organization in this fractal model, the resistance reactivity is equal to the shear stress reactivity
= 
of Eq. A15.
 | (A19) |
Figure 1 illustrates the impact of hemodynamic reactivity with branching on normalized values of Rcum and
for the simplest case of a symmetrically bifurcating network where
b = 2 and
d = 21/D0.
 |
FOOTNOTES
|
|---|
Address for reprint requests and other correspondence: S. H. Bennett, Neonatology Surge I Rm. 1121, Univ. of California, Davis, CA 95616 (e-mail: shbennett{at}ucdavis.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
|
|---|
- Barnard J, Wilson P, Moore T, Thompson W, Taylor A. Effect of nitric oxide and cyclooxygenase products on vascular resistance in dog and rat lungs. J Appl Physiol 74: 29402948, 1993.[Abstract/Free Full Text]
- Belik J, Keeley FW, Baldwin F, Rabinovitch M. Pulmonary hypertension and vascular remodeling in fetal sheep. Am J Physiol Heart Circ Physiol 266: H2303H2309, 1994.[Abstract/Free Full Text]
- Bennett S, Eldridge M, Zaghi D, Zaghi S, Milstein J, Goetzman B. Form and function of fetal and neonatal pulmonary arterial bifurcations. Am J Physiol Heart Circ Physiol 279: H3047H3057, 2000.[Abstract/Free Full Text]
- Bennett SH, Goetzman BW, Milstein JM, Pannu JS. Role of arterial design on pulse wave reflection in a fractal pulmonary network. J Appl Physiol 80: 10331056, 1996.[Abstract/Free Full Text]
- Black SM, Bekker JM, Johengen MJ, Parry AJ, Soifer SJ, Fineman JR. Altered regulation of the ET-1 cascade in lambs with increased pulmonary blood flow and pulmonary hypertension. Pediatr Res 47: 97106, 2000.[Web of Science][Medline]
- Bousamra M, Rossi R, Jacobs E, Parviz M, Busch C, Nelin L, Haworth S, Dawson C. Systemic lobar shunting induces advanced pulmonary vasculopathy. J Thorac Cardiovasc Surg 120: 8898, 2000.[Abstract/Free Full Text]
- Budhiraja R, Tuder R, Hassoun PM. Endothelial dysfunction in pulmonary hypertension. Circulation 109: 159165, 2004.[Free Full Text]
- Cassin S, Dawes G, Mott J, Ross B, Strang L. The vascular resistance of the foetal and newly ventilated lung of the lamb. J Physiol 171: 6179, 1963.[Web of Science]
- Cornfield D, Resnik E, Herron J, Reinhartz O, Fineman J. Pulmonary vascular K+ channel expression and vasoreactivity in a model of congenital heart disease. Am J Physiol Lung Cell Mol Physiol 283: L1210L1219, 2002.[Abstract/Free Full Text]
- Dakshinamurti S. Pathophysiologic mechanisms of persistent pulmonary hypertension of the newborn. Pediatr Pulmonol 39: 492503, 2005.[CrossRef][Web of Science][Medline]
- Davies PF. Flow-mediated endothelial mechanotransduction. Physiol Rev 75: 519560, 1995.[Abstract/Free Full Text]
- Dawson C. Role of pulmonary vasomotion in physiology of the lung. Physiol Rev 64: 544616, 1984.[Free Full Text]
- Dawson CA, Krenz GS, Karau KL, Haworth ST, Hanger CC, Linehan JH. Structure-function relationships in the pulmonary artery tree. J Appl Physiol 86: 569583, 1999.[Abstract/Free Full Text]
- Dekker R, Thienen JV, Rohlena J, Jager SD, Elderkamp Y, Seppen J, deVries C, Blessen E, Berkel TV, Pannekeoek H, Horrevoets A. Endothelial KLF2 links local arterial shear stress levels to the expression of vascular tone-regulating genes. Am J Pathol 167: 609618, 2005.[Abstract/Free Full Text]
- Emery C, Bee D, Barer G. Mechanical properties and reactivity of vessels in isolated perfused lungs of chronically hypoxic rats. Clin Sci (Lond) 61: 569580, 1981.[Medline]
- Engelberg J, Dubois A. Mechanics of pulmonary circulation in isolated rabbit lungs. Am J Physiol 196: 401414, 1959.[Abstract/Free Full Text]
- Fineman JR, Soifer SJ, Heyman MA. Regulation of pulmonary vascular tone in the perinatal period. Annu Rev Physiol 57: 115134, 1995.[CrossRef][Web of Science][Medline]
- Fineman JR, Soifer SJ, Heymann MA. The role of pulmonary vascular endothelium in perinatal pulmonary circulatory regulation. Semin Perinatol 15: 5862, 1991.[Web of Science][Medline]
- Fry D. Acute vascular endothelial changes associated with increased blood velocity gradients. Circ Res 22: 165197, 1968.[Abstract/Free Full Text]
- Galie N, Torbicki A, Barst R, Dartevelle P, Haworth S, Higenbottam T, Olschewski H, Peacock A, Pietra G, Rubin LJ, Simonneau G, Priori SG, Garcia MA, Blanc JJ, Budaj A, Cowie M, Dean V, Deckers J, Burgos EF, Lekakis J, Lindahl B, Mazzotta G, McGregor K, Morais J, Oto A, Smiseth OA, Barbera JA, Gibbs S, Hoeper M, Humbert M, Naeije R, Pepke-Zaba J, Force T. Guidelines on diagnosis and treatment of pulmonary arterial hypertension. The Task Force on Diagnosis and Treatment of Pulmonary Arterial Hypertension of the European Society of Cardiology. Eur Heart J 25: 22432278, 2004.[Free Full Text]
- Garcia-Cardena G, Comander J, Anderson KR, Blackman BR, Gimbrone MA Jr. Biomechanical activation of vascular endothelium as a determinant of its functional phenotype. Proc Natl Acad Sci USA 98: 44784485, 2001.[Abstract/Free Full Text]
- Gregory T, Ellsworth M, Newell J. Hypoxic pulmonary vasoconstriction during steady and pulsatile flow in ferrets. J Appl Physiol 57: 205212, 1984.[Abstract/Free Full Text]
- Griffith T, Edwards D, Randall M. Blood flow and optimal vascular topography: role of the endothelium. Basic Res Cardiol 86: 8996, 1991.[Web of Science][Medline]
- Haworth S, Linehan J, Bronikowski T, Dawson C. A hemodynamic model representation of the dog lung. J Appl Physiol 79: 1526, 1991.
- Heath D and Edwards J. The pathology of hypertensive pulmonary vascular disease. Circulation 18: 533547, 1958.[Web of Science][Medline]
- Hoffman JIE, Rudolph AM, Heyman MA. Pulmonary vascular disease with congenital heart lesions: pathologic features and causes. Circulation 64: 873877, 1981.[Abstract/Free Full Text]
- Horsfield K. Some mathematical properties of branching trees with application to the respiratory system. Bull Math Biol 38: 305315, 1976.[Web of Science][Medline]
- Horsfield K, Cumming G. Morphometry of the bronchial tree in man. J Appl Physiol 24: 373383, 1968.[Free Full Text]
- Horsfield K, Woldenberg MJ. Diameters and cross-sectional areas of branches in the human pulmonary arterial tree. Anat Rec 223: 245251, 1989.[CrossRef][Medline]
- Kleiber M. The Fire of Life: an Introduction to Animal Energetics. New York: Wiley, 1961.
- Krenz GS, Linehan JH, Dawson CA. A fractal continuum model of the pulmonary arterial tree. J Appl Physiol 72: 22252237, 1992.[Abstract/Free Full Text]
- Krishnan A, Linehan JH, Rickaby DA, Dawson CA. Cat lung hemodynamics: comparison of experimental results and model predictions. J Appl Physiol 61: 20232034, 1986.[Abstract/Free Full Text]
- Kurz H, Sandau KS. Modeling of blood vessel development-bifurcation pattern and hemodynamics, optimality and allometry. Comments Theor Biol 4: 261291, 1997.
- LaBarbera M. Principles of design of fluid transport systems in biology. Science 249: 9921000, 1990.[Abstract/Free Full Text]
- Langille B, O'Donnell F. Reductions in arterial diameter produced by chronic decreases in blood flow are endothelium-dependent. Science 231: 405407, 1986.[Abstract/Free Full Text]
- Langille L. Integrity of arterial endothelium following acute exposure to high shear stress. Biorheology 21: 333346, 1984.[Web of Science][Medline]
- Laskey W, Ferrari VA, Palevsky HI, Kussmaul WG. Pulmonary artery hemodynamics in primary pulmonary hypertension. J Am Coll Cardiol 21: 406412, 1993.[Abstract]
- Limpert E, Stahel W, Abbt M. Log-normal distributions across the sciences: keys and clues. Bioscience 51: 341352, 2001.[CrossRef][Web of Science]
- Maggiorini M, Brimioulle S, De Canniere D, Delcroix M, Wauthy P, Naeije R. Pulmonary vascular impedance response to hypoxia in dogs and minipigs: effects of inhaled nitric oxide. J Appl Physiol 79: 11561162, 1995.[Abstract/Free Full Text]
- Maggiorini M, Melot C, Gilbert E, Vermeulen F, Naeije R. Pulmonary vascular resistance in dogs and minipigseffects of hypoxia and inhaled nitric oxide. Respir Physiol 111: 213222, 1998.[CrossRef][Web of Science][Medline]
- Mata-Greenwood E, Meyrick B, Steinhorn RH, Fineman JR, Black SM. Alterations in TGF-
1 expression in lambs with increased pulmonary blood flow and pulmonary hypertension. Am J Physiol Lung Cell Mol Physiol 285: L209L221, 2003.[Abstract/Free Full Text] - Murphy JD, Rabinovitch M, Goldstein JD, Reid LM. The structural basis of persistent pulmonary hypertension of the newborn infant. J Pediatr 98: 962967, 1981.[CrossRef][Web of Science][Medline]
- Peake M, Harabin A, Brennan N, Sylvester J. Steady-state vascular responses to graded hypoxia in isolated lungs of five species. J Appl Physiol 51: 12141219, 1981.[Abstract/Free Full Text]
- Rabinovitch M, Haworth S, Castenada A, Nadas A, and Reid L. Lung biopsy in congenital heart disease: a morphometric approach to vascular disease. Circulation 58: 11071122, 1978.[Abstract/Free Full Text]
- Rabinovitch M, Keane J, Norwood W, Castenada A, Reid L. Vascular structure in lung biopsy tissue correlated with pulmonary hemodynamic findings after repair of congenital heart defects. Circulation 69: 655667, 1984.[Abstract/Free Full Text]
- Reddy VM, Meyrick B, Wong J, Khoor A, Liddicoat JR, Hanley FL, Fineman JR. In utero placement of aortopulmonary shuntsa model of postnatal pulmonary hypertension with increased pulmonary blood flow in lambs. Circulation 92: 606613, 1995.[Abstract/Free Full Text]
- Reeves J, Taylor A. Pulmonary hemodynamics and fluid exchange in the lungs during exercise. In: Handbook of Physiology. Exercise Regulation and Integration of Multiple Systems. Bethesda, MD: Am. Physiol. Soc., 1996, sect. 12, chapt. 13, p. 585613.
- Reid DL, Thornberg KL. Pulmonary pressure-flow relationships in the fetal lamb during in utero ventilation. J Appl Physiol 69: 16301636, 1990.[Abstract/Free Full Text]
- Reid LM. Vascular remodeling. In: The Pulmonary Circulation, edited by Fishman A. Philadelphia: Univ. of Pennsylvania Press, 1990, p. 259282.
- Rendas A, Lennox S, Reid L. Aortopulmonary shunts in growing pigs. J Thorac Cardiovasc Surg 77: 109118, 1979.[Abstract]
- Roy AG, Woldenberg MJ. A generalization of the optimal models of arterial branching. Bull Math Biol 44: 349360, 1982.[Web of Science][Medline]
- Rudolph AM. The changes in the circulation after birth: their importance in congenital heart disease. Circulation 41: 342359, 1970.
- Sherman TF. On connecting large vessels to small: the meaning of Murray's law. J Gen Physiol 78: 431453, 1981.[Abstract/Free Full Text]
- Singhal S, Henderson R, Horsfield K, Harding K, Cumming G. Morphometry of the human pulmonary arterial tree. Circ Res 33: 190197, 1973.[Abstract/Free Full Text]
- Stenmark K, Mecham R. Cellular and molecular mechanisms of pulmonary vascular remodeling. Annu Rev Physiol 59: 89144, 1997.[CrossRef][Web of Science][Medline]
- Vasinav S, Patel D, Atabek H, Deshpande M, Plowman F, Vossoughi J. Determination of the local erosion stress of the canine endothelium using a jet impingement method. J Biomech Eng 105: 7783, 1983.[Web of Science][Medline]
- Weibel E. Design of biological organisms and fractal geometry. In: Fractals in Biology and Medicine, edited by Nonnenmacher T, Losa G, and Weibel E. Basel: Birkhauser, 1994, p. 6885.
- West G, Brown J, Enquist B. A general model for ontogenic growth. Nature 413: 628631, 2001.[CrossRef]
- West G, Brown J, Enquist B. A general model for the origin of allometric scaling laws in biology. Science 276: 122126, 1997.[Abstract/Free Full Text]
- Zar H. Biostatistical Analysis. Englewood Cliffs, NJ: Prentice-Hall, 1984.
Copyright © 2007 by the American Physiological Society.