Reproducibility analysis of diffusion tensor indices and fiber architecture of human calf muscles in vivo at 1.5 Tesla in neutral and plantarflexed ankle positions at rest

Consistency of muscle fibers directionality in human thigh derived from diffusion-weighted MRI

Objective.Diffusion-weighted MR imaging (DW-MRI) is known to quantify muscle fiber directionality and thus may be useful for radiotherapy target definition in sarcomas. Here, we investigate the variability of tissue anisotropy derived from diffusion tensor (DT) in the human thigh to establish the baseline parameters and protocols for DW-MRI acquisition for future studies in sarcoma patients.Approach.We recruited ten healthy volunteers to acquire diffusion-weighted MR images of the left and right thigh. DW-MRI data were used to reconstruct DT eigenvectors within each individual thigh muscle. Deviations of the principal eigenvector from its mean were calculated for different experimental conditions.Main results.Within the majority of muscles in most subjects, the mode of the histogram of the angular deviation of the principal eigenvector of the water DT from its muscle-averaged value did not exceed 20°. On average for all subjects, the mode ranged from 15° to 24°. Deviations much larger than 20° were observed in muscles far from the RF coil, including cases with significant amounts of subcutaneous fat and muscle deformation under its own weight.Significance.Our study is a robust characterization of angular deviations of muscle fiber directionality in the thigh as determined by DW-MRI. We show that an appropriate choice of experimental conditions reduces the variability of the observed directionality. Precise determination of tissue directionality will enable reproducible models of microscopic tumor spread, with future application in defining the clinical target volume for soft tissue sarcoma.

Assessment of myofiber microstructure changes due to atrophy and recovery with time-dependent diffusion MRI

Current clinical MRI evaluation of musculature largely focuses on nonquantitative assessments (including T1-, T2- and PD-weighted images), which may vary greatly between imaging systems and readers. This work aims to determine the efficacy of a quantitative approach to study the microstructure of muscles at the cellular level with the random permeable barrier model (RPBM) applied to time-dependent diffusion tensor imaging (DTI) for varying diffusion time. Patients (N = 15, eight males and seven females) with atrophied calf muscles due to immobilization of one leg in a nonweight-bearing cast, were enrolled after providing informed consent. Their calf muscles were imaged with stimulated echo diffusion for DTI, T1-mapping and RPBM modeling. Specifically, After cast removal, both calf muscles (atrophied and contralateral control leg) were imaged with MRI for all patients, with follow-up scans to monitor recovery of the atrophied leg for six patients after 4 and 8 weeks. We compare RPBM-derived microstructural metrics: myofiber diameter, a, and sarcolemma permeability, κ, along with macroscopic anatomical parameters (muscle cross-sectional area, fiber orientation, <θ>, and T1 relaxation). ROC analysis was used to compare parameters between control and atrophied muscle, while the Friedman test was used to evaluate the atrophied muscle longitudinally. We found that the RPBM framework enables measurement of microstructural parameters from diffusion time-dependent DTI, of which the myofiber diameter is a stronger predictor of intramuscular morphological changes than either macroscopic (anatomical) measurements or empirical diffusion parameters. This work demonstrates the potential of RPBM to assess pathological changes in musculature that seem undetectable with standard diffusion and anatomical MRI.

Diffusion Tensor Imaging of Skeletal Muscle Contraction Using Oscillating Gradient Spin Echo

