The impact of diffusion tensor imaging tractography settings on muscle fascicle architecture and diffusion parameter estimates: Tract length, completion, and curvature are most sensitive to tractography settings.

Diffusion-tensor (DT)-MRI tractography provides information about properties relevant to muscle health and function, including estimates of architectural properties such as fascicle length, pennation angle, and curvature and diffusion properties such as mean diffusivity (MD) and fractional anisotropy (FA). Tractography settings, including integration algorithms, thresholds for early tract termination, and tract smoothing approaches, impact the accuracy of the muscle property estimates. However, muscle DT-MRI tractography is performed using a variety of these settings, complicating comparisons between different studies. The effects of different tractography settings on muscle architecture estimates have not been fully explored, and optimized settings for muscle tractography have not yet been determined. We examined the influence of integration algorithm and termination check settings combined with a range of step sizes, termination criteria, and smoothing polynomial orders on tract characteristics, completion/reason for termination, and goodness of fit between fiber tracts and smoothing polynomials using 3-T DT-MR images of the lower leg muscles of seven healthy adults. We found that tract length and completion were highly sensitive to strict FA and intersegment angle thresholds (25%-69% reduction in complete fiber tracts from lowest to highest minimum FA threshold and 11%-36% reduction from highest to lowest intersegment angle threshold). Higher order polynomials (third and fourth order vs. second order) better fit the muscle fiber trajectories, but curvature estimates were highly sensitive to smoothing polynomial order (3.9-6.6 m-1 increase for second- vs. fourth-order fitting polynomials). Step size impacted curvature estimates, albeit to a lesser degree. Integration algorithm had little impact, and mean pennation angle, and tract-based FA and MD, were relatively insensitive to all parameters. The results demonstrate which muscle diffusion measures and architectural estimates are most sensitive to varying tractography settings and support the need for consistent reporting of tractography details to aid interpretation and comparison of results between studies.

Exercise MR of Skeletal Muscles, the Heart, and the Brain.

Magnetic resonance (MR) imaging (MRI) is routinely used to evaluate organ morphology and pathology in the human body at rest or in combination with pharmacological stress as an exercise surrogate. With MR during actual physical exercise, we can assess functional characteristics of tissues and organs under real-life stress conditions. This is particularly relevant in patients with limited exercise capacity or exercise intolerance, and where complaints typically present only during physical activity, such as in neuromuscular disorders, inherited metabolic diseases, and heart failure. This review describes practical and physiological aspects of exercise MR of skeletal muscles, the heart, and the brain. The acute effects of physical exercise on these organs are addressed in the light of various dynamic quantitative MR readouts, including phosphorus-31 MR spectroscopy (31P-MRS) of tissue energy metabolism, phase-contrast MRI of blood flow and muscle contraction, real-time cine MRI of cardiac performance, and arterial spin labeling MRI of muscle and brain perfusion. Exercise MR will help advancing our understanding of underlying mechanisms that contribute to exercise intolerance, which often proceed structural and anatomical changes in disease. Its potential to detect disease-driven alterations in organ function, perfusion, and metabolism under physiological stress renders exercise MR stress testing a powerful noninvasive imaging modality to aid in disease diagnosis and risk stratification. Although not yet integrated in most clinical workflows, and while some applications still require thorough validation, exercise MR has established itself as a comprehensive and versatile modality for characterizing physiology in health and disease in a noninvasive and quantitative way. EVIDENCE LEVEL: 5 TECHNICAL EFFICACY: Stage 1.

A deep learning approach for fast muscle water T2 mapping with subject specific fat T2 calibration from multi-spin-echo acquisitions.

