Imaging two-dimensional displacements and strains in skeletal muscle during joint motion by cine DENSE MR

Compositional and Functional MRI of Skeletal Muscle: A Review

Due to its exceptional sensitivity to soft tissues, MRI has been extensively utilized to assess anatomical muscle parameters such as muscle volume and cross-sectional area. Quantitative Magnetic Resonance Imaging (qMRI) adds to the capabilities of MRI, by providing information on muscle composition such as fat content, water content, microstructure, hypertrophy, atrophy, as well as muscle architecture. In addition to compositional changes, qMRI can also be used to assess function for example by measuring muscle quality or through characterization of muscle deformation during passive lengthening/shortening and active contractions. The overall aim of this review is to provide an updated overview of qMRI techniques that can quantitatively evaluate muscle structure and composition, provide insights into the underlying biological basis of the qMRI signal, and illustrate how qMRI biomarkers of muscle health relate to function in healthy and diseased/injured muscles. While some applications still require systematic clinical validation, qMRI is now established as a comprehensive technique, that can be used to characterize a wide variety of structural and compositional changes in healthy and diseased skeletal muscle. Taken together, multiparametric muscle MRI holds great potential in the diagnosis and monitoring of muscle conditions in research and clinical applications. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 2.

The need for speed - Does the force-velocity property significantly alter strain distributions within skeletal muscle?

Skeletal muscles are complex structures with nonlinear constitutive properties. This complexity often requires finite element (FE) modeling to better understand muscle behavior and response to activation, especially the fiber strain distributions that can be difficult to measure in vivo. However, many FE muscle models designed to study fiber strain do not include force-velocity behavior. To investigate force-velocity property impact on strain distributions within skeletal muscle, we modified a muscle constitutive model with active and passive force-length properties to include force-velocity properties. We implemented the new constitutive model as a plugin for the FE software FEBio and applied it to four geometries: 1) a single element, 2) a multiple-element model representing a single fiber, 3) a model of tapering fibers, and 4) a model representing the bicep femoris long head (BFLH) morphology. Maximum fiber velocity and boundary conditions of the finite element models were varied to test their influence on fiber strain distribution. We found that force-velocity properties in the constitutive model behaved as expected for the single element and multi-element conditions. In the tapered fiber models, fiber strain distributions were impacted by changes in maximum fiber velocity; the range of strains increased with maximum fiber velocity, which was most noted in isometric contraction simulations. In the BFLH model, maximum fiber velocity had minimal impact on strain distributions, even in the context of sprinting. Taken together, the combination of muscle model geometry, activation, and displacement parameters play a critical part in determining the magnitude of impact of force-velocity on strain distribution.

Muscle fiber strain rates in the lower leg during ankle dorsi-/plantarflexion exercise

Static quantitative magnetic resonance imaging (MRI) provides readouts of structural changes in diseased muscle, but current approaches lack the ability to fully explain the loss of contractile function. Muscle contractile function can be assessed using various techniques including phase-contrast MRI (PC-MRI), where strain rates are quantified. However, current two-dimensional implementations are limited in capturing the complex motion of contracting muscle in the context of its three-dimensional (3D) fiber architecture. The MR acquisitions (chemical shift-encoded water-fat separation scan, spin echo-echoplanar imaging with diffusion weighting, and two time-resolved 3D PC-MRI) wereperformed at 3 T. PC-MRI acquisitions and performed with and without load at 7.5% of the maximum voluntary dorsiflexion contraction force. Acquisitions (3 T, chemical shift-encoded water-fat separation scan, spin echo-echo planar imaging with diffusion weighting, and two time-resolved 3D PC-MRI) were performed with and without load at 7.5% of the maximum voluntary dorsiflexion contraction force. Strain rates and diffusion tensors were calculated and combined to obtain strain rates along and perpendicular to the muscle fibers in seven lower leg muscles during the dynamic dorsi-/plantarflexion movement cycle. To evaluate strain rates along the proximodistal muscle axis, muscles were divided into five equal segments. t-tests were used to test if cyclic strain rate patterns (amplitude > 0) were present along and perpendicular to the muscle fibers. The effects of proximal-distal location and load were evaluated using repeated measures ANOVAs. Cyclic temporal strain rate patterns along and perpendicular to the fiber were found in all muscles involved in dorsi-/plantarflexion movement (p < 0.0017). Strain rates along and perpendicular to the fiber were heterogeneously distributed over the length of most muscles (p < 0.003). Additional loading reduced strain rates of the extensor digitorum longus and gastrocnemius lateralis muscle (p < 0.001). In conclusion, the lower leg muscles involved in cyclic dorsi-/plantarflexion exercise showed cyclic fiber strain rate patterns with amplitudes that varied between muscles and between the proximodistal segments within the majority of muscles.

In vivo imaging of skeletal muscle form and function: 50 years of insight

Skeletal muscle form and function has fascinated scientists for centuries. Our understanding of muscle function has long been driven by advancements in imaging techniques. For example, the sliding filament theory of muscle, which is now widely leveraged in biomechanics research, stemmed from observations made possible by scanning electron microscopy. Over the last 50 years, advancing in medical imaging, combined with ingenuity and creativity of biomechanists, have provide a wealth of new and important insights into in vivo human muscle function. Incorporation of in vivo imaging has also advanced computational modeling and allowed our research to have an impact in many clinical populations. While this review does not provide a comprehensive or meta-analysis of the all the in vivo muscle imaging work over the last five decades, it provides a narrative about the past, present, and future of in vivo muscle imaging.