Diffusion tensor imaging (DTI) measures water diffusion in skeletal muscle tissue and allows for muscle assessment in a broad range of neuromuscular diseases. However, current DTI measurements, typically performed using pulsed gradient spin echo (PGSE) diffusion encoding, are limited to the assessment of non-contracted musculature, therefore providing limited insight into muscle contraction mechanisms and contraction abnormalities. In this study, we propose the use of an oscillating gradient spin echo (OGSE) diffusion encoding strategy for DTI measurements to mitigate the effect of signal voids in contracted muscle and to obtain reliable diffusivity values. Two OGSE sequences with encoding frequencies of 25 and 50 Hz were tested in the lower leg of five healthy volunteers with relaxed musculature and during active dorsiflexion and plantarflexion, and compared with a conventional PGSE approach. A significant reduction of areas of signal voids using OGSE compared with PGSE was observed in the tibialis anterior for the scans obtained in active dorsiflexion and in the soleus during active plantarflexion. The use of PGSE sequences led to unrealistically elevated axial diffusivity values in the tibialis anterior during dorsiflexion and in the soleus during plantarflexion, while the corresponding values obtained using the OGSE sequences were significantly reduced. Similar findings were seen for radial diffusivity, with significantly higher diffusivity measured in plantarflexion in the soleus muscle using the PGSE sequence. Our preliminary results indicate that DTI with OGSE diffusion encoding is feasible in human musculature and allows to quantitatively assess diffusion properties in actively contracting skeletal muscle. OGSE holds great potential to assess microstructural changes occurring in the skeletal muscle during contraction, and for non-invasive assessment of contraction abnormalities in patients with muscle diseases.

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 (λ₁, λ₂, and λ₃), 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 λ₂ (4.9% IVIM; 5.5% standard), λ₁ (5.3% IVIM; 7.5% standard), and λ₃ (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.

Electronic supplementary material

The online version of this article (10.1007/s00330-019-06403-5) contains supplementary material, which is available to authorized users.

Multi-center evaluation of stability and reproducibility of quantitative MRI measures in healthy calf muscles

The purpose of this study was to evaluate temporal stability, multi-center reproducibility and the influence of covariates on a multimodal MR protocol for quantitative muscle imaging and to facilitate its use as a standardized protocol for evaluation of pathology in skeletal muscle. Quantitative T2, quantitative diffusion and four-point Dixon acquisitions of the calf muscles of both legs were repeated within one hour. Sixty-five healthy volunteers (31 females) were included in one of eight 3-T MR systems. Five traveling subjects were examined in six MR scanners. Average values over all slices of water-T2 relaxation time, proton density fat fraction (PDFF) and diffusion metrics were determined for seven muscles. Temporal stability was tested with repeated measured ANOVA and two-way random intraclass correlation coefficient (ICC). Multi-center reproducibility of traveling volunteers was assessed by a two-way mixed ICC. The factors age, body mass index, gender and muscle were tested for covariance. ICCs of temporal stability were between 0.963 and 0.999 for all parameters. Water-T2 relaxation decreased significantly (P < 10^-3 ) within one hour by ~ 1 ms. Multi-center reproducibility showed ICCs within 0.879-0.917 with the lowest ICC for mean diffusivity. Different muscles showed the highest covariance, explaining 20-40% of variance for observed parameters. Standardized acquisition and processing of quantitative muscle MRI data resulted in high comparability among centers. The imaging protocol exhibited high temporal stability over one hour except for water T2 relaxation times. These results show that data pooling is feasible and enables assembling data from patients with neuromuscular diseases, paving the way towards larger studies of rare muscle disorders.

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.

Muscle Architecture Assessment: Strengths, Shortcomings and New Frontiers of in Vivo Imaging Techniques

Skeletal muscle structural assembly (and its remodeling in response to loading-unloading states) can be investigated macroscopically by assessing muscle architecture, described as fascicle geometric disposition within the muscle. Over recent decades, various medical imaging techniques have been developed to facilitate the in vivo assessment of muscle architecture. However, the main advantages and limitations of these methodologies have been fragmentally discussed. In the present article, the main techniques used for the evaluation of muscle architecture are presented: conventional B-mode ultrasonography, extended-field-of-view ultrasound, 3-D ultrasound and magnetic resonance imaging-based diffusion tensor imaging. By critically discussing potentials and shortcomings of each methodology, we aim to provide readers with an overview of both established and new techniques for the in vivo assessment of muscle architecture. This review may serve as decision guidance facilitating selection of the appropriate technique to be applied in biomedical research or clinical routine.

Multiparametric quantitative MRI assessment of thigh muscles in limb-girdle muscular dystrophy 2A and 2B

INTRODUCTION

The aim of this study was to apply quantitative MRI (qMRI) to assess structural modifications in thigh muscles of subjects with limb girdle muscular dystrophy (LGMD) 2A and 2B with long disease duration.

METHODS

Eleven LGMD2A, 9 LGMD2B patients and 11 healthy controls underwent a multi-parametric 3T MRI examination of the thigh. The protocol included structural T1-weighted images, DIXON sequences for fat fraction calculation, T2 values quantification and diffusion MRI. Region of interest analysis was performed on 4 different compartments (anterior compartment, posterior compartment, gracilis, sartorius).

RESULTS

Patients showed high levels of fat infiltration as measured by DIXON sequences. Sartorius and anterior compartment were more infiltrated in LGMD2B than LGMD2A patients. T2 values were mildly reduced in both disorders. Correlations between clinical scores and qMRI were found.

CONCLUSIONS

qMRI measures may help to quantify muscular degeneration, but careful interpretation is needed when fat infiltration is massive. Muscle Nerve 58: 550-558, 2018.

Skeletal muscle diffusion tensor-MRI fiber tracking: rationale, data acquisition and analysis methods, applications and future directions

The mechanical functions of muscles involve the generation of force and the actuation of movement by shortening or lengthening under load. These functions are influenced, in part, by the internal arrangement of muscle fibers with respect to the muscle's mechanical line of action. This property is known as muscle architecture. In this review, we describe the use of diffusion tensor (DT)-MRI muscle fiber tracking for the study of muscle architecture. In the first section, the importance of skeletal muscle architecture to function is discussed. In addition, traditional and complementary methods for the assessment of muscle architecture (brightness-mode ultrasound imaging and cadaver analysis) are presented. Next, DT-MRI is introduced and the structural basis for the reduced and anisotropic diffusion of water in muscle is discussed. The third section discusses issues related to the acquisition of skeletal muscle DT-MRI data and presents recommendations for optimal strategies. The fourth section discusses methods for the pre-processing of DT-MRI data, the available approaches for the calculation of the diffusion tensor and the seeding and propagating of fiber tracts, and the analysis of the tracking results to measure structural properties pertinent to muscle biomechanics. Lastly, examples are presented of how DT-MRI fiber tracking has been used to provide new insights into how muscles function, and important future research directions are highlighted. Copyright © 2016 John Wiley & Sons, Ltd.

How does passive lengthening change the architecture of the human medial gastrocnemius muscle?

There are few comprehensive investigations of the changes in muscle architecture that accompany muscle contraction or change in muscle length in vivo. For this study, we measured changes in the three-dimensional architecture of the human medial gastrocnemius at the whole muscle level, the fascicle level and the fiber level using anatomical MRI and diffusion tensor imaging (DTI). Data were obtained from eight subjects under relaxed conditions at three muscle lengths. At the whole muscle level, a 5.1% increase in muscle belly length resulted in a reduction in both muscle width (mean change -2.5%) and depth (-4.8%). At the fascicle level, muscle architecture measurements obtained at 3,000 locations per muscle showed that for every millimeter increase in muscle-tendon length above the slack length, average fascicle length increased by 0.46 mm, pennation angle decreased by 0.27° (0.17° in the superficial part and 0.37° in the deep part), and fascicle curvature decreased by 0.18 m^-1 There was no evidence of systematic variation in architecture along the muscle's long axis at any muscle length. At the fiber level, analysis of the diffusion signal showed that passive lengthening of the muscle increased diffusion along fibers and decreased diffusion across fibers. Using these measurements across scales, we show that the complex shape changes that muscle fibers, whole muscles, and aponeuroses of the medial gastrocnemius undergo in vivo cannot be captured by simple geometrical models. This justifies the need for more complex models that link microstructural changes in muscle fibers to macroscopic changes in architecture.NEW & NOTEWORTHY Novel MRI and DTI techniques revealed changes in three-dimensional architecture of the human medial gastrocnemius during passive lengthening. Whole muscle belly width and depth decreased when the muscle lengthened. Fascicle length, pennation, and curvature changed uniformly or near uniformly along the muscle during passive lengthening. Diffusion of water molecules in muscle changes in the same direction as fascicle strains.

Assessment of passive muscle elongation using Diffusion Tensor MRI: Correlation between fiber length and diffusion coefficients

In this study we investigated the changes in fiber length and diffusion parameters as a consequence of passive lengthening and stretching of the calf muscles. We hypothesized that changes in radial diffusivity (RD) are caused by changes in the muscle fiber cross sectional area (CSA) as a consequence of lengthening and shortening of the muscle. Diffusion Tensor MRI (DT-MRI) measurements were made twice in five healthy volunteers, with the foot in three different positions (30° plantarflexion, neutral position and 15° dorsiflexion). The muscles of the calf were manually segmented on co-registered high resolution anatomical scans, and maps of RD and axial diffusivity (AD) were reconstructed from the DT-MRI data. Fiber tractography was performed and mean fiber length was calculated for each muscle group. Significant negative correlations were found between the changes in RD and changes in fiber length in the dorsiflexed and plantarflexed positions, compared with the neutral foot position. Changes in AD did not correlate with changes in fiber length. Assuming a simple cylindrical model with constant volume for the muscle fiber, the changes in the muscle fiber CSA were calculated from the changes in fiber length. In line with our hypothesis, we observed a significant positive correlation of the CSA with the measured changes in RD. In conclusion, we showed that changes in diffusion coefficients induced by passive muscle stretching and lengthening can be explained by changes in muscle CSA, advancing the physiological interpretation of parameters derived from skeletal muscle DT-MRI.

Transplanted Endothelial Progenitor Cells Improve Ischemia Muscle Regeneration in Mice by Diffusion Tensor MR Imaging

Endothelial progenitor cells (EPCs) play an important role in repairing ischemia tissues. Diffusion tensor imaging (DTI) was applied to detect the architectural organization of skeletal muscle. This study investigated the feasibility and accuracy of using the DTI to evaluate effectiveness of EPCs treatment. Mouse bone marrow-derived EPCs were isolated, cultured, characterized, and transplanted to hindlimb ischemia mice model. DTI was performed on the hindlimb at postischemia time points. The edema regions of diffusion restriction (high signal in diffusion weighted imaging) were decreased in the ischemic muscle of EPCs treated mice after 14 days compared with the controls. These results from DTI show the lower apparent diffusion coefficient and eigenvalues (λ1, λ2, and λ3) and the higher fractional anisotropy and fiber counts of ischemic muscle on 7 and 14 days after EPCs treatment compared to the controls. There was a significant correlation between fiber counts calculated by DTI and survival fibers evaluated by histological section (r = 0.873, P < 0.01). Our study demonstrated that the time frame for muscle fiber regeneration after EPCs transplantation was significantly shortened in vivo. DTI could be a useful tool for noninvasive evaluation of muscle tissue damage and repair in animal models and patient with ischemic diseases.

Techniques and applications of skeletal muscle diffusion tensor imaging: A review

Diffusion tensor imaging (DTI) is increasingly applied to study skeletal muscle physiology, anatomy, and pathology. The reason for this growing interest is that DTI offers unique, noninvasive, and potentially diagnostically relevant imaging readouts of skeletal muscle structure that are difficult or impossible to obtain otherwise. DTI has been shown to be feasible within most skeletal muscles. DTI parameters are highly sensitive to patient-specific properties such as age, body mass index (BMI), and gender, but also to more transient factors such as exercise, rest, pressure, temperature, and relative joint position. However, when designing a DTI study one should not only be aware of sensitivity to the above-mentioned factors but also the fact that the DTI parameters are dependent on several acquisition parameters such as echo time, b-value, and diffusion mixing time. The purpose of this review is to provide an overview of DTI studies covering the technical, demographic, and clinical aspects of DTI in skeletal muscles. First we will focus on the critical aspects of the acquisition protocol. Second, we will cover the reported normal variance in skeletal muscle diffusion parameters, and finally we provide an overview of clinical studies and reported parameter changes due to several (patho-)physiological conditions.

Fiber orientation measurements by diffusion tensor imaging improve hydrogen-1 magnetic resonance spectroscopy of intramyocellular lipids in human leg muscles

Twelve healthy subjects underwent hydrogen-1 magnetic resonance spectroscopy ([Formula: see text]) acquisition ([Formula: see text]), diffusion tensor imaging (DTI) with a [Formula: see text]-value of [Formula: see text], and fat-water magnetic resonance imaging (MRI) using the Dixon method. Subject-specific muscle fiber orientation, derived from DTI, was used to estimate the lipid proton spectral chemical shift. Pennation angles were measured as 23.78 deg in vastus lateralis (VL), 17.06 deg in soleus (SO), and 8.49 deg in tibialis anterior (TA) resulting in a chemical shift between extramyocellular lipids (EMCL) and intramyocellular lipids (IMCL) of 0.15, 0.17, and 0.19 ppm, respectively. IMCL concentrations were [Formula: see text], [Formula: see text], and [Formula: see text] in SO, VL, and TA, respectively. Significant differences were observed in IMCL and EMCL pairwise comparisons in SO, VL, and TA ([Formula: see text]). Strong correlations were observed between total fat fractions from [Formula: see text] and Dixon MRI for VL ([Formula: see text]), SO ([Formula: see text]), and TA ([Formula: see text]). Bland-Altman analysis between fat fractions (FFMRS and FFMRI) showed good agreement with small limits of agreement (LoA): [Formula: see text] (LoA: [Formula: see text] to 0.69%) in VL, [Formula: see text] (LoA: [Formula: see text] to 1.33%) in SO, and [Formula: see text] (LoA: [Formula: see text] to 0.47%) in TA. The results of this study demonstrate the variation in muscle fiber orientation and lipid concentrations in these three skeletal muscle types.

Dynamic diffusion-tensor measurements in muscle tissue using the single-line multiple-echo diffusion-tensor acquisition technique at 3T

When diffusion biomarkers display transient changes, i.e. in muscle following exercise, traditional diffusion-tensor imaging (DTI) methods lack the temporal resolution to resolve the dynamics. This article presents an MRI method for dynamic diffusion-tensor acquisitions on a clinical 3T scanner. This method, the Single-Line Multiple-Echo Diffusion-Tensor Acquisition Technique (SL-MEDITATE), achieves a high temporal resolution (4 s) by rapid diffusion encoding through the acquisition of multiple echoes with unique diffusion sensitization and limiting the readout to a single line volume. The method is demonstrated in a rotating anisotropic phantom, a flow phantom with adjustable flow speed and in vivo skeletal calf muscle of healthy volunteers following a plantar flexion exercise. The rotating and flow-varying phantom experiments show that SL-MEDITATE correctly identifies the rotation of the first diffusion eigenvector and the changes in diffusion-tensor parameter magnitudes, respectively. Immediately following exercise, the in vivo mean diffusivity (MD) time courses show, before the well-known increase, an initial decrease that is not typically observed in traditional DTI. In conclusion, SL-MEDITATE can be used to capture transient changes in tissue anisotropy in a single line. Future progress might allow for dynamic DTI when combined with appropriate k-space trajectories and compressed sensing reconstruction.

Diffusion tensor imaging and T2 mapping in early denervated skeletal muscle in rats

BACKGROUND

To evaluate the temporal changes of diffusion tensor imaging (DTI) indices, T2 values, and visual signal intensity on various fat suppression techniques in the early state of denervated skeletal muscle in a rat model.

METHODS

Institutional Animal Care and Use Committee approval was obtained. Sciatic nerves of eight rats were transected for irreversible neurotmesis model. We examined normal lower leg and denervated muscles at 3 days, 1 week, and 2 weeks on a 3 Tesla MR. fractional anisotropy (FA), mean apparent diffusion coefficient (mADC), and T2 values were measured by using DTI and T2 mapping scan. We subjectively classified the signal intensity change on various fat suppression images into the following three grades: negative, suspicious, and definite change. Wilcoxon-sign rank test and Kruskal-Wallis test were used for the comparison of FA, mADC, T2 values. McNemar's test was used for comparing signal intensity change among fat suppression techniques.

RESULTS

FA values of denervated muscles at 3 days (0.35 ± 0.06), 1 week (0.29 ± 0.04), and 2 weeks (0.34 ± 0.05) were significantly (P < 0.05) lower than that in the control group (0.54 ± 0.17). mADC of denervated muscles decreased without statistically significant (P > 0.05) change. T2 values were significantly increased at 1 week (38.11 ± 6.42 ms, P = 0.017) and markedly increased at 2 weeks (46.53 ± 5.17 ms, P = 0.012). The grade of visual signal intensity change on chemical shift selective fat saturation, STIR and IDEAL images were identical in all cases (P = 1.000).

CONCLUSION

FA and T2 values can demonstrate the early temporal changes in denervated rat skeletal muscle.

Muscle changes detected with diffusion-tensor imaging after long-distance running

PURPOSE

To develop a protocol for diffusion-tensor imaging (DTI) of the complete upper legs and to demonstrate feasibility of detection of subclinical sports-related muscle changes in athletes after strenuous exercise, which remain undetected by using conventional T2-weighted magnetic resonance (MR) imaging with fat suppression.

MATERIALS AND METHODS

The research was approved by the institutional ethics committee review board, and the volunteers provided written consent before the study. Five male amateur long-distance runners underwent an MR examination (DTI, T1-weighted MR imaging, and T2-weighted MR imaging with fat suppression) of both upper legs 1 week before, 2 days after, and 3 weeks after they participated in a marathon. The tensor eigenvalues (λ1, λ2, and λ3), the mean diffusivity, and the fractional anisotropy (FA) were derived from the DTI data. Data per muscle from the three time-points were compared by using a two-way mixed-design analysis of variance with a Bonferroni posthoc test.

RESULTS

The DTI protocol allowed imaging of both complete upper legs with adequate signal-to-noise ratio and within a 20-minute imaging time. After the marathon, T2-weighted MR imaging revealed grade 1 muscle strains in nine of the 180 investigated muscles. The three eigenvalues, mean diffusivity, and FA were significantly increased (P < .05) in the biceps femoris muscle 2 days after running. Mean diffusivity and eigenvalues λ1 and λ2 were significantly (P < .05) increased in the semitendinosus and gracilis muscles 2 days after the marathon.

CONCLUSION

A feasible method for DTI measurements of the upper legs was developed that fully included frequently injured muscles, such as hamstrings, in one single imaging session. This study also revealed changes in DTI parameters that over time were not revealed by qualitative T2-weighted MR imaging with fat suppression.

Multi-parametric MRI characterization of healthy human thigh muscles at 3.0 T - relaxation, magnetization transfer, fat/water, and diffusion tensor imaging

Muscle diseases commonly have clinical presentations of inflammation, fat infiltration, fibrosis, and atrophy. However, the results of existing laboratory tests and clinical presentations are not well correlated. Advanced quantitative MRI techniques may allow the assessment of myo-pathological changes in a sensitive and objective manner. To progress towards this goal, an array of quantitative MRI protocols was implemented for human thigh muscles; their reproducibility was assessed; and the statistical relationships among parameters were determined. These quantitative methods included fat/water imaging, multiple spin-echo T2 imaging (with and without fat signal suppression, FS), selective inversion recovery for T1 and quantitative magnetization transfer (qMT) imaging (with and without FS), and diffusion tensor imaging. Data were acquired at 3.0 T from nine healthy subjects. To assess the repeatability of each method, the subjects were re-imaged an average of 35 days later. Pre-testing lifestyle restrictions were applied to standardize physiological conditions across scans. Strong between-day intra-class correlations were observed in all quantitative indices except for the macromolecular-to-free water pool size ratio (PSR) with FS, a metric derived from qMT data. Two-way analysis of variance revealed no significant between-day differences in the mean values for any parameter estimate. The repeatability was further assessed with Bland-Altman plots, and low repeatability coefficients were obtained for all parameters. Among-muscle differences in the quantitative MRI indices and inter-class correlations among the parameters were identified. There were inverse relationships between fractional anisotropy (FA) and the second eigenvalue, the third eigenvalue, and the standard deviation of the first eigenvector. The FA was positively related to the PSR, while the other diffusion indices were inversely related to the PSR. These findings support the use of these T1 , T2 , fat/water, and DTI protocols for characterizing skeletal muscle using MRI. Moreover, the data support the existence of a common biophysical mechanism, water content, as a source of variation in these parameters.

Myoarchitectural and functional alterations in rabbit external anal sphincter muscle following experimental surgical trauma

Obstetrical trauma to external anal sphincter (EAS) is extremely common; however, its role in the development of anal incontinence is not clear. We examined the regenerative process and functional impact of experimental surgical trauma to EAS muscle in an animal model. Surgical myotomy, a craniocaudal incision extending along the entire length and thickness of the EAS, was performed in rabbits. Animals were allowed to recover, and anal pressures were recorded at weekly intervals for 12 wk using a custom-designed probe system to determine the length-tension property of EAS muscle. Animals were killed at predetermined time intervals, and the anal canal was harvested for histochemical studies (for determination of muscle/connective tissue/collagen) and sarcomere length measurement. In addition, magnetic resonance diffusion tensor imaging (MR-DTI) and fiber tracking was performed to determine myoarchitectural changes in the EAS. Myotomy of the EAS muscle resulted in significant impairment of its length-tension property that showed only partial recovery during the 12-wk study period. Histology revealed marked increase in the fibrosis (connective tissue = 69% following myotomy vs. 28% in controls) at 3 wk, which persisted at 12 wk. Immunostaining studies confirmed deposition of collagen in the fibrotic tissue. There was no change in the sarcomere length following myotomy. MR-DTI studies revealed disorganized muscle fiber orientation in the regenerating muscle. We conclude that, following experimental injury, the EAS muscle heals with an increase in the collagen content and loss of normal myoarchitecture, which we suspect is the cause of impaired EAS function.

Time-dependent diffusion in skeletal muscle with the random permeable barrier model (RPBM): application to normal controls and chronic exertional compartment syndrome patients

The purpose of this work was to carry out diffusion tensor imaging (DTI) at multiple diffusion times Td in skeletal muscle in normal subjects and chronic exertional compartment syndrome (CECS) patients and analyze the data with the random permeable barrier model (RPBM) for biophysical specificity. Using an institutional review board approved HIPAA-compliant protocol, seven patients with clinical suspicion of CECS and eight healthy volunteers underwent DTI of the calf muscle in a Siemens MAGNETOM Verio 3 T scanner at rest and after treadmill exertion at four different T(d) values. Radial diffusion values λ(rad) were computed for each of seven different muscle compartments and analyzed with RPBM to produce estimates of free diffusivity D(0), fiber diameter a, and permeability κ. Fiber diameter estimates were compared with measurements from literature autopsy reference for several compartments. Response factors (post/pre-exercise ratios) were computed and compared between normal controls and CECS patients using a mixed-model two-way analysis of variance. All subjects and muscle compartments showed nearly time-independent diffusion along and strongly time-dependent diffusion transverse to the muscle fibers. RPBM estimates of fiber diameter correlated well with corresponding autopsy reference. D(0) showed significant (p < 0.05) increases with exercise for volunteers, and a increased significantly (p < 0.05) in volunteers. At the group level, response factors of all three parameters showed trends differentiating controls from CECS patients, with patients showing smaller diameter changes (p = 0.07), and larger permeability increases (p = 0.07) than controls. Time-dependent diffusion measurements combined with appropriate tissue modeling can provide enhanced microstructural specificity for in vivo tissue characterization. In CECS patients, our results suggest that high-pressure interfiber edema elevates free diffusion and restricts exercise-induced fiber dilation. Such specificity may be useful in differentiating CECS from other disorders or in predicting its response to either physical therapy or fasciotomy.

Fast low-angle shot diffusion tensor imaging with stimulated echo encoding in the muscle of rabbit shank

In the past, spin-echo (SE) echo planar imaging(EPI)-based diffusion tensor imaging (DTI) has been widely used to study the fiber structure of skeletal muscles in vivo. However, this sequence has several shortcomings when measuring restricted diffusion in small animals, such as its sensitivity to susceptibility-related distortions and a relatively short applicable diffusion time. To address these limitations, in the current work, a stimulated echo acquisition mode (STEAM) MRI technique, in combination with fast low-angle shot (FLASH) readout (turbo-STEAM MRI), was implemented and adjusted for DTI in skeletal muscles. Signal preparation using stimulated echoes enables longer effective diffusion times, and thus the detection of restricted diffusion within muscular tissue with intracellular distances up to 100 µm. Furthermore, it has a reduced penalty for fast T2 muscle signal decay, but at the expense of 50% signal loss compared with a SE preparation. Turbo-STEAM MRI facilitates high-resolution DTI of skeletal muscle without introducing susceptibility-related distortions. To demonstrate its applicability, we carried out rabbit in vivo measurements on a human whole-body 3 T scanner. DTI parameters of the shank muscles were extracted, including the apparent diffusion coefficient, fractional anisotropy, eigenvalues and eigenvectors. Eigenvectors were used to calculate maps of structural parameters, such as the planar index and the polar coordinates θ and ϕ of the largest eigenvector. These parameters were compared between three muscles. θ and ϕ showed clear differences between the three muscles, reflecting different pennation angles of the underlying fiber structures. Fiber tractography was performed to visualize and analyze the architecture of skeletal pennate muscles. Optimization of tracking parameters and utilization of T2 -weighted images for improved muscle boundary detection enabled the determination of additional parameters, such as the mean fiber length. The presented results support the applicability of turbo-STEAM MRI as a promising method for quantitative DTI analysis and fiber tractography in skeletal muscles.

Heterogeneity of muscle fat infiltration in children with spina bifida

Children with spina bifida have well recognized functional deficits of muscle, but little is known about the associated changes in muscle anatomy and composition. This study used water-fat magnetic resonance imaging (MRI) to measure fat infiltration in the lower extremity muscles of 11 children with myelomeningocele, the most severe form of spina bifida. MRI measurements of muscle fat fraction (FF) were compared against manual muscle test (MMT) scores for muscle strength. The FF measurements were objective and reliable with mean inter-rater differences of <2% and intraclass correlation coefficients>0.98. There was a significant inverse relationship between muscle FF and MMT scores (P ≤ 0.001). Surprisingly, however, muscles with negligible strength (MMT 0-1) exhibited a bimodal distribution of FF with one group having FF>70% and another group having FF<20%. The MRI also revealed striking heterogeneity amongst individual muscles in the same muscle group (e.g., 4% fat in one participant's lateral gastrocnemius vs. 88% in her medial gastrocnemius), as well as significant asymmetry in FF in one participant with asymmetric strength and sensation. These results suggest that quantitative water-fat MRI may serve as a biomarker for muscle degeneration which may reveal subclinical changes useful for predicting functional potential and prognosis.

Multiple-echo diffusion tensor acquisition technique (MEDITATE) on a 3T clinical scanner

This article describes the concepts and implementation of an MRI method, the multiple-echo diffusion tensor acquisition technique (MEDITATE), which is capable of acquiring apparent diffusion tensor maps in two scans on a 3T clinical scanner. In each MEDITATE scan, a set of RF pulses generates multiple echoes, the amplitudes of which are diffusion weighted in both magnitude and direction by a pattern of diffusion gradients. As a result, two scans acquired with different diffusion weighting strengths suffice for accurate estimation of diffusion tensor imaging (DTI) parameters. The MEDITATE variation presented here expands previous MEDITATE approaches to adapt to the clinical scanner platform, such as exploiting longitudinal magnetization storage to reduce T2 weighting. Fully segmented multi-shot Cartesian encoding is used for image encoding. MEDITATE was tested on isotropic (agar gel), anisotropic diffusion phantoms (asparagus) and in vivo skeletal muscle in healthy volunteers with cardiac gating. Comparisons of accuracy were performed with standard twice-refocused spin echo (TRSE) DTI in each case and good quantitative agreement was found between diffusion eigenvalues, mean diffusivity and fractional anisotropy derived from TRSE DTI and from the MEDITATE sequence. Orientation patterns were correctly reproduced in both isotropic and anisotropic phantoms, and approximately for in vivo imaging. This illustrates that the MEDITATE method of compressed diffusion encoding is feasible on the clinical scanner platform. With future development and employment of appropriate view-sharing image encoding, this technique may be used in clinical applications requiring time-sensitive acquisition of DTI parameters such as dynamical DTI in muscle.

Stimulated echo diffusion tensor imaging and SPAIR T2 -weighted imaging in chronic exertional compartment syndrome of the lower leg muscles

PURPOSE

To evaluate the performance of diffusion tensor imaging (DTI) in the evaluation of chronic exertional compartment syndrome (CECS) as compared to T2 -weighted (T2w) imaging.

MATERIALS AND METHODS

Using an Institutional Review Board (IRB)-approved, Health Insurance Portability and Accountability Act (HIPAA)-compliant protocol, spectral adiabatic inversion recovery (SPAIR) T2w imaging and stimulated echo DTI were applied to eight healthy volunteers and 14 suspected CECS patients before and after exertion. Longitudinal and transverse diffusion eigenvalues, mean diffusivity (MD), and fractional anisotropy (FA) were measured in seven calf muscle compartments, which in patients were classified by their response on T2w: normal (<20% change), and CECS (>20% change). Mixed model analysis of variance compared subject groups and compartments in terms of response factors (post/pre-exercise ratios) of DTI parameters.

RESULTS

All diffusivities significantly increased (P < 0.0001) and FA decreased (P = 0.0014) with exercise. Longitudinal diffusion responses were significantly smaller than transversal diffusion responses (P < 0.0001). Nineteen of 98 patient compartments were classified as CECS on T2w. MD increased by 3.8 ± 3.4% (volunteer), 7.4 ± 4.2% (normal), and 9.1 ± 7.0% (CECS) with exercise.

CONCLUSION

DTI shows promise as an ancillary imaging method in the diagnosis and understanding of the pathophysiology in CECS. Future studies may explore its utility in predicting response to treatment.

Polynomial fitting of DT-MRI fiber tracts allows accurate estimation of muscle architectural parameters

Fiber curvature is a functionally significant muscle structural property, but its estimation from diffusion-tensor magnetic resonance imaging fiber tracking data may be confounded by noise. The purpose of this study was to investigate the use of polynomial fitting of fiber tracts for improving the accuracy and precision of fiber curvature (κ) measurements. Simulated image data sets were created in order to provide data with known values for κ and pennation angle (θ). Simulations were designed to test the effects of increasing inherent fiber curvature (3.8, 7.9, 11.8 and 15.3 m(-1)), signal-to-noise ratio (50, 75, 100 and 150) and voxel geometry (13.8- and 27.0-mm(3) voxel volume with isotropic resolution; 13.5-mm(3) volume with an aspect ratio of 4.0) on κ and θ measurements. In the originally reconstructed tracts, θ was estimated accurately under most curvature and all imaging conditions studied; however, the estimates of κ were imprecise and inaccurate. Fitting the tracts to second-order polynomial functions provided accurate and precise estimates of κ for all conditions except very high curvature (κ=15.3 m(-1)), while preserving the accuracy of the θ estimates. Similarly, polynomial fitting of in vivo fiber tracking data reduced the κ values of fitted tracts from those of unfitted tracts and did not change the θ values. Polynomial fitting of fiber tracts allows accurate estimation of physiologically reasonable values of κ, while preserving the accuracy of θ estimation.

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.

Diffusion-Tensor MRI Based Skeletal Muscle Fiber Tracking

A skeletal muscle's function is strongly influenced by the internal organization and geometric properties of its fibers, a property known as muscle architecture. Diffusion-tensor magnetic resonance imaging-based fiber tracking provides a powerful tool for non-invasive muscle architecture studies, has three-dimensional sensitivity, and uses a fixed frame of reference. Significant advances have been made in muscle fiber tracking technology, including defining seed points for fiber tracking, quantitatively characterizing muscle architecture, implementing denoising procedures, and testing validity and repeatability. Some examples exist of how these data can be integrated with those from other advanced MRI and computational methods to provide novel insights into muscle function. Perspectives are offered regarding future directions in muscle diffusion-tensor imaging, including needs to develop an improved understanding for the microstructural basis for reduced and anisotropic diffusion, establish the best practices for data acquisition and analysis, and integrate fiber tracking with other physiological data.

Simultaneous temperature and regional blood volume measurements in human muscle using an MRI fast diffusion technique

The thermal dependence of the translational diffusion coefficient and of the regional blood volume was investigated in vivo by using a special MR pulsed gradient technique with reduced sensitivity to bulk tissue motion. Measurements were done at 0.5 T, using a small gradient insert. The diffusion coefficient of muscle water was calibrated against thermocouple-measured temperature in vitro, both with the muscle fibers parallel and perpendicular to the diffusion gradient. The maximum muscle temperature variation obtained by percutaneous conduction was -8.8 +/- 1.6 degrees C under cooling and +3.7 +/- 1.6 degrees C under heating, from basal state. Simultaneously the fractional regional blood volume decreased by a factor of 3.5 under cooling and increased by a factor of 2.7 under heating. Due to the interdependence of microcirculation and tissue temperature, this technique may be used to follow heat production or deposition in living tissue (muscle exercise, electromagnetic irradiation, etc.).