This work presents a deep learning approach for rapid and accurate muscle water T2 with subject-specific fat T2 calibration using multi-spin-echo acquisitions. This method addresses the computational limitations of conventional bi-component Extended Phase Graph fitting methods (nonlinear-least-squares and dictionary-based) by leveraging fully connected neural networks for fast processing with minimal computational resources. We validated the approach through in vivo experiments using two different MRI vendors. The results showed strong agreement of our deep learning approach with reference methods, summarized by Lin's concordance correlation coefficients ranging from 0.89 to 0.97. Further, the deep learning method achieved a significant computational time improvement, processing data 116 and 33 times faster than the nonlinear least squares and dictionary methods, respectively. In conclusion, the proposed approach demonstrated significant time and resource efficiency improvements over conventional methods while maintaining similar accuracy. This methodology makes the processing of water T2 data faster and easier for the user and will facilitate the utilization of the use of a quantitative water T2 map of muscle in clinical and research studies.

A registration strategy to characterize DTI-observed changes in skeletal muscle architecture due to passive shortening.

Skeletal muscle architecture is a key determinant of muscle function. Architectural properties such as fascicle length, pennation angle, and curvature can be characterized using Diffusion Tensor Imaging (DTI), but acquiring these data during a contraction is not currently feasible. However, an image registration-based strategy may be able to convert muscle architectural properties observed at rest to their contracted state. As an initial step toward this long-term objective, the aim of this study was to determine if an image registration strategy could be used to convert the whole-muscle average architectural properties observed in the extended joint position to those of a flexed position, following passive rotation. DTI and high-resolution fat/water scans were acquired in the lower leg of seven healthy participants on a 3T MR system in +20° (plantarflexion) and -10° (dorsiflexion) foot positions. The diffusion and anatomical images from the two positions were used to propagate DTI fiber-tracts from seed points along a mesh representation of the aponeurosis of fiber insertion. The -10° and +20° anatomical images were registered and the displacement fields were used to transform the mesh and fiber-tracts from the +20° to the -10° position. Student's paired t-tests were used to compare the mean architectural parameters between the original and transformed fiber-tracts. The whole-muscle average fiber-tract length, pennation angle, curvature, and physiological cross-sectional areas estimates did not differ significantly. DTI fiber-tracts in plantarflexion can be transformed to dorsiflexion position without significantly affecting the average architectural characteristics of the fiber-tracts. In the future, a similar approach could be used to evaluate muscle architecture in a contracted state.

Accelerated IVIM-corrected DTI in acute hamstring injury: towards a clinically feasible acquisition time.

BACKGROUND: Intravoxel incoherent motion (IVIM)-corrected diffusion tensor imaging (DTI) potentially enhances return-to-play (RTP) prediction after hamstring injuries. However, the long scan times hamper clinical implementation. We assessed accelerated IVIM-corrected DTI approaches in acute hamstring injuries and explore the sensitivity of the perfusion fraction (f) to acute muscle damage. METHODS: Athletes with acute hamstring injury received DTI scans of both thighs < 7 days after injury and at RTP. For a subset, DTI scans were repeated with multiband (MB) acceleration. Data from standard and MB-accelerated scans were fitted with standard and accelerated IVIM-corrected DTI approach using high b-values only. Segmentations of the injury and contralateral healthy muscles were contoured. The fitting methods as well as the standard and MB-accelerated scan were compared using linear regression analysis. For sensitivity to injury, Δ(injured minus healthy) DTI parameters between the methods and the differences between injured and healthy muscles were compared (Wilcoxon signed-rank test). RESULTS: The baseline dataset consisted of 109 athletes (16 with MB acceleration); 64 of them received an RTP scan (8 with MB acceleration). Linear regression of the standard and high-b DTI fitting showed excellent agreement. With both fitting methods, standard and MB-accelerated scans were comparable. Δ(injured minus healthy) was similar between standard and accelerated methods. For all methods, all IVIM-DTI parameters except f were significantly different between injured and healthy muscles. CONCLUSIONS: High-b DTI fitting with MB acceleration reduced the scan time from 11:08 to 3:40 min:s while maintaining sensitivity to hamstring injuries; f was not different between healthy and injured muscles. RELEVANCE STATEMENT: The accelerated IVIM-corrected DTI protocol, using fewer b-values and MB acceleration, reduced the scan time to under 4 min without affecting the sensitivity of the quantitative outcome parameters to hamstring injuries. This allows for routine clinical monitoring of hamstring injuries, which could directly benefit injury treatment and monitoring. KEY POINTS: • Combining high-b DTI-fitting and multiband-acceleration dramatically reduced by two thirds the scan time. • The accelerated IVIM-corrected DTI approaches maintained the sensitivity to hamstring injuries. • The IVIM-derived perfusion fraction was not sensitive to hamstring injuries.