A 3D model of the soleus reveals effects of aponeuroses morphology and material properties on complex muscle fascicle behavior

The soleus is an important plantarflexor muscle with complex fascicle and connective tissue arrangement. In this study we created an image-based finite element model representing the 3D structure of the soleus muscle and its aponeurosis connective tissue, including distinct fascicle architecture of the posterior and anterior compartments. The model was used to simulate passive and active soleus lengthening during ankle motion to predict tissue displacements and fascicle architecture changes. Both the model's initial architecture and changes incurred during passive lengthening were consistent with prior in vivo data from diffusion tensor imaging. Model predictions of active lengthening were consistent with axial plane muscle displacements that we measured in eight subjects' lower legs using cine DENSE (Displacement Encoding with Stimulated Echoes) MRI during eccentric dorsiflexion. Regional strains were variable and nonuniform in the model, but average fascicle strains were similar between the compartments for both passive (anterior: 0.18 ± 0.06, posterior: 0.19 ± 0.05) and active (anterior: 0.12 ± 0.05, posterior: 0.13 ± 0.06) lengthening and were two- to three-times greater than muscle belly strain (0.06). We used additional model simulations to investigate the effects of aponeurosis material properties on muscle deformation, by independently varying the longitudinal or transverse stiffness of the posterior or anterior aponeurosis. Results of model variations elucidate how properties of soleus aponeuroses contribute to fascicle architecture changes. Greater longitudinal stiffness of posterior compared to anterior aponeurosis promoted more uniform spatial distribution of muscle tissue deformation. Reduced transverse stiffness in both aponeuroses resulted in larger differences between passive and active soleus lengthening.

Suppression of artifact-generating echoes in cine DENSE using deep learning

PURPOSE

To use deep learning for suppression of the artifact-generating T1 -relaxation echo in cine displacement encoding with stimulated echoes (DENSE) for the purpose of reducing the scan time.

METHODS

A U-Net was trained to suppress the artifact-generating T1 -relaxation echo using complementary phase-cycled data as the ground truth. A data-augmentation method was developed that generates synthetic DENSE images with arbitrary displacement-encoding frequencies to suppress the T1 -relaxation echo modulated for a range of frequencies. The resulting U-Net (DAS-Net) was compared with k-space zero-filling as an alternative method. Non-phase-cycled DENSE images acquired in shorter breath-holds were processed by DAS-Net and compared with DENSE images acquired with phase cycling for the quantification of myocardial strain.

RESULTS

The DAS-Net method effectively suppressed the T1 -relaxation echo and its artifacts, and achieved root Mean Square(RMS) error = 5.5 ± 0.8 and structural similarity index = 0.85 ± 0.02 for DENSE images acquired with a displacement encoding frequency of 0.10 cycles/mm. The DAS-Net method outperformed zero-filling (root Mean Square error = 5.8 ± 1.5 vs 13.5 ± 1.5, DAS-Net vs zero-filling, P < .01; and structural similarity index = 0.83 ± 0.04 vs 0.66 ± 0.03, DAS-Net vs zero-filling, P < .01). Strain data for non-phase-cycled DENSE images with DAS-Net showed close agreement with strain from phase-cycled DENSE.

CONCLUSION

The DAS-Net method provides an effective alternative approach for suppression of the artifact-generating T1 -relaxation echo in DENSE MRI, enabling a 42% reduction in scan time compared to DENSE with phase-cycling.

Quantitative Magnetic Resonance Imaging Assessment of the Quadriceps Changes during an Extreme Mountain Ultramarathon

INTRODUCTION/PURPOSE

Extreme ultra-endurance races are growing in popularity, but their effects on skeletal muscles remain mostly unexplored. This longitudinal study explores physiological changes in mountain ultramarathon athletes' quadriceps using quantitative magnetic resonance imaging (MRI) coupled with serological biomarkers. The study aimed to monitor the longitudinal effect of the race and recovery and to identify local inflammatory and metabolic muscle responses by codetection of biological markers.

METHODS

An automatic image processing framework was designed to extract imaging-based biomarkers from quantitative MRI acquisitions of the upper legs of 20 finishers at three time points. The longitudinal effect of the race was demonstrated by analyzing the image markers with dedicated biostatistical analysis.

RESULTS

Our framework allows for a reliable calculation of statistical data not only inside the whole quadriceps volume but also within each individual muscle head. Local changes in MRI parameters extracted from quantitative maps were described and found to be significantly correlated with principal serological biomarkers of interest. A decrease in the PDFF after the race and a stable paramagnetic susceptibility value were found. Pairwise post hoc tests suggested that the recovery process differs among the muscle heads.

CONCLUSIONS

This longitudinal study conducted during a prolonged and extreme mechanical stress showed that quantitative MRI-based markers of inflammation and metabolic response can detect local changes related to the prolonged exercise, with differentiated involvement of each head of the quadriceps muscle as expected in such eccentric load. Consistent and efficient extraction of the local biomarkers enables to highlight the interplay/interactions between blood and MRI biomarkers. This work indeed proposes an automatized analytic framework to tackle the time-consuming and mentally exhausting segmentation task of muscle heads in large multi-time-point cohorts.

