# Influences of Desmin and Keratin 19 on Passive Biomechanical Properties of Mouse Skeletal Muscle

Authors: Sameer B. Shah; James M. Love; Andrea O'Neill; Richard M. Lovering; Robert J. Bloch
Journal: Journal of Biomedicine and Biotechnology
 (2012)

---

## Abstract

In skeletal muscle fibers, forces must be transmitted between the plasma membrane and the intracellular contractile lattice, and within this lattice between adjacent myofibrils. Based on their prevalence, biomechanical properties and localization, desmin and keratin intermediate filaments (IFs) are likely to participate in structural connectivity and force transmission. We examined the passive load-bearing response of single fibers from the extensor digitorum longus (EDL) muscles of young (3 months) and aged (10 months) wild-type, desmin-null, K19-null, and desmin/K19 double-null mice. Though fibers are more compliant in all mutant genotypes compared to wild-type, the structural response of each genotype is distinct, suggesting multiple mechanisms by which desmin and keratin influence the biomechanical properties of myofibers. This work provides additional insight into the influences of IFs on structure-function relationships in skeletal muscle. It may also have implications for understanding the progression of desminopathies and other IF-related myopathies.


---

## Body

## 1. Introduction

Lateral force transmission, both across the sarcolemma and within a single myofiber, is believed to be an important component of the biomechanical function of muscle [1–4]. Forces must be transmitted both between the plasma membrane and the intracellular myofibrillar lattice, and also within this lattice, from myofibril to myofibril. Based on their prevalence, biomechanical properties and localization, intermediate filaments (IFs) are strong candidates to regulate such lateral transmission. IFs composed of desmin and keratins 8 and 19 have been identified as key components of mature skeletal muscle [5–7]. Considerable evidence links these IFs to roles in structural connectivity and mechanical function.Desmin, a type III IF expressed exclusively in striated muscle, localizes around Z-disks in the extramyofibrillar space, around mitochondria, at the periphery of cell nuclei, and at costameres, which are periodic sarcolemmal protein complexes lying over the Z-disks and M-bands of peripheral myofibrils [8–14]. Structurally, the loss of desmin results in misaligned myofibrils and disrupted costameres [15, 16]. Mechanically, muscles lacking desmin display reduced maximal isometric force, increased fatigability, poor transverse and longitudinal coupling, and slight increases in passive stiffness, particularly in older animals [17–19]. Although there are discrepancies regarding the susceptibility of desmin-null muscles to injury, these may be attributable to differences in the strains of desmin-null mice probed, the identity of muscles tested, or injury protocols imposed [17, 18].Keratin 19 (K19) and keratin 8 (K8), type I and type II IF proteins, respectively, are the best characterized of the keratins expressed in mature striated muscle. At the sarcolemma, K8 and K19 colocalize with and specifically interact with costameric proteins, including dystrophin, overlying both the Z-disks and M-bands of adjacent myofibrils [20, 21]. Within the myofibrillar lattice, K19 and K8 concentrate primarily around Z-disks [21]. Structurally, the absence of K19 results in the partial disruption of costameres and the dramatic separation of the sarcolemma from adjacent myofibrils by a large gap within which mitochondria accumulate [18, 22]. A small but significant shift in the lateral spacing between myofibrils is also observed in the absence of K19; however, myofibrillar desmin localization is undisturbed [18, 22]. Physiologically, the absence of K19 results in decreased maximal isometric force and increased plasma creatine kinase levels [18, 22]. Double-null muscles, lacking both K19 and desmin, display increased susceptibility to injury and significant increases in membrane instability. Decreases in specific tension are greater than that of K19-null muscles, but equal to those in desmin-null muscles [18].Given the possibility that extracellular matrix remodeling could influence the mechanical properties of whole muscle [23–25], experiments on single fibers have been performed to examine structural connectivity during passive mechanical stretch [26, 27]. Analysis of desmin-null fibers confirmed the necessity of desmin for myofibrillar and nuclear alignment during mechanical loading. Unlike observations in whole muscle, however, isolated desmin-null fibers do not appear to display misaligned costameres, even in the presence of considerable shear force. In addition, desmin-null fibers are more compliant than their wild-type counterparts.In this study, we have extended the examination of biomechanical roles for IFs in single fibers. We investigated the structural changes and tensile load-bearing response of passively loaded fibers harvested from the extensor digitorum longus (EDL) muscles of young (3 months) and middle-aged (10 months) wild-type, desmin-null, K19-null, and desmin/K19 double-null FVB mice. We observed that, although fibers are more compliant in both desmin and K19-null knockout mice, the structural changes that occur within the two types of fibers during loading are different. Our data suggest that IFs influence force transmission in myofibers through multiple biomechanical pathways.


## 2. Methods

### 2.1. Animals

Female FVB mice homozygous for the desmin-null (des −/−), K19-null (K19 −/−), and desmin/K19 double-null genotypes were bred and genotyped as described previously [18] and compared to age- and sex-matched wild-type FVB mice (Taconic, Hudson, NY, USA). Mice were studied at two ages, 3 months and 10 months. They were euthanized by carbon dioxide inhalation immediately prior to tissue harvest. Institutional Animal Care and Usage (IACUC) Committees of the University of Maryland School of Medicine and the University of Maryland, College Park, approved all animal use protocols.


### 2.2. Fiber Dissection and Attachment to Mechanical Apparatus

Tissue harvest, muscle storage, and fiber dissection protocols have been described previously [27]. Briefly, EDL muscles were dissected from mouse hindlimbs, bathed in relaxing solution for 60 min and placed in a glycerol-based storage solution for a minimum of 24 h and a maximum of one week at −20°C. All solutions were made in the presence of the protease inhibitor leupeptin. Segments of single intact fibers were carefully dissected and transferred in relaxing solution to a chamber housed in a custom-made mechanical apparatus. This apparatus is similar to that previously used [26, 27], with the only difference being the modification of the structural frame and stage insert to be mounted on an inverted wide-field fluorescence microscope (Nikon TE-2000PFS, Melville, NY, USA). Within the chamber, the segment was securely tied with 10-0 monofilament suture to pins projecting from a force transducer (Aurora Scientific 405A, Aurora, Ontario, Canada) and a rotational bearing (Newport MTRS, Irvine, Calif, USA). This configuration enabled the high-resolution imaging of fiber and sarcomere geometry during passive loading. To minimize any experimental artifacts associated with glycerol-mediated skinning of fibers during storage, segments displaying abnormal discoloration, localized swelling, or otherwise irregular geometry or appearance were discarded.


### 2.3. Passive Mechanical Protocol and Imaging

The fiber segment was brought to its resting length, determined as the knot-to-knot segment length at which passive tension was just measurable above the noise level of the force transducer. The segment was preconditioned by three loading cycles of 15% fiber strain, resulting in a change in resting length of <3% from the resting length before preconditioning. The segment was then loaded in increments of ~10% of the segment length at a rate of 0.05 mm/s. At each length, after two minutes of stress relaxation, images of the fiber were captured in a middle region of the fiber, and steady-state force (within 5% of the final force value) was recorded. During the 2 min of stress relaxation, the imaging field was adjusted to visualize the same region and focal plane of the fiber at each strain. Phase-contrast imaging was performed on an inverted microscope (Nikon TE-2000PFS, Melville, NY), using 10x and 20x CFI Plan Apo objectives (Nikon). Exposure times were 5–20 ms; image acquisition was controlled by Elements software (Nikon, Melville, NY, USA). A custom built environmental chamber (Precision Plastics, Beltsville, Md, USA) maintained temperature (25°C) and humidity during imaging.


