Abstract
Recent studies showing significant changes in valvular matrix composition with age offer design criteria for age-specific tissue-engineered heart valves. However, knowledge regarding aging-related changes in valvular material properties is limited. Therefore, 6-week, 6-month, and 6-year-old porcine aortic valves (AV) and mitral valves (MV) were subjected to uniaxial tensile testing. In addition to standard material parameters, the radius of transition curvature (RTC) was measured to assess the acuteness of the transition region of the tension–strain curve. Radially, the MV had greater stiffness and a smaller RTC compared with the AV. Circumferentially, the center of the MV anterior leaflet (MVAC) had the highest stiffness (MVAC > AV > MV free edge [MVF]), greater stress relaxation (MVAC > MVF/AV), lowest extensibility (MVAC < AV < MVF), and smaller RTC compared with MVF (AV < MVAC < MVF). AV and MV radial strips had a larger RTC compared with circumferential strips. Aging elevated stiffness for MV and AV radial and circumferential strips, elevated stress relaxation in AV and MVF circumferential strips, and increased RTC for MV radial and MVF circumferential strips. In conclusion, there are significant age-related differences in the material properties of heart valves, which parallel differences in tissue composition and structure, likely impact valve function, and highlight the need for age-specific design goals for tissue-engineered heart valves.
Introduction
Valves are considered to be quasi-viscoelastic,13,14 highly anisotropic tissues15–20 that demonstrate stress relaxation14,17,20–24 but low hysteresis13,14 and no creep.14,23 Tensile testing of valves results in a characteristic load-elongation curve with a low slope pre-transition region in which the collagen fibers are crimped and the elastic fibers are bearing load, followed by a transition region during which loading is transferred to the collagen fibers, and then a post-transition region whose high slope reflects load bearing by the collagen fibers.21,22,25 Valve anisotropy, in which valves are stiffer and less extensible in the circumferential direction compared with the radial direction,15–20 is attributed to collagen fiber alignment in the circumferential direction.15,26 Stiffness also varies regionally within the mitral valve (MV) (i.e., between the anterior center [MVAC] and the free edge [MVF]), 15 as does ECM composition in these regions. 27
Therefore, the primary objective of this study was to assess age-related differences in the material properties of porcine MV and AV and to consider these differences in the context of valve microstructure and ECM composition. In addition to standard material parameters, we developed a novel parameter termed “radius of transition curvature” (RTC) to quantify the acuteness of the transition from the pre-transition to post-transition region of the tension–strain curve.
Materials and Methods
Tissue sample procurement
Based on comparable anatomy between porcine and human valves,28–30 valves from 6-week, 6-month, and 6-year-old pigs were used (corresponding to child, young adult, and older adult ages in humans 11 ). Porcine hearts were obtained from abattoirs (from Fisher Ham and Meat [Spring, TX] for 6-week-old and 6-month-old pigs; from Animal Technologies [Tyler, TX] for 6-year-old pigs) and measured to normalize leaflet thickness to heart size (maximum circumference and longitudinal length from the origin of the pulmonary artery to heart apex). MVs and AVs were dissected from the hearts, chordae were removed from the MV, and circumferential and radial strips of 5 mm diameter were cut from both MV and AV (Table 1). Given known differences between the material properties of the MVAC (defined as the central region without chordal attachments) and MVF, 15 circumferential strips were taken from both of these regions (Fig. 1A). MV radial strips spanned the entire leaflet height from annulus to free edge. Tissues were stored in phosphate-buffered saline (PBS) at 4°C until testing and tissues were tested within 4 days of harvesting (based on studies showing material properties of valves do not change in that time period31,32).