Dynamic MRI of plantar flexion: A comprehensive repeatability study of electrical stimulation-gated muscle contraction standardized on evoked force

Quantification of skeletal muscle contraction in Magnetic Resonance Imaging (MRI) is a non-invasive method for studying muscle motion and deformation. The aim of this study was to evaluate the repeatability of quantitative measures such as strain, based on single slice dynamic MRI synchronized with neuromuscular electrical stimulation (NMES) and standardized to a similar relative force level across various individuals. Unilateral electrical stimulation of the triceps surae muscles was applied in eight volunteers during single-slice, three-directional phase contrast MRI acquisition at a 3T MRI scanner. To assess repeatability, the same process was executed on two different days by standardizing the stimulation aiming at evoking a fixed percentage of their maximal voluntary force in the same position. Except from the force, the effect of using the current as reference was evaluated on day two as a secondary acquisition. Finally, the presence of fatigue induced by NMES was assessed (on day one) by examining the difference between consecutive measurements. Strain maps were derived from the acquired slice at every time point; distribution of strain in the muscle and peak strain over the muscle of interest were evaluated for repeatability. It was found that fatigue did not have an appreciable effect on the results. The stimulation settings based on evoked force produced more repeatable results with respect to using the current as the only reference, with an intraclass correlation coefficient between different days of 0.95 for the former versus 0.88 for the latter. In conclusion, for repeatable strain imaging it is advisable to record the force output of the evoked contraction and use that for the standardization of the NMES setup rather than the current.

Role of the Extracellular Matrix in Loss of Muscle Force With Age and Unloading Using Magnetic Resonance Imaging, Biochemical Analysis, and Computational Models

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.

OpenArm 2.0: Automated Segmentation of 3D Tissue Structures for Multi-Subject Study of Muscle Deformation Dynamics

We present a novel neural-network-based pipeline for segmentation of 3D muscle and bone structures from localized 2D ultrasound data of the human arm. Building from the U-Net [1] neural network framework, we examine various data augmentation techniques and training data sets to both optimize the network's performance on our data set and hypothesize strategies to better select training data, minimizing manual annotation time while maximizing performance. We then employ this pipeline to generate the OpenArm 2.0 data set, the first factorial set of multi-subject, multi-angle, multi-force scans of the arm with full volumetric annotation of the biceps and humerus. This data set has been made available on SimTK (https://simtk.org/projects/openarm) to enable future exploration of muscle force modeling, improved musculoskeletal graphics, and assistive device control.

Exploration of New Contrasts, Targets, and MR Imaging and Spectroscopy Techniques for Neuromuscular Disease - A Workshop Report of Working Group 3 of the Biomedicine and Molecular Biosciences COST Action BM1304 MYO-MRI

Neuromuscular diseases are characterized by progressive muscle degeneration and muscle weakness resulting in functional disabilities. While each of these diseases is individually rare, they are common as a group, and a large majority lacks effective treatment with fully market approved drugs. Magnetic resonance imaging and spectroscopy techniques (MRI and MRS) are showing increasing promise as an outcome measure in clinical trials for these diseases. In 2013, the European Union funded the COST (co-operation in science and technology) action BM1304 called MYO-MRI (www.myo-mri.eu), with the overall aim to advance novel MRI and MRS techniques for both diagnosis and quantitative monitoring of neuromuscular diseases through sharing of expertise and data, joint development of protocols, opportunities for young researchers and creation of an online atlas of muscle MRI and MRS. In this report, the topics that were discussed in the framework of working group 3, which had the objective to: Explore new contrasts, new targets and new imaging techniques for NMD are described. The report is written by the scientists who attended the meetings and presented their data. An overview is given on the different contrasts that MRI can generate and their application, clinical needs and desired readouts, and emerging methods.

Contemporary image-based methods for measuring passive mechanical properties of skeletal muscles in vivo

Skeletal muscles' primary function in the body is mechanical: to move and stabilize the skeleton. As such, their mechanical behavior is a key aspect of their physiology. Recent developments in medical imaging technology have enabled quantitative studies of passive muscle mechanics, ranging from measurements of intrinsic muscle mechanical properties, such as elasticity and viscosity, to three-dimensional muscle architecture and dynamic muscle deformation and kinematics. In this review we summarize the principles and applications of contemporary imaging methods that have been used to study the passive mechanical behavior of skeletal muscles. Elastography measurements can provide in vivo maps of passive muscle mechanical parameters, and both MRI and ultrasound methods are available (magnetic resonance elastography and ultrasound shear wave elastography, respectively). Both have been shown to differentiate between healthy muscle and muscles affected by a broad range of clinical conditions. Detailed muscle architecture can now be depicted using diffusion tensor imaging, which not only is particularly useful for computational modeling of muscle but also has potential in assessing architectural changes in muscle disorders. More dynamic information about muscle mechanics can be obtained using a range of dynamic MRI methods, which characterize the detailed internal muscle deformations during motion. There are several MRI techniques available (e.g., phase-contrast MRI, displacement-encoded MRI, and "tagged" MRI), each of which can be collected in synchrony with muscle motion and postprocessed to quantify muscle deformation. Together, these modern imaging techniques can characterize muscle motion, deformation, mechanical properties, and architecture, providing complementary insights into skeletal muscle function.

