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.</p

Model-based reconstructions for intravoxel incoherent motion and diffusion-tensor imaging parameter map estimations

PURPOSE: Intravoxel incoherent motion (IVIM) imaging and diffusion-tensor imaging (DTI) facilitate non-invasive quantification of tissue perfusion and diffusion. Both are promising biomarkers in various diseases and a combined acquisition is therefore desirable. This comes with challenges, including noisy parameter maps and long scan times, especially for the perfusion fraction f and pseudo-diffusion coefficient D*. A model-based reconstruction has the potential to overcome these challenges. As a first step, our goal was to develop a model-based reconstruction framework for IVIM and combined IVIM-DTI parameter estimation. METHODS: The IVIM and IVIM-DTI models were implemented in the PyQMRI model-based reconstruction framework and validated with simulations and in vivo data. Commonly used voxel-wise nonlinear least-squares fitting was used as reference. Simulations with the IVIM and IVIM-DTI models were performed with 100 noise realizations to assess accuracy and precision. Diffusion-weighted data were acquired for IVIM-reconstruction in the liver (n=5) as well as for IVIM-DTI in the kidneys (n=5) and lower-leg muscles (n=6) of healthy volunteers. The median and interquartile range (IQR) values of the IVIM and IVIM-DTI parameters were compared to assess bias and precision. RESULTS: With model-based reconstruction, the parameter maps exhibit less noise, which was most pronounced in the f and D* maps both in the simulations and in vivo. The bias values in the simulations were comparable between model-based reconstruction and the reference method. The IQR was lower with model-based reconstruction compared to the reference for all parameters. CONCLUSION: Model-based reconstruction is feasible for IVIM and IVIM-DTI and improves the precision of the parameter estimates, particularly for f and D* maps

Multi-parametric MR in Becker muscular dystrophy patients

Quantitative MRI and MRS of muscle are increasingly being used to measure individual pathophysiological processes in Becker muscular dystrophy (BMD). In particular, muscle fat fraction was shown to be highly associated with functional tests in BMD. However, the muscle strength per unit of contractile cross-sectional area is lower in patients with BMD compared with healthy controls. This suggests that the quality of the non-fat-replaced (NFR) muscle tissue is lower than in healthy controls. Consequently, a measure that reflects changes in muscle tissue itself is needed. Here, we explore the potential of water T2 relaxation times, diffusion parameters and phosphorus metabolic indices as early disease markers in patients with BMD. For this purpose, we examined these measures in fat-replaced (FR) and NFR lower leg muscles in patients with BMD and compared these values with those in healthy controls. Quantitative proton MRI (three-point Dixon, multi-spin-echo and diffusion-weighted spin-echo echo planar imaging) and 2D chemical shift imaging 31P MRS data were acquired in 24 patients with BMD (age 18.8-66.2 years) and 13 healthy controls (age 21.3-63.6 years). Muscle fat fractions, phosphorus metabolic indices, and averages and standard deviations (SDs) of the water T2 relaxation times and diffusion tensor imaging (DTI) parameters were assessed in six individual leg muscles. Phosphodiester levels were increased in the NFR and FR tibialis anterior, FR peroneus and FR gastrocnemius lateralis muscles. No clear pattern was visible for the other metabolic indices. Increased T2 SD was found in the majority of FR muscles compared with NFR and healthy control muscles. No differences in average water T2 relaxation times or DTI indices were found between groups. Overall, our results indicate that primarily muscles that are further along in the disease process showed increases in T2 heterogeneity and changes in some metabolic indices. No clear differences were found for the DTI indices between groups

A diffusion tensor-based method facilitating volumetric assessment of fiber orientations in skeletal muscle

Background The purpose of this study was to develop a DTI-based method to quantitatively assess fiber angles and changes therein in leg muscles in order to facilitate longitudinal studies on muscle fiber architectural adaptations in healthy subjects. Methods The upper legs of five volunteers were scanned twice on the same day. The right lower legs of five volunteers were scanned twice with the ankle in three positions, i.e. -15° dorsiflexion, 0° neutral position, and 30° plantarflexion. The MRI protocols consisted of a noise scan, a 3- point mDixon scan and a DTI scan. Fiber-angle color maps were generated for four muscles in the upper legs and two muscles in the lower leg. Voxel-wise fiber angles (θ) were calculated from the angle between the principal eigenvector of the diffusion tensor and a reference line defined between the origo and insertion points of each muscle. Bland-Altman analysis, intraclass correlation coefficient (ICC), coefficient of variation (CV%), minimal detectable change (MDC), standard error (SE) and Friedman test were used for assessing the feasibility of this method and in order to have an indication of the repeatability and the sensitivity. Results Bland-Altman analysis showed good repeatability (CV%<10 and 0.7≤ICC≤0.9) with exception of the Tibialis Anterior (TA) muscle in dorsiflexion position(CV%: 12.2) and the Semitendinosus (ST) muscle (left leg) (CV%: 11.4). The best repeatability metrics were found for the SOL muscle in neutral position (CV%: 2.6). Changes in average θ in TA and SOL with ankle positions were observed in accordance with expected agonist and antagonist functions of both muscles. For example, for the anterior left compartment the change in fiber angle Δθ with respect to the neutral position Δθ = -1.6° ± 0.8° and 2.2° ± 2.8° (p = 0.008), for dorsiflexion and plantarflexion, respectively. Conclusion Our method facilitates fast inspection and quantification of muscle fiber angles in the lower and upper leg muscles in rest and detection of changes in lower-leg muscle fiber angles with varying ankle angles