Predicted effects of image acquisition and analysis conditions on DTMRI tractography-based muscle architecture estimates.

PURPOSE: To quantify the effects of the intrinsic signal pattern, image acquisition conditions, and data analysis conditions on diffusion-tensor MRI (DTMRI) tractography-based muscle architecture estimates using a sampling-reconstruction assessment framework. METHODS: Numerical models of muscles were constructed with realistic architectural properties. DTMRI signals were computed at signal-to-noise ratio (SNR) of 24-96 and common voxel sizes. Fiber tracking was performed, and the results were compared with the known architectural properties. RESULTS: SNR exerted the most significant impact on the outcome. The outcome variables approached asymptotes at SNR ≈ 54. Large in-plane voxel dimensions reduced the similarity between reconstructed fibers and the known architectural properties. Higher order polynomials helped reconstruct fibers with more complicated geometry but overfit noise for less complex geometries. The intrinsic fiber curvature also affected the robustness of polynomial smoothing to SNR. Other conditions, such as the fiber dimensionality, voxel aspect ratio, and slice thickness, did not affect the outcomes. CONCLUSION: SNR ≥ 54 is recommended for accurate muscle architecture characterization using DTMRI. Averaged across all simulated conditions, the greatest percent errors under SNR = 54 were -5.6% and -4.0% for the pennation angle and fiber-tract length estimates, respectively. For fiber tracts with intermediate intrinsic curvature, the greatest percent error for the curvature estimate was 9.8% for SNR = 54. Smaller in-plane voxel size (≤1.5 mm) is preferred to minimize the estimation error in architectural properties. If necessary, slice thickness may be adjusted within typical ranges to achieve sufficient SNR when slices are aligned near the fiber direction. Third-order polynomial fitting is appropriate for smoothing fiber tracts.

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.

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.

Diffusion-tensor magnetic resonance imaging captures increased skeletal muscle fibre diameters in Becker muscular dystrophy.

BACKGROUND: Becker muscular dystrophy (BMD) is an X-linked disorder characterized by slow, progressive muscle damage and muscle weakness. Hallmarks include fibre-size variation and replacement of skeletal muscle with fibrous and adipose tissues, after repeated cycles of regeneration. Muscle histology can detect these features, but the required biopsies are invasive, are difficult to repeat and capture only small muscle volumes. Diffusion-tensor magnetic resonance imaging (DT-MRI) is a potential non-invasive alternative that can calculate muscle fibre diameters when applied with the novel random permeable barrier model (RPBM). In this study, we assessed muscle fibre diameters using DT-MRI in BMD patients and healthy controls and compared these with histology. METHODS: We included 13 BMD patients and 9 age-matched controls, who underwent water-fat MRI and DT-MRI at multiple diffusion times, allowing RPBM parameter estimation in the lower leg muscles. Tibialis anterior muscle biopsies were taken from the contralateral leg in 6 BMD patients who underwent DT-MRI and from an additional 32 BMD patients and 15 healthy controls. Laminin and Sirius-red stainings were performed to evaluate muscle fibre morphology and fibrosis. Twelve ambulant patients from the MRI cohort underwent the North Star ambulatory assessment, and 6-min walk, rise-from-floor and 10-m run/walk functional tests. RESULTS: RPBM fibre diameter was significantly larger in BMD patients (P = 0.015): mean (SD) = 68.0 (25.3) µm versus 59.4 (19.2) µm in controls. Inter-muscle differences were also observed (P ≤ 0.002). Both inter- and intra-individual RPBM fibre diameter variability were similar between groups. Laminin staining agreed with the RPBM, showing larger median fibre diameters in patients than in controls: 72.5 (7.9) versus 63.2 (6.9) µm, P = 0.006. However, despite showing similar inter-individual variation, patients showed more intra-individual fibre diameter variability than controls-mean variance (SD) = 34.2 (7.9) versus 21.4 (6.9) µm, P < 0.001-and larger fibrosis areas: median (interquartile range) = 21.7 (5.6)% versus 14.9 (3.4)%, P < 0.001. Despite good overall agreement of RPBM and laminin fibre diameters, they were not associated in patients who underwent DT-MRI and muscle biopsy, perhaps due to lack of colocalization of DT-MRI with biopsy samples. CONCLUSIONS: DT-MRI RPBM metrics agree with histology and can quantify changes in muscle fibre size that are associated with regeneration without the need for biopsies. They therefore show promise as imaging biomarkers for muscular dystrophies.

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