Spatially resolved kinetics of skeletal muscle exercise response and recovery with multiple echo diffusion tensor imaging (MEDITI): a feasibility study

OBJECTIVES

We describe measurement of skeletal muscle kinetics with multiple echo diffusion tensor imaging (MEDITI). This approach allows characterization of the microstructural dynamics in healthy and pathologic muscle.

MATERIALS AND METHODS

In a Siemens 3-T Skyra scanner, MEDITI was used to collect dynamic DTI with a combination of rapid diffusion encoding, radial imaging, and compressed sensing reconstruction in a multi-compartment agarose gel rotation phantom and within in vivo calf muscle. An MR-compatible ergometer (Ergospect Trispect) was employed to enable in-scanner plantar flexion exercise. In a HIPAA-compliant study with written informed consent, post-exercise recovery of DTI metrics was quantified in eight volunteers. Exercise response of DTI metrics was compared with that of T2-weighted imaging and characterized by a gamma variate model.

RESULTS

Phantom results show quantification of diffusivities in each compartment over its full dynamic rotation. In vivo calf imaging results indicate larger radial than axial exercise response and recovery in the plantar flexion-challenged gastrocnemius medialis (fractional response: nT2w = 0.385 ± 0.244, nMD = 0.163 ± 0.130, nλ₁ = 0.110 ± 0.093, nλ_rad = 0.303 ± 0.185). Diffusion and T2-weighted response magnitudes were correlated (e.g., r = 0.792, p = 0.019 for nMD vs. nT2w).

CONCLUSION

We have demonstrated the feasibility of MEDITI for capturing spatially resolved diffusion tensor data in dynamic systems including post-exercise skeletal muscle recovery following in-scanner plantar flexion.

Ultrasound Imaging of Muscle Contraction of the Tibialis Anterior in Patients with Facioscapulohumeral Dystrophy

A need exists for biomarkers to diagnose, quantify and longitudinally follow facioscapulohumeral muscular dystrophy (FSHD) and many other neuromuscular disorders. Furthermore, the pathophysiological mechanisms leading to muscle weakness in most neuromuscular disorders are not completely understood. Dynamic ultrasound imaging (B-mode image sequences) in combination with speckle tracking is an easy, applicable and patient-friendly imaging tool to visualize and quantify muscle deformation. This dynamic information provides insight in the pathophysiological mechanisms and may help to distinguish the various stages of diseased muscle in FSHD. In this proof-of-principle study, we applied a speckle tracking technique to 2-D ultrasound image sequences to quantify the deformation of the tibialis anterior muscle in patients with FSHD and in healthy controls. The resulting deformation patterns were compared with muscle ultrasound echo intensity analysis (a measure of fat infiltration and dystrophy) and clinical outcome measures. Of the four FSHD patients, two patients had severe peroneal weakness and two patients had mild peroneal weakness on clinical examination. We found a markedly varied muscle deformation pattern between these groups: patients with severe peroneal weakness showed a different motion pattern of the tibialis anterior, with overall less displacement of the central tendon region, while healthy patients showed a non-uniform displacement pattern, with the central aponeurosis showing the largest displacement. Hence, dynamic muscle ultrasound of the tibialis anterior muscle in patients with FSHD revealed a distinctively different tissue deformation pattern among persons with and without tibialis anterior weakness. These findings could clarify the understanding of the pathophysiology of muscle weakness in FSHD patients. In addition, the change in muscle deformation shows good correlation with clinical measures and quantitative muscle ultrasound measurements. In conclusion, dynamic ultrasound in combination with speckle tracking allows the study of the effects of muscle pathology in relation to strength, force transmission and movement generation. Although further research is required, this technique can develop into a biomarker to quantify muscle disease severity.

Relationship of changes in strain rate indices estimated from velocity-encoded MR imaging to loss of muscle force following disuse atrophy

PURPOSE

This study explores changes in strain rate (SR) (rate of regional deformation) parameters extracted from velocity-encoded MRI and their relationship to muscle force loss following 4-week unilateral lower limb suspension in healthy humans.

METHODS

Two-dimensional SR maps were derived from three directional velocity-encoded MR phase-contrast images of the medial gastrocnemius in seven subjects. Atrophy-related and regional differences in the SR eigenvalues, angle between the SR and muscle fiber (SR-fiber angle), and strain rates in the fiber basis were statistically analyzed using analysis of variance and linear regression.

RESULTS

During isometric contraction, SR in the fiber cross section (SR_in-plane ) was significantly lower, and the SR-fiber angle was significantly higher postsuspension (P < 0.05). On multiple variable regression analysis, the volume of medial gastrocnemius, SR_in-plane , and SR-fiber angle were significantly associated with force and changes in the, and the SR eigenvalues and shear SR were significantly associated with change in force with disuse.

CONCLUSIONS

Changes in SR-fiber angle, SR_in-plane , and shear SR as well as their ability to predict force and force changes may reflect the role of remodeling of the extracellular matrix in disuse atrophy and its functional consequence in reducing lateral transmission of force. Magn Reson Med 79:912-922, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Characterization of three dimensional volumetric strain distribution during passive tension of the human tibialis anterior using Cine Phase Contrast MRI