### 2.4. Calculation of Biomechanical Parameters

Image analysis was performed using Elements software and Micron software (Westover Scientific, Mill Creek, Wash, USA). Fiber cross-sectional area was calculated using a circular approximation, based on diameters measured from a phase contrast image of the fiber at each strain. Area was multiplied by fiber length to estimate fiber volume, assuming a cylindrical fiber. Cauchy (true) stress (σC) was calculated by dividing force by the current cross-sectional area at each strain, and LaGrangian stress (σLG) was calculated by dividing force by the cross-sectional area of the undeformed fiber. Fiber strain (εF) was calculated as the ratio of the imposed deformation to the resting knot-to-knot fiber length. Sarcomere length (SL) was calculated based on the average spacing between 20 consecutive Z-disks, as measured from phase images of fibers at each fiber strain. SL in different regions of the fiber varied by less than 0.1 μm, and therefore SLs were measured at three locations per image and averaged. Sarcomere strain (εSL) was calculated as the ratio of the change in SL to the SL at resting fiber length. Strain in the sarcolemmal plasma membrane in the vicinity of the imaged sarcomeres, or local membrane strain, was calculated based on changes in the spacing between membrane markers, or pieces of residual debris on the external surface of the fiber; local membrane strain (εM) was calculated as the ratio of the change in spacing between membrane markers to marker spacing at resting fiber length. Radial strain (εR) was calculated as the ratio of the change in fiber diameter to the fiber diameter at resting fiber length. Poisson’s ratio (ν) was calculated based on the power-law relationship between longitudinal and radial strain. Biomechanical parameters and their relationships with each other are summarized in Table 1.Table 1
Summary and explanation of biomechanical parameters.
Parameter (units)SymbolFormulaLay description of parameter and relevanceForce (N)F—Nonnormalized forceUndeformed fiber length (m)Lo—Initial knot-knot fiber lengthUndeformed diameter (m)do—Initial fiber diameterUndeformed area (m2)Ao—Initial fiber cross-sectional areaDeformed fiber length (m)Lf—Fiber length following longitudinal deformationDeformed diameter (m)df—Fiber diameter following longitudinal deformationDeformed area (m2)Af—Fiber cross-sectional area following longitudinal deformationLongitudinal fiber strain (dimensionless)εL(Lf-Lo)/LoNormalized longitudinal deformation—percentage fiber extensionFiber volumeVfLf*AdFiber volume after deformation, assuming a cylindrical geometryRadial strain (dimensionless)εRdf-do/doNormalized lateral deformation—percentage reduction in fiber diameter following longitudinal deformation; increased radial strain represents stronger coupling between longitudinal and lateral responses.Cauchy (true) stress (Pa)σCF/AdForce normalized to deformed area—couples force generation to changes in 3D-fiber geometryLaGrangian stress (Pa)σLGF/AoForce normalized to undeformed area—measure of inherent load-bearing ability, decoupled from changes in 3D-fiber geometryPoisson’s ratio (dimensionless)ν(1+εL)=(1+εL)-νMeasure of relation between longitudinal strain and radial strain—lateral strain transmission


### 2.5. Statistics

Regression analysis was performed using Excel (Microsoft, Seattle, Wash, USA). One-way ANOVAs examining the effect of genotype were performed independently on mechanical data from young and aged mice, since an additional double-null group was evaluated in aged mice. Data from double-null mice were not considered for subsequent two-way ANOVA analysis, which examined effects and interactions of age and genotype on mechanical parameters. ANOVA and post hoc two-tailedt-tests assuming unequal variance were performed using Prism (Graph Pad Software, Inc., La Jolla, Calif, USA). Values are given as mean ± SEM.


## 2.1. Animals

Female FVB mice homozygous for the desmin-null (des −/−), K19-null (K19 −/−), and desmin/K19 double-null genotypes were bred and genotyped as described previously [18] and compared to age- and sex-matched wild-type FVB mice (Taconic, Hudson, NY, USA). Mice were studied at two ages, 3 months and 10 months. They were euthanized by carbon dioxide inhalation immediately prior to tissue harvest. Institutional Animal Care and Usage (IACUC) Committees of the University of Maryland School of Medicine and the University of Maryland, College Park, approved all animal use protocols.


## 2.2. Fiber Dissection and Attachment to Mechanical Apparatus

Tissue harvest, muscle storage, and fiber dissection protocols have been described previously [27]. Briefly, EDL muscles were dissected from mouse hindlimbs, bathed in relaxing solution for 60 min and placed in a glycerol-based storage solution for a minimum of 24 h and a maximum of one week at −20°C. All solutions were made in the presence of the protease inhibitor leupeptin. Segments of single intact fibers were carefully dissected and transferred in relaxing solution to a chamber housed in a custom-made mechanical apparatus. This apparatus is similar to that previously used [26, 27], with the only difference being the modification of the structural frame and stage insert to be mounted on an inverted wide-field fluorescence microscope (Nikon TE-2000PFS, Melville, NY, USA). Within the chamber, the segment was securely tied with 10-0 monofilament suture to pins projecting from a force transducer (Aurora Scientific 405A, Aurora, Ontario, Canada) and a rotational bearing (Newport MTRS, Irvine, Calif, USA). This configuration enabled the high-resolution imaging of fiber and sarcomere geometry during passive loading. To minimize any experimental artifacts associated with glycerol-mediated skinning of fibers during storage, segments displaying abnormal discoloration, localized swelling, or otherwise irregular geometry or appearance were discarded.


## 2.3. Passive Mechanical Protocol and Imaging

The fiber segment was brought to its resting length, determined as the knot-to-knot segment length at which passive tension was just measurable above the noise level of the force transducer. The segment was preconditioned by three loading cycles of 15% fiber strain, resulting in a change in resting length of <3% from the resting length before preconditioning. The segment was then loaded in increments of ~10% of the segment length at a rate of 0.05 mm/s. At each length, after two minutes of stress relaxation, images of the fiber were captured in a middle region of the fiber, and steady-state force (within 5% of the final force value) was recorded. During the 2 min of stress relaxation, the imaging field was adjusted to visualize the same region and focal plane of the fiber at each strain. Phase-contrast imaging was performed on an inverted microscope (Nikon TE-2000PFS, Melville, NY), using 10x and 20x CFI Plan Apo objectives (Nikon). Exposure times were 5–20 ms; image acquisition was controlled by Elements software (Nikon, Melville, NY, USA). A custom built environmental chamber (Precision Plastics, Beltsville, Md, USA) maintained temperature (25°C) and humidity during imaging.


## 2.4. Calculation of Biomechanical Parameters