(
MV, mitral valve; AV, aortic valve; NC, noncoronary; MVAC, MV between the anterior center; and MVF, MV between the free edge; LC, left coronary; RC, right coronary.
Mechanical testing
The thickness of the tissue was measured at four locations using a displacement gauge (Mitutoyo, Kawasaki, Japan). Tissue samples were glued to balsa wood with cyanoacrylate and clamped within an EnduraTec ELF 3200 (Bose, Minnetonka, MN) such that the edge of the balsa wood aligned with the end of the grips (Fig. 1B). The balsa wood improved traction of the tissue within the grips, thereby preventing slippage during testing. Mechanical testing was performed in a 37°C PBS bath as described previously. 33 Briefly, mechanical testing began with the tissue in an unloaded state and consisted of 25 cycles of preconditioning at 1 Hz followed by one load-elongation cycle (0.5 Hz). The tissue was then allowed to rest for 15 s before another three cycles of preconditioning and a stress relaxation test (duration 100 s).
Data analysis
Gauge length, defined as the length at which the tissue first starts bearing load, was calculated by fitting a cubic function to the pre-transition region of the load-elongation curve and then determining the tissue length corresponding to the local minimum, similar to the method used by Carew and Vesely.
34
Recorded displacement values were then divided by gauge length to calculate strains and load converted to tension by dividing by tissue width. The slope of the post-transition portion of the tension–strain curve (post-transition “stiffness”) was calculated by first fitting a straight line to the maximum load and the preceding four points. Additional preceding points were added and the slope of the linear least-squares fit line was recalculated until its value changed by more than 1%. The pre-transition stiffness was determined as the slope of the linear least-squares fit between zero strain to one half of “end-transition” strain (where the end-transition point was the closest point to the origin that was used in the post-transition slope calculation). Extensibility was defined as the intersection of the linear least-squares fit of the post-transition region with the x-axis (strain). The stress relaxation data were fit to the following two-phase decay equation:
The Kfast and Kslow values were calculated iteratively using Graphpad Prism (GraphPad Software, La Jolla, CA). The percentage relaxation (%SR) was calculated as the difference between the initial load and the load remaining at 100 s as a percentage of the initial load. The y0 value refers to the plateau stress in the tissue (where t is approaching 100 s). A0 and A1 are best-fit values from a nonlinear regression referring to the span of the fast and slow portions of the curve. The RTC was calculated by first dividing the y-axis of the tension–strain curve by 100 N/m to make both axes dimensionless and then rotating the tension–strain curve clockwise so that pre-transition and post-transition lines were at equal angles from the vertical (which converted the vertex of the transition region to a global minimum). A hyperbola was then fit to the points in the transition region (from 0 to 1.5 × end-transition strain) and the RTC was calculated as the inverse of the second derivative of this global minimum.
Histology and immunohistochemistry
To evaluate the valvular material behavior in context of the ECM composition, representative tissue strips equivalent to those used for mechanical testing were fixed in 10% formalin overnight, dehydrated, paraffin embedded, and sectioned for histological examination. Movat pentachrome was used to visualize the different ECM components and leaflet layers. Immunohistochemistry was also performed as previously described11,35 to identify and localize the proteoglycans (PGs) decorin, biglycan, and versican (decorin [LF122] and biglycan [LF104] from Dr. Larry Fisher, NIH, Bethesda, MD; 36 versican [2B1] from Associates of Cape Cod, East Falmouth, MA), the glycosaminoglycan (GAG) hyaluronan (HA, using HA-binding protein; Associates of Cape Cod), and collagen type III (Col III, which has been shown to change with age, 11 LF69; from Dr. Larry Fisher 37 ). Based on studies showing that valve cell contractility can affect valvular material behavior, 38 staining was also performed for the myofibroblast marker nonmuscle myosin heavy chain-IIB (NMM; Covance, Berkeley, CA) as previously described. 39 Briefly, immunohistochemical staining consisted of antigen retrieval (citrate buffer; Biocare Medical, Concord, CA), except for PGs in which chondroitinase ABC was used 35 (Associates of Cape Cod), followed by quenching of endogenous peroxidases. Nonspecific staining was blocked using goat serum buffer, followed by application of primary antibody. Rinses in PBS were performed before application of the appropriate species-specific secondary antibody. Positive staining was visualized using Vectastain Elite ABC and diaminobenzidine kits (Vector Laboratories, Burlingame, CA), followed by hematoxylin counterstaining. Radially oriented strips of the MV and AV, which were previously stained in this same manner as described in a previous publication,10,11 were also used to provide context.
Statistical analysis
Multifactorial analysis of variance (ANOVA) was performed using SigmaStat (SPSS, Chicago, IL). When the data were normally distributed, an ANOVA test was used. When the data set was not normally distributed, a rank transform was performed before the ANOVA. In both cases the level of significance was set at 0.05.
Results
Thickness
All test strips cut from the MV were thicker than in the AV (p < 0.001), whether in the circumferential direction (Fig. 2) or radial direction. For all valve strips, thickness increased with age (each p < 0.001), as reported previously for human valves. 40 The normalized thickness of circumferential and radial strips (whether normalized to the circumferential or longitudinal dimension of the heart) of 6-month olds was greater than for other age groups (p ≤ 0.035 for each valve or leaflet region; thickness normalized to circumference).

(
Post-transition and pre-transition stiffnesses
In the post-transition region of the tension–strain curves, MV and AV radial stiffness was less than circumferential (p < 0.001), consistent with previously published studies.15–20 The radial stiffness of MV was greater than for AV (p < 0.001; Fig. 3). Among circumferential strips, the stiffness of MVAC was greatest, followed by AV, and then MVF (p = 0.005). With increasing age, MV and AV radial stiffness increased (each p ≤ 0.003), as did MV and AV circumferential stiffness (each p ≤ 0.010).

Slope of post-transition region of tension–strain curves for (
In the pre-transition region of the curves, there were no significant differences between circumferential and radial stiffnesses in either the AV or MV. Amongst MV and AV circumferential samples, MVF pre-transition stiffness was greater than MVAC, followed by AV (Fig. 4, p < 0.001). With increasing age, the pre-transition stiffness of MV samples (both circumferential and radial) and AV radial samples increased (each p < 0.001).

Slope of the pre-transition region of tension–strain curves for (
Stress relaxation
In the AV, radial sections showed a higher Kfast but a lower Kslow than circumferential strips (Table 2, each p < 0.001). In the MV, however, the Kfast and Kslow of MVAC circumferential strips were higher than for MV radial (each p < 0.005). Among radially oriented sections, the AV Kfast and Kslow were both greater than for the MV (each p ≤ 0.002). Among circumferential sections, Kfast of MVAC was greater than AV, followed by MVF (p < 0.001), but Kslow was greatest in AV (AV > MVAC > MVF, p < 0.001). With increasing age, AV circumferential Kfast was reduced (p = 0.006), as was MVF circumferential Kslow (p = 0.046). With respect to the relative magnitudes of relaxation, %SR was greater for AV circumferential than for AV radial (Fig. 5, p < 0.001) and greater for MVAC circumferential, followed by MVF circumferential, then MV radial (p < 0.001). Among circumferential strips, %SR of MVAC was greater than for either AV or MVF (p = 0.014). With increasing age, %SR was increased in AV circumferential strips (p = 0.009) and in MVF circumferential sections (p = 0.005).

Total percentage of relaxation (SR) for (
Units for Kfast and Kslow are N/s. Data represent mean ± standard error of the mean. Samples with the same numeric superscripts indicate significant differences between Kfast as well as Kslow values for those samples across ages, except four in which Kfast was only significantly different between samples and six in which only Kslow was significantly different between samples. All differences marked represent p < 0.05.
Significant difference between indicated ages within a given valve region/direction.
Rad, radial; circ, circumferential.
Extensibility
In both the AV and MV, circumferential strips were less extensible than radial strips (p < 0.001; Fig. 6). There were no differences between extensibilities of AV and MV strips. In both MV and AV radial strips, as well as in MVAC circumferential strips, extensibility was reduced with age (each p ≤ 0.010). There was no effect of age on the extensibility of MVF circumferential strips. Among AV circumferential strips, the 6-month-old samples were the least extensible (p < 0.001).

Extensibility for (
Radius of Transition Curvature (RTC)
AV circumferential RTC was less than AV radial RTC (p < 0.001; Fig. 7). MVAC circumferential RTC was less than MVF circumferential RTC and then MV radial RTC (p < 0.001). Among circumferential strips, AV RTC was less than MVAC RTC, which was less than MVF RTC (p < 0.001). With age, MVF circumferential and MV radial RTC increased (p = 0.003 and p = 0.033), but no age-associated changes were evident in AV strips.

Radius of transition curvature: (
Differences between AV leaflets
In 6-year-old AVs, the noncoronary (NC) leaflet was thicker than the left or right coronary leaflets (p = 0.007), as previously reported for human valves.40,41 However, there was no significant difference in thickness between AV leaflets in the other age groups. For both radially and circumferentially oriented strips, Kfast of the NC leaflet was greater than for the other two leaflets (each p ≤ 0.037). However, the circumferential %SR for the NC leaflet was less than for the other two leaflets (p = 0.039). RTC of AV NC leaflets was less than the coronary leaflets (p = 0.031).
Differences in valve matrix composition
By Movat staining, it was evident that the annulus portion of MV radial strips (corresponding to the MVAC) contained a higher proportion of the collagen-rich fibrosa layer compared with AV radial strips (Fig. 8). MVAC circumferential strips contained a much thicker fibrosa than found in MVF (Fig. 9). With age, the fibrosa thickened and collagen content increased throughout the valve layers, but this increase was particularly evident in the fibrosa layer and in the ventricularis layer, where collagen surrounded elastic fibers. Interestingly, there was also a marbling of PG/GAGs (see inset in Fig. 8) in the fibrosa of the 6-year-old MV radial annulus and MVAC circumferential section, which was not evident in the younger ages. Immunohistochemistry revealed that the PGs versican, biglycan, and decorin, as well as the GAG HA, were present within this marbling. Circumferential strips also contained more Col III with age (Fig. 10), as found previously for radial strips. 11 The myofibroblast marker NMM was more strongly expressed in 6-year-old strips of MVAC than in younger valves.

Representative Movat-stained tissue sample sections for 6-week-old and 6-year-old radial strips. In Movat pentachrome-stained tissue, yellow = aligned collagen, black = elastic fibers, and green/blue = proteoglycans (PGs)/glycosaminoglycans (GAGs). To facilitate comparison between AV and MV, valves were all oriented with the fibrosa at the bottom of the image and sized to show the complete leaflet thickness allowing visualization of all leaflet layers. Because of differences in leaflet thickness (see Fig. 2), magnifications differ between images. Scale bars equal 200 μm. Note that the MV radial annulus is almost exclusively fibrosa (corresponding to the MVAC), whereas AV radial annulus shows a considerably smaller proportion of the leaflet thickness composed of fibrosa. Also of note is the marbling of PG/GAG with collagen in the fibrosa of 6-year-old MV radial annular strips (see close-up view in inset) and MVAC circumferential strips compared with 6-week-old strips. This marbling may, in part, account for the increase in SR seen with age (see Discussion section). Color images available online at www.liebertonline.com/ten.

Representative Movat-stained tissue sample sections for 6-week-old, 6-month-old, and 6-year-old circumferential strips. In Movat pentachrome-stained tissue, yellow = aligned collagen, black = elastic fibers, and green/blue = PGs/GAGs. To facilitate comparison between AV and MV, valves were all oriented with the fibrosa at the bottom of the image. Images in column 1 show high magnification of 6-week-old valves (all images in column 1 are the same magnification). Images in columns 2–4 are lower magnification to demonstrate the differences in leaflet thickness (all images in columns 2–4 are the same magnification). Scale bars equal 200 μm. Compared with MVF, MVAC circumferential sections contain much higher proportion of fibrosa. For all valve regions an increase in collagen was noted throughout the valve layers with age, but most particularly in the fibrosa and ventricularis. Color images available online at www.liebertonline.com/ten.

MVAC circumferential sections showing increased Col III and NMM in 6-year-old sections compared with 6-month-old sections. Similar differences were noted for the staining of Col III for other valve regions between the 6-year-old aged sections and the younger ages. Scale bar equals 200 μm. Insets in NMM images are 1.5 × higher magnification for improved visualization of cells. Col III, collagen type III; NMM, nonmuscle myosin heavy chain-IIB. Color images available online at www.liebertonline.com/ten.
Discussion
This study showed for the first time that there are significant age-related differences in numerous material properties of MV and AV. At each age, the valves were more extensible and less stiff radially than circumferentially, as is well known.15–20 All valves, however, demonstrated age-associated increases in post-transition stiffness and most also showed increased pre-transition stiffness. In general, extensibility was reduced with age, especially radially. Stress relaxation parameters also varied with age, particularly circumferentially. The overall shape of the tension–strain curves also changed with age, as demonstrated by the novel parameter RTC, which generally increased with age in the MV and was larger radially than circumferentially. These age-related changes in material properties paralleled an increase in collagen content, particularly in the fibrosa and ventricularis (which is normally rich in elastic fibers), and a marbling appearance of the fibrosa, in which PG/GAGs were interspersed throughout the collagen. Because this study tested MV and AV from the same hearts, it was also possible to demonstrate differences in material properties between MV and AV.
With increasing age, the leaflets showed greater post-transition stiffness both radially and circumferentially. This finding is consistent with age-related stiffening of other connective tissues such as cartilage 42 and tendon, 43 which are mainly attributed to increased collagen crosslinking. 43 Leaflet composition explains much of the differences in stiffness among sample groups. For example, the preponderance of the collagenous fibrosa layer within the MV annulus likely contributes to the greater radial stiffness of MV compared with AV and the greater circumferential stiffness of MVAC compared with MVF. Age-related increases in fibrosa thickness and collagen content shown here for numerous valve regions, as well as previously reported increases in collagen1,2 and in crosslinking2,4 in aging valves, could explain the increased stiffness with age. A greater proportion of cells in the 6-year-old MVAC expressed the activated myofibroblast phenotype, which has been linked with increased leaflet stiffness. 38 In addition, collagen crosslinking has been directly correlated with the stiffness of circumferential strips of human AV, although no such correlation was found for radial strips. 26 Combined with evidence suggesting that the ventricularis is a key contributor to valvular material behavior in the radial direction,22,44 the age-related increase in radial stiffness could be due to increased collagen in the ventricularis, as observed in the Movat-stained sections.
Interestingly, the pre-transition stiffness also significantly increased with age. This parameter may be due in part to the stiffness of the elastic fibers, 45 but is not often reported15,16,18,46 and has never previously been studied with respect to age. The increase in pre-transition stiffness with age could reflect the stiffening of aging elastic fibers, 43 as well as increased collagenous reinforcement of the elastic fibers, which was observed in Movat-stained sections. It was also noteworthy that the age-related increases in the thickness of the valve leaflets 40 were not proportional to the increase in heart size, possibly due to the ventricular hypertrophy found in older animals.47,48
Accompanying the age-related increase in the post-transition stiffness was a pronounced reduction in radial extensibility of older valves. Extensibility assesses the amount of tissue that can stretch before the collagen fibers are fully uncrimped and able to bear load; 21 its magnitude reflects both collagen fiber crimp and alignment. The finding that extensibility is greater radially than circumferentially confirms previous reports;16–18,20 this anisotropy permits the leaflets to stretch radially during valve closure. These findings are also consistent with a previous report on an age-related reduction in maximal radial stretch of human cryopreserved valves. 9 A reduction in extensibility with age could be attributed to factors that impact collagen uncrimping, such as increased collagen crosslinking2,4 or elastic fiber fragmentation, 43 which can lead to gradual permanent tissue stretch and less-crimped collagen fibers. Indeed, in tendon collagen, crimp amplitude decreases and crimp wavelength increases with age; 49 these extensibility changes motivate future study of collagen crimp in aging heart valves. Similarly, heart valves subjected to glutaraldehyde fixation (crosslinking) under pressure demonstrated reduced collagen crimp and extensibility. 50
The novel finding of an increase in the percentage of circumferential stress relaxation (in AV and MVF) with age could be due to many factors. First, there is more HA with age 10 as well as marbling of PGs/GAGs in the collagenous fibrosa, which would allow more sliding of collagen bundles relative to one another. 51 Alternatively, the greater numbers of the PGs decorin and biglycan, 10 which bind to the surface of collagen fibrils, 52 could be transferring more load from collagen to other matrix components, resulting in greater %SR. Indeed, age-related changes in %SR were most evident in the circumferential direction, the predominant collagen fiber direction. 53 Reductions in the functionality 43 and content 10 of elastic fibers could also increase stress relaxation, according to studies on pericardium. 54 Lastly, the greater abundance of Col III, which dissipates more energy than does collagen type I, 55 could have contributed to the greater %SR with age. There is some disagreement in the heart valve literature regarding anisotropy in stress relaxation.14,17,20,23,24 Our finding of greater %SR in circumferential strips compared with radial strips corroborates the work by Lee et al., who measured the %SR of fresh leaflets over 1000 s using uniaxial testing. 20
The viscoelasticity of valve tissues has been examined using a wide variety of constitutive models,17,19,33,56–59 but the effects of age have not been previously investigated. In particular, the quasi-linear viscoelastic (QLV) model developed by Fung 60 has been applied to fresh and fixed valve tissues;19,33,56 analysis using this technique, however, is computationally intense. In this study, a Maxwell–Weichert model with three elements was used to examine biphasic decay behavior of the stress relaxation curve.61,62 The age-related decreases in Kfast of AV circumferential samples and in Kslow of MVF samples both indicate more viscous behavior, which could relate to the increased abundance of PGs and GAGs. 10 The Kfast and Kslow values reported in this study also serve as baseline parameters for evaluating mechanical properties of TEHVs. While current TEHV studies have focused primarily on linear material properties, understanding the viscoelastic properties will continue to be of significant value, especially as more advanced materials are developed.
The novel parameter RTC, which quantifies the acuteness of the transition region of the tension–strain curve, can be compared across samples for any soft tissues that display bilinear tension–strain or stress–strain curves. Collagen fibers that are highly aligned and uniformly crimped will have a very acute transfer of load bearing from elastic fibers to the straightened collagen fibers, and hence a very small RTC. If collagen fibers are not uniformly crimped, their transition to load bearing would be more gradual, resulting in a greater RTC. A larger RTC would also result from a random fiber alignment, in which some fibers would be initially aligned with the applied load and begin bearing load at low strains, but the remaining collagen bundles would require varying amounts of additional strain and/or rotational realignment before bearing load. 21 Indeed, the angular distribution of collagen fibers significantly changes within this transition region. 21 The smaller RTC of MVAC circumferential strips compared with MVF likely indicates greater alignment 53 and homogeneity of collagen fiber crimping in MVAC compared with the other valve regions. In the valve, collagen bundles are predominantly oriented circumferentially, 53 and therefore, RTC would be lower circumferentially than radially, as found here. The age-related increase in RTC of MVF and MV radial strips could be due to greater collagen crosslinking2,4 and a greater abundance of network-forming Col III. In future, it will be important to relate this parameter to changes in the soft tissue microstructure during this transition region of loading.
This study is also the first to show differences in material properties between NC and coronary AV leaflets, such as lower %SR and RTC in the NC leaflet. These differences may be due to the lack of coronary blood flow within the NC sinus; simulations of human AVs have shown that NC experiences higher stress 63 and greater diastolic pressure loading than the other leaflets. 64 The NC leaflet also contains smaller diameter collagen fiber bundles, 65 which may permit easier rotation of collagen fibers as they align in the direction of loading, and hence contribute to a lower RTC. An additional subset of valves was frozen at −20°C for 2 months before testing; as previously demonstrated by Clark, 18 significant differences were found between fresh AV and previously frozen AV (data not shown).
Conclusion
In this study, AV and MV from the same porcine hearts were subjected to tensile testing. In the analysis of the tensile testing data, we developed a novel parameter, RTC, which may be useful in quantifying the dynamic stress–strain behavior of biological tissues. Overall, our results showed that between AV and MV, and within each valve, there are profound heterogeneities and age-related changes that reflect the ECM and microstructural composition. In addition, the effects of age on the heart valve tissues differ between valves, valve regions, and material testing orientations. Given the documented similarities between porcine and human valves, these results provide further justification (as well as baseline data) for the development of age-specific TEHVs. These age-related changes in valve material properties likely contribute to proper physiologic function of valves at different ages. For instance, the age-related increase in leaflet stiffness likely allows the AV and MV to withstand the increase in cardiac pressures and concomitant decrease in aortic compliance with age. 12 Clearly, considerable work remains in understanding the structural and mechanical complexity of heart valves, their ECM and cell composition, and the contributions of these characteristics to valve function.
Footnotes
Acknowledgments
The authors thank Janet Barzilla, Ph.D., Joyce J. Kuo, B.S., and Daniel Laucirica in the Grande-Allen lab, and Anita Mol, Ph.D., and Martijn Cox, M.Sc., from the Eindhoven University of Technology exchange program. Special thanks also to Scott Baggett, Ph.D., for statistical expertise; Benjamin K. Stephens, Ph.D., for mathematical consultation regarding curve fitting and radius of curvature analysis; Roman M. Natoli, Ph.D., for consultation on constitutive modeling; and Larry Fisher, Ph.D., for the generous gift of multiple antibodies.
Disclosure Statement
No competing financial interests exist.