Intramuscular pressure correlates strongly with muscle tension and is a promising tool for quantifying individual muscle force. However, clinical application is impeded by measurement variability that is not fully understood. Previous studies point to regional differences in IMP, specifically increasing pressure with muscle depth. Based on conservation of mass, intramuscular pressure and volumetric strain distributions may be inversely related. Therefore, we hypothesized volumetric strain would decrease with muscle depth. To test this we quantified 3D volumetric strain in the tibialis anterior of 12 healthy subjects using Cine Phase Contrast Magnetic Resonance Imaging. Cine Phase Contrast data were collected while a custom apparatus rotated the subjects' ankle continuously between neutral and plantarflexion. A T2-weighted image stack was used to define the resting tibials anterior position. Custom and commercial post-processing software were used to quantify the volumetric strain distribution. To characterize regional strain changes, the muscle was divided into superior-inferior sections and either medial-lateral or anterior-posterior slices. Mean volumetric strain was compared across the sections and slices. As hypothesized, volumetric strain demonstrated regional differences with a decreasing trend from the anterior (superficial) to the posterior (deep) muscle regions. Statistical tests showed significant main effects and interactions of superior-inferior and anterior-posterior position as well as superior-inferior and medial-lateral position on regional strain. These data support our hypothesis and imply a potential relationship between regional volumetric strain and intramuscular pressure. This finding may advance our understanding of intramuscular pressure variability sources and lead to more reliable measurement solutions in the future.

Method of quantifying 3D strain distribution in skeletal muscle using cine phase contrast MRI

Intramuscular pressure (IMP), a correlate of muscle tension, may fill an important clinical testing void. A barrier to implementing this measure clinically is its non-uniform distribution, which is not fully understood. Pressure is generated by changes in fluid mass and volume, therefore 3D volumetric strain distribution may affect IMP distribution. The purpose of this study was to develop a method for quantifying 3D volumetric strain distribution in the human tibialis anterior (TA) during passive tension using cine phase contrast (CPC) MRI and to assess its accuracy and precision.Five healthy subjects each participated in three data collections. A custom MRI-compatible apparatus repeatedly rotated a subject's ankle between 0° and 26° plantarflexion while CPC MRI data were collected. Additionally, T2-weighted images of the lower leg were collected both before and after the CPC data collection with the ankle stationary at both 0° and 26° plantarflexion for TA muscle segmentation. A 3D hexahedral mesh was generated based on the TA surface before CPC data collection with the ankle at 0° plantarflexion and the node trajectories were tracked using the CPC data. The volumetric strain of each element was quantified.Three tests were employed to assess the measure accuracy and precision. First, to quantify leg position drift, the TA segmentations were compared before and after CPC data collection. The Hawsdorff distance measure (error) was 1.5  ±  0.7 mm. Second, to assess the surface node trajectory accuracy, the deformed mesh surface was compared to the TA segmented at 26° of ankle plantarflexion. This error was 0.6  ±  0.2 mm. Third, the standard deviation of volumetric strain across the three data collections was calculated for each element and subject. The median between-day variability across subjects and mesh elements was 0.06 mm3 mm(-3) (95% confidence interval 0.01 to 0.18 mm3 mm(-3)). Overall the results demonstrated excellent accuracy and precision.

Computational models predict larger muscle tissue strains at faster sprinting speeds

INTRODUCTION

Proximal biceps femoris musculotendon strain injury has been well established as a common injury among athletes participating in sports that require sprinting near or at maximum speed; however, little is known about the mechanisms that make this muscle tissue more susceptible to injury at faster speeds.

PURPOSE

This study aimed to quantify localized tissue strain during sprinting at a range of speeds.

METHODS

Biceps femoris long head (BFlh) musculotendon dimensions of 14 athletes were measured on magnetic resonance (MR) images and used to generate a finite-element computational model. The model was first validated through comparison with previous dynamic MR experiments. After validation, muscle activation and muscle-tendon unit length change were derived from forward dynamic simulations of sprinting at 70%, 85%, and 100% maximum speed and used as input to the computational model simulations. Simulations ran from midswing to foot contact.

RESULTS

The model predictions of local muscle tissue strain magnitude compared favorably with in vivo tissue strain measurements determined from dynamic MR experiments of the BFlh. For simulations of sprinting, local fiber strain was nonuniform at all speeds, with the highest muscle tissue strain where injury is often observed (proximal myotendinous junction). At faster sprinting speeds, increases were observed in fiber strain nonuniformity and peak local fiber strain (0.56, 0.67, and 0.72 for sprinting at 70%, 85%, and 100% maximum speed). A histogram of local fiber strains showed that more of the BFlh reached larger local fiber strains at faster speeds.

CONCLUSIONS

At faster sprinting speeds, peak local fiber strain, fiber strain nonuniformity, and the amount of muscle undergoing larger strains are predicted to increase, likely contributing to the BFlh muscle's higher injury susceptibility at faster speeds.

Dynamic ultrasound imaging applications to quantify musculoskeletal function