Intravoxel incoherent motion (IVIM) imaging and diffusion tensor imaging (DTI) facilitate noninvasive 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. 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 the 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. With model-based reconstruction, the parameter maps exhibited 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 with the reference for all parameters. In 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.

Diffusion tensor imaging and quantitative T2 mapping to monitor muscle recovery following hamstring injury.

MRI examinations are accurate for diagnosing sports-related acute hamstring injuries. However, sensitive imaging methods for assessing recovery of these injuries are lacking. Diffusion tensor imaging (DTI) and quantitative T2 (qT2) mapping have both shown promise for assessing recovery of muscle micro trauma and exercise effects. The purpose of this study was to explore the potential of DTI and qT2 mapping for monitoring the muscle recovery processes after acute hamstring injury. In this prospective study, athletes with an acute hamstring injury underwent a 3-T MRI examination of the injured and contralateral hamstrings including DTI and qT2 measurements at three time points: (1) within 1 week after sustaining the injury, (2) 2 weeks after time point 1, and (3) return to play (RTP). A linear mixed model was used for time-effect analysis and paired t-tests for the detection of differences between injured and uninjured muscles. Forty-one athletes (age 27.8 ± 7 years; two females and 39 males) were included. Mean RTP time was 50 (range 12-169) days. A significant time effect was found for mean diffusivity, radial diffusivity, and the second and third eigenvalues (p ≤ 0.001) in the injured muscles. Fractional anisotropy (p = 0.40), first eigenvalue (p = 0.02), and qT2 (p = 0.61) showed no significant time effect. All DTI indices, except for fractional anisotropy, were significantly elevated compared with control muscles right after the injury (p < 0.001). Values normalized during the recovery period, with no significant differences between control and injured muscles at RTP (p values ranged from 0.08 to 0.51). Mean qT2 relaxation times in injured muscles were not significantly elevated compared with control muscles at any time point (p > 0.04). In conclusion, DTI can be used to monitor recovery after an acute hamstring injury. Future work should explore the potential of DTI indices to predict RTP and recovery times in athletes after an acute strain injury.

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.

Multi-parametric quantitative magnetic resonance imaging of the upper arm muscles of patients with spinal muscular atrophy.

Quantitative magnetic resonance imaging (qMRI) is frequently used to map the disease state and disease progression in the lower extremity muscles of patients with spinal muscular atrophy (SMA). This is in stark contrast to the almost complete lack of data on the upper extremity muscles, which are essential for carrying out daily activities. The aim of this study was therefore to assess the disease state in the upper arm muscles of patients with SMA in comparison with healthy controls by quantitative assessment of fat fraction, diffusion indices, and water T2 relaxation times, and to relate these measures to muscle force. We evaluated 13 patients with SMA and 15 healthy controls with a 3-T MRI protocol consisting of DIXON, diffusion tensor imaging, and T2 sequences. qMRI measures were compared between groups and related to muscle force measured with quantitative myometry. Fat fraction was significantly increased in all upper arm muscles of the patients with SMA compared with healthy controls and correlated negatively with muscle force. Additionally, fat fraction was heterogeneously distributed within the triceps brachii (TB) and brachialis muscle, but not in the biceps brachii muscle. Diffusion indices and water T2 relaxation times were similar between patients with SMA and healthy controls, but we did find a slightly reduced mean diffusivity (MD), λ1, and λ3 in the TB of patients with SMA. Furthermore, MD was positively correlated with muscle force in the TB of patients with SMA. The variation in fat fraction further substantiates the selective vulnerability of muscles. The reduced diffusion tensor imaging indices, along with the positive correlation of MD with muscle force, point to myofiber atrophy. Our results show the feasibility of qMRI to map the disease state in the upper arm muscles of patients with SMA. Longitudinal data in a larger cohort are needed to further explore qMRI to map disease progression and to capture the possible effects of therapeutic interventions.