Image analysis was performed using Elements software and Micron software (Westover Scientific, Mill Creek, Wash, USA). Fiber cross-sectional area was calculated using a circular approximation, based on diameters measured from a phase contrast image of the fiber at each strain. Area was multiplied by fiber length to estimate fiber volume, assuming a cylindrical fiber. Cauchy (true) stress (σC) was calculated by dividing force by the current cross-sectional area at each strain, and LaGrangian stress (σLG) was calculated by dividing force by the cross-sectional area of the undeformed fiber. Fiber strain (εF) was calculated as the ratio of the imposed deformation to the resting knot-to-knot fiber length. Sarcomere length (SL) was calculated based on the average spacing between 20 consecutive Z-disks, as measured from phase images of fibers at each fiber strain. SL in different regions of the fiber varied by less than 0.1 μm, and therefore SLs were measured at three locations per image and averaged. Sarcomere strain (εSL) was calculated as the ratio of the change in SL to the SL at resting fiber length. Strain in the sarcolemmal plasma membrane in the vicinity of the imaged sarcomeres, or local membrane strain, was calculated based on changes in the spacing between membrane markers, or pieces of residual debris on the external surface of the fiber; local membrane strain (εM) was calculated as the ratio of the change in spacing between membrane markers to marker spacing at resting fiber length. Radial strain (εR) was calculated as the ratio of the change in fiber diameter to the fiber diameter at resting fiber length. Poisson’s ratio (ν) was calculated based on the power-law relationship between longitudinal and radial strain. Biomechanical parameters and their relationships with each other are summarized in Table 1.Table 1
Summary and explanation of biomechanical parameters.
Parameter (units)SymbolFormulaLay description of parameter and relevanceForce (N)F—Nonnormalized forceUndeformed fiber length (m)Lo—Initial knot-knot fiber lengthUndeformed diameter (m)do—Initial fiber diameterUndeformed area (m2)Ao—Initial fiber cross-sectional areaDeformed fiber length (m)Lf—Fiber length following longitudinal deformationDeformed diameter (m)df—Fiber diameter following longitudinal deformationDeformed area (m2)Af—Fiber cross-sectional area following longitudinal deformationLongitudinal fiber strain (dimensionless)εL(Lf-Lo)/LoNormalized longitudinal deformation—percentage fiber extensionFiber volumeVfLf*AdFiber volume after deformation, assuming a cylindrical geometryRadial strain (dimensionless)εRdf-do/doNormalized lateral deformation—percentage reduction in fiber diameter following longitudinal deformation; increased radial strain represents stronger coupling between longitudinal and lateral responses.Cauchy (true) stress (Pa)σCF/AdForce normalized to deformed area—couples force generation to changes in 3D-fiber geometryLaGrangian stress (Pa)σLGF/AoForce normalized to undeformed area—measure of inherent load-bearing ability, decoupled from changes in 3D-fiber geometryPoisson’s ratio (dimensionless)ν(1+εL)=(1+εL)-νMeasure of relation between longitudinal strain and radial strain—lateral strain transmission


## 2.5. Statistics

Regression analysis was performed using Excel (Microsoft, Seattle, Wash, USA). One-way ANOVAs examining the effect of genotype were performed independently on mechanical data from young and aged mice, since an additional double-null group was evaluated in aged mice. Data from double-null mice were not considered for subsequent two-way ANOVA analysis, which examined effects and interactions of age and genotype on mechanical parameters. ANOVA and post hoc two-tailedt-tests assuming unequal variance were performed using Prism (Graph Pad Software, Inc., La Jolla, Calif, USA). Values are given as mean ± SEM.


## 3. Results

### 3.1. Intermediate Filament Influences on Passive Properties of Fibers from Aged Animals

We examined the response of EDL fibers harvested from 10-month-old wild-type, desmin-null, K19-null, and desmin/K19 double-null mice to tensile loading to examine the relative contributions of each intermediate filament to fiber passive mechanical properties. A summary and a brief description of measured and calculated parameters are provided in Table1.For stress-strain curves, regression coefficients were used to compare the relative compliance of fiber populations. Cauchy stress was plotted against fiber strain for each experiment and fit using quadratic regression (Figure1(a), cf. [26]). Fibers from all genotypes displayed a strong quadratic relationship between stress and fiber strain (r2>0.78 for all experiments), indicating an elastic modulus that varies linearly with strain. ANOVA revealed a significant effect of genotype on the quadratic regression coefficient (P<0.0001). Post hoc comparison of the quadratic coefficients revealed that wild-type fibers were significantly stiffer than desmin-null (P<0.002), K19-null (P<0.006), and double-null (P<0.002) fibers (Figure 1(b)). There were no significant differences in the quadratic regression coefficients of K19-null and desmin-null or desmin-null and double-null fibers. Interestingly, however, a significant difference (P<0.03) was observed between K19-null and double-null fibers for Cauchy stress-strain curves.Mechanical response of fibers from aged animals. Raw data (a, c, e) and regression coefficients (b, d, f) for (a-b) Cauchy stress-fiber strain (quadratic regression); (c-d) LaGrangian stress-fiber strain (linear regression); (e-f) radial strain-fiber strain (exponential regression) for fibers from 10-month-old wild-type (blue), desmin-null (red), K19-null (orange), and double-null (green) animals (n=7 fibers/genotype). Different symbols within a color group represent data from individual experiments. Displayed regression lines in (a), (c), and (e) represent averaged coefficients for individual experiments reported in (b), (d), and (f). * Indicates significant difference from wild-type; § indicates significant difference from K19-null; ¥ indicates significant difference from corresponding genotype in younger animals. P values are reported in the text.
(a)(b)(c)(d)(e)(f)To determine whether these differences in fiber populations were strictly a result of changes in fiber geometry or also inherent load-bearing ability, LaGrangian stress was plotted against fiber strain (Figure1(c)). Use of the undeformed area to calculate stress had the effect of flattening the curves, so linear regression was used (r2>0.81 for all experiments) to describe the relationship between stress and fiber strain. ANOVA again revealed a significant effect of genotype on the regression coefficient (P<0.0001). Comparison of linear coefficients confirmed that wild-type fibers were significantly stiffer than desmin-null (P<0.005), K19-null (P<0.007), and double-null (P<0.004) fibers (Figure 1(d)). There were no significant differences in the linear regression coefficients of K19-null and desmin-null, desmin-null and double-null, or K19-null and double-null fibers. Consistent with results from Cauchy stress-strain relationships, though, a strong but nonsignificant trend (P<0.09) was observed for the latter.Analysis of radial strain revealed additional differences in the mechanical response of fibers from different genotypes (Figure1(e)). First-order exponential functions were used to describe the relationship between radial and longitudinal strains (r2>0.83 for all experiments). Although, as expected, ANOVA revealed a significant effect of genotype on the exponential decay constant (P<0.04), post hoc comparison (Figure 1(f)) revealed that wild-type fibers exhibited tighter coupling between longitudinal and radial strains than desmin-null (P<0.04), but not K19-null (P<0.24) fibers. Wild-type fibers also trended towards a higher time constant (tighter coupling) than double-null fibers, but were not significantly different (P<0.06). As was the case for stress-strain relationships, there was no difference between the relationship between radial and longitudinal strain for desmin-null and double-null fibers.We further examined the coupling of longitudinal and radial strains in the context of fiber anisotropy, by calculating Poisson’s ratio. A power-law was used to fit the relationship between longitudinal and radial strain, with the exponent denoting Poisson’s ratio for each fiber.R2 values for each experiment were reasonably high (r2>0.62), though not as strong as an exponential regression fit. For all genotypes, Poisson’s ratios exceeded 0.5, indicating significant anisotropy. There were no significant differences in Poisson’s ratio between any of the genotypes (Mean ± SEM for wild-type: 0.81 ± 0.10; des −/−: 0.63 ± 0.08; ker −/−: 0.71 ± 0.05; double KO: 0.61 ± 0.04), though a trend towards a slightly higher Poisson’s ratio was observed for wild-type fibers compared to double-null fibers (P<0.09).Consistent with the high Poisson’s ratios, when changes in fiber volume were plotted against longitudinal strain, each genotype displayed a decrease in volume with increasing strain, indicating that loading was not isovolumic. This decrease was usually linear, and therefore linear regression was used to fit this data (r2>0.68). Wild-type fibers displayed the most dramatic decrease in volume (Mean slope ± SEM: −0.38 ± 0.11), though there were no significant differences in first-order regression slopes between any of the genotypes (des −/−: −0.16 ± 0.05; ker −/−: −0.22 ± 0.04; double KO: −0.13 ± 0.04). As for Poisson’s ratio, a trend towards a slightly higher magnitude of regression slope was observed for wild-type fibers compared to double-null fibers (P<0.08). Significant differences between parameters are summarized in Table 2.Table 2
Summary of parameters that display significant post hoc differences (P<0.05) between genotypes for fibers harvested from aged mice. Raw data, means, and standard errors are provided in Figure 1. Parameters that display non-significant trends (P<0.09) are provided in parentheses. NS: No significant difference.
Genotypedes −/−K19 −/−Double −/−+/+σLG, σC, εRσLG, σCσLG, σC, (εR), (ν), (Vf)des −/−NSNSK19 −/−σC, (σLG), (εR)