Advances in imaging methods have led to new capability to study muscle and tendon motion in vivo. Direct measurements of muscle and tendon kinematics using imaging may lead to improved understanding of musculoskeletal function. This review presents quantitative ultrasound methods for muscle dynamics that can be used to assess in vivo musculoskeletal function when integrated with other conventional biomechanical measurements.

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

PURPOSE

The strain rate (SR) tensor measures the principal directions and magnitude of the instantaneous deformation; this study aims to track age-related changes in the 2D SR tensor in the medial gastrocnemius during passive joint rotation and active isometric contraction.

METHODS

SR tensors were derived from velocity encoded magnetic resonance phase-contrast images in nine young (28 years) and eight senior (78 years) women. Strain rates along and in the cross-section of the fiber were calculated from the SR tensor and used to derive the out-plane SR. Age-related and regional differences in the SR eigenvalues, orientation, and the angle between the SR and muscle fiber (SR-fiber angle) were statistically analyzed.

RESULTS

SR along the fiber was significantly different between the cohorts during isometric contraction with higher values in the young (P < 0.05). The SR-fiber angle was larger in the young for both motion types but this difference was not statistically significant. Significant regional differences in the SR indices was seen in passive joint rotation (P < 0.05) for both cohorts.

CONCLUSION

SR mapping reflects age-related and regional differences during active and passive motion respectively; this may arise from differences in contractility (active motion) and elastic properties (active and passive motion).

Accurate high-resolution measurements of 3-D tissue dynamics with registration-enhanced displacement encoded MRI

Displacement fields are important to analyze deformation, which is associated with functional and material tissue properties often used as indicators of health. Magnetic resonance imaging (MRI) techniques like DENSE and image registration methods like Hyperelastic Warping have been used to produce pixel-level deformation fields that are desirable in high-resolution analysis. However, DENSE can be complicated by challenges associated with image phase unwrapping, in particular offset determination. On the other hand, Hyperelastic Warping can be hampered by low local image contrast. The current work proposes a novel approach for measuring tissue displacement with both DENSE and Hyperelastic Warping, incorporating physically accurate displacements obtained by the latter to improve phase characterization in DENSE. The validity of the proposed technique is demonstrated using numerical and physical phantoms, and in vivo small animal cardiac MRI.

Magnetic resonance imaging assessment of mechanical interactions between human lower leg muscles in vivo

Evidence on epimuscular myofascial force transmission (EMFT) was shown for undissected muscle in situ. We hypothesize that global length changes of gastrocnemius muscle-tendon complex in vivo will cause sizable and heterogeneous local strains within all muscles of the human lower leg. Our goal is to test this hypothesis. A method was developed and validated using high-resolution 3D magnetic resonance image sets and Demons nonrigid registration algorithm for performing large deformation analyses. Calculation of strain tensors per voxel in human muscles in vivo allowed quantifying local heterogeneous tissue deformations and volume changes. After hip and knee movement (Δ knee angle ≈ 25 deg) but without any ankle movement, local lengthening within m. gastrocnemius was shown to occur simultaneously with local shortening (maximally by +34.2% and -32.6%, respectively) at different locations. Moreover, similar local strains occur also within other muscles, despite being kept at constant muscle-tendon complex length. This is shown for synergistic m. soleus and deep flexors, as well as for antagonistic anterior crural and peroneal muscle groups: minimum peak lengthening and shortening equaled 23.3% and 25.54%, respectively despite global isometric conditions. These findings confirm our hypothesis and show that in vivo, muscles are in principle not independent mechanically.

Automated motion estimation for 2-D cine DENSE MRI

Cine displacement encoding with stimulated echoes (DENSE) is a magnetic resonance (MR) method that directly encodes tissue displacement into MR phase images. This technique has successfully interrogated many forms of tissue motion, but is most commonly used to evaluate cardiac mechanics. Currently, motion analysis from cine DENSE images requires manually delineated anatomical structures. An automated analysis would improve measurement throughput, simplify data interpretation, and potentially access important physiological information during the MR exam. In this paper, we present the first fully automated solution for the estimation of tissue motion and strain from 2-D cine DENSE data. Results using both simulated and human cardiac cine DENSE data indicate good agreement between the automated algorithm and the standard semi-manual analysis method.

MRI of weight bearing and movement

Conventional, static magnetic resonance imaging (MRI) is able to provide a vast amount of information regarding the anatomy and pathology of the musculoskeletal system. However, patients, especially those whose pain is position dependent or elucidated by movement, may benefit from more advanced imaging techniques that allow for the acquisition of functional information. This manuscript reviews a variety of advancements in MRI techniques that are used to image the musculoskeletal system dynamically, while in motion or under load. The methodologies, advantages and drawbacks of stress MRI, cine-phase contrast MRI and real-time MRI are discussed as each has helped to advance the field by providing a scientific basis for understanding normal and pathological musculoskeletal anatomy and function. Advancements in dynamic MR imaging will certainly lead to improvements in the understanding, prevention, diagnosis and treatment of musculoskeletal disorders. It is difficult to anticipate that dynamic MRI will replace conventional MRI, however, dynamic MRI may provide additional valuable information to findings of conventional MRI.

Activation and aponeurosis morphology affect in vivo muscle tissue strains near the myotendinous junction

