| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Integrative Physiology |
From the Department of Cell Biology and Anatomy (J.T.B., T.C.M., D.S., D.T., R.R.M.) and the Department of Pediatrics (Cardiology) (J.T.B., T.C.M.), Cardiovascular Developmental Biology Center, Childrens Research Institute, Medical University of South Carolina, Charleston; and the Institute of Animal Physiology and Genetics (D.S.), Videnska, and the Institute of Anatomy, 1st Faculty of Medicine (D.S.), Charles University, Prague, Czech Republic.
Correspondence to Jonathan T. Butcher, Department of Biomedical Engineering, 270 Olin Hall, Cornell University, Ithaca NY. Email jtb47{at}cornell.edu
| Abstract |
|---|
|
|
|---|
Key Words: chick development modeling ultrasound flow aspiration
| Introduction |
|---|
|
|
|---|
Various mutant models demonstrate that genetic defects compromising valve structural maturation result in severe regurgitation, dilatation, and lethality, suggesting that the appropriate maturation of valve matrix is essential for continued function.6,7 These models add weight to the argument that structural properties of AV cushion tissue may be critical determinants of AV cushion function, and that abnormal AV cushion tissue properties (other than mass, which has been well documented) may result in physiological abnormalities in blood flow and cardiac mechanical function before gross morphological disturbances are seen. A link between myocardial elastic and viscous material properties and structural development exists,8,9 and changes in these properties from altered hemodynamic loading argue that mechanical forces do shape embryonic cardiac structural development.1,911 The relationship between mechanics and structure in adult valves1012 has been studied, but not in embryonic cushions. Therefore, connections between cushion structure, hemodynamics, and biomechanics may form a basis for quantifying normal and abnormal embryonic valve function. In this paper we have combined state of the art ultrasonography, optical mapping, and a custom mechanical testing system to quantify and relate these parameters. We found a remarkable fine tuning of the AV cushion biomechanics that coincides with a transition in cardiac pumping mechanisms.
| Materials and Methods |
|---|
|
|
|---|
Mesomechanical Testing
Pipette aspiration systems have been developed and used for measuring mechanical properties of cells and tissues that are too small to be tested by conventional techniques.17 We developed a similar system to measure material properties of embryonic cushions (supplemental Figure III). The superior and inferior AV cushions at different stages of development were isolated from avian hearts at HH17, HH21, and H25, representing early, mid, and late cushion formation stages with myocardial wall intact. Cushions were placed in an isotonic bath supplemented with BSA and other amino acids (RB1 medium, kind gift of S. Kubalak Medical University of South Carolina, Charleston) and positioned to the tip of smoothed tip glass micropipette, and adhered at the central portion of the endocardial surface by a small vacuum pressure not capable of distending the tissue (P
0.1 Pa). Aspirated tissue length was then measured simultaneously with applied vacuum pressure. Increased pressure resulted in incrementally less aspirated length, suggesting a nonlinear material response. Previous computational and experimental studies by Ohashi et al and Aoki et al investigated the ability of pipette aspiration to measure nonlinear finite elasticity of soft tissues.17,18 Geometric influences were negligible if the tissue sample was at least 5 times the radius of the pipette in diameter and 4 times the radius of the pipette in thickness, which was the case in all of our tissue samples (supplemental Figure III). They then used strain energy based pseudoelasticity theory to model changes in local tissue mechanics,18,19 and found that the principal tissue stress was equivalent to the applied pressure. To apply this theory to embryonic cushions, it was therefore assumed that the cushion material response was homogeneous, isotropic, and nonlinear hyperelastic. Billiar and Sacks postulated a pseudoelasticity-based constitutive model for adult valve segments that incorporated an exponential strain energy formulation.20 As a prerequisite for pseudoelastic theory, cushions were preconditioned with
20 loading cycles at low pressure (0<P<1.0 Pa) before quasi-static loading to remove previous strain history.21 4 to 10 cushions were tested per anatomical position and stage. The nonlinear cushion loading curve (pressure versus stretch ratio) was then modeled using a similar theory as previously described21 and curve fit by Newton-Gaussian iteration. Material coefficients and effective modulus were compared between cushion location and developmental stage using ANOVA with P<0.05 considered significant, and the data were also assessed for curve fitting by Normalized Standard Estimation of the Error (NSE).22 Additional details are provided in the online supplement.
Structural Constituent Modeling
Avian superior atrioventricular cushions were isolated as previously described from HH21 hearts and enzymatically digested by collagenase 2 (Case; Worthington, 300 U/mL) to remove collagen, hyaluronidase (HAD; Sigma, 100 U/mL) to remove glycosaminoglycans (GAGs), or cytochalasin D (CD; Sigma, 1 µmol/L) to remove cellular traction forces by inhibiting actin polymerization. Each treatment was incubated at 37°C 5% CO2 for 6 hours under gentle rocking, with untreated cushions serving as controls. 5 to 10 cushions were used per treatment. Mechanical testing was conducted as before, and the resulting modeling curves compared as before. The collagen, cell, and glycosaminoglycan (GAG) specific loading curves were derived from the individual treatment results by appropriating strain energy. Mixture theory was then used to predict material responses based on tissue composition. This assumes the total tissue material response is equivalent to the sum of the individual constituent contributions multiplied by their volume fractions. Serial sections through the AV of 4% paraformaldehyde fixed, paraffin embedded chick embryos from HH1725 were stained with Movat pentachrome to identify volume fractions of cells, collagen, and glycosaminoglycans using color thresholding. These volume fractions were combined with the derived component material curves using mixture theory to model whole cushion responses, which were compared with the actual loading curves using NSE. Additional details are available in the online supplement.
| Results |
|---|
|
|
|---|
|
|
|
HH21 endocardial and myocardial AV contours (Figure 2, middle panels and supplemental Movie II) showed persistent thickness throughout the cardiac cycle, ranging from 140 to 260 µm, rather than longitudinal treadmilling. The cushions rocked caudally such that their coaptation occurred at a point (or line in 3D) that translated along the AV canal. Also unlike HH17, changes in tissue thickness were approximately 0.20 cycle (72 degrees) out of phase with the contraction of the myocardium. Peak cushion thickness occurred within the range of peak AV myocardial contraction (again 30%), but minimum thickness occurred while the myocardium is rapidly contracting. Blood flow at this stage was also markedly different. Though still biphasic, the initial phase had much slower blood velocity than the second, which now carried most of the flow energy. AV canal flow pathways were still mostly open during the cardiac cycle (Figure 3).
At HH25, as can be seen in Figure 2, bottom panels (and supplemental Movie III), there was no longer any longitudinal tissue motion through the AV canal and the cushions coapted along their entire length simultaneously. The cushions came together and parted rapidly, and the principle flow was a single atrially propelled jet of blood. Cushion thickness varied in phase with the myocardial contraction, becoming thinnest (260 µm) at peak myocardial relaxation and thickest (450 µm) at peak contraction (now only 16%). Average peak blood velocity during this phase was 17.1 cm/s (Figure 3). Our reported AV Doppler velocity profiles compare well with previously published studies.23,24
Comparison of tissue motions and blood flow from each stage suggested that the flow regulation mechanics of the AV canal were tied directly to the material properties of the cushions (Figure 4). A significant inverse linear correlation is apparent when comparing the tissue wave speed and peak blood velocity between embryos over stages 17 to 25 (R=0.899). The nearly acellular gelatinous AV cushions of HH17 heart undulated as a wave with myocardial contraction. The tissue wave speed was approximately 8 mm/s, slower than the peak blood flow. Myocardial tissue wave velocity was diminished at HH21 (
6 mm/s), whereas peak blood flow velocity was increased (
11 cm/s). There was persistent though changing thickness of the cushions throughout the cardiac cycle during compression and stretching. At HH25, there was almost no tissue wave (
2 mm/s) yet the peak blood velocity was much higher (
20 cm/s). The cushions at this stage change less in thickness over the cardiac cycle than the other stages, suggesting that they were more rigid. These findings highlight changes of AV canal cushion mechanics as being critical to the normal function of the AV canal in regulating unidirectional blood flow.
|
Alterations in Atrioventricular Conduction Patterns
The propagation of the myocardial depolarization signal changes significantly between HH17 and HH25. As shown in Figure 5, depolarization at HH17 was activated at the superior aspect of the primitive atrial segment (*) and progressed in a pseudo-linear manner as indicated by the isochrone lines through the ventricular and outflow segments, compatible with peristaltic/suction pumping. At HH21 the conduction pattern began to transition to a binodal configuration. Depolarization progressed radially through the atrium to the atrioventricular canal, followed by earliest ventricular activation of a portion of the superior ventricle near the inner curvature. By HH25, depolarization of the atria was followed by ventricular activation that was earliest at a zone clearly distant from the AV canal and which then spread independently through the right and left ventricular myocardial tissues (*), supporting function as independent piston-like pumps. To determine the velocity of conduction propagation through the atrioventricular canal at these different stages, the length of the canal was measured using the B-Mode ultrasound images, this length being between the initial and final cushion extremes. This length was then divided by the conduction delay (approximately 49, 40, and 46 ms, respectively), which resulted in HH17, HH21, HH25 conduction velocities being 20.1, 18.4, and 10.3 mm/s,. respectively. These results show that AV conduction velocity is slowing concomitant with the mesenchymal growth of the cushions, as has been previously shown.25
|
Nonlinear Pseudoelastic Material Modeling of Developing AV Valve Cushions
To quantify changes in AV cushion material properties during these stages, a mesomechanical test system was developed to apply tensile tests to these small tissues while limiting gripping artifacts. The stress-strain loading curves of isolated valve cushions all showed a monotonically increasing nonlinear mechanical response. As evidenced by low NSE values, our data were modeled well by pseudoelasticity theory, but there were distinct differences between developmental stage (Figure 6 and supplemental Table I). Generally, both AV cushions were extremely pliable at HH17 (EEff
0.15 Pa) but became successively more rigid at HH21 (EEff
0.85 Pa) and HH25 (EEff
3.6 Pa). Effective moduli were statistically significant between stages (P<0.05), but not between cushions. Statistical differences between linear and nonlinear coefficients (supplemental Table I) suggested that the mechanical response of HH17 inferior and superior AV cushions were different, but similar at HH21 and HH25.
|
Composite Modeling Predicts Cushion Material Properties
Specific enzymatic digestion treatments were applied to HH21 AV cushions to ascertain the contribution of collagen, GAGs, and cells to the material properties. As shown in Figure 7A and 7B, hyaluronidase digestion resulted in a very nonlinear rigid cushion, whereas collagenase digestion resulted in an extremely fragile, more linear elastic tissue. Cytochalasin D treatment of cushions resulted in minimal differences compared with control tissues, and only at large stretch ratios (also supplemental Table I) and without significant differences in effective modulus, suggesting that cell traction forces are not significant contributors to the material properties of the AV cushions at these stages of development.
|
Histological staining of AV valve cushions at HH17 through HH25 shows dramatic changes in cushion morphology and composition (supplemental Figure IV). It is important to note that from a biomechanical perspective, the structural composition of the material is the driving factor in the material properties of the tissue rather than the size of the tissue or absolute amounts of the constituents. HH17 AV cushions are comprised mostly of GAGs (hyaluronan), but some invading cells are present. Between HH21 and HH25, AV cushions increase in collagen and cell proportion, and image analysis shows that collagen and cell content approach 30% each of the total cushion volume (Figure 8A and 8B).
|
Interestingly, the derived individual component loading curves (Figure 7C) were similar in trend to the stage specific cushion loading curves, suggesting that cushion structure may regulate tissue biomechanics. To this end, composite modeling of the cushions at HH17 and HH25 was achieved by combining the component strain energies for GAGs, collagen, and cell traction fractions determined from HH21 cushions with the HH17 and HH25 histologically determined volume fractions. Because CD treatments resulted in little differences in material response and cushion mesenchymal cells were surrounded by collagen fibers at HH21 and HH25, we postulated that the collagen volume fractions at these stages were better represented mechanically by adding to it the cell volume fraction. The component fractions are shown in supplemental Table II, and the resulting curves are shown in Figure 8C and 8D. The mixture model generally predicted the measured material response to tensile stress at HH17 and HH25 for both superior and inferior cushions. The low NSE values for these curves indicate that these curves are still relatively accurate given the number of data points and small number of empirical coefficients. These results show that good prediction of the mechanical properties of valvular cushions were possible by composite modeling of the component volume fractions.
| Discussion |
|---|
|
|
|---|
A seminal study by Foruohar et al determined that the early zebrafish heart functions as a suction pump through the several criteria, most notably that (1) maximum cushion thickness occurs at peak myocardial contraction, (2) blood velocity exceeds that of tissue velocity, and (3) wave propagation is initiated by a single myocardial source.5 Our HH17 data are consistent with these criteria, suggesting that this stage chick heart may also behave like a suction pump, in contrast to the peristaltic mechanism previously suggested. A thermodynamic characterization of piston pumping is a volume changedriven propulsion of fluid, with negligible contribution to pumping by the orifice except to throttle the outlet flow. Our data at HH25 is again consistent with this notion, and therefore appears that the HH21 AV may be a transitional stage between 2 pumping styles.
Our data show that material properties of embryonic prevalvular tissues change in concert with changes in their observed mechanics. Using the enzymatic digestions and Cytochalasin D treatments, the individual contributions of collagen, GAGs, and cell traction forces were approximated. Collagen contributed to the majority of the cushion tissue nonlinear rigidity, as has been shown for numerous soft tissues. GAG (mostly hyaluronan in these tissues) material properties were found to be extremely weak, extensible, and mostly linear. Previous studies measuring the material properties of cartilage, which also contains glycosaminoglycans, showed that hydration and swelling also contribute to biomechanical response.26 GAGs are complex, highly coiled chains that exhibit entropic elasticity: they resist uncoiling from stretching to maintain disorder. The mechanical contribution of GAGs in valve cushions is likely attributable to a combination of entropic elasticity and altered extracellular matrix hydration properties.
We used mixture theory to combine the stress-strain behavior attributable to individual components and found that material properties could be reasonably predicted across developmental stage and anatomical position. Similar formulations were applied to predict material behavior in the aorta as a function of collagen, elastin, and cell contractions, and could predict functional consequences of alterations in component configurations.22 Deviations in the model at large strains of HH25 AV cushions may be caused by unmodeled interaction effects between constituents or from additional cell/matrix components not considered. This is most clearly demonstrated by the apparent negative strain energy induced by HAD treatment in Figure 7C, resulting in a more rigid tissue but thermodynamically impossible. We believe that the HAD effect is the result of the removal of an important but as yet unknown interaction effect between GAGs and collagen. Indeed, it is likely that all of the assumptions posed in this model will need to be revisited in more complex experimental and computational models,20 but nevertheless these results demonstrate the utility of this simplified model for predicting material response to changes in matrix.
Hemodynamics may be an important factor in cushion tissue strengthening required for embryonic development to progress. Temporal increases in blood pressure and blood flow velocity must be supported by changes in cushion tissue structure to ensure proper structural integrity. The increasing cushion tissue elastic rigidity with developmental stage reported here may help to maintain morphological integrity under increasing mechanical stress, as has been postulated in adult valves.27 The nonlinear nature of the stress-strain curve suggests the tissue is pliable under lighter loads but has an innate resistance to deformation by increased loads. This suggests that cushions can preserve their structure under a wide range of hemodynamic forces, and in that way inhibit regurgitation under temporary variations from ideal hemodynamic conditions.
Appropriate transitions in AV cushion properties are likely to include signaling through the endocardial cells that line the lumenal surface of the cushions. It is already known that these cells are unique in their ability to undergo EMT, but they may also play a role in regulating post-EMT cushion morphogenesis through integration of hemodynamic stimuli and signaling of underlying mesenchyme. Vascular endothelial cells are the primary mechanosensitive agents in normal vascular function. Germane to this discussion, it was recently shown that adult valvular endothelial cells possess mechanosensitive functions not shared by vascular endothelial cells,28 and can stimulate valvular interstitial cells in coculture in response to shear stress.29 This raises the question whether similar mechanosensitivity mechanisms are factors in normal and abnormal cushion and valvular development. Several studies have shown that obstructing or altering blood flow patterns results in defective heart development.1,30,31 One recent study demonstrated changes in embryonic endocardial expression of mechanosensitive genes in vivo by altered hemodynamic loading.32 These and the results of this study point to endocardial cells as excellent candidates for the mechanosensitive signaling in normal morphological and physiologic cushion morphogenesis.
How endothelial signaling corresponds with myocardial contractile function is currently unknown, but many of the signals implicated in EMT, such as bone morphogenetic proteins, transforming growth factor beta, and vascular endothelial growth factor, require coordination between the myocardium and endocardium.7,3335 The role of hemodynamics in EMT is still controversial. One study suggests that disruption of myocardial contraction (and likely flow) is sufficient to inhibit EMT.36 The fact that EMT can be initiated in vitro through the expression or inhibition of a variety of signals (reviewed by Person et al37) without hemodynamics does not rule out a necessary role for hemodynamic processes in vivo.
Our article therefore adds to the conclusion that biomechanics play an important signaling role in valvulogenesis. There are many examples where genetic deficiencies result in altered biomechanical properties, which over time lead to early valvular degeneration, as in some palliated congenital heart defects38 as well as functional defects like bicuspid aortic valve, myxomatous valves, and mitral valves in patients with Marfan syndrome. It may be that biomechanical abnormalities arising from a number of genetic causes may contribute to the development of abnormal morphological phenotypes such as common AV canal as well as normal structures that subsequently become dysfunctional prematurely. New biomechanical tools and analytical techniques such as those presented in this study will help to characterize more subtle phenotypes with greater predictive power and, ultimately, clinical potential.
| Acknowledgments |
|---|
This work was funded in part by an NIH training grant in pediatric cardiology (HL007710 to J.T.B.), RR16434 (to R.R.M., T.C.M.), HL/HD67135 (to T.C.M.), VZ 206100-3 from the Ministry of Education of the Czech Republic, IRP IAPG AVOZ50450515 from the Academy of Sciences of the Czech Republic, and Purkinje Fellowship from the Academy of Sciences of the Czech Republic (to D.S.). J.T.B., T.C.M., D.T., and R.R.M. conducted their work in a facility constructed with support from the National Institutes of Health, Grant Number C06 RR018823 from the Extramural Research Facilities Program of the National Center for Research Resources.
Disclosures
None.
| Footnotes |
|---|
| References |
|---|
|
|
|---|
Related Article:
This article has been cited by other articles:
![]() |
M. Pho, W. Lee, D. R. Watt, C. Laschinger, C. A. Simmons, and C. A. McCulloch Cofilin is a marker of myofibroblast differentiation in cells from porcine aortic cardiac valves Am J Physiol Heart Circ Physiol, April 1, 2008; 294(4): H1767 - H1778. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. C. W. Groenendijk, K. Van der Heiden, B. P. Hierck, and R. E. Poelmann The Role of Shear Stress on ET-1, KLF2, and NOS-3 Expression in the Developing Cardiovascular System of Chicken Embryos in a Venous Ligation Model Physiology, December 1, 2007; 22(6): 380 - 389. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. B. Keller New Insights Into the Developmental Biomechanics of the Atrioventricular Valves Circ. Res., May 25, 2007; 100(10): 1399 - 1401. [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Circulation Research Home | Subscriptions | Archives | Feedback | Authors | Help | AHA Journals Home | Search Copyright © 2007 American Heart Association, Inc. All rights reserved. Unauthorized use prohibited. |