### 3.2. Intermediate Filament Influences on Passive Properties of Fibers from Young Animals

Passive biomechanical analysis was then performed on EDL fibers harvested from 3-month-old wild-type, desmin-null, and K19-null mice, to probe whether phenotypes observed in the older 10-month-old animals originate at earlier time points. As was observed for older animals, fibers from all genotypes displayed a strong quadratic relationship between Cauchy stress and fiber strain (r2>0.82 for all experiments; Figure 2(a)) and a linear relationship between LaGrangian stress and fiber strain (r2>0.86 for all experiments; Figure 2(c)). For both Cauchy and LaGrangian stress-strain, ANOVA revealed a significant effect of genotype on regression coefficients (quadratic: P<0.002; linear: P<0.0002, resp.). Post hoc comparison of the quadratic coefficients revealed that wild-type fibers were significantly stiffer than desmin-null (P<0.01) and K19-null (P<0.01) fibers (Figure 2(b)). Similarly, comparison of linear coefficients confirmed that wild-type fibers were significantly stiffer than desmin-null (P<0.01) and K19-null (P<0.01) fibers (Figure 2(d)).Mechanical response of fibers from young animals. Raw data (a, c, e) and regression coefficients (b, d, f) for (a-b) Cauchy stress-fiber strain (quadratic regression); (c-d) LaGrangian stress-fiber strain (linear regression); (e-f) radial strain-fiber strain (exponential regression) for fibers from 3-month-old wild-type (blue), desmin-null (red), and K19-null (orange) animals (n=5 fibers/genotype). Different symbols within a color group represent data from individual experiments. Displayed regression lines in (a), (c), and (e) represent averaged coefficients for individual experiments reported in (b), (d), and (f). * indicates significant difference from wild-type. P values are reported in the text.
(a)(b)(c)(d)(e)(f)Analysis of radial strain in fibers from younger animals, in contrast to observations in fibers from older animals, also did not reveal any differences between desmin-null and K19-null genotypes (Figure2(e)). First-order exponential functions were again used to describe the relationship between radial and longitudinal strain (r2>0.78 for all experiments). ANOVA revealed a significant effect of genotype on the exponential decay constant (P<0.0003). Post hoc comparison of the time constants revealed that wild-type fibers exhibited a stronger inverse relationship, and thus tighter coupling, between longitudinal and radial strains than desmin-null (P<0.01) and K19-null (P<0.01) fibers (Figure 1(f)).This finding was consistent with our analysis of Poisson’s ratio. As for fibers from older animals, regression revealed a power-law relationship between longitudinal and radial strains (r2>0.61), and Poisson’s ratios for all genotypes far exceeded 0.5, indicating considerable anisotropy. There were no significant differences in Poisson’s ratio between desmin-null and K19-null genotypes (Mean ± SEM for des −/−: 1.08 ± 0.12; ker −/−: 1.19 ± 0.08), but the Poisson’s ratio of wild-type fibers (2.43 ± 0.33) was significantly higher than those of both desmin-null and K19-null fibers (P<0.02). With respect to changes in volume in response to loading, as for older animals, all genotypes decreased in volume with increased deformation, and wild-type fibers displayed the most dramatic decrease in volume (Mean slope ± SEM: −1.08 ± 0.19). This decrease was significantly greater (P<0.01) than that in both desmin-null and K19-null fibers (des −/−: 0.24 ± 0.05; ker −/−: 0.24 ± 0.02). Significant differences between parameters are summarized in Table 3.Table 3
Summary of parameters that display significant post hoc differences (P<0.05) between genotypes for fibers harvested from young mice. Raw data, means, and standard errors are provided in Figure 2. NS: no significant difference.
Genotypedes −/−K19 −/−+/+σLG, σC, εR, ν, VfσLG, σC, εR, ν, VfDes −/−NSThough we have not yet compared the properties of fibers in 10-month-old to 3-month-old double-null mice, we confirmed differences in passive mechanical parameters between young and aged wild-type, desmin-null, and K19-null mice using two-way ANOVA. Effects of genotype (P<0.0001), but no effect of age and no interaction were found on stress-strain and volume-strain relationships. However, a significant effect of genotype (P<0.0001) and age (P<0.0001) and an interaction between genotype and age (P<0.0001) were observed for the coupling between radial strain and longitudinal strain as indicated by both exponential and power-law relationships. Post hoc analysis revealed a significant reduction in the exponential time constant and Poisson’s ratio of longitudinal and radial strain for wild-type (exp: P<0.008; Poisson’s ratio: P<0.006), desmin-null (exp: P<0.005; Poisson’s ratio: P<0.02), and K19-null (exp: P<0.001; Poisson’s ratio: P<0.001) fibers from aged mice compared to their younger counterparts (Figures 1(e) and 1(f) versus Figures 2(e) and 2(f)).


### 3.3. Structural Bases for Differences in Biomechanical Properties