Hamstring strain injury is one of the most common injuries in athletes, particularly for sports that involve high speed running. The aims of this study were to determine whether muscle activation and internal morphology influence in vivo muscle behavior and strain injury susceptibility. We measured tissue displacement and strains in the hamstring muscle injured most often, the biceps femoris long head muscle (BFLH), using cine DENSE dynamic magnetic resonance imaging. Strain measurements were used to test whether strain magnitudes are (i) larger during active lengthening than during passive lengthening and (ii) larger for subjects with a relatively narrow proximal aponeurosis than a wide proximal aponeurosis. Displacement color maps showed higher tissue displacement with increasing lateral distance from the proximal aponeurosis for both active lengthening and passive lengthening, and higher tissue displacement for active lengthening than passive lengthening. First principal strain magnitudes were averaged in a 1cm region near the myotendinous junction, where injury is most frequently observed. It was found that strains are significantly larger during active lengthening (0.19 SD 0.09) than passive lengthening (0.13 SD 0.06) (p<0.05), which suggests that elevated localized strains may be a mechanism for increased injury risk during active as opposed to passive lengthening. First principal strains were higher for subjects with a relatively narrow aponeurosis width (0.26 SD 0.15) than wide (0.14 SD 0.04) (p<0.05). This result suggests that athletes who have BFLH muscles with narrow proximal aponeuroses may have an increased risk for BFLH strain injuries.

Clinical application of diffusion tensor magnetic resonance imaging in skeletal muscle

Diffusion tensor magnetic resonance imaging (DTI) is increasingly applied in the detection and characterization of skeletal muscle. This promising technique has aroused much enthusiasm and generated high expectations, because it is able to provide some specific information of skeletal muscle that is not available from other imaging modalities. Compared with conventional MRI, DTI could reconstruct the trajectories of skeletal muscle fibers. It makes it possible to non-invasively detect several physiological values (diffusion values), like fractional anisotropy (FA) and apparent diffusion coefficient (ADC), which have a great association with the muscle physiology and pathology. Furthermore, other advantages of DTI are the capability of investigating the muscle biomechanics and also investigate the pathological condition of skeletal muscle. Finally, several challenges, which limit the wide application of DTI in skeletal muscle, were discussed. It is believed that this review may arouse in-depth studies on the clinical application of DTI in skeletal muscle in future.

Combined diffusion and strain tensor MRI reveals a heterogeneous, planar pattern of strain development during isometric muscle contraction

The purposes of this study were to create a three-dimensional representation of strain during isometric contraction in vivo and to interpret it with respect to the muscle fiber direction. Diffusion tensor MRI was used to measure the muscle fiber direction of the tibialis anterior (TA) muscle of seven healthy volunteers. Spatial-tagging MRI was used to measure linear strains in six directions during separate 50% maximal isometric contractions of the TA. The strain tensor (E) was computed in the TA's deep and superficial compartments and compared with the respective diffusion tensors. Diagonalization of E revealed a planar strain pattern, with one nonzero negative strain (ε(N)) and one nonzero positive strain (ε(P)); both strains were larger in magnitude (P < 0.05) in the deep compartment [ε(N) = -40.4 ± 4.3%, ε(P) = 35.1 ± 3.5% (means ± SE)] than in the superficial compartment (ε(N) = -24.3 ± 3.9%, ε(P) = 6.3 ± 4.9%). The principal shortening direction deviated from the fiber direction by 24.0 ± 1.3° and 39.8 ± 6.1° in the deep and superficial compartments, respectively (P < 0.05, deep vs. superficial). The deviation of the shortening direction from the fiber direction was due primarily to the lower angle of elevation of the shortening direction over the axial plane than that of the fiber direction. It is concluded that three-dimensional analyses of strain interpreted with respect to the fiber architecture are necessary to characterize skeletal muscle contraction in vivo. The deviation of the principal shortening direction from the fiber direction may relate to intramuscle variations in fiber length and pennation angle.

Dynamic imaging of skeletal muscle contraction in three orthogonal directions

In this study, a multidimensional strain estimation method using biplane ultrasound is presented to assess local relative deformation (i.e., local strain) in three orthogonal directions in skeletal muscles during induced and voluntary contractions. The method was tested in the musculus biceps brachii of five healthy subjects for three different types of muscle contraction: 1) excitation of the muscle with a single electrical pulse via the musculocutaneous nerve, resulting in a so-called “twitch” contraction; 2) a train of five pulses at 10 Hz and 20 Hz, respectively, to obtain a submaximum tetanic contraction; and 3) voluntary contractions at 30, 60, and 100% of maximum contraction force. Results show that biplane ultrasound strain imaging is feasible. The method yielded adequate performance using the radio frequency data in tracking the tissue motion and enabled the measurement of local deformation in both the vertical direction (orthogonal to the arm) and in the horizontal directions (parallel and perpendicular to direction of the arm) in two orthogonal cross sections of the muscle. The twitch experiments appeared to be reproducible in all three directions, and high strains in vertical (25 to 30%) and horizontal (−20% to −10%) directions were measured. Visual inspection of both the ultrasound data, as well as the strain data, revealed a relaxation that was significantly slower than the force decay. The pulse train experiments nicely illustrated the performance of our technique: 1) similar patterns of force and strain waveforms were found; and 2) each stimulation frequency yielded a different strain pattern, e.g., peak vertical strain was 40% during 10-Hz stimulation and 60% during 20-Hz stimulation. The voluntary contraction patterns were found to be both practically feasible and reproducible, which will enable muscles and more natural contraction patterns to be examined without the need of electrical stimulation.