Baseline fat fraction is a strong predictor of disease progression in Becker muscular dystrophy.

In Becker muscular dystrophy (BMD), muscle weakness progresses relatively slowly, with a highly variable rate among patients. This complicates clinical trials, as clinically relevant changes are difficult to capture within the typical duration of a trial. Therefore, predictors for disease progression are needed. We assessed if temporal increase of fat fraction (FF) in BMD follows a sigmoidal trajectory and whether fat fraction at baseline (FFbase) could therefore predict FF increase after 2 years (ΔFF). Thereafter, for two different MR-based parameters, we tested the additional predictive value to FFbase. We used 3-T Dixon data from the upper and lower leg, and multiecho spin-echo MRI and 7-T 31 P MRS datasets from the lower leg, acquired in 24 BMD patients (age: 41.4 [SD 12.8] years). We assessed the pattern of increase in FF using mixed-effects modelling. Subsequently, we tested if indicators of muscle damage like standard deviation in water T2 (stdT2 ) and the phosphodiester (PDE) over ATP ratio at baseline had additional value to FFbase for predicting ∆FF. The association between FFbase and ΔFF was described by the derivative of a sigmoid function and resulted in a peak ΔFF around 0.45 FFbase (fourth-order polynomial term: t = 3.7, p < .001). StdT2 and PDE/ATP were not significantly associated with ∆FF if FFbase was included in the model. The relationship between FFbase and ∆FF suggests a sigmoidal trajectory of the increase in FF over time in BMD, similar to that described for Duchenne muscular dystrophy. Our results can be used to identify muscles (or patients) that are in the fast progressing stage of the disease, thereby facilitating the conduct of clinical trials.

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.

Magnetic resonance reveals mitochondrial dysfunction and muscle remodelling in spinal muscular atrophy.

Genetic therapy has changed the prognosis of hereditary proximal spinal muscular atrophy, although treatment efficacy has been variable. There is a clear need for deeper understanding of underlying causes of muscle weakness and exercise intolerance in patients with this disease to further optimize treatment strategies. Animal models suggest that in addition to motor neuron and associated musculature degeneration, intrinsic abnormalities of muscle itself including mitochondrial dysfunction contribute to the disease aetiology. To test this hypothesis in patients, we conducted the first in vivo clinical investigation of muscle bioenergetics. We recruited 15 patients and 15 healthy age and gender-matched control subjects in this cross-sectional clinico-radiological study. MRI and 31P magnetic resonance spectroscopy, the modality of choice to interrogate muscle energetics and phenotypic fibre-type makeup, was performed of the proximal arm musculature in combination with fatiguing arm-cycling exercise and blood lactate testing. We derived bioenergetic parameter estimates including: blood lactate, intramuscular pH and inorganic phosphate accumulation during exercise, and muscle dynamic recovery constants. A linear correlation was used to test for associations between muscle morphological and bioenergetic parameters and clinico-functional measures of muscle weakness. MRI showed significant atrophy of triceps but not biceps muscles in patients. Maximal voluntary contraction force normalized to muscle cross-sectional area for both arm muscles was 1.4-fold lower in patients than in controls, indicating altered intrinsic muscle properties other than atrophy contributed to muscle weakness in this cohort. In vivo31P magnetic resonance spectroscopy identified white-to-red remodelling of residual proximal arm musculature in patients on the basis of altered intramuscular inorganic phosphate accumulation during arm-cycling in red versus white and intermediate myofibres. Blood lactate rise during arm-cycling was blunted in patients and correlated with muscle weakness and phenotypic muscle makeup. Post-exercise metabolic recovery was slower in residual intramuscular white myofibres in patients demonstrating mitochondrial ATP synthetic dysfunction in this particular fibre type. This study provides the first in vivo evidence in patients that degeneration of motor neurons and associated musculature causing atrophy and muscle weakness in 5q spinal muscular atrophy type 3 and 4 is aggravated by disproportionate depletion of myofibres that contract fastest and strongest. Our finding of decreased mitochondrial ATP synthetic function selectively in residual white myofibres provides both a possible clue to understanding the apparent vulnerability of this particular fibre type in 5q spinal muscular atrophy types 3 and 4 as well as a new biomarker and target for therapy.