A diffusion tensor-based method facilitating volumetric assessment of fiber orientations in skeletal muscle

Background The purpose of this study was to develop a DTI-based method to quantitatively assess fiber angles and changes therein in leg muscles in order to facilitate longitudinal studies on muscle fiber architectural adaptations in healthy subjects. Methods The upper legs of five volunteers were scanned twice on the same day. The right lower legs of five volunteers were scanned twice with the ankle in three positions, i.e. -15° dorsiflexion, 0° neutral position, and 30° plantarflexion. The MRI protocols consisted of a noise scan, a 3- point mDixon scan and a DTI scan. Fiber-angle color maps were generated for four muscles in the upper legs and two muscles in the lower leg. Voxel-wise fiber angles (θ) were calculated from the angle between the principal eigenvector of the diffusion tensor and a reference line defined between the origo and insertion points of each muscle. Bland-Altman analysis, intraclass correlation coefficient (ICC), coefficient of variation (CV%), minimal detectable change (MDC), standard error (SE) and Friedman test were used for assessing the feasibility of this method and in order to have an indication of the repeatability and the sensitivity. Results Bland-Altman analysis showed good repeatability (CV%<10 and 0.7≤ICC≤0.9) with exception of the Tibialis Anterior (TA) muscle in dorsiflexion position(CV%: 12.2) and the Semitendinosus (ST) muscle (left leg) (CV%: 11.4). The best repeatability metrics were found for the SOL muscle in neutral position (CV%: 2.6). Changes in average θ in TA and SOL with ankle positions were observed in accordance with expected agonist and antagonist functions of both muscles. For example, for the anterior left compartment the change in fiber angle Δθ with respect to the neutral position Δθ = -1.6° ± 0.8° and 2.2° ± 2.8° (p = 0.008), for dorsiflexion and plantarflexion, respectively. Conclusion Our method facilitates fast inspection and quantification of muscle fiber angles in the lower and upper leg muscles in rest and detection of changes in lower-leg muscle fiber angles with varying ankle angles

The repeatability of bilateral diffusion tensor imaging (DTI) in the upper leg muscles of healthy adults

Objectives: Assessment of the repeatability of diffusion parameter estimations in the upper leg muscles of healthy adults over the time course of 2 weeks, from a simultaneous bilateral upper leg DTI measurement. Methods: SE-EPI DTI datasets were acquired at 3 T in the upper legs of 15 active adults at a time interval of 2 weeks. ROIs were manually drawn for four quadriceps and three hamstring muscles of both legs. The following DTI parameters were analyzed: 1st, 2nd, and 3rd eigenvalue (λ1, λ2, and λ3), mean diffusivity (MD), and fractional anisotropy (FA). DTI parameters per muscle were calculated with and without intravoxel incoherent motion (IVIM) correction together with SNR levels per muscle. Bland-Altman plots and within-subject coefficient of variation (wsCV) were calculated. Left-right differences between muscles were assessed. Results: The Bland-Altman analysis showed good repeatability of all DTI parameters except FA for both the IVIM-corrected and standard data. wsCV values show that MD has the highest repeatability (4.5% IVIM; 5.6% standard), followed by λ2 (4.9% IVIM; 5.5% standard), λ1 (5.3% IVIM; 7.5% standard), and λ3 (5.7% IVIM; 5.7% standard). wsCV values of FA were 15.2% for the IVIM-corrected data and 13.9% for the standard analysis. The SNR (41.8 ± 16.0 right leg, 41.7 ± 17.1 left leg) and wsCV values were similar for the left and right leg and no left-right bias was detected. Conclusions: Repeatability was good for standard DTI data and slightly better for IVIM-corrected DTI data. Our protocol is suitable for DTI of the upper legs with overall good SNR. Key Points: • The presented DTI protocol is repeatable and therefore suitable for bilateral DT imaging of the upper legs. • Additional B1+calibrations improve SNR and repeatability. • Correcting for perfusion effects improves repeatability

Multi‐parametric MR in Becker muscular dystrophy patients