The influence of prior hamstring injury on lengthening muscle tissue mechanics

Hamstring strain injuries often occur near the proximal musculotendon junction (MTJ) of the biceps femoris. Post-injury remodeling can involve scar tissue formation, which may alter contraction mechanics and influence re-injury risk. The purpose of this study was to assess the affect of prior hamstring strain injury on muscle tissue displacements and strains during active lengthening contractions. Eleven healthy and eight subjects with prior biceps femoris injuries were tested. All previously injured subjects had since returned to sport and exhibited evidence of residual scarring along the proximal aponeurosis. Subjects performed cyclic knee flexion-extension on an MRI-compatible device using elastic and inertial loads, which induced active shortening and lengthening contractions, respectively. CINE phase-contrast imaging was used to measure tissue velocities within the biceps femoris during these tasks. Numerical integration of the velocity information was used to estimate two-dimensional tissue displacement and strain fields during muscle lengthening. The largest tissue motion was observed along the distal MTJ, with the active lengthening muscle exhibiting significantly greater and more homogeneous tissue displacements. First principal strain magnitudes were largest along the proximal MTJ for both loading conditions. The previously injured subjects exhibited less tissue motion and significantly greater strains near the proximal MTJ. We conclude that localized regions of high tissue strains during active lengthening contractions may predispose the proximal biceps femoris to injury. Furthermore, post-injury remodeling may alter the in-series stiffness seen by muscle tissue and contribute to the relatively larger localized tissue strains near the proximal MTJ, as was observed in this study.

A Magnetic Resonance-Compatible Loading Device for Dynamically Imaging Shortening and Lengthening Muscle Contraction Mechanics

The purpose of this study was to design and test a magnetic resonance (MR)-compatible device to induce either shortening or lengthening muscle contractions for use during dynamic MR imaging. The proposed device guides the knee through cyclic flexion-extension, while either elastic or inertial loads are imposed on the hamstrings. Ten subjects were tested in a motion capture laboratory to evaluate the repeatability of limb motion and imposed loads. Image data were subsequently obtained for all ten subjects using cine phase contrast imaging. Subjects achieved ~30 deg of knee joint motion, with individual subjects remaining within ~1 deg of their average motion across 56 repeated cycles. The maximum hamstring activity and loading occurred when the knee was flexed for the elastic loading condition (shortening contraction), and extended for the inertial loading condition (lengthening contraction). Repeat MR image acquisitions of the same loading condition resulted in similar tissue velocities, while spatial variations in velocity data were clearly different between loading conditions. The proposed device can enable dynamic imaging of the muscle under different types of loads, which has the potential to improve our understanding of basic muscle mechanics, identify potential causes of muscle injury, and provide a basis for quantitatively assessing injury effects at the tissue level. Slight modifications to the device design and/or subject positioning could allow for imaging of the quadriceps or the knee.

Wave attenuation as a measure of muscle quality as measured by magnetic resonance elastography: initial results

Advances in imaging technologies such as magnetic resonance elastography (MRE) have allowed researchers to gain insights into muscle function in vivo. MRE has been used to examine healthy and diseased muscle by calculating shear modulus. However, additional information can be measured from visualizing a mechanical wave as it passes through a tissue. One such measurable quantity is wave attenuation. The purpose of this study was to determine if a simple measure of wave attenuation could be used to distinguish between healthy and diseased muscle. Twenty seven subjects (14 healthy controls, 7 hyperthyroid myopathy patients, 6 myositis patients) participated in this study. Wave amplitude was determined along a linear profile through the center of the muscle, and an exponential decay curve was fit to the data. This measure was able to find significant differences in attenuation between healthy and diseased muscle. Furthermore, four hyperthyroid myopathy subjects who were tested following treatment all showed improvement by this measure. A likely reason for patients with hyperthyroid myopathy and myositis behaving similarly is that this measurement may reflect similar changes in the muscle extracellular matrix. In addition to modulus, attenuation seems to be an important parameter to measure in skeletal muscle. Further research is needed to investigate other potential measures of attenuation as well as examining other potential measures that can be found from visualizing wave propagation. Future studies should also include muscle biopsies to confirm that the changes seen are as a result of changes in extracellular matrix structure.

Longitudinal strain estimation in incompressible cylindrical tissues from magnetic resonance imaging

In this paper, we present a simple approach for estimating the average longitudinal strains from models reconstructed from medical images. It can be used for many incompressible generalized cylindrical tissues, such as tendons, ligaments, and fusiform muscles; the major deformation directions of these soft tissues are along the longitudinal axes. The method is especially useful when pre- and post-deformation tissue correspondences are difficult to establish directly from images for various reasons, such as insufficient image resolution, homogenous image intensity, and noise. Incompressibility, which is accepted as a good approximation for soft tissues, is exploited as a constraint on the tissue deformation. Experiments with Magnetic Resonance Imaging (MRI) of tissue phantoms and computer simulations show that the method is accurate and practical even in the presence of noise. Finally, we demonstrate the usefulness of our approach on studying extraocular muscle deformation.