Usha Sinha1
Vadim Malis2
Jiun-Shyan Chen3
Robert Csapo4
Ryuta Kinugasa5,6
Marco Vincenzo Narici7
Shantanu Sinha8*- 1Department of Physics, San Diego State University, San Diego, CA, United States
- 2Department of Physics, University of California, San Diego, San Diego, CA, United States
- 3Department of Structural Engineering, University of California, San Diego, San Diego, CA, United States
- 4Research Unit for Orthopaediic Sports Medicine and Injury Prevention, ISAG, Private University for Health Sciences, Medical Informatics and Technology, Hall in Tirol, Austria
- 5Department of Human Sciences, Kanagawa University, Yokohama, Japan
- 6Computational Engineering Applications Unit, Advanced Center for Computing and Communication, RIKEN, Saitama, Japan
- 7Institute of Physiology, Department of Biomedical Sciences, University of Padua, Padua, Italy
- 8Department of Radiology, University of California, San Diego, San Diego, CA, United States
The focus of this review is the application of advanced MRI to study the effect of aging and disuse related remodeling of the extracellular matrix (ECM) on force transmission in the human musculoskeletal system. Structural MRI includes (i) ultra-low echo times (UTE) maps to visualize and quantify the connective tissue, (ii) diffusion tensor imaging (DTI) modeling to estimate changes in muscle and ECM microstructure, and (iii) magnetization transfer contrast imaging to quantify the macromolecular fraction in muscle. Functional MRI includes dynamic acquisitions during contraction cycles enabling computation of the strain tensor to monitor muscle deformation. Further, shear strain extracted from the strain tensor may be a potential surrogate marker of lateral transmission of force. Biochemical and histological analysis of muscle biopsy samples can provide “gold-standard” validation of some of the MR findings. The review summarizes biochemical studies of ECM adaptations with age and with disuse. A brief summary of animal models is included as they provide experimental confirmation of longitudinal and lateral force transmission pathways. Computational muscle models enable exploration of force generation and force pathways and elucidate the link between structural adaptations and functional consequences. MR image findings integrated in a computational model can explain and predict subject specific functional changes to structural adaptations. Future work includes development and validation of MRI biomarkers using biochemical analysis of muscle tissue as a reference standard and potential translation of the imaging markers to the clinic to noninvasively monitor musculoskeletal disease conditions and changes consequent to rehabilitative interventions.
Introduction
The loss of muscle mass with aging and with unloading (e.g., disuse) has been studied extensively. Accompanying the loss of muscle mass (termed sarcopenia) is a disproportionately greater loss of muscle strength as we age, termed dynapenia (Goodpaster et al., 2006). The role of muscular and neural determinants of muscle force loss with age has been investigated in several human and animal studies (Morse et al., 2005; Faulkner et al., 2007; Degens et al., 2009; Aagaard et al., 2010; Ballak et al., 2014). In addition to these determinants, the role of the extracellular matrix (ECM) in muscle force loss is also being increasingly recognized (Ramaswamy et al., 2011; Zhang and Gao, 2014).
The ECM of the muscle is arranged hierarchically with the endomysium around the muscle fiber, the perimysium surrounding a group of muscle fibers, and the epimysium around the whole muscle. Scanning electron microscopy (SEM) images of the connective tissue reveal that the endomysium is a highly ordered network that surrounds individual muscle fibers (Purslow and Trotter, 1994) while the perimysium includes extended cables organized as bundles of collagen fibers that terminates on muscle cells (Gillies and Lieber, 2011). Collagen is the main structural protein in ECM and may be present up to 10% by weight. Older muscles typically demonstrate fibrotic morphology (Lieber and Ward, 2013) that is characterized by an increase in collagen content and non-enzymatic cross-linking of collagen fibers (Haus et al., 2007), and results in a stiffening of the ECM (Wood et al., 2014). Measurements of the passive mechanical properties of the ECM reveal that it is inherently stiffer than muscle fibers and that the aging ECM has a higher modulus than the young ECM; the latter was attributed to an increase of densely packed extensively cross-linked collagen (Meyer and Lieber, 2011). Mechanical properties are measured on dissected, skinned single fibers placed in an activating solution using a device to measure force and length (Wood et al., 2014).
Aging and unloading significantly alter the concentration of force transfer proteins (e.g., decrease in dystrophin protein with aging) that transmit force both inside the muscle fibers (Hughes et al., 2017) and in the ECM (Kragstrup et al., 2011), changing their content, orientation, and composition. Age-related ECM remodeling process is accompanied by changes in the biochemical composition of the ECM, which features a higher ratio of type I to type III collagen fibers in older individuals (Hindle et al., 2009). In addition, augmented concentrations of type IV collagens but lower laminin contents have been observed in the basal lamina of slow twitch muscles of rat soleus and rectus femoris muscles (Kovanen et al., 1988). Transcriptional profiling studies performed in rat soleus muscles have shown that a large proportion of genes exerting a biological role in the ECM and cell adhesion is decreased in older muscles (Pattison et al., 2003). This observation suggests that the age-associated fibrogenic processes are driven by a decreased degradation capacity, rather than by increased formation of collagenous structures. In particular, the activity of matrix metalloproteinases, i.e., enzymes responsible for collagen degradation were reported to decrease at older age in rat gastrocnemius, digitorum profundus, and soleus muscles (de Sousa Neto et al., 2018). The interplay between age-associated ECM remodeling and the development of fibrosis in mouse gastrocnemius muscle is the subject of current investigations (Stearns-Reider et al., 2017; Forcina et al., 2019). As compared to the increasing number of studies investigating age-associated ECM remodeling, scant evidence exists on changes promoted by prolonged disuse. Early studies performed by Karpakka et al. (1990, 1991) showed that immobilization of rats led to reduced hydroxylase activity and hydroxyproline (an amino acid constituting collagens) in the soleus and tibialis anterior muscles. More recent studies, similarly performed in rodents, suggest that a disuse-induced increase in the activity of matrix metalloproteinase 9 (gelatinase B) may be responsible for decreased expressions of collagens I and IV (Giannelli et al., 2005; Liu et al., 2010). Changes in collagen content (in rat gastrocnemius, soleus, and tibialis anterior as well as in human vastus lateralis muscles) in response to short-term immobilization are usually small (Savolainen et al., 1988; Haus et al., 2007), which may be explained by the relatively slow turnover rate of collagenous structures. Prolonged disuse, however, may cause a proteolytic imbalance that is responsible for the breakdown of basal lamina structures. Such ECM remodeling may lead to altered permeability and, ultimately, promote disuse atrophy. It is important to note that aging and disuse may induce changes in the structure and biochemical composition of the ECM, which are expected to be associated with impaired transmission of contractile force through lateral pathways [i.e., through shearing between adjacent muscle fibers across costameres and the network of intramuscular connective tissue (IMCT)] (Csapo et al., 2020).
The current accepted structural models of muscle tissue include muscle fibers embedded in a matrix of connective tissue ECM. Force generated in the muscle fiber is transmitted to the ECM at multiple focal adhesions and subsequently to the tendon via the ECM (Huijing, 2003; Purslow, 2010). Myofascial pathways of force have been observed at various levels including between adjacent fibers (Huijing, 1999), multi-tendon muscle (Zhang and Gao, 2014), and muscle synergists (Bernabei et al., 2015). In multi-tendon muscles, force pathways exist through the myotendinous junction as well as through the ECM of the muscle bellies. In the multi-tendon extensor digitorum longus (EDL), it was shown that tenotomy severing the myotendinous junction had a smaller effect on the measured force than myotomy (severing the ECM between muscle heads) (Zhang and Gao, 2014). This effect was more pronounced in younger than aging rats, implying that lateral transmission of force (LTF) is reduced in the aging muscle. Architectural gear ratio (AGR, defined as whole muscle velocity to muscle fiber velocity), hypothesized to arise from the interplay of contractile forces and connective tissue constraints, has been identified as one of the determinants of mechanical performance in pennate muscles (Brainerd and Azizi, 2005, Eng et al., 2018). Pennate muscles can achieve a gear ratio greater than one since they can shorten through a combination of fiber shortening and rotation. Further, the gearing is variable across contractions with muscles operating at high gear during low force contractions and at low gear during high force contractions (Eng et al., 2018). Variable gearing seen in younger cohorts is lost in the muscle of aging rats presumably from the disruption in the interaction between contractile and ECM tissue leading to the decline in muscle performance (Azizi et al., 2008; Holt et al., 2016; Dick and Wakeling, 2017). AGR has been determined from in situ studies with sonomicrometry crystals embedded along with fibers to measure fiber length changes (Azizi et al., 2008), and from ultrasound (US) and MRI on human subjects (Shin et al., 2009b; Dick and Wakeling, 2017). Studies also show that the restriction of radial expansion by connective tissues results in decreased muscle fiber shortening and generation of less mechanical work (Azizi et al., 2017). Further, a stiff aging ECM may limit muscle performance possibly from both a smaller longitudinal contraction as well as a decrease in myofascial force transmission pathways via LTF (Ramaswamy et al., 2011, Zhang and Gao, 2014). In young rodent muscles, at least 80% of force transmission occurs laterally through the ECM (Huijing et al., 1998), whereas in frail old muscle, LTF through the ECM is reduced ~60% (Ramaswamy et al., 2011; Zhang and Gao, 2014).
The above summary clearly highlights that most of the tools used to study ECM are invasive, limiting comprehensive investigations of the contribution of the ECM to dynapenia in humans. This review focuses on structural and functional MRI to map ECM adaptations and consequence on muscle performance in healthy, aging, and disused muscle. Subject-specific anatomical (muscle compartments and morphology), voxel level compositional (fat, connective tissue, and muscle fractions), and fiber morphology (pennation angle and fiber length) models using MRI data will more accurately predict force output as well as force transmission pathways. Modeling-derived force transmission pathways can be used to validate MR derived biomarkers currently hypothesized to be surrogate markers of lateral transmission of force. It should be noted that the proposed MR indices are novel and validation with biopsy analysis, animal models, or computational models are either ongoing or need to be performed in future studies.
Structural Imaging of the Ecm
The distribution and content of the ECM in humans have so far proven challenging to investigate with in vivo imaging, since its physical properties (very low spin-spin relaxation time, T2) renders it “invisible” in routine MRI. Consequently, there are few imaging-based studies published to date that have reported quantitative data of the ECM/fibrosis within skeletal muscles. Imaging these short T2 tissues requires special pulse sequences with ultra-low echo times (UTEs).
Ultra-Low TE Imaging of Short T2 Intra-Muscular Connective Tissue
UTE imaging has been implemented by several groups and primarily used to study cortical bone, tendons, and cartilage (Sinha et al., 2007; Du et al., 2007a,b). Very few studies have investigated UTEs for muscle connective tissue quantification (Csapo et al., 2014; Ugarte et al., 2017a). Since adipose tissue is also present intramuscularly and a voxel can contain fat, water, and connective tissue fractions, and as these fractions change with age, it is important to be able to quantify all three tissues in muscle and investigate their variations with age. Fat quantification is an area of intense research and several recent excellent reviews on skeletal muscle adipose quantification in normal and diseased muscle are available (Grimm et al., 2018; Hu et al., 2019; Marty and Carlier, 2019). The IDEAL-IQ imaging sequence that accounts for the multiple fat resonance peaks has been validated for accurate fat fraction (Eskreis-Winkler et al., 2018).
Intramuscular adipose tissue (IMAT) and IMCT in the triceps surae (TS) muscles have been quantified in a cohort of young and aging subjects, using quantitative fat/water and UTE sequences (Csapo et al., 2014). The latter study found an 18% greater TS muscle volume and a 39% higher muscle force in the younger cohort compared to the aging cohort (differences calculated from the average of young and aging subjects normalized to the average of all subjects), implying that 45% of force loss was not due to muscle size (ratio of volume to force differences). Further, in the gastrocnemius muscles (computed together for lateral and medial), IMCT increased by 64% and IMAT increased by 40% in the aging cohort, reflecting a 36% increase in the total amount of noncontractile tissue in the aging muscle (Figures 1A–C). A multichannel fuzzy cluster mean algorithm was implemented to automatically label IMAT/IMCT from UTEs images of calf images (Figure 1A; Ugarte et al., 2017a,b). Recent advances in image analysis enable accurate quantification of the CT, adipose, and water fraction in muscle from UTEs images (Araujo et al., 2017). US imaging also provides information on fat infiltration and fibrosis but is largely limited to 2D imaging and is semi-quantitative; however, US and MRI may provide complimentary information (Nijboer-Oosterveld et al., 2011; Arts et al., 2012; Mul et al., 2018).
Figure 1. 3D muscle volumes with the noncontractile tissue shown segmented (automated MCFCM algorithm) in green for a young (left) and aging (right) subject (A) reproduced with permission from Wiley publisher of Ugarte et al. (2017a). Box plots reproduced with permission from Springer publisher of Csapo et al. (2014) showing adipose (B) and connective tissue (C) content in young and aging subjects in the calf muscles (MG, medial gastrocnemius; LG, lateral gastrocnemius; SOL, soleus). Significant differences are marked by an asterisk. Baseline and eigenvalue (D|| and D⊥) images shown at two diffusion times (short: 30 ms and long: 700 ms) for young and aging subjects (D). Tertiary eigenvalues decrease with diffusion time and are significantly different between young and aging subjects at longer diffusion times indicating that at these longer times when the average diffusion path length approaches fiber cell dimensions, the tertiary eigenvalue representing diffusion in the fiber cross-section is more sensitive to the fiber dimensions. The significant difference may potentially indicate that aging muscle fiber has a smaller diameter than young muscle fiber (D) reproduced with permission from ISMRM publisher of Malis et al. (2019a). Microstructural parameters of the model (D0, free diffusion coefficient; S/V, surface to volume ratio of fibers; a, muscle fiber diameter; and τD, residence time, which is the average time spent within a muscle fiber) for young and aging subjects (D). Histological image of muscle tissue (top), segmented cross-section (middle), and 3D model produced by extruding the 2D model along the z-direction (blue: muscle fiber, red: connective tissue) (E) reproduced with permission from Wiley publisher of Zhang et al. (2019). The micromechanical model prediction of force generated by three muscle types (young, adult, and old) is shown for isometric, concentric, and eccentric contractions. Young, adult, and old muscle tissue were modeled with passive stiffness as low, medium, and high and connective tissue thickness as thin, medium thick, and thick, respectively; λ is the stretch ratio (F) reproduced with permission from Wiley publisher of Zhang et al. (2019).
Magnetization Transfer Contrast Imaging
Quantification of the macromolecular fraction (such as collagen in connective tissue) can also be performed by quantitative magnetization transfer (qMT) imaging (Ramani et al., 2002). In the qMT method, macromolecules are selectively saturated with an off-resonance pulse and the effect of this saturation is seen as a decrease in signal intensity of the free protons. qMT models the signal intensity as a function of the magnetization transfer (MT) pulse offset and power to yield several parameters including the macromolecular fraction (Ramani et al., 2002). Recently, qMT has been applied to skeletal muscle and collagen has been proposed as the most likely macromolecule measured though it has not been validated with biopsy analysis (Sinclair et al., 2010; Li et al., 2014). A recent animal study to monitor muscle tissue formation found correlations between magnetization transfer ratio (MTR) with specific skeletal markers and to contractility; this study indicated the potential of MTR to serve as a surrogate marker of muscle fiber formation (rather than collagen) (Rottmar et al., 2016). Age-related changes in MTsat, a semi-quantitative index of MT showed that MTsat decreases with age (White, 2019). The latter finding is surprising since if the measured macromolecule fraction reflects collagen, and collagen is known to increase with age, it would be anticipated that the MTsat would increase with age (Alnaqeeb et al., 1984). These studies emphasize that validation based on biochemical analysis of biopsy samples is critical to identifying the macromolecules detected by qMT, MTsat, or UTEs techniques and is currently underway in the author’s laboratories.
Diffusion Tensor Imaging and Modeling
Diffusion tensor imaging (DTI) is an established technique to probe the architecture of anisotropic structures such as skeletal and cardiac muscle fibers. Several excellent reviews discuss the imaging techniques, diffusion eigenvalues, and fractional anisotropy, and muscle fiber tracking including fiber length and pennation angles (Froeling et al., 2013; Oudeman et al., 2016). Diffusion eigenvalues reflect diffusion along the muscle fiber axis (primary eigenvalue) while secondary and tertiary eigenvalues maybe proportional to the fiber diameter (Galbán et al., 2004, 2005; Karampinos et al., 2007). However, the resolution of DTI-MRI precludes direct observations at the microscopic scale and modeling is required to link tissue microarchitecture to DTI indices.
Models for diffusion in muscle tissue were only proposed recently (Karampinos et al., 2009; Fieremans et al., 2017). A two-compartmental model, with elliptical muscle fibers as first compartment embedded in the second compartment, the ECM has been proposed (Karampinos et al., 2009). Model derived muscle fiber diameters of calf muscles (Karampinos et al., 2009) were in the range of muscle fiber diameters obtained from prior histological studies of excised muscle (Behan et al., 2002). Another model (Fieremans et al., 2017), termed as the random permeable barrier model (RPBM) fits the time dependence of D⊥ (average of the secondary and tertiary eigenvalues) to extract microstructural parameters. The time dependence of D⊥ is obtained from a stimulated echo sequence at diffusion times of the order of several hundred milliseconds (Figure 1D). Estimates of the fiber diameter in the calf muscle using RPBM were comparable to autopsy reports from literature (Polgar et al., 1973) though DTI modeling underestimated fiber diameters (Sigmund et al., 2014).
The two-compartment model applied to DTI data pre‐ and post-unloading showed that muscle fiber, permeability, and intracellular volume decreased and collagen increased while muscle ellipticity decreased post-suspension; the latter could potentially be explained by the lack of a mechanical stimulus during unloading (Malis et al., 2018b). The RPBM was applied to time dependent DTI data on old and young cohorts (Malis et al., 2019a). Significant differences were found in the tertiary eigenvalue at the longest diffusion times; the latter is most sensitive to the fiber diameter (Figure 1D). The RPBM model yielded a lower volume fraction for the aging muscle; the latter is a measure of the membrane’s ability to hinder diffusion. It is possible that aging muscle may have compromised sarcolemma integrity that makes it more permeable. While the DTI modeling with age and disuse are preliminary analysis, they confirm that modeling enables one to monitor fiber and ECM adaptations at the microstructural level (Malis et al., 2018b, 2019a).
Strain and Strain Rate Mapping Using Dynamic Mri
Muscle tissue deformation is quantified using MRI sequences that map voxel displacement (displacement encoded imaging) or voxel velocity [velocity encoded phase contrast imaging, (VE-PC)] or MR tagging to track voxel displacements [spatial modulation of magnetization, (SPAMM)] (Tan et al., 2018). Strain heterogeneity along the muscle axis as well as the muscle fiber has been reported in the calf muscle and in the biceps brachii (Zhong et al., 2008; Shin et al., 2009b). Using the SPAMM technique, 3D strain and diffusion tensors were evaluated in the anterior tibialis and revealed a planar strain pattern where the principal shortening direction deviated from the muscle fiber direction indicating the presence of shear strain in the tissue (Englund et al., 2011). The deviation was attributed to the variations in fiber length leading to a spatial variability in the absolute length change that in turn results in a deviation of the average strain away from the fiber direction toward the longest fibers; the inter-fiber connections are mediated by the ECM (Englund et al., 2011). Heterogeneous length changes along the fascicles of human medial gastrocnemius muscle measured by combining magnetic resonance imaging deformation analyses and DTI were reported (Pamuk et al., 2016; Karakuzu et al., 2017). Computational models reveal that the strain heterogeneity arises primarily from a variation in fascicle lengths and the curvature of the fascicles. Nonuniform strain is hypothesized to affect force production and length range of force exertion (Yucesoy et al., 2002) and is an important determinant for the outcome of orthopedic surgery (Yucesoy et al., 2006; Yucesoy, 2010).
Recent studies focus on the extraction of strain rate (SR) indices from VE-PC images and the changes of these indices with age and with disuse (Sinha et al., 2015a, 2018; Malis et al., 2018a). The SR tensor was calculated from the spatial gradient of the velocity images (acquired during isometric contraction) and diagonalized to obtain the eigenvalues (SRfiber and SRin-plane) and eigenvectors. Figure 2A is a schematic of muscle fiber and ECM deformation during contraction specifying the SR components derived in the principal axis (SRfiber and SRin-plane) and in the muscle fiber axis (SRff and SRcc) reference frames. Maximum shear strain (SRmax) was calculated by rotating the SR tensor by 45° from the principal axis frame and the potential origin of the shear strain is schematically shown in Figure 2A (Sinha et al., 2015b). It was hypothesized that the MR measured shear strain arises from the shear of the endomysium as the muscle fiber contracts (Figure 2A). This hypothesis is in conformance with computational models (Gao et al., 2008a; Sharafi and Blemker, 2011). SRfiber, SRin-plane, and SRmax were significantly lower in the aging cohort during isometric contraction (Figures 2B,D; Sinha et al., 2015a). On multiple variable regressions, maximum shear strain rate (SRmax) and normal strain rate in the fiber cross-section (SRcc) were significantly associated with force (young and old cohorts) (Sinha et al., 2015a). In disuse atrophy, SRin-plane and SRfiber angle were significantly changed with unloading (Figures 2C,E). Changes in SR-fiber angle, SRin-plane, and shear SR as well as their ability to predict force and force changes may reflect the role of remodeling of the ECM in disuse atrophy and its functional consequence in reducing LTF (Malis et al., 2018a).
Figure 2. Schematic of a muscle fiber and endomysium with the principal basis and fiber basis defined with respect to the muscle fiber. The muscle fiber is shown contracting (green arrows) while the endomysium experiences a shear strain (end attached to the muscle contracts with the muscle fiber while the other end stays fixed resulting in a shear strain). The direction of the strain rate components in the two basis sets (SRfiber and SRin-plane in the principal basis and SRff and SRcc in the muscle fiber basis) as well as the angle between these two basis sets, SR-fiber angle are shown. The thick (unfilled) arrows show the lateral transmission of force pathways: from the cross-bridges of the myofilaments through the protein network on the sarcolemma to the extracellular matrix. Arrows indicate normal and shear strains (A) reproduced with permission from Wiley publisher of Sinha et al. (2018). Strain rate maps of the normal strain rates (SRfiber and SRin-plane) and shear strain rates (SRmax) in a young and aging subject at the same force level. The color maps are color-coded according to the legend provided with the figure. Color maps are overlaid on the magnitude image at the corresponding frame; maps show lower values in both the normal strain rates in the aging cohort (B) reproduced with permission from Wiley publisher of Sinha et al. (2018). Streamlines of eigenvectors corresponding to SRfiber overlaid on fiber directions (in black, tracked from fascicle orientations) for one slice of the same subject pre‐ and post-suspension at the peak of the contraction phase. This shows the difference in angle enclosed by SRfiber and muscle fiber (SR-fiber angle) (C) reproduced with permission from Wiley publisher of Malis et al. (2018a). Box plots of SR indices evaluated in the MG that showed significant difference between young and aging subjects (D). Box plots of SR-fiber angle showing significant differences in all three regions of the MG on unloading (E). Significant differences in (D) and in (E) are denoted by ‘*’.
Computational models have identified that endomysium shear is the mechanism by which force is transmitted laterally (Gao et al., 2008b; Sharafi and Blemker, 2010). The experimentally measured MR shear strain was potentially attributed to the shear in the endomysium (Sinha et al., 2015a; Malis et al., 2018a). Further, shear SR was a significant predictor of force in the age and disuse studies highlighting that if shear strain is a measure of LTF, these results indicate the decreasing magnitude of LTF with age and with disuse (Sinha et al., 2015a, 2018; Malis et al., 2018a). The studies highlight that, while a direct noninvasive measurement of LTF is not possible, shear strain rate may potentially serve as a surrogate measure of LTF (Sinha et al., 2015a; Malis et al., 2018a). The increased stiffness of the aging muscle may result in lower shear in the ECM and this is confirmed in the experimental studies (Sinha et al., 2018). Recent advances in fast sequences integrating compressed sensing with VE-PC acquisitions have enabled the study of the full 3D strain tensor at different %MVC at significantly lesser number of contraction cycles (Malis et al., 2019b, 2020). It is anticipated that the ratio of shear to normal strains will decrease with increasing %MVC for aging subjects (compared to younger subjects) reflecting the decreased contribution of LTF in aging muscle.
US strain imaging based on speckle tracking has been used to measure strain and strain rate during skeletal muscle dynamic imaging as well as muscle stiffness (Lopata et al., 2010; Frich et al., 2019). High temporal resolution and real time capabilities are advantages while restriction to 2D-planes and limited field of view are disadvantages of US compared to dynamic MRI.
Computational Modeling
Subject specific computational models using imaging data will more accurately predict force output as well as force transmission pathways (Blemker et al., 2007). MRI can provide anatomical (3D morphology of muscles including size and shape), voxel level compositional (fat, connective tissue, and muscle fractions), and fiber morphology (pennation angle, fiber length, and curvature) information enabling the generation of subject specific models.
The computational modeling of skeletal muscles has advanced from 1D lumped-parameter model (Bahler, 1968; McMahon, 1984; Zajac, 1989), 3D sophisticated macro-scale muscle models (Martins et al., 1998; Johansson et al., 2000; Lemos et al., 2004; Böl and Reese, 2008; Chi et al., 2010), to multi-scale models (Sharafi and Blemker, 2010; Virgilio et al., 2015; Spyrou et al., 2017, 2019; Zhang et al., 2019). Force transmission between the matrix and fibers using an interconnected fiber-matrix mesh discretization has been investigated (Yucesoy et al., 2002, 2006). Imaging advancements now allow measurements of the deformations of contracting skeletal muscles (Drace and Pelc, 1994; Maganaris et al., 1998; Pappas et al., 2002; Finni et al., 2003; Kinugasa et al., 2008; Shin et al., 2009a), and they also offer model parameters for computational investigation of deformation variations in muscle fibers and aponeuroses (van der Linden et al., 1998; Johansson et al., 2000; Yucesoy et al., 2002; Oomens et al., 2003; Blemker and Delp, 2005; Chi et al., 2010). Positive to negative strain variations in aponeuroses, nonuniform strain changes in sarcomeres, and asymmetrical muscle fiber cross-sectional deformations, have also been reported in a nonlinear FE muscle model (Chi et al., 2010). Several specially designed FE models have demonstrated how muscle force output can be strongly influenced by passive material properties (Chi et al., 2012). The reproducing kernel particle method (RKPM) (Chen et al., 1996, 2002) allows the model to be constructed directly from MR image voxels avoiding the complexity of mesh-based models and has been applied to muscle modeling (Basava et al., 2014).
A micromechanical model based on the concept of repeating unit cells extracted from histological cross-sections of rabbit muscles was proposed to study the effects of microstructural geometry in both fiber and fascicle levels on the macroscopic along-fiber shear modulus of skeletal muscle in normal and disease conditions (Sharafi and Blemker, 2010; Virgilio et al., 2015). Recently, a micromechanical computational model (Figure 1E) demonstrated that the increased stiffness and thickness of the passive materials of aging muscle tissue reduces the force generation capability under concentric contraction while maintaining the force generation capability under eccentric contraction (Figure 1F; Zhang et al., 2019).
While it has been acknowledged that the wealth of data from MRI can be leveraged to create subject specific models (Blemker et al., 2007), this has not been fully exploited. Subject specific models created for MRI data (including fiber direction from DTI) on a young and old subject predicted a force decrease with age that matched the experimental data (Basava et al., 2014). Further, the availability of functional MRI (strain patterns) provides experimental data for validation and fine tuning of computational models. It is also anticipated that subject specific endomysial thickness (MRI diffusion modeling) and stiffness (dynamic MRI) can be incorporated into micromechanical models (Zhang et al., 2019).
Summary and Future Prospects
Advanced MRI techniques provide a number of indices to probe the structure and function of muscle fiber and the ECM. In order to establish these indices as imaging markers, MRI derived structural/functional indices await validation from biopsy studies. Structural MRI data can be used to generate subject specific computational models allowing predictions of force and strain distributions; the models in turn can be validated using the functional MR data. Future simulation studies using computational models may provide the causal link between shear in the ECM, shear strain from MRI, and lateral force transmission pathways allowing one to establish shear strain from MRI as a surrogate marker of LTF. In addition, a direct way to establish the link between in vivo MR strain indices and LTF will require animal model studies. It is acknowledged that MRI is only one of many tools including US imaging and biopsy analysis in the arsenal of tools for the comprehensive assessment of the ECM. The focus is to identify and validate MRI biomarkers, which then can be combined with markers from complementary modalities to advance the diagnosis, tracking of natural progression, and response to therapy in different musculoskeletal conditions ranging from sarcopenia, disuse atrophy to dystrophies.
Author Contributions
US was the primary writer of the draft. VM, J-SC, RC, RK, MN, and SS contributed to sub-sections of the draft. US, VM, and SS contributed to the editing and proof-reading. VM produced the figures.
Funding
National Institute on Aging. Grant Number: R01AG056999.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Aagaard, P., Suetta, C., Caserotti, P., Magnusson, S. P., and Kjaer, M. (2010). Role of the nervous system in sarcopenia and muscle atrophy with aging: strength training as a countermeasure. Scand. J. Med. Sci. Sports 20, 49–64. doi: 10.1111/j.1600-0838.2009.01084.x
Alnaqeeb, M. A., Al Zaid, N. S., and Goldspink, G. (1984). Connective tissue changes and physical properties of developing and ageing skeletal muscle. J. Anat. 139, 677–689.
Araujo, E. C. A., Azzabou, N., Vignaud, A., Guillot, G., and Carlier, P. G. (2017). Quantitative ultrashort TE imaging of the short-T2 components in skeletal muscle using an extended echo-subtraction method. Magn. Reson. Med. 78, 997–1008. doi: 10.1002/mrm.26489
Arts, I. M. P., Schelhaas, H. J., Verrijp, K. C. N., Zwarts, M. J., Overeem, S., van der Laak, J. A. W. M., et al. (2012). Intramuscular fibrous tissue determines muscle echo intensity in amyotrophic lateral sclerosis. Muscle Nerve 45, 449–450. doi: 10.1002/mus.22254
Azizi, E., Brainerd, E. L., and Roberts, T. J. (2008). Variable gearing in pennate muscles. Proc. Natl. Acad. Sci. U. S. A. 105, 1745–1750. doi: 10.1073/pnas.0709212105
Azizi, E., Deslauriers, A. R., Holt, N. C., and Eaton, C. E. (2017). Resistance to radial expansion limits muscle strain and work. Biomech. Model. Mechanobiol. 16, 1633–1643. doi: 10.1007/s10237-017-0909-3
Bahler, A. S. (1968). Modeling of mammalian skeletal muscle. I.E.E.E. Trans. Biomed. Eng. 15, 249–257. doi: 10.1109/tbme.1968.4502575
Ballak, S. B., Degens, H., de Haan, A., and Jaspers, R. T. (2014). Aging related changes in determinants of muscle force generating capacity: a comparison of muscle aging in men and male rodents. Ageing Res. Rev. 14, 43–55. doi: 10.1016/j.arr.2014.01.005
Basava, R. R., Chen, J. -S., Zhang, Y., Sinha, S., Sinha, U., Hodgson, J., et al. (2014). “Pixel based meshfree modeling of skeletal muscles.” in Medical image computing and computer-assisted intervention—MICCAI’98 lecture notes in computer science. eds. W. M. Wells, A. Colchester, and S. Delp (Cham: Springer International Publishing), 316–327.
Behan, W. M. H., Cossar, D. W., Madden, H. A., and McKay, I. C. (2002). Validation of a simple, rapid, and economical technique for distinguishing type 1 and 2 fibres in fixed and frozen skeletal muscle. J. Clin. Pathol. 55, 375–380. doi: 10.1136/jcp.55.5.375
Bernabei, M., van Dieën, J. H., Baan, G. C., and Maas, H. (2015). Significant mechanical interactions at physiological lengths and relative positions of rat plantar flexors. J. Appl. Physiol. 118, 427–436. doi: 10.1152/japplphysiol.00703.2014
Blemker, S. S., Asakawa, D. S., Gold, G. E., and Delp, S. L. (2007). Image-based musculoskeletal modeling: applications, advances, and future opportunities. J. Magn. Reson. Imaging 25, 441–451. doi: 10.1002/jmri.20805
Blemker, S. S., and Delp, S. L. (2005). A 3D model of muscle reveals the causes of nonuniform strains in the biceps brachii. J. Biomech. 38, 657–665. doi: 10.1016/j.jbiomech.2004.04.009
Böl, M., and Reese, S. (2008). Micromechanical modeling of skeletal muscles based on the finite element method. Comput. Methods Biomech. Biomed. Engin. 11, 489–504. doi: 10.1080/10255840701771750
Brainerd, E. L., and Azizi, E. (2005). Muscle fiber angle segment bulging and architectural gear ratio in segmented musculature. J. Exp. Biol. 208, 3249–3261. doi: 10.1242/jeb.01770
Chen, J. -S., Pan, C., Wu, C. T., and Liu, W. K. (1996). Reproducing kernel particle methods for large deformation analysis of non-linear structures. Comput. Methods Appl. Mech. Eng. 139, 195–227. doi: 10.1016/S0045-7825(96)01083-3
Chen, J. -S., Yoon, S., and Wu, C. T. (2002). Non-linear version of stabilized conforming nodal integration for galerkin mesh-free methods. Int. J. Numer. Methods Eng. 53, 2587–2615. doi: 10.1002/nme.338
Chi, S. -W., Chen, J. -S., Edgerton, V. R., and Sinha, S. (2012). Finite element modeling of passive material influence on the deformation and force output of skeletal muscle. J. Mech. Behav. Biomed. Mater. 9, 163–183. doi: 10.1016/j.jmbbm.2012.01.010
Chi, S. -W., Hodgson, J., Chen, J. -S., Reggie Edgerton, V., Shin, D. D., Roiz, R. A., et al. (2010). Finite element modeling reveals complex strain mechanics in the aponeuroses of contracting skeletal muscle. J. Biomech. 43, 1243–1250. doi: 10.1016/j.jbiomech.2010.01.005
Csapo, R., Gumpenberger, M., and Wessner, B. (2020). Skeletal muscle extracellular matrix—what do we know about its composition, regulation, and physiological roles? A narrative review. Front. Physiol. 11:253. doi: 10.3389/fphys.2020.00253
Csapo, R., Malis, V., Sinha, U., Du, J., and Sinha, S. (2014). Age-associated differences in triceps surae muscle composition and strength—an MRI-based cross-sectional comparison of contractile, adipose, and connective tissue. BMC Musculoskelet. Disord. 15:209. doi: 10.1186/1471-2474-15-209
de Sousa Neto, I. V., Durigan, J. L. Q., Guzzoni, V., Tibana, R. A., Prestes, J., de Araujo, H. S. S., et al. (2018). Effects of resistance training on matrix metalloproteinase activity in skeletal muscles and blood circulation during aging. Front. Physiol. 9:190. doi: 10.3389/fphys.2018.00190
Degens, H., Erskine, R. M., and Morse, C. I. (2009). Disproportionate changes in skeletal muscle strength and size with resistance training and ageing. J. Musculoskelet. Neuronal Interact. 9, 123–129.
Dick, T. J. M., and Wakeling, J. M. (2017). Shifting gears: dynamic muscle shape changes and force-velocity behavior in the medial gastrocnemius. J. Appl. Physiol. 123, 1433–1442. doi: 10.1152/japplphysiol.01050.2016
Drace, J. E., and Pelc, N. J. (1994). Measurement of skeletal muscle motion in vivo with phase-contrast MR imaging. J. Magn. Reson. Imaging 4, 157–163. doi: 10.1002/jmri.1880040211
Du, J., Hamilton, G., Takahashi, A., Bydder, M., and Chung, C. B. (2007a). Ultrashort echo time spectroscopic imaging (UTESI) of cortical bone. Magn. Reson. Med. 58, 1001–1009. doi: 10.1002/mrm.21397
Du, J., Sinha, S., Chung, C. B., and Bydder, G. M. (2007b). Magnetic resonance imaging of the menisci of the knee using a 3D ultrashort TE pulse sequence. Chicago, IL.
Eng, C. M., Azizi, E., and Roberts, T. J. (2018). Structural determinants of muscle gearing during dynamic contractions. Integr. Comp. Biol. 58, 207–218. doi: 10.1093/icb/icy054
Englund, E. K., Elder, C. P., Xu, Q., Ding, Z., and Damon, B. M. (2011). Combined diffusion and strain tensor MRI reveals a heterogeneous, planar pattern of strain development during isometric muscle contraction. Am. J. Physiol. Regul. Integr. Comp. Physiol. 300, R1079–R1090. doi: 10.1152/ajpregu.00474.2010
Eskreis-Winkler, S., Corrias, G., Monti, S., Zheng, J., Capanu, M., Krebs, S., et al. (2018). IDEAL-IQ in an oncologic population: meeting the challenge of concomitant liver fat and liver iron. Cancer Imaging 18:1409. doi: 10.1186/s40644-018-0167-3
Faulkner, J. A., Larkin, L. M., Claflin, D. R., and Brooks, S. V. (2007). Age-related changes in the structure and function of skeletal muscles. Clin. Exp. Pharmacol. Physiol. 34, 1091–1096. doi: 10.1111/j.1440-1681.2007.04752.x
Fieremans, E., Lemberskiy, G., Veraart, J., Sigmund, E. E., Gyftopoulos, S., and Novikov, D. S. (2017). In vivo measurement of membrane permeability and myofiber size in human muscle using time-dependent diffusion tensor imaging and the random permeable barrier model. NMR Biomed. 30:e3612. doi: 10.1002/nbm.3612
Finni, T., Hodgson, J. A., Lai, A. M., Edgerton, V. R., and Sinha, S. (2003). Nonuniform strain of human soleus aponeurosis-tendon complex during submaximal voluntary contractions in vivo. J. Appl. Physiol. 95, 829–837. doi: 10.1152/japplphysiol.00775.2002
Forcina, L., Miano, C., and Scicchitano, B. M. (2019). Signals from the niche: insights into the role of IGF-1 and IL-6 in modulating skeletal muscle fibrosis. Cell 8:232. doi: 10.3390/cells8030232
Frich, L. H., Lambertsen, K. L., Hjarbaek, J., Dahl, J. S., and Holsgaard-Larsen, A. (2019). Musculoskeletal application and validation of speckle-tracking ultrasonography. BMC Musculoskelet. Disord. 20:192. doi: 10.1186/s12891-019-2562-8
Froeling, M., Nederveen, A. J., Nicolay, K., and Strijkers, G. J. (2013). DTI of human skeletal muscle: the effects of diffusion encoding parameters, signal-to-noise ratio and T2 on tensor indices and fiber tracts. NMR Biomed. 26, 1339–1352. doi: 10.1002/nbm.2959
Galbán, C. J., Maderwald, S., Uffmann, K., de Greiff, A., and Ladd, M. E. (2004). Diffusive sensitivity to muscle architecture: a magnetic resonance diffusion tensor imaging study of the human calf. Eur. J. Appl. Physiol. 93, 253–262. doi: 10.1007/s00421-004-1186-2
Galbán, C. J., Maderwald, S., Uffmann, K., and Ladd, M. E. (2005). A diffusion tensor imaging analysis of gender differences in water diffusivity within human skeletal muscle. NMR Biomed. 18, 489–498. doi: 10.1002/nbm.975
Gao, Y., Waas, A. M., Faulkner, J. A., Kostrominova, T. Y., and Wineman, A. S. (2008a). Micromechanical modeling of the epimysium of the skeletal muscles. J. Biomech. 41, 1–10. doi: 10.1016/j.jbiomech.2007.08.008
Gao, Y., Wineman, A. S., and Waas, A. M. (2008b). Mechanics of muscle injury induced by lengthening contraction. Ann. Biomed. Eng. 36, 1615–1623. doi: 10.1007/s10439-008-9547-3
Giannelli, G., De Marzo, A., Marinosci, F., and Antonaci, S. (2005). Matrix metalloproteinase imbalance in muscle disuse atrophy. Histol. Histopathol. 20, 99–106. doi: 10.14670/HH-20.99
Gillies, A. R., and Lieber, R. L. (2011). Structure and function of the skeletal muscle extracellular matrix. Muscle Nerve 44, 318–331. doi: 10.1002/mus.22094
Goodpaster, B. H., Park, S. W., Harris, T. B., Kritchevsky, S. B., Nevitt, M., Schwartz, A. V., et al. (2006). The loss of skeletal muscle strength, mass, and quality in older adults: the health aging and body composition study. J. Gerontol. A Biol. Sci. Med. Sci. 61, 1059–1064. doi: 10.1093/gerona/61.10.1059
Grimm, A., Meyer, H., Nickel, M. D., Nittka, M., Raithel, E., Chaudry, O., et al. (2018). Repeatability of dixon magnetic resonance imaging and magnetic resonance spectroscopy for quantitative muscle fat assessments in the thigh. J. Cachexia. Sarcopenia Muscle 9, 1093–1100. doi: 10.1002/jcsm.12343
Haus, J. M., Carrithers, J. A., Trappe, S. W., and Trappe, T. A. (2007). Collagen, cross-linking, and advanced glycation end products in aging human skeletal muscle. J. Appl. Physiol. 103, 2068–2076. doi: 10.1152/japplphysiol.00670.2007
Hindle, A. G., Horning, M., Mellish, J. -A. E., and Lawler, J. M. (2009). Diving into old age: muscular senescence in a large-bodied, long-lived mammal, the weddell seal (Leptonychotes weddellii). J. Exp. Biol. 212, 790–796. doi: 10.1242/jeb.025387
Holt, N. C., Danos, N., Roberts, T. J., and Azizi, E. (2016). Stuck in gear: age-related loss of variable gearing in skeletal muscle. J. Exp. Biol. 219, 998–1003. doi: 10.1242/jeb.133009
Hu, H. H., Branca, R. T., Hernando, D., Karampinos, D. C., Machann, J., McKenzie, C. A., et al. (2019). Magnetic resonance imaging of obesity and metabolic disorders: summary from the 2019 ISMRM workshop. Magn. Reson. Med. 83, 1565–1576. doi: 10.1002/mrm.28103
Hughes, D. C., Marcotte, G. R., Marshall, A. G., West, D. W. D., Baehr, L. M., Wallace, M. A., et al. (2017). Age-related differences in dystrophin: impact on force transfer proteins, membrane integrity, and neuromuscular junction stability. J. Gerontol. A Biol. Sci. Med. Sci. 72, 640–648. doi: 10.1093/gerona/glw109
Huijing, P. A. (1999). Muscle as a collagen fiber reinforced composite: a review of force transmission in muscle and whole limb. J. Biomech. 32, 329–345. doi: 10.1016/s0021-9290(98)00186-9
Huijing, P. A. (2003). Muscular force transmission necessitates a multilevel integrative approach to the analysis of function of skeletal muscle. Exerc. Sport Sci. Rev. 31, 167–175. doi: 10.1097/00003677-200310000-00003
Huijing, P. A., Baan, G. C., and Rebel, G. T. (1998). Non-nyotendinous force transmission in rat extensor digitorum longus muscle. J. Exp. Biol. 201, 683–691.
Johansson, T., Meier, P., and Blickhan, R. (2000). A finite-element model for the mechanical analysis of skeletal muscles. J. Theor. Biol. 206, 131–149. doi: 10.1006/jtbi.2000.2109
Karakuzu, A., Pamuk, U., Ozturk, C., Acar, B., and Yucesoy, C. A. (2017). Magnetic resonance and diffusion tensor imaging analyses indicate heterogeneous strains along human medial gastrocnemius fascicles caused by submaximal plantar-flexion activity. J. Biomech. 57, 69–78. doi: 10.1016/j.jbiomech.2017.03.028
Karampinos, D. C., King, K. F., Sutton, B. P., and Georgiadis, J. G. (2007). In vivo study of cross-sectional skeletal muscle fiber asymmetry with diffusion-weighted MRI. Conf. Proc. IEEE Eng. Med. Biol. Soc. 2007, 327–330. doi: 10.1109/IEMBS.2007.4352290
Karampinos, D. C., King, K. F., Sutton, B. P., and Georgiadis, J. G. (2009). Myofiber ellipticity as an explanation for transverse asymmetry of skeletal muscle diffusion MRI in vivo signal. Ann. Biomed. Eng. 37, 2532–2546. doi: 10.1007/s10439-009-9783-1
Karpakka, J., Väänänen, K., Orava, S., and Takala, T. E. (1990). The effects of preimmobilization training and immobilization on collagen synthesis in rat skeletal muscle. Int. J. Sports Med. 11, 484–488. doi: 10.1055/s-2007-1024842
Karpakka, J., Virtanen, P., Väänänen, K., Orava, S., and Takala, T. E. (1991). Collagen synthesis in rat skeletal muscle during immobilization and remobilization. J. Appl. Physiol. 70, 1775–1780. doi: 10.1152/jappl.1991.70.4.1775
Kinugasa, R., Shin, D., Yamauchi, J., Mishra, C., Hodgson, J. A., Edgerton, V. R., et al. (2008). Phase-contrast MRI reveals mechanical behavior of superficial and deep aponeuroses in human medial gastrocnemius during isometric contraction. J. Appl. Physiol. 105, 1312–1320. doi: 10.1152/japplphysiol.90440.2008
Kovanen, V., Suominen, H., Risteli, J., and Risteli, L. (1988). Type IV collagen and laminin in slow and fast skeletal muscle in rats—effects of age and life-time endurance training. Coll. Relat. Res. 8, 145–153. doi: 10.1016/S0174-173X(88)80026-8
Kragstrup, T. W., Kjaer, M., and Mackey, A. L. (2011). Structural, biochemical, cellular, and functional changes in skeletal muscle extracellular matrix with aging. Scand. J. Med. Sci. Sports 21, 749–757. doi: 10.1111/j.1600-0838.2011.01377.x
Lemos, R. R., Epstein, M., Herzog, W., and Wyvill, B. (2004). A framework for structured modeling of skeletal muscle. Comput. Methods Biomech. Biomed. Engin. 7, 305–317. doi: 10.1080/10255840412331317398
Li, K., Dortch, R. D., Welch, E. B., Bryant, N. D., Buck, A. K. W., Towse, T. F., et al. (2014). Multi-parametric MRI characterization of healthy human thigh muscles at 3.0 T—relaxation, magnetization transfer, fat/water, and diffusion tensor imaging. NMR Biomed. 27, 1070–1084. doi: 10.1002/nbm.3159
Lieber, R. L., and Ward, S. R. (2013). Cellular mechanisms of tissue fibrosis. 4. Structural and functional consequences of skeletal muscle fibrosis. Am. J. Physiol. Cell Physiol. 305, C241–C252. doi: 10.1152/ajpcell.00173.2013
Liu, X., Lee, D. J., Skittone, L. K., Natsuhara, K., and Kim, H. T. (2010). Role of gelatinases in disuse-induced skeletal muscle atrophy. Muscle Nerve 41, 174–178. doi: 10.1002/mus.21463
Lopata, R. G. P., van Dijk, J. P., Pillen, S., Nillesen, M. M., Maas, H., Thijssen, J. M., et al. (2010). Dynamic imaging of skeletal muscle contraction in three orthogonal directions. J. Appl. Physiol. 109, 906–915. doi: 10.1152/japplphysiol.00092.2010
Maganaris, C. N., Baltzopoulos, V., and Sargeant, A. J. (1998). In vivo measurements of the triceps surae complex architecture in man: implications for muscle function. J. Physiol. 512, 603–614. doi: 10.1111/j.1469-7793.1998.603be.x
Malis, V., Sinha, U., Csapo, R., Narici, M., and Sinha, S. (2018a). Relationship of changes in strain rate indices estimated from velocity-encoded MR imaging to loss of muscle force following disuse atrophy. Magn. Reson. Med. 79, 912–922. doi: 10.1002/mrm.26759
Malis, V., Sinha, U., Csapo, R., Narici, M., Smitaman, E., and Sinha, S. (2018b). Diffusion tensor imaging and diffusion modeling: application to monitoring changes in the medial gastrocnemius in disuse atrophy induced by unilateral limb suspension. J. Magn. Reson. Imaging 49, 1655–1664. doi: 10.1002/jmri.26295
Malis, V., Sinha, U., and Sinha, S. (2019b). Compressed sensing velocity encoded phase contrast imaging: monitoring skeletal muscle kinematics. Proceedings of the 27th Annual Meeting of ISMRM; Montreal.
Malis, V., Sinha, U., and Sinha, S. (2020). Compressed sensing velocity encoded phase contrast imaging: monitoring skeletal muscle kinematics. Magn. Reson. Med. 84, 142–156. doi: 10.1002/mrm.28100
Malis, V., Sinha, S., Smitaman, E., and Sinha, U. (2019a). Permeable barrier modeling of age induced changes in the time dependent diffusion eigenvalues. Proceedings of the 27th Annual Meeting of ISMRM; Montreal.
Martins, J. A. C., Pires, E. B., Salvado, R., and Dinis, P. B. (1998). A numerical model of passive and active behavior of skeletal muscles. Comput. Methods Appl. Mech. Eng. 151, 419–433. doi: 10.1016/S0045-7825(97)00162-X
Marty, B., and Carlier, P. G. (2019). MR fingerprinting for water T1 and fat fraction quantification in fat infiltrated skeletal muscles. Magn. Reson. Med. 83, 621–634. doi: 10.1002/mrm.27960
Meyer, G. A., and Lieber, R. L. (2011). Elucidation of extracellular matrix mechanics from muscle fibers and fiber bundles. J. Biomech. 44, 771–773. doi: 10.1016/j.jbiomech.2010.10.044
Morse, C. I., Thom, J. M., Reeves, N. D., Birch, K. M., and Narici, M. V. (2005). In vivo physiological cross-sectional area and specific force are reduced in the gastrocnemius of elderly men. J. Appl. Physiol. 99, 1050–1055. doi: 10.1152/japplphysiol.01186.2004
Mul, K., Horlings, C. G. C., Vincenten, S. C. C., Voermans, N. C., van Engelen, B. G. M., and van Alfen, N. (2018). Quantitative muscle MRI and ultrasound for facioscapulohumeral muscular dystrophy: complementary imaging biomarkers. J. Neurol. 265, 2646–2655. doi: 10.1007/s00415-018-9037-y
Nijboer-Oosterveld, J., van Alfen, N., and Pillen, S. (2011). New normal values for quantitative muscle ultrasound: obesity increases muscle echo intensity. Muscle Nerve 43, 142–143. doi: 10.1002/mus.21866
Oomens, C. W. J., Maenhout, M., van Oijen, C. H., Drost, M. R., and Baaijens, F. P. (2003). Finite element modeling of contracting skeletal muscle. Philos. Trans. R. Soc. Lond. B Biol. Sci. 358, 1453–1460. doi: 10.1098/rstb.2003.1345
Oudeman, J., Nederveen, A. J., Strijkers, G. J., Maas, M., Luijten, P. R., and Froeling, M. (2016). Techniques and applications of skeletal muscle diffusion tensor imaging: a review. J. Magn. Reson. Imaging 43, 773–788. doi: 10.1002/jmri.25016
Pamuk, U., Karakuzu, A., Ozturk, C., Acar, B., and Yucesoy, C. A. (2016). Combined magnetic resonance and diffusion tensor imaging analyses provide a powerful tool for in vivo assessment of deformation along human muscle fibers. J. Mech. Behav. Biomed. Mater. 63, 207–219. doi: 10.1016/j.jmbbm.2016.06.031
Pappas, G. P., Asakawa, D. S., Delp, S. L., Zajac, F. E., and Drace, J. E. (2002). Nonuniform shortening in the biceps brachii during elbow flexion. J. Appl. Physiol. 92, 2381–2389. doi: 10.1152/japplphysiol.00843.2001
Pattison, J. S., Folk, L. C., Madsen, R. W., Childs, T. E., and Booth, F. W. (2003). Transcriptional profiling identifies extensive downregulation of extracellular matrix gene expression in sarcopenic rat soleus muscle. Physiol. Genomics 15, 34–43. doi: 10.1152/physiolgenomics.00040.2003
Polgar, J., Johnson, M. A., Weightman, D., and Appleton, D. (1973). Data on fiber size in thirty-six human muscles an autopsy study. J. Neurol. Sci. 19, 307–318. doi: 10.1016/0022-510x(73)90094-4
Purslow, P. P. (2010). Muscle fascia and force transmission. J. Bodyw. Mov. Ther. 14, 411–417. doi: 10.1016/j.jbmt.2010.01.005
Purslow, P. P., and Trotter, J. A. (1994). The morphology and mechanical properties of endomysium in series-fibred muscles: variations with muscle length. J. Muscle Res. Cell Motil. 15, 299–308. doi: 10.1007/BF00123482
Ramani, A., Dalton, C., Miller, D. H., Tofts, P. S., and Barker, G. J. (2002). Precise estimate of fundamental in-vivo MT parameters in human brain in clinically feasible times. Magn. Reson. Imaging 20, 721–731. doi: 10.1016/S0730-725X(02)00598-2
Ramaswamy, K. S., Palmer, M. L., van der Meulen, J. H., Renoux, A., Kostrominova, T. Y., Michele, D. E., et al. (2011). Lateral transmission of force is impaired in skeletal muscles of dystrophic mice and very old rats. J. Physiol. 589, 1195–1208. doi: 10.1113/jphysiol.2010.201921
Rottmar, M., Haralampieva, D., Salemi, S., Eberhardt, C., Wurnig, M. C., Boss, A., et al. (2016). Magnetization transfer MR imaging to monitor muscle tissue formation during myogenic in vivo differentiation of muscle precursor cells. Radiology 281, 436–443. doi: 10.1148/radiol.2016152330
Savolainen, J., Väänänen, K., Puranen, J., Takala, T. E., Komulainen, J., and Vihko, V. (1988). Collagen synthesis and proteolytic activities in rat skeletal muscles: effect of cast-immobilization in the lengthened and shortened positions. Arch. Phys. Med. Rehabil. 69, 964–969.
Sharafi, B., and Blemker, S. S. (2010). A micromechanical model of skeletal muscle to explore the effects of fiber and fascicle geometry. J. Biomech. 43, 3207–3213. doi: 10.1016/j.jbiomech.2010.07.020
Sharafi, B., and Blemker, S. S. (2011). A mathematical model of force transmission from intrafascicularly terminating muscle fibers. J. Biomech. 44, 2031–2039. doi: 10.1016/j.jbiomech.2011.04.038
Shin, D. D., Hodgson, J. A., Edgerton, V. R., and Sinha, S. (2009a). In vivo intramuscular fascicle-aponeuroses dynamics of the human medial gastrocnemius during plantarflexion and dorsiflexion of the foot. J. Appl. Physiol. 107, 1276–1284. doi: 10.1152/japplphysiol.91598.2008
Sigmund, E. E., Sigmund, E. E., Novikov, D. S., Novikov, D. S., Sui, D., Ukpebor, O., et al. (2014). Time-dependent diffusion in skeletal muscle with the random permeable barrier model (RPBM): application to normal controls and chronic exertional compartment syndrome patients. NMR Biomed. 27, 519–528. doi: 10.1002/nbm.3087
Sinclair, C. D. J., Samson, R. S., Thomas, D. L., Weiskopf, N., Lutti, A., Thornton, J. S., et al. (2010). Quantitative magnetization transfer in in vivo healthy human skeletal muscle at 3 T. Magn. Reson. Med. 64, 1739–1748. doi: 10.1002/mrm.22562
Sinha, S., Du, J., Takahashi, A. M., Chung, C. B., and Bydder, G. M. (2007). Estimation of T1 and T2 of deep radial and calcified layers in human patella from ultra short echo time (UTE) imaging on a 3T clinical system. Berlin.
Sinha, U., Malis, V., Csapo, R., Moghadasi, A., Kinugasa, R., and Sinha, S. (2015a). Age-related differences in strain rate tensor of the medial gastrocnemius muscle during passive plantarflexion and active isometric contraction using velocity encoded MR imaging: potential index of lateral force transmission. Magn. Reson. Med. 73, 1852–1863. doi: 10.1002/mrm.25312
Sinha, U., Malis, V., Csapo, R., Narici, M., and Sinha, S. (2018). Shear strain rate from phase contrast velocity encoded MRI: application to study effects of aging in the medial gastrocnemius muscle. J. Magn. Reson. Imaging 48, 1351–1357. doi: 10.1002/jmri.26030
Sinha, U., Malis, V., Xue, Y., and Sinha, S. (2015b). Age-related differences in diffusion tensor indices and fiber architecture in the medial and lateral gastrocnemius. J. Magn. Reson. Imaging 41, 941–953. doi: 10.1002/jmri.24641
Spyrou, L. A., Agoras, M., and Danas, K. (2017). A homogenization model of the voigt type for skeletal muscle. J. Theor. Biol. 414, 50–61. doi: 10.1016/j.jtbi.2016.11.018
Spyrou, L. A., Brisard, S., and Danas, K. (2019). Multiscale modeling of skeletal muscle tissues based on analytical and numerical homogenization. J. Mech. Behav. Biomed. Mater. 92, 97–117. doi: 10.1016/j.jmbbm.2018.12.030
Stearns-Reider, K. M., D’Amore, A., Beezhold, K., Rothrauff, B., Cavalli, L., Wagner, W. R., et al. (2017). Aging of the skeletal muscle extracellular matrix drives a stem cell fibrogenic conversion. Aging Cell 16, 518–528. doi: 10.1111/acel.12578
Tan, K., Jugé, L., Hatt, A., Cheng, S., and Bilston, L. E. (2018). Measurement of large strain properties in calf muscles in vivo using magnetic resonance elastography and spatial modulation of magnetization. NMR Biomed. 31:e3925. doi: 10.1002/nbm.3925
Ugarte, V., Sinha, U., Malis, V., Csapo, R., and Sinha, S. (2017b). 3D multimodal spatial fuzzy segmentation of intramuscular connective and adipose tissue from ultrashort TE MR images of calf muscle. Magn. Reson. Med. 77, 870–883. doi: 10.1002/mrm.26156
Ugarte, V., Sinha, U., Malis, V., and Sinha, S. (2017a). Automated segmentation of intramuscular connective tissue (IMCT) from skeletal muscle in presence of artifacts: application to changes in IMCT in an unilateral limb suspension induced acute atrophy model in the Plantarflexors. Honolulu, HI: John Wiley & Sons, Ltd.
van der Linden, B. J. J. J., Koopman, H. F. J. M., Grootenboer, H. J., and Huijing, P. A. (1998). Modeling functional effects of muscle geometry. J. Electromyogr. Kinesiol. 8, 101–109. doi: 10.1016/S1050-6411(97)00026-6
Virgilio, K. M., Martin, K. S., Peirce, S. M., and Blemker, S. S. (2015). Multiscale models of skeletal muscle reveal the complex effects of muscular dystrophy on tissue mechanics and damage susceptibility. Interface Focus 5:20140080. doi: 10.1098/rsfs.2014.0080
White, J. (2019). Quantitative magnetization transfer to monitor age related effects in skeletal muscle.
Wood, L. K., Kayupov, E., Gumucio, J. P., Mendias, C. L., Claflin, D. R., and Brooks, S. V. (2014). Intrinsic stiffness of extracellular matrix increases with age in skeletal muscles of mice. J. Appl. Physiol. 117, 363–369. doi: 10.1152/japplphysiol.00256.2014
Yucesoy, C. A. (2010). Epimuscular myofascial force transmission implies novel principles for muscular mechanics. Exerc. Sport Sci. Rev. 38, 128–134. doi: 10.1097/JES.0b013e3181e372ef
Yucesoy, C. A., Koopman, B. H. F. J. M., Huijing, P. A., and Grootenboer, H. J. (2002). Three-dimensional finite element modeling of skeletal muscle using a two-domain approach: linked fiber-matrix mesh model. J. Biomech. 35, 1253–1262. doi: 10.1016/s0021-9290(02)00069-6
Yucesoy, C. A., Maas, H., Koopman, B. H. F. J. M., Grootenboer, H. J., and Huijing, P. A. (2006). Mechanisms causing effects of muscle position on proximo-distal muscle force differences in extra-muscular myofascial force transmission. Med. Eng. Phys. 28, 214–226. doi: 10.1016/j.medengphy.2005.06.004
Zajac, F. E. (1989). Muscle and tendon: properties, models, scaling, and application to biomechanics and motor control. Crit. Rev. Biomed. Eng. 17, 359–411.
Zhang, Y., Chen, J. -S., He, Q., He, X., Basava, R. R., Hodgson, J., et al. (2019). Microstructural analysis of skeletal muscle force generation during aging. Int. J. Numer. Method. Biomed. Eng. 52:B125. doi: 10.1002/cnm.3295
Zhang, C., and Gao, Y. (2014). Effects of aging on the lateral transmission of force in rat skeletal muscle. J. Biomech. 47, 944–948. doi: 10.1016/j.jbiomech.2014.01.026
Keywords: extracellular matrix, age and disuse related muscle force loss, structural muscle MRI, strain imaging, lateral transmission of force
Citation: Sinha U, Malis V, Chen J-S, Csapo R, Kinugasa R, Narici MV and Sinha S (2020) Role of the Extracellular Matrix in Loss of Muscle Force With Age and Unloading Using Magnetic Resonance Imaging, Biochemical Analysis, and Computational Models. Front. Physiol. 11:626. doi: 10.3389/fphys.2020.00626
Edited by:
Can A. Yucesoy, Boğaziçi University, TurkeyReviewed by:
Carolyn Eng, Yale University, United StatesDavid C. Hughes, The University of Iowa, United States
Copyright © 2020 Sinha, Malis, Chen, Csapo, Kinugasa, Narici and Sinha. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shantanu Sinha, shsinha@ucsd.edu