Muscle architecture is associated with muscle fat replacement in Duchenne and Becker muscular dystrophies.

INTRODUCTION/AIMS: Duchenne and Becker muscular dystrophies (DMD and BMD, respectively) are characterized by fat replacement of different skeletal muscles in a specific temporal order. Given the structural role of dystrophin in skeletal muscle mechanics, muscle architecture could be important in the progressive pathophysiology of muscle degeneration. Therefore, the aim of this study was to assess the role of muscle architecture in the progression of fat replacement in DMD and BMD. METHODS: We assessed the association between literature-based leg muscle architectural characteristics and muscle fat fraction from 22 DMD and 24 BMD patients. Dixon-based magnetic resonance imaging estimates of fat fractions at baseline and 12 (only DMD) and 24 months were related to fiber length and physiological cross-sectional area (PCSA) using age-controlled linear mixed modeling. RESULTS: DMD and BMD muscles with long fibers and BMD muscles with large PCSAs were associated with increased fat fraction. The effect of fiber length was stronger in muscles with larger PCSA. DISCUSSION: Muscle architecture may explain the pathophysiology of muscle degeneration in dystrophinopathies, in which proximal muscles with a larger mass (fiber length × PCSA) are more susceptible, confirming the clinical observation of a temporal proximal-to-distal progression. These results give more insight into the mechanical role in the pathophysiology of muscular dystrophies. Ultimately, this new information can be used to help support the selection of current and the development of future therapies.

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.

Selection Approach to Identify the Optimal Biomarker Using Quantitative Muscle MRI and Functional Assessments in Becker Muscular Dystrophy.

OBJECTIVE: To identify the best quantitative fat-water MRI biomarker for disease progression of leg muscles in Becker muscular dystrophy (BMD) by applying a stepwise approach based on standardized response mean (SRM) over 24 months, correlations with baseline ambulatory tests, and reproducibility. METHODS: Dixon fat-water imaging was performed at baseline (n = 24) and 24 months (n = 20). Fat fractions (FF) were calculated for 3 center slices and the whole muscles for 19 muscles and 6 muscle groups. Contractile cross-sectional area (cCSA) was obtained from the center slice. Functional assessments included knee extension and flexion force and 3 ambulatory tests (North Star Ambulatory Assessment [NSAA], 10-meter run, 6-minute walking test). MRI measures were selected using SRM (≥0.8) and correlation with all ambulatory tests (ρ ≤ -0.8). Measures were evaluated based on intraclass correlation coefficient (ICC) and SD of the difference. Sample sizes were calculated assuming 50% reduction in disease progression over 24 months in a clinical trial with 1:1 randomization. RESULTS: Median whole muscle FF increased between 0.2% and 2.6% without consistent cCSA changes. High SRMs and strong functional correlations were found for 8 FF but no cCSA measures. All measures showed excellent ICC (≥0.999) and similar SD of the interrater difference. Whole thigh 3 center slices FF was the best biomarker (SRM 1.04, correlations ρ ≤ -0.81, ICC 1.00, SD 0.23%, sample size 59) based on low SD and acquisition and analysis time. CONCLUSION: In BMD, median FF of all muscles increased over 24 months. Whole thigh 3 center slices FF reduced the sample size by approximately 40% compared to NSAA.