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 e