Quantitative MRI and MRS of muscle are increasingly being used to measure individual pathophysiological processes in Becker muscular dystrophy (BMD). In particular, muscle fat fraction was shown to be highly associated with functional tests in BMD. However, the muscle strength per unit of contractile cross-sectional area is lower in patients with BMD compared with healthy controls. This suggests that the quality of the non-fat-replaced (NFR) muscle tissue is lower than in healthy controls. Consequently, a measure that reflects changes in muscle tissue itself is needed. Here, we explore the potential of waterT(2)relaxation times, diffusion parameters and phosphorus metabolic indices as early disease markers in patients with BMD. For this purpose, we examined these measures in fat-replaced (FR) and NFR lower leg muscles in patients with BMD and compared these values with those in healthy controls. Quantitative proton MRI (three-point Dixon, multi-spin-echo and diffusion-weighted spin-echo echo planar imaging) and 2D chemical shift imaging(31)P MRS data were acquired in 24 patients with BMD (age 18.8-66.2 years) and 13 healthy controls (age 21.3-63.6 years). Muscle fat fractions, phosphorus metabolic indices, and averages and standard deviations (SDs) of the waterT(2)relaxation times and diffusion tensor imaging (DTI) parameters were assessed in six individual leg muscles. Phosphodiester levels were increased in the NFR and FR tibialis anterior, FR peroneus and FR gastrocnemius lateralis muscles. No clear pattern was visible for the other metabolic indices. IncreasedT(2)SD was found in the majority of FR muscles compared with NFR and healthy control muscles. No differences in average waterT(2)relaxation times or DTI indices were found between groups. Overall, our results indicate that primarily muscles that are further along in the disease process showed increases inT(2)heterogeneity and changes in some metabolic indices. No clear differences were found for the DTI indices between groups.Neuro Imaging Researc

A MATLAB toolbox for muscle diffusion-tensor MRI tractography

Diffusion-tensor MRI fiber tractography has been used to reconstruct skeletal muscle architecture, but remains a specialized technique using custom-written data processing routines. In this work, we describe the public release of a software toolbox having the following design objectives: accomplish the pre-processing tasks of file input, image registration, denoising, and diffusion-tensor calculation; allow muscle-specific methods for defining seed points; make fiber-tract architectural measurements referenced to tendinous structures; visualize fiber tracts and other muscle structures of interest; analyze the goodness of outcomes; and provide a programming structure that allows the addition of new capabilities in future versions. The proper function of the code was verified using simulated datasets. The toolbox capabilities for characterizing human muscle structure in vivo were demonstrated in a case study. These capabilities included measurements of muscle morphology; contractile and non-contractile tissue volumes; fiber-tract length, pennation angle, curvature; and the physiological cross-sectional area,. The free public release of this software is a first step in creating of a community of users who use these tools in studies of muscle physiology and biomechanics. Users may further contribute to code development. Along with simulated and actual datasets for benchmarking, these tools will further create mechanisms for enhancing scientific rigor and developing and validating new code features. Planned future developments include additional options for image pre-processing, development of a graphical user interface, analysis of architectural patterns during muscle contraction, and integration of functional imaging data

Effect of two eccentric hamstring exercises on muscle architectural characteristics assessed with diffusion tensor MRI

Objectives: To evaluate the effect of a Nordic hamstring exercise or Diver hamstring exercise intervention on biceps femoris long head, semitendinosus and semimembranosus muscle's fascicle length and orientation through diffusion tensor imaging (DTI) with magnetic resonance imaging. Methods: In this three-arm, single-center, randomized controlled trial, injury-free male basketball players were randomly assigned to a Nordic, Diver hamstring exercise intervention or control group. The primary outcome was the DTI-derived fascicle length and orientation of muscles over 12 weeks. Results: Fifty-three participants were included for analysis (mean age 22 ± 7 years). Fascicle length in the semitendinosus over 12 weeks significantly increased in the Nordic-group (mean [M]: 20.8 mm, 95% confidence interval [95% CI]: 7.8 to 33.8) compared with the Control-group (M: 0.9 mm, 95% CI: −7.1 to 8.9), mean between-groups difference: 19.9 mm, 95% CI: 1.9 to 37.9, p = 0.026. Fascicle orientation in the biceps femoris long head over 12 weeks significantly decreased in the Diver-group (mean: -2.6°, 95% CI: −4.1 to −1.0) compared with the Control-group (mean: −0.2°, 95% CI: −1.4 to 1.0), mean between-groups difference: -2.4°, 95% CI: −4.7 to −0.1, p = 0.039. Conclusion: The Nordic hamstring exercise intervention did significantly increase the fascicle length of the semitendinosus and the Diver hamstring exercise intervention did significantly change the orientation of fascicles of the biceps femoris long head. As both exercises are complementary to each other, the combination is relevant for preventing hamstring injuries