In an effort to elucidate potential structural bases for differences in passive mechanical function between genotypes, we used phase-contrast microscopy to image the response of fibers from 10-month-old animals of each genotype at higher resolution (Figures3–7). Wild-type fibers displayed well aligned Z-disks, even at high fiber strains (Figure 4). In contrast, desmin-null fibers displayed the characteristic “basket-weave” pattern of misaligned myofibrils (Figure 5). Despite the misalignment of myofibrils in desmin-null mice, local sarcomere strain scaled consistently with fiber strain in both wild-type and desmin-null fibers (Figures 3(a) and 3(b), 4, and 5). Additionally, membrane strain was also tightly correlated with sarcomere strain in both wild-type and desmin-null fibers, indicating continuity between the sarcolemma and myofibrillar lattice (Figures 3(c) and 3(d)).Relationships between global fiber strain and local membrane and sarcomere strain. Raw data (a, c) and regression coefficients (b, d) for (a-b) fiber strain-sarcomere strain; (c-d) membrane strain-sarcomere strain for fibers from 10-month-old wild-type (blue), desmin-null (red), K19-null (orange), and double-null (green) animals (n=3 fibers/genotype). Different symbols within a color group represent data from individual experiments. Displayed regression lines in (a) and (c) represent averaged coefficients for individual experiments reported in (b) and (d).
(a)(b)(c)(d)Figure 4
Geometric and structural response of a 10-month-old wild-type fiber to tensile loading. The same region of the fiber, imaged using phase-contrast microscopy, is shown at five different deformations. Fiber strains, radial strains, and sarcomere lengths are reported to the right of the fiber for each deformation. The fiber displays well aligned Z-lines (dark vertical stripes), consistently increasing sarcomere strain, and decreasing fiber diameter with increased fiber deformation. A sample region of the fiber that displays particularly well aligned Z-lines that are experiencing some shear, inferred through the slanting of Z-lines, is marked with a white box. The boxed region has been computationally enlarged 2x for clarity. At high strains (εF>0.8), artifactual thickening of  Z-lines occurs due to shear-induced offsets of Z-lines in deeper regions of the fiber. Bar = 50 μm.Figure 5
Geometric and structural response of a 10-month-old desmin-null fiber to tensile loading. The same region of the fiber, imaged using phase-contrast microscopy, is shown at five different deformations. Fiber strains, radial strains, and sarcomere lengths are reported to the right of the fiber for each deformation. The fiber displays misaligned Z-lines, but consistently increasing sarcomere strain and decreasing fiber diameter with increased fiber deformation with increased fiber deformation. Misalignment is characterized by a striking basketweave pattern where no clear Z-line orientation may be inferred. A sample region of the fiber that displays this pattern is marked with a white box, and may be compared to the well aligned region shown in Figure4. The boxed region has been computationally enlarged 2x for clarity. Desmin-null fibers display dramatically more debris on their external membrane compared to wild-type fibers. This observation is consistent with the increased “stickiness” experienced during fiber dissection, possibly resulting from differences in extracellular matrix composition. Bar = 50 μm.Geometric and structural response of a 10-month-old K19-null fiber to tensile loading. (a) The same region of the fiber, imaged using phase-contrast microscopy, is shown at five different deformations. Fiber strains, radial strains, and sarcomere lengths are reported to the right of the fiber for each deformation. The fiber displays aligned Z-lines, but discontinuity between fiber deformation and sarcomere strain. This fiber also displays considerable shear, as seen by the slanting of Z-lines (dark lines). A sample region of the fiber that displays particularly well aligned Z-lines that are experiencing some shear is marked with a white box. The boxed region has been computationally enlarged 2x for clarity. There is no slippage or damage of the fiber at the (b) left edge or (c) right edge, imaged in the vicinity of the knot and pins (starred black regions in the periphery of the image) at a moderate fiber strain of 0.5. This may be concluded from the gradual reduction in of the slant of Z-disks away from the knot, which is located to the left in (b) and to the right in (c). Bar = 50μm.
(a)(b)(c)Figure 7
Geometric and structural response of a 10-month-old double-null fiber to tensile loading. The same region of the fiber, imaged using phase-contrast microscopy, is shown at five different deformations. Fiber strains, radial strains, and sarcomere lengths are reported to the right of the fiber for each deformation. Fiber displays both misaligned Z-lines as well as discontinuity between fiber, sarcomere, and radial strains. A sample region of the fiber that displays the basketweave pattern characteristic of misalignment is marked with a white box. The boxed region has been also computationally enlarged 2x for clarity. As for desmin-null fibers, double-null fibers display dramatically more debris on their external membrane compared to wild-type fibers. Bar = 50μm.In contrast, K19-null fibers appeared to have well aligned myofibrils, even in the presence of considerable shear (Figure6). Though differences did not reach statistical significance due to the small sample size, there was considerable discontinuity between fiber strain and sarcomere strain for K19-null fibers compared to wild-type or desmin-null fibers (variability of each genotype in Figures 3(a) and 3(b); examples in Figures 4–6). Examination of the entire length of the fiber at multiple time points during experimentation suggested that this was not simply an artifact of fiber slippage at the knots or visible damage outside of the imaging window (Figures 6(b) and 6(c)). Additionally, similar to findings in desmin-null fibers, the absence of K19 did not appear to disrupt coupling between the sarcolemma and the contractile apparatus significantly, based on the strong correlation between local membrane strain and sarcomere strain (Figures 3(c) and 3(d)).Double-null fibers, as expected, displayed phenotypes from both desmin-null and K19-null fibers. Visually, a lack of Z-disk alignment similar to that in desmin-null fibers was observed (Figure7). As in K19-null fibers, however, there was also considerable discontinuity between fiber strain and sarcomere strain in double-null fibers (Figures 3(a), 3(b), and 7). As for all genotypes, there was no apparent decoupling between the sarcolemma and underlying myofibrils (Figures 3(c) and 3(d)).


## 3.1. Intermediate Filament Influences on Passive Properties of Fibers from Aged Animals

We examined the response of EDL fibers harvested from 10-month-old wild-type, desmin-null, K19-null, and desmin/K19 double-null mice to tensile loading to examine the relative contributions of each intermediate filament to fiber passive mechanical properties. A summary and a brief description of measured and calculated parameters are provided in Table1.For stress-strain curves, regression coefficients were used to compare the relative compliance of fiber populations. Cauchy stress was plotted against fiber strain for each experiment and fit using quadratic regression (Figure1(a), cf. [26]). Fibers from all genotypes displayed a strong quadratic relationship between stress and fiber strain (r2>0.78 for all experiments), indicating an elastic modulus that varies linearly with strain. ANOVA revealed a significant effect of genotype on the quadratic regression coefficient (P<0.0001). Post hoc comparison of the quadratic coefficients revealed that wild-type fibers were significantly stiffer than desmin-null (P<0.002), K19-null (P<0.006), and double-null (P<0.002) fibers (Figure 1(b)). There were no significant differences in the quadratic regression coefficients of K19-null and desmin-null or desmin-null and double-null fibers. Interestingly, however, a significant difference (P<0.03) was observed between K19-null and double-null fibers for Cauchy stress-strain curves.Mechanical response of fibers from aged animals. Raw data (a, c, e) and regression coefficients (b, d, f) for (a-b) Cauchy stress-fiber strain (quadratic regression); (c-d) LaGrangian stress-fiber strain (linear regression); (e-f) radial strain-fiber strain (exponential regression) for fibers from 10-month-old wild-type (blue), desmin-null (red), K19-null (orange), and double-null (green) animals (n=7 fibers/genotype). Different symbols within a color group represent data from individual experiments. Displayed regression lines in (a), (c), and (e) represent averaged coefficients for individual experiments reported in (b), (d), and (f). * Indicates significant difference from wild-type; § indicates significant difference from K19-null; ¥ indicates significant difference from corresponding genotype in younger animals. P values are reported in the text.
(a)(b)(c)(d)(e)(f)To determine whether these differences in fiber populations were strictly a result of changes in fiber geometry or also inherent load-bearing ability, LaGrangian stress was plotted against fiber strain (Figure1(c)). Use of the undeformed area to calculate stress had the effect of flattening the curves, so linear regression was used (r2>0.81 for all experiments) to describe the relationship between stress and fiber strain. ANOVA again revealed a significant effect of genotype on the regression coefficient (P<0.0001). Comparison of linear coefficients confirmed that wild-type fibers were significantly stiffer than desmin-null (P<0.005), K19-null (P<0.007), and double-null (P<0.004) fibers (Figure 1(d)). There were no significant differences in the linear regression coefficients of K19-null and desmin-null, desmin-null and double-null, or K19-null and double-null fibers. Consistent with results from Cauchy stress-strain relationships, though, a strong but nonsignificant trend (P<0.09) was observed for the latter.Analysis of radial strain revealed additional differences in the mechanical response of fibers from different genotypes (Figure1(e)). First-order exponential functions were used to describe the relationship between radial and longitudinal strains (r2>0.83 for all experiments). Although, as expected, ANOVA revealed a significant effect of genotype on the exponential decay constant (P<0.04), post hoc comparison (Figure 1(f)) revealed that wild-type fibers exhibited tighter coupling between longitudinal and radial strains than desmin-null (P<0.04), but not K19-null (P<0.24) fibers. Wild-type fibers also trended towards a higher time constant (tighter coupling) than double-null fibers, but were not significantly different (P<0.06). As was the case for stress-strain relationships, there was no difference between the relationship between radial and longitudinal strain for desmin-null and double-null fibers.We further examined the coupling of longitudinal and radial strains in the context of fiber anisotropy, by calculating Poisson’s ratio. A power-law was used to fit the relationship between longitudinal and radial strain, with the exponent denoting Poisson’s ratio for each fiber.R2 values for each experiment were reasonably high (r2>0.62), though not as strong as an exponential regression fit. For all genotypes, Poisson’s ratios exceeded 0.5, indicating significant anisotropy. There were no significant differences in Poisson’s ratio between any of the genotypes (Mean ± SEM for wild-type: 0.81 ± 0.10; des −/−: 0.63 ± 0.08; ker −/−: 0.71 ± 0.05; double KO: 0.61 ± 0.04), though a trend towards a slightly higher Poisson’s ratio was observed for wild-type fibers compared to double-null fibers (P<0.09).Consistent with the high Poisson’s ratios, when changes in fiber volume were plotted against longitudinal strain, each genotype displayed a decrease in volume with increasing strain, indicating that loading was not isovolumic. This decrease was usually linear, and therefore linear regression was used to fit this data (r2>0.68). Wild-type fibers displayed the most dramatic decrease in volume (Mean slope ± SEM: −0.38 ± 0.11), though there were no significant differences in first-order regression slopes between any of the genotypes (des −/−: −0.16 ± 0.05; ker −/−: −0.22 ± 0.04; double KO: −0.13 ± 0.04). As for Poisson’s ratio, a trend towards a slightly higher magnitude of regression slope was observed for wild-type fibers compared to double-null fibers (P<0.08). Significant differences between parameters are summarized in Table 2.Table 2
Summary of parameters that display significant post hoc differences (P<0.05) between genotypes for fibers harvested from aged mice. Raw data, means, and standard errors are provided in Figure 1. Parameters that display non-significant trends (P<0.09) are provided in parentheses. NS: No significant difference.
Genotypedes −/−K19 −/−Double −/−+/+σLG, σC, εRσLG, σCσLG, σC, (εR), (ν), (Vf)des −/−NSNSK19 −/−σC, (σLG), (εR)


## 3.2. Intermediate Filament Influences on Passive Properties of Fibers from Young Animals

Passive biomechanical analysis was then performed on EDL fibers harvested from 3-month-old wild-type, desmin-null, and K19-null mice, to probe whether phenotypes observed in the older 10-month-old animals originate at earlier time points. As was observed for older animals, fibers from all genotypes displayed a strong quadratic relationship between Cauchy stress and fiber strain (r2>0.82 for all experiments; Figure 2(a)) and a linear relationship between LaGrangian stress and fiber strain (r2>0.86 for all experiments; Figure 2(c)). For both Cauchy and LaGrangian stress-strain, ANOVA revealed a significant effect of genotype on regression coefficients (quadratic: P<0.002; linear: P<0.0002, resp.). Post hoc comparison of the quadratic coefficients revealed that wild-type fibers were significantly stiffer than desmin-null (P<0.01) and K19-null (P<0.01) fibers (Figure 2(b)). Similarly, comparison of linear coefficients confirmed that wild-type fibers were significantly stiffer than desmin-null (P<0.01) and K19-null (P<0.01) fibers (Figure 2(d)).Mechanical response of fibers from young animals. Raw data (a, c, e) and regression coefficients (b, d, f) for (a-b) Cauchy stress-fiber strain (quadratic regression); (c-d) LaGrangian stress-fiber strain (linear regression); (e-f) radial strain-fiber strain (exponential regression) for fibers from 3-month-old wild-type (blue), desmin-null (red), and K19-null (orange) animals (n=5 fibers/genotype). Different symbols within a color group represent data from individual experiments. Displayed regression lines in (a), (c), and (e) represent averaged coefficients for individual experiments reported in (b), (d), and (f). * indicates significant difference from wild-type. P values are reported in the text.
(a)(b)(c)(d)(e)(f)Analysis of radial strain in fibers from younger animals, in contrast to observations in fibers from older animals, also did not reveal any differences between desmin-null and K19-null genotypes (Figure2(e)). First-order exponential functions were again used to describe the relationship between radial and longitudinal strain (r2>0.78 for all experiments). ANOVA revealed a significant effect of genotype on the exponential decay constant (P<0.0003). Post hoc comparison of the time constants revealed that wild-type fibers exhibited a stronger inverse relationship, and thus tighter coupling, between longitudinal and radial strains than desmin-null (P<0.01) and K19-null (P<0.01) fibers (Figure 1(f)).This finding was consistent with our analysis of Poisson’s ratio. As for fibers from older animals, regression revealed a power-law relationship between longitudinal and radial strains (r2>0.61), and Poisson’s ratios for all genotypes far exceeded 0.5, indicating considerable anisotropy. There were no significant differences in Poisson’s ratio between desmin-null and K19-null genotypes (Mean ± SEM for des −/−: 1.08 ± 0.12; ker −/−: 1.19 ± 0.08), but the Poisson’s ratio of wild-type fibers (2.43 ± 0.33) was significantly higher than those of both desmin-null and K19-null fibers (P<0.02). With respect to changes in volume in response to loading, as for older animals, all genotypes decreased in volume with increased deformation, and wild-type fibers displayed the most dramatic decrease in volume (Mean slope ± SEM: −1.08 ± 0.19). This decrease was significantly greater (P<0.01) than that in both desmin-null and K19-null fibers (des −/−: 0.24 ± 0.05; ker −/−: 0.24 ± 0.02). Significant differences between parameters are summarized in Table 3.Table 3
Summary of parameters that display significant post hoc differences (P<0.05) between genotypes for fibers harvested from young mice. Raw data, means, and standard errors are provided in Figure 2. NS: no significant difference.
Genotypedes −/−K19 −/−+/+σLG, σC, εR, ν, VfσLG, σC, εR, ν, VfDes −/−NSThough we have not yet compared the properties of fibers in 10-month-old to 3-month-old double-null mice, we confirmed differences in passive mechanical parameters between young and aged wild-type, desmin-null, and K19-null mice using two-way ANOVA. Effects of genotype (P<0.0001), but no effect of age and no interaction were found on stress-strain and volume-strain relationships. However, a significant effect of genotype (P<0.0001) and age (P<0.0001) and an interaction between genotype and age (P<0.0001) were observed for the coupling between radial strain and longitudinal strain as indicated by both exponential and power-law relationships. Post hoc analysis revealed a significant reduction in the exponential time constant and Poisson’s ratio of longitudinal and radial strain for wild-type (exp: P<0.008; Poisson’s ratio: P<0.006), desmin-null (exp: P<0.005; Poisson’s ratio: P<0.02), and K19-null (exp: P<0.001; Poisson’s ratio: P<0.001) fibers from aged mice compared to their younger counterparts (Figures 1(e) and 1(f) versus Figures 2(e) and 2(f)).


## 3.3. Structural Bases for Differences in Biomechanical Properties

In an effort to elucidate potential structural bases for differences in passive mechanical function between genotypes, we used phase-contrast microscopy to image the response of fibers from 10-month-old animals of each genotype at higher resolution (Figures3–7). Wild-type fibers displayed well aligned Z-disks, even at high fiber strains (Figure 4). In contrast, desmin-null fibers displayed the characteristic “basket-weave” pattern of misaligned myofibrils (Figure 5). Despite the misalignment of myofibrils in desmin-null mice, local sarcomere strain scaled consistently with fiber strain in both wild-type and desmin-null fibers (Figures 3(a) and 3(b), 4, and 5). Additionally, membrane strain was also tightly correlated with sarcomere strain in both wild-type and desmin-null fibers, indicating continuity between the sarcolemma and myofibrillar lattice (Figures 3(c) and 3(d)).Relationships between global fiber strain and local membrane and sarcomere strain. Raw data (a, c) and regression coefficients (b, d) for (a-b) fiber strain-sarcomere strain; (c-d) membrane strain-sarcomere strain for fibers from 10-month-old wild-type (blue), desmin-null (red), K19-null (orange), and double-null (green) animals (n=3 fibers/genotype). Different symbols within a color group represent data from individual experiments. Displayed regression lines in (a) and (c) represent averaged coefficients for individual experiments reported in (b) and (d).
(a)(b)(c)(d)Figure 4
Geometric and structural response of a 10-month-old wild-type fiber to tensile loading. The same region of the fiber, imaged using phase-contrast microscopy, is shown at five different deformations. Fiber strains, radial strains, and sarcomere lengths are reported to the right of the fiber for each deformation. The fiber displays well aligned Z-lines (dark vertical stripes), consistently increasing sarcomere strain, and decreasing fiber diameter with increased fiber deformation. A sample region of the fiber that displays particularly well aligned Z-lines that are experiencing some shear, inferred through the slanting of Z-lines, is marked with a white box. The boxed region has been computationally enlarged 2x for clarity. At high strains (εF>0.8), artifactual thickening of  Z-lines occurs due to shear-induced offsets of Z-lines in deeper regions of the fiber. Bar = 50 μm.Figure 5
Geometric and structural response of a 10-month-old desmin-null fiber to tensile loading. The same region of the fiber, imaged using phase-contrast microscopy, is shown at five different deformations. Fiber strains, radial strains, and sarcomere lengths are reported to the right of the fiber for each deformation. The fiber displays misaligned Z-lines, but consistently increasing sarcomere strain and decreasing fiber diameter with increased fiber deformation with increased fiber deformation. Misalignment is characterized by a striking basketweave pattern where no clear Z-line orientation may be inferred. A sample region of the fiber that displays this pattern is marked with a white box, and may be compared to the well aligned region shown in Figure4. The boxed region has been computationally enlarged 2x for clarity. Desmin-null fibers display dramatically more debris on their external membrane compared to wild-type fibers. This observation is consistent with the increased “stickiness” experienced during fiber dissection, possibly resulting from differences in extracellular matrix composition. Bar = 50 μm.Geometric and structural response of a 10-month-old K19-null fiber to tensile loading. (a) The same region of the fiber, imaged using phase-contrast microscopy, is shown at five different deformations. Fiber strains, radial strains, and sarcomere lengths are reported to the right of the fiber for each deformation. The fiber displays aligned Z-lines, but discontinuity between fiber deformation and sarcomere strain. This fiber also displays considerable shear, as seen by the slanting of Z-lines (dark lines). A sample region of the fiber that displays particularly well aligned Z-lines that are experiencing some shear is marked with a white box. The boxed region has been computationally enlarged 2x for clarity. There is no slippage or damage of the fiber at the (b) left edge or (c) right edge, imaged in the vicinity of the knot and pins (starred black regions in the periphery of the image) at a moderate fiber strain of 0.5. This may be concluded from the gradual reduction in of the slant of Z-disks away from the knot, which is located to the left in (b) and to the right in (c). Bar = 50μm.
(a)(b)(c)Figure 7
Geometric and structural response of a 10-month-old double-null fiber to tensile loading. The same region of the fiber, imaged using phase-contrast microscopy, is shown at five different deformations. Fiber strains, radial strains, and sarcomere lengths are reported to the right of the fiber for each deformation. Fiber displays both misaligned Z-lines as well as discontinuity between fiber, sarcomere, and radial strains. A sample region of the fiber that displays the basketweave pattern characteristic of misalignment is marked with a white box. The boxed region has been also computationally enlarged 2x for clarity. As for desmin-null fibers, double-null fibers display dramatically more debris on their external membrane compared to wild-type fibers. Bar = 50μm.In contrast, K19-null fibers appeared to have well aligned myofibrils, even in the presence of considerable shear (Figure6). Though differences did not reach statistical significance due to the small sample size, there was considerable discontinuity between fiber strain and sarcomere strain for K19-null fibers compared to wild-type or desmin-null fibers (variability of each genotype in Figures 3(a) and 3(b); examples in Figures 4–6). Examination of the entire length of the fiber at multiple time points during experimentation suggested that this was not simply an artifact of fiber slippage at the knots or visible damage outside of the imaging window (Figures 6(b) and 6(c)). Additionally, similar to findings in desmin-null fibers, the absence of K19 did not appear to disrupt coupling between the sarcolemma and the contractile apparatus significantly, based on the strong correlation between local membrane strain and sarcomere strain (Figures 3(c) and 3(d)).Double-null fibers, as expected, displayed phenotypes from both desmin-null and K19-null fibers. Visually, a lack of Z-disk alignment similar to that in desmin-null fibers was observed (Figure7). As in K19-null fibers, however, there was also considerable discontinuity between fiber strain and sarcomere strain in double-null fibers (Figures 3(a), 3(b), and 7). As for all genotypes, there was no apparent decoupling between the sarcolemma and underlying myofibrils (Figures 3(c) and 3(d)).


## 4. Discussion

We characterized the mechanical response of skeletal muscle fibers from wild-type, desmin-null, K19-null, and desmin/K19 double-null mice to passive tensile loading. Using a mechanical testing device that also enabled high-resolution imaging, we also characterized structural changes in fibers from each genotype during deformation. Our study extends previous work describing the roles of IFs in skeletal muscle morphology and function [15, 17–22, 26] by showing that IFs comprised of desmin and of keratins each contribute to the biomechanical integrity of myofibrils in the sarcoplasm and their links to the sarcolemma. Our work also demonstrates that the weakening of these links in myofibers null for either of these proteins becomes more pronounced in older mice.Experiments on wild-type and desmin-null fibers from 3-month-old and 10-month-old FVB mice confirmed and extended our previous characterization of fibers of these genotypes from 129/SVJ mice [26]. Reported regression coefficients and magnitudes for Cauchy stresses, LaGrangian stresses, and radial strains are consistent with those previously reported, suggesting that the murine strain does not influence passive mechanical properties in skeletal muscle. Differences in elastic moduli between wild-type and desmin-null fibers, reported for both LaGrangian and Cauchy stress-strain curves, indicate a role for desmin in both defining inherent load-bearing capabilities within the muscle fiber and regulating changes in fiber geometry that further influence force transmission (Figures 1(a)–1(d), and 2(a)–2(d)). The latter finding is also consistent with significantly decreased coupling between longitudinal and radial strains and reduced changes in volume following loading in desmin-null fibers compared to wild-type (Figures 1(e), 1(f), 2(e), and 2(f), Poisson’s ratios).Given the misalignment of Z-disks observed in single fibers and in whole muscle (Figure5; [15, 16, 18, 26]), it is likely that a major site of desmin’s influence on force transmission is within the myofibrillar lattice. Within this region, desmin enhances the effectiveness in transmission of shear loads and coordinates tensile load-bearing between myofibrils. However, it remains an open question whether desmin also influences force transmission at the sarcolemma. On one hand, this study reports a strong correlation of membrane strain and sarcomere strain (Figures 3(c) and 3(d)). This would suggest limited influence of desmin at the sarcolemma, consistent with the tight coupling between Z-disks and costameric proteins observed in single desmin-null EDL fibers [26], and intact costameres in quadriceps muscle from desmin-null mice [15]. On the other hand, significant costamere disruption is seen in tibialis anterior and EDL muscle [15], suggesting that desmin does, in fact, influence costamere, and thus membrane, coupling to the underlying myofibrillar lattice.Examination of K19-null fibers from three-month-old FVB mice indicated, for the first time, a role for keratins in regulating the passive mechanical properties of muscle. Similar to desmin-null fibers, K19-null fibers are more compliant than wild-type fibers, based on LaGrangian and Cauchy stress-strain curves (Figures1(a)–1(d), and 2(a)–2(d)). In addition, K19-null fibers from younger mice exhibit poor coupling between longitudinal and radial strain (Figures 2(e) and 2(f), Poisson’s ratios). Interestingly, coupling between longitudinal and radial strain, and thus perhaps the efficiency of lateral force transmission, decreases for all genotypes with aging (Figures 1(e), 1(f), 2(e), and 2(f), Poisson’s ratios). This may be a result of myofibril loss and a consequent increase in intramyofibrillar space [28, 29]. In contrast to desmin-null fibers, though, coupling between longitudinal and radial strain was not significantly decreased in 10-month-old K19-null fibers compared to wild-type fibers (Figures 1(e) and 1(f)), and Z-disks appeared to remain well aligned, suggesting different mechanisms by which K19 and desmin influence passive biomechanical properties. Taken together, these results raise the possibility that genotypic differences in longitudinal and radial coupling may in part be a result of differences in Z-disk connectivity, both longitudinally and laterally [16, 18, 22, 26].Also unlike desmin-null fibers, K19-null fibers reveal a poor correlation between fiber strain and sarcomere strain (Figures3(a), 3(b), 5, and 6). The reason for this lack of correlation remains unclear. Based on disrupted costameres and the dramatic accumulation of mitochondria in some subsarcolemmal regions of K19-null muscle [22], decoupling between the membrane and myofibrillar response to deformation was predicted; however, coupling between the plasma membrane and the contractile apparatus does not appear to be dramatically disrupted, at least over a length scale of 100–200 μm (Figures 3(c) and 3(d)). There are at least two reasonable explanations for this paradox, which are also applicable to understanding desmin-mediated connectivity at the costameres. First, when disrupted costameres have been observed in whole muscle sections [15, 22], there appear to be at least a few costameres in a given fiber that are intact. This limited connectivity may be sufficient to tether the membrane to the underlying cytoskeleton, though not necessarily strongly or effectively. The fraction of intact costameres is likely to be even higher in successfully teased (i.e., more robust) single fibers from IF-null muscles (this study, [26]), further enhancing apparent connectivity. The second explanation arises from the fact that multiple structural proteins, including other intermediate filaments such as K8, localize to costameres and thus likely contribute to the mechanical properties of the connection between individual costameres and underlying Z-disks (cf. [30, 31]). In this case, costameres may still be weakly coupled to myofibrils via remaining or compensatory structural proteins. Examination of more subtle differences in connectivity at the scale of individual or a few sarcomeres under conditions of radial tensile loading (e.g., [30]) may reveal differences in local compliance that could then underlie differences in the efficiency of lateral force transmission.Additional insight into the different roles of desmin and keratin emerged from the examination of desmin/K19 double-null fibers in 10-month-old mice. The similar fiber compliance of double-null fibers and desmin-null fibers, but increased compliance compared to K19-null fibers (Figures1(a)–1(d)) parallels decreased myofibrillar organization within desmin-null and double-null fibers compared to K19-null fibers (Figures 5–7). Our findings in single myofibers are also consistent with biomechanical results in whole muscle, where deficits in specific tension are similar and more severe in desmin-null and double-null tibialis anterior muscle than in K19-null muscle [18]. These data support the hypothesis that myofibrils must be engaged and coordinated to maximally support mechanical loads. It will be interesting to continue the characterization of double-null fibers in younger mice, where phenotypic differences between desmin-null and K19-null fibers were not significant, to address developmental differences in IF function.Collectively, integration of mechanical and structural data from this work and previous studies support a model in which membrane-myofibril connectivity, guided by two sets of IFs, composed of keratins and of desmin, dictates the radial transmission of tensile and shear loads into and out of the cell. In an outside-in model, this connectivity would enable the recruitment of longitudinal load-bearing structures, predominantly coordinated by desmin, within the myofibrillar lattice. Conversely, in an inside-out model, connectivity between myofibrils and the membrane would enable the reliable transmission of myofibrillar forces, again coordinated by desmin, to the extracellular matrix.Finally, this work illustrates the importance of studying single muscle fibers to the interpretation of biomechanical data from knockout models. Increased passive compliance measured in desmin-null fibers (Figures1 and 2; [26]) is contrary to the increased passive stiffness reported in desmin-null muscles [17, 32]. These differences are most likely due to secondary remodeling of extracellular matrix, though this remains to be tested directly. Given the increasing relevance of extracellular matrix remodeling to our understanding of myopathies and muscular dystrophies [23–25], comparative studies between single fibers and whole muscle or fiber bundles may reveal interesting and important details regarding roles of intermediate filaments in initiating physiologically relevant extracellular signaling pathways.


## 5. Conclusions

We have characterized the structural and mechanical response of single fibers from wild-type, desmin-null, K19-null, and desmin/K19 double-null mice to passive tensile loading. Through this work, we provide the first direct evidence, in an experimental system free from influences of the extracellular matrix, that desmin and keratin 19 differentially influence passive force transmission in single muscle fibers. Our work has strong implications for our understanding of the role of intermediate filaments on structure-function relationships in skeletal muscle. In addition, our work is important for interpreting structural and biomechanical changes resulting from desminopathies and other myopathies linked to intermediate filaments and their associated proteins [33–37].



---

*Source: 704061-2012-01-12.xml*