DTI of human skeletal muscle: the effects of diffusion encoding parameters, signal-to-noise ratio and T2 on tensor indices and fiber tracts

From Voxels to Physiology: A Review of Diffusion Magnetic Resonance Imaging Applications in Skeletal Muscle

Skeletal muscle has a classic structure function relationship; both skeletal muscle microstructure and architecture are directly related to force generating capacity. Biopsy, the gold standard for evaluating muscle microstructure, is highly invasive, destructive to muscle, and provides only a small amount of information about the entire volume of a muscle. Similarly, muscle fiber lengths and pennation angles, key features of muscle architecture predictive of muscle function, are traditionally studied via cadaveric dissection. Noninvasive techniques such as diffusion magnetic resonance imaging (dMRI) offer quantitative approaches to study skeletal muscle microstructure and architecture. Despite its prevalence in applications for musculoskeletal research, clinical adoption is hindered by a lack of understanding regarding its sensitivity to clinically important biomarkers such as muscle fiber cross-sectional area. This review aims to elucidate how dMRI has been utilized to study skeletal muscle, covering fundamentals of muscle physiology, dMRI acquisition techniques, dMRI modeling, and applications where dMRI has been leveraged to noninvasively study skeletal muscle changes in response to disease, aging, injury, and human performance. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 2.

Stimulated echo acquisition mode (STEAM) diffusion tensor imaging with different diffusion encoding times in the supraspinatus muscle: Test-retest reliability and comparison to spin echo diffusion tensor imaging

Diffusion tensor imaging (DTI) provides insight into the skeletal muscle microstructure and can be acquired using a stimulated echo acquisition mode (STEAM)-based approach to quantify time-dependent tissue diffusion. This study examined diffusion metrics and signal-to-noise ratio (SNR) in the supraspinatus muscle obtained with a STEAM-DTI sequence with different diffusion encoding times (Δ) and compared them to measures from a spin echo (SE) sequence. Ten healthy subjects (mean age 31.5 ± 4.7 years; five females) underwent 3-Tesla STEAM and SE-DTI of the shoulder in three sessions. STEAM was acquired with Δ of 100/200/400/600 ms. The diffusion encoding time in SE scans was 19 ms (b = 500 s/mm2). Region of interest-based measurement of fractional anisotropy (FA), mean diffusivity (MD), and SNR was performed. Intraclass correlation coefficients (ICCs) were computed to assess test-retest reliability. ANOVA with post-hoc pairwise tests was used to compare measures between different Δ of STEAM as well as STEAM and SE, respectively. FA was significantly higher (FASTEAM: 0.38-0.46 vs. FASE: 0.26) and MD significantly lower (MDSTEAM: 1.20-1.33 vs. MDSE: 1.62 × 10-3 mm2/s) in STEAM compared to SE (p < 0.001, respectively). SNR was significantly higher for SE (72.3 ± 8.7) than for STEAM (p < 0.001). ICCs were excellent for FA in STEAM (≥0.911) and SE (0.960). For MD, ICCs were good for STEAM100ms-600ms (≥0.759) and SE (0.752). STEAM and SE exhibited excellent reliability for FA and good reliability for MD in the supraspinatus muscle. SNR was significantly higher in SE compared to STEAM.

End-divergent architecture diversifies within-muscle mechanical action in human gluteus maximus in vivo

A muscle's mechanical action is affected by its architecture. However, less is known about the architecture of muscles with broad attachments: "end-divergent" muscles. Potential regional variation of fascicle orientation in end-divergent muscles suggests that their mechanical action varies by region. Here, we comprehensively examined 3D architecture and potential action of the human gluteus maximus (typical end-divergent muscle) in vivo. The gluteus maximus fascicles were three-dimensionally reconstructed over the whole muscle belly using diffusion tensor imaging and tractography. We calculated the force fraction and moment-arm length about the hip joint for individual muscle fascicles, and their product (specific torque, an estimate of torque-generating capacity for a given cross-sectional area). We found that the specific torque for hip extension and external rotation tended to be greater in the distal than the other regions, whereas that for hip abduction appeared to be greater in the proximal than the other regions. Notably, the distal-lateral region exhibited a negative specific torque for hip abduction, indicating that fascicles in this region act for hip "adduction". These findings indicate that end-divergent architecture diversifies within-muscle mechanical action in terms of directions as well as magnitudes in vivo.

Optimal Setup and Parameters of Diffusion-Weighted Magnetic Resonance Imaging for Translational Evaluation of a Tumor Progression Model for Soft Tissue Sarcomas

Purpose

Defining a microscopic tumor infiltration boundary is critical to the success of radiation therapy. Currently, radiation oncologists use margins to geometrically expand the visible tumor for radiation treatment planning in soft tissue sarcomas (STS). Image-based models of tumor progression would be critical to personalize the treatment radiation field to the pattern of sarcoma spread. Evaluation of these models is necessary to demonstrate feasibility in the clinical setting. This study presents an imaging protocol for the preclinical evaluation of a tumor progression model in extremity STS.

Methods and Materials

We recruited 7 healthy volunteers and acquired diffusion-weighted magnetic resonance imaging (DW-MRI) images of the thigh on a magnetic resonance imaging scanner used for imaging cancer patients in a radiation oncology department. We developed a protocol that includes positioning the patient, configuring the radiofrequency coils, and setting the DW-MRI sequence parameters. To find the optimal parameter configuration, the image signal-to-noise ratio (SNR) and the directional variability (DV) of the principal eigenvector of the diffusion tensor were calculated.

Results

The mean SNR across all trials and 12 thigh muscles was 41, with a range of 12 to 72. The mean DV was 13° and ranged from 11° to 23°. The longest scan time was 22 minutes and 58 seconds, and the shortest was 11 minutes and 46 seconds. For the high-resolution image with a voxel volume of 1.3 × 1.3 × 6 mm3 and 38 slices, the optimal parameters were found to be a repetition time of 8000 ms, 12 signal averages, and 6 gradient directions. This configuration resulted in a scan time of 11 minutes and 46 seconds, an SNR of 34, and a DV of 13°.

Conclusions

A DW-MRI scan duration acceptable for imaging cancer patients was achieved with an image quality suitable for reproducible modeling of tumor infiltration. The developed protocol can be used for preclinical evaluation in STS patients.

Investigating skeletal muscle micro-trauma with time-dependent diffusion and the random permeable barrier model

Repeated muscle micro-trauma may cause severe muscle damage. Diffusion tensor imaging (DTI) exhibits sensitivity to microstructural changes in skeletal muscle. We hypothesize that longer diffusion times enhance sensitivity to micro-trauma and that membrane permeability increases with micro-trauma. We obtained DTI scans of the thighs in ten male runners 1 week before (TP-0), 24-48 h after (TP-1), and 2 weeks after (TP-2) they completed a marathon. Diffusion times were 28.1, 116.7, and 316.7 ms. The random permeable barrier model (RPBM) was fitted to the radial diffusivities to obtain estimates for fiber diameter and membrane permeability. Hamstring and quadriceps muscles were manually segmented. A linear mixed model assessed variations across time points and diffusion times within the DTI and RPBM parameters and assessed sensitivity to micro-trauma by comparing %-changes in DTI parameters at TP-1 and TP-2 to TP-0. All DTI parameters except FA significantly changed between TP-0 and TP-1, and between TP-1 and TP-2. The %-change did not differ between diffusion times. The permeability increased at TP-1 and TP-2 compared to TP-0. In conclusion, longer diffusion times did not improve sensitivity to micro-trauma. The increased permeability post-marathon underscores the potential of RPBM-derived parameters as a biomarker for micro-trauma and muscle injuries.

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.

A Comparison of Skeletal Muscle Diffusion Tensor Imaging Tractography Seeding Methods

The internal arrangement of a muscle's fibers with respect to its mechanical line of action (muscle architecture) is a major determinant of muscle function. Muscle architecture can be quantified using diffusion tensor magnetic resonance imaging-based tractography, which propagates streamlines from a set of seed points by integrating vectors that represent the direction of greatest water diffusion (and by inference, the local fiber orientation). Previous work has demonstrated that tractography outcomes are sensitive to the method for defining seed points, but this sensitivity has not been fully examined. To do so, we developed a realistic simulated muscle architecture and implemented four novel methods for tract seeding: seeding along the muscle-aponeurosis boundary with an updated procedure for rounding seed points prior to lookup in the muscle boundary mask and diffusion tensor matrix (APO-3); voxel-based seeding throughout the muscle volume at a user-specified spatial frequency (VXL-1); voxel-based seeding throughout the muscle volume at a variable spatial frequency (VXL-2), and seeding near external and internal muscle boundaries (VXL-3). We then implemented these methods in an example human dataset. The updated aponeurosis seeding procedures allow more accurate and robust tract propagation from seed points. The voxel-based seeding methods had quantification outcomes that closely matched the updated aponeurosis seeding method. Further, the voxel-based methods can accelerate the overall workflow and may be beneficial in high throughput analysis of multi-muscle datasets. Continued evaluation of these methods in a wider range of muscle architectures is warranted.

Diffusion tensor imaging of hamstring muscles after acute strain injury and throughout recovery in collegiate athletes

OBJECTIVE

To identify the region of interest (ROI) to represent injury and observe between-limb diffusion tensor imaging (DTI) microstructural differences in muscle following hamstring strain injury.

MATERIALS AND METHODS

Participants who sustained a hamstring strain injury prospectively underwent 3T-MRI of bilateral thighs using T1, T2, and diffusion-weighted imaging at time of injury (TOI), return to sport (RTS), and 12 weeks after RTS (12wks). ROIs were using the hyperintense region on a T2-weighted sequence: edema, focused edema, and primary muscle injured excluding edema (no edema). Linear mixed-effects models were used to compare diffusion parameters between ROIs and timepoints and limbs and timepoints.

RESULTS

Twenty-four participants (29 injuries) were included. A significant ROI-by-timepoint interaction was detected for all diffusivity measures. The edema and focused edema ROIs demonstrated increased diffusion at TOI compared to RTS for all diffusivity measures (p-values < 0.006), except λ1 (p-values = 0.058-0.12), and compared to 12wks (p-values < 0.02). In the no edema ROI, differences in diffusivity measures were not observed (p-values > 0.82). At TOI, no edema ROI diffusivity measures were lower than the edema ROI (p-values < 0.001) but not at RTS or 12wks (p-values > 0.69). A significant limb-by-timepoint interaction was detected for all diffusivity measures with increased diffusion in the involved limb at TOI (p-values < 0.001) but not at RTS or 12wks (p-values > 0.42). Significant differences in fractional anisotropy over time or between limbs were not detected.

CONCLUSION

Hyperintensity on T2-weighted imaging used to define the injured region holds promise in describing muscle microstructure following hamstring strain injury by demonstrating between-limb differences at TOI but not at follow-up timepoints.

Diffusion tensor imaging reveals subclinical alterations in muscles of patients with Becker muscular dystrophy

OBJECTIVES

Becker muscular dystrophy (BMD) is a relatively less investigated neuromuscular disease, partially overlapping the phenotype of Duchenne dystrophy (DMD). Physiopathological and anatomical patterns are still not comprehensively known, despite recent effort in the search of early biomarkers. Aim of this study was to selectively compare normal appearing muscles of BMD with healthy controls.

METHODS

Among a pool of 40 BMD patients and 20 healthy controls, Sartorius and gracilis muscles were selected on the basis of a blinded clinical quantitative/qualitative evaluation, if classified as normal (0 or 1 on Mercuri scale) and subsequently segmented on diffusion tensor MRI scans with a tractographic approach. Diffusion derived parameters were extracted.

RESULTS

Non-parametric testing revealed significant differences between normal and normal appearing BMD derived parameters in both muscles, the difference being more evident in sartorius. Bonferroni-corrected P-values (<.05) of Mann-Whitney test could discriminate between BMD and controls for standard deviation of all diffusion parameters (mean diffusivity, fractional anisotropy, axial and radial diffusivity) in both sartorius and gracilis, while in sartorius the significant difference was found also in the average values of the same parameters (with exception of RD).

CONCLUSIONS

This method could identify microstructural alterations in BMD normal appearing sartorius and gracilis.

ADVANCES IN KNOWLEDGE

Diffusion based MRI could be able to identify possible early or subclinical microstructural alterations in dystrophic patients with BMD.

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.

Impact of Deep Learning Denoising Algorithm on Diffusion Tensor Imaging of the Growth Plate on Different Spatial Resolutions

To assess the impact of a deep learning (DL) denoising reconstruction algorithm applied to identical patient scans acquired with two different voxel dimensions, representing distinct spatial resolutions, this IRB-approved prospective study was conducted at a tertiary pediatric center in compliance with the Health Insurance Portability and Accountability Act. A General Electric Signa Premier unit (GE Medical Systems, Milwaukee, WI) was employed to acquire two DTI (diffusion tensor imaging) sequences of the left knee on each child at 3T: an in-plane 2.0 × 2.0 mm2 with section thickness of 3.0 mm and a 2 mm3 isovolumetric voxel; neither had an intersection gap. For image acquisition, a multi-band DTI with a fat-suppressed single-shot spin-echo echo-planar sequence (20 non-collinear directions; b-values of 0 and 600 s/mm2) was utilized. The MR vendor-provided a commercially available DL model which was applied with 75% noise reduction settings to the same subject DTI sequences at different spatial resolutions. We compared DTI tract metrics from both DL-reconstructed scans and non-denoised scans for the femur and tibia at each spatial resolution. Differences were evaluated using Wilcoxon-signed ranked test and Bland-Altman plots. When comparing DL versus non-denoised diffusion metrics in femur and tibia using the 2 mm × 2 mm × 3 mm voxel dimension, there were no significant differences between tract count (p = 0.1, p = 0.14) tract volume (p = 0.1, p = 0.29) or tibial tract length (p = 0.16); femur tract length exhibited a significant difference (p < 0.01). All diffusion metrics (tract count, volume, length, and fractional anisotropy (FA)) derived from the DL-reconstructed scans, were significantly different from the non-denoised scan DTI metrics in both the femur and tibial physes using the 2 mm3 voxel size (p < 0.001). DL reconstruction resulted in a significant decrease in femorotibial FA for both voxel dimensions (p < 0.01). Leveraging denoising algorithms could address the drawbacks of lower signal-to-noise ratios (SNRs) associated with smaller voxel volumes and capitalize on their better spatial resolutions, allowing for more accurate quantification of diffusion metrics.

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.

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.

Sensitivity analyses of probabilistic and deterministic DTI tractography methodologies for studying arm muscle architecture

PURPOSE

To determine the sensitivity profiles of probabilistic and deterministic DTI tractography methods in estimating geometric properties in arm muscle anatomy.

METHODS

Spin-echo diffusion-weighted MR images were acquired in the dominant arm of 10 participants. Both deterministic and probabilistic tractography were performed in two different muscle architectures of the parallel-structured biceps brachii (and the pennate-structured flexor carpi ulnaris. Muscle fascicle geometry estimates and number of fascicles were evaluated with respect to tractography turning angle, polynomial fitting order, and SNR. The DTI tractography estimated fascicle lengths were compared with measurements obtained from conventional cadaveric dissection and ultrasound modalities.

RESULTS

The probabilistic method generally estimated fascicle lengths closer to ranges reported by conventional methods than the deterministic method, most evident in the biceps brachii (p > 0.05), consisting of longer, arc-like fascicles. For both methods, a wide turning angle (50º-90°) generated fascicle lengths that were in close agreement with conventional methods, most evident in the flexor carpi ulnaris (p > 0.05), consisting of shorter, feather-like fascicles. The probabilistic approach produced at least two times more fascicles than the deterministic approach. For both approaches, second-order fitting yielded about double the complete tracts as third-order fitting. In both muscles, as SNR decreased, deterministic tractography produced less fascicles but consistent geometry (p > 0.05), whereas probabilistic tractography produced a consistent number but altered geometry of fascicles (p < 0.001).

CONCLUSION

Findings from this study provide best practice recommendations for implementing DTI tractography in skeletal muscle and will inform future in vivo studies of healthy and pathological muscle structure.

Double pulsed field gradient diffusion MRI to assess skeletal muscle microstructure

PURPOSE

Preliminary study to determine whether double pulsed field gradient (PFG) diffusion MRI is sensitive to key features of muscle microstructure related to function.

METHODS

The restricted diffusion profile of molecules in models of muscle microstructure derived from histology were systematically simulated using a numerical simulation approach. Diffusion tensor subspace imaging analysis of the diffusion signal was performed, and spherical anisotropy (SA) was calculated for each model. Linear regression was used to determine the predictive capacity of SA on the fiber area, fiber diameter, and surface area to volume ratio of the models. Additionally, a rat model of muscle hypertrophy was scanned using a single PFG and a double PFG pulse sequence, and the restricted diffusion measurements were compared with histological measurements of microstructure.

RESULTS

Excellent agreement between SA and muscle fiber area (r2  = 0.71; p < 0.0001), fiber diameter (r2  = 0.83; p < 0.0001), and surface area to volume ratio (r2  = 0.97; p < 0.0001) in simulated models was found. In a scanned rat leg, the distribution of these microstructural features measured from histology was broad and demonstrated that there is a wide variance in the microstructural features observed, similar to the SA distributions. However, the distribution of fractional anisotropy measurements in the same tissue was narrow.

CONCLUSIONS

This study demonstrates that SA-a scalar value from diffusion tensor subspace imaging analysis-is highly sensitive to muscle microstructural features predictive of function. Furthermore, these techniques and analysis tools can be translated to real experiments in skeletal muscle. The increased dynamic range of SA compared with fractional anisotropy in the same tissue suggests increased sensitivity to detecting changes in tissue microstructure.

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.

Inter-vendor and inter-observer reliability of diffusion tensor imaging in the musculoskeletal system: a multiscanner MR study

BACKGROUND

To evaluate the inter-observer and inter-vendor reliability of diffusion tensor imaging parameters in the musculoskeletal system.

METHODS

This prospective study included six healthy volunteers three men (mean age: 42; range: 31-52 years) and three women (mean age: 36; range: 30-44 years). Each subject was scanned using different 3 Tesla magnetic resonance scanners from three different vendors at three different sites bilaterally. First, the intra-class correlation coefficient was used to determine between-observers agreement for overall measurements and clinical sites. Next, between-group comparisons were made through the nonparametric Friedman's test. Finally, the Bland-Altman method was used to determine agreement among the three scanner measurements, comparing them two by two.

RESULTS

A total of 792 measurement were calculated. ICC reported high levels of agreement between the two observers. ICC related to MD, FA, and RD measurements ranged from 0.88 (95% CI 0.85-0.90) to 0.95 (95% CI 0.94-0.96), from 0.85 (95% CI 0.81-0.88) to 0.95 (95% CI 0.93-0.96), and from 0.89 (0.85-0.90) to 0.92 (0.90-0.94). No statistically significant inter-vendor differences were observed. The Bland-Altmann method confirmed a high correlation between parameter values.

CONCLUSION

An excellent inter-observer and inter-vendor reliability was found in our study.

A handbook for beginners in skeletal muscle diffusion tensor imaging: physical basis and technical adjustments

Magnetic resonance imaging (MRI) of skeletal muscle is routinely performed using morphological sequences to acquire anatomical information. Recently, there is an increasing interest in applying advanced MRI techniques that provide pathophysiologic information for skeletal muscle evaluation to complement standard morphologic information. Among these advanced techniques, diffusion tensor imaging (DTI) has emerged as a potential tool to explore muscle microstructure. DTI can noninvasively assess the movement of water molecules in well-organized tissues with anisotropic diffusion, such as skeletal muscle. The acquisition of DTI studies for skeletal muscle assessment requires specific technical adjustments. Besides, knowledge of DTI physical basis and skeletal muscle physiopathology facilitates the evaluation of this advanced sequence and both image and parameter interpretation. Parameters derived from DTI provide a quantitative assessment of muscle microstructure with potential to become imaging biomarkers of normal and pathological skeletal muscle. KEY POINTS: • Diffusion tensor imaging (DTI) allows to evaluate the three-dimensional movement of water molecules inside biological tissues. • The skeletal muscle structure makes it suitable for being evaluated with DTI. • Several technical adjustments have to be considered for obtaining robust and reproducible DTI studies for skeletal muscle assessment, minimizing potential artifacts.

Measuring healthy female nulliparous pubovisceral muscle from diffusion kurtosis imaging

This study explores the feasibility of using diffusion kurtosis imaging (DKI) in the pelvic floor region and assesses the water diffusivity of the pubovisceral muscle. Twenty-seven healthy young nulliparous females underwent DKI at 3.0 T that included 15 gradient directions and three b values (0, 750, and 1500 s/mm² ). The diffusion tensor and diffusion kurtosis metrics values of the pubovisceral muscle were measured after image processing. Two independent sample t-tests, a paired-samples t-test, and a nonparametric hypothesis test were performed as appropriate to compare the differences among different metrics. Twenty-six subjects (mean ± standard deviation age, 25 ± 2 years) were successfully analyzed by measuring the diffusion tensor and diffusion kurtosis metrics of the bilateral pubovisceral muscles. The metrics included mean kurtosis, axial kurtosis, radial kurtosis, fractional anisotropy, mean diffusivity, axial diffusivity, and radial diffusivity. We found no statistically significant differences for these measurement values between the left and right pubovisceral muscles (p = 0.271-0.931). However, radial kurtosis was greater than axial kurtosis in both pubovisceral muscles (p < 0.001) and axial diffusivity was lower than radial diffusivity in both pubovisceral muscles (p < 0.001). We deem the application of DKI technology to the pelvic floor region to be feasible.

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.

Morphological comparison of masseter muscle fibers in the mandibular rest and open positions using diffusion tensor imaging

BACKGROUND

The masseter muscle has a complicated multipennate internal structure and exhibits functional differentiation when performing various stomatognathic functions. It is important to understand the internal structural changes of the muscle during functioning to elucidate characteristic muscle disorders such as local myalgia. Diffusion tensor imaging (DTI) may be useful for investigating the internal structural features of muscle.

OBJECTIVES

To evaluate the features of masseter muscle fibers in human participants using DTI fiber tractography, and to elucidate the structural differences in the masseter muscle between the mandibular rest and open positions.

METHODS

Five healthy men (age 31±7 years) underwent DTI and T1-weighted MRI of the right masseter muscle in the mandibular rest and open positions. MR images were used as a reference for muscle layer segmentation (superficial, intermediate, and deep). DTI fiber tractography of the masseter muscle was performed and the orientation of the DTI fibers was analyzed in each layer using coordinates based on the Frankfurt horizontal plane.

RESULTS

The DTI fiber orientation of the deep layer significantly changed between the mandibular rest and open positions in the frontal plane (p<0.05, Wilcoxon rank sum test). However, no significant change was found in the superficial and intermediate layers.

CONCLUSION

DTI fiber tractography confirmed regional differences in the orientation change of the masseter muscle fibers between different mandibular positions. The results may support the existence of functional partitioning inside the masseter muscle and suggest that DTI may be useful for the evaluation of muscle fibers in multipennate muscles.

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.

Potential Role of MRI Imaging for Myofascial Pain: A Scoping Review for the Clinicians and Theoretical Considerations

The most common cause of chronic musculoskeletal pain is chronic myofascial pain syndrome (MPS). MPS often presents with increased muscle stiffness, and the myofascial trigger point (MTrP). Imaging modalities have been used to identify the MTrP, but their role in the detection and diagnosis of MPS remains unclear. The purpose of this review was to identify evidence in literature for the use of imaging in the role of classifying and explaining the physiology of MTrPs. Since few imaging techniques have been performed on MTrPs, we explored the imaging techniques that can effectively image complex skeletal muscle microstructure, and how they could be used. As part of a scoping review, we conducted a systematic search from three medical databases (CINAHL, EMBASE and MEDLINE) from year to year to analyze past MTrP imaging, as well as analyzing imaging techniques performed on the microstructure of muscle. Previously, ultrasound has been used to differentiate active, latent MTrPs, but these studies do not adequately address their underlying anatomical structure. MRI remains the standard method of imaging skeletal muscle. The existing MRI literature suggests that the DTI technique can quantify muscle injury, strain, and structure. However, theoretically, HARDI and DKI techniques seem to provide more information for complex structural areas, although these modalities have a disadvantage of longer scan times and have not been widely used on skeletal muscle. Our review suggests that DTI is the most effective imaging modality that has been used to define the microstructure of muscle and hence, could be optimal to image the MTrP. HARDI and DKI are techniques with theoretical potential for analysis of muscle, which may provide more detailed information representative of finer muscle structural features. Future research utilizing MRI techniques to image muscle are necessary to provide a more robust means of imaging skeletal muscle and the MTrP.

Varying diffusion time to discriminate between simulated skeletal muscle injury models using stimulated echo diffusion tensor imaging

PURPOSE

Evaluate the relationship between muscle microstructure, diffusion time (Δ), and the diffusion tensor (DT) to identify the optimal Δ where changes in muscle fiber size may be detected.

METHODS

The DT was simulated in models with histology informed geometry over a range of Δ with a stimulated echo DT imaging (DTI) sequence using the numerical simulation application DifSim. The difference in the DT at each Δ between healthy and injured skeletal muscle models was calculated, to identify the optimal Δ at which changes in muscle fiber size may be detected. The random permeable barrier model (RPBM) was used to estimate muscle microstructure from the simulated DT measurements, which were compared to the ground truth.

RESULTS

Across all models, fractional anisotropy provided greater contrast between injured and control models than diffusivity measurements. Compared to control models, in atrophic injury models, the greatest difference in the DT was found between 90 ms and 250 ms. In models with acute edema, the contrast between injured and control muscle increased with increasing diffusion time, although these models had smaller mean fiber areas. RPBM systematically underestimated fiber size but accurately estimated surface area-to-volume ratio of simulated models.

CONCLUSION

These findings may better inform pulse sequence parameter selection when performing DTI experiments in vivo. If only a single diffusion experiment can be performed, the selected Δ should be ~170 ms to maximize the ability to discriminate between different injury models. Ideally several diffusion times between 90 ms and 500 ms should be sampled in order to maximize diffusion contrast, particularly when the disease process is unknown.

Quantification of disease progression in spinal muscular atrophy with muscle MRI-a pilot study

OBJECTIVES

Quantitative MRI (qMRI) of muscles is a promising tool to measure disease progression or to assess therapeutic effects in neuromuscular diseases. Longitudinal imaging studies are needed to show sensitivity of qMRI in detecting disease progression in spinal muscular atrophy (SMA). In this pilot study we therefore studied one-year changes in quantitative MR parameters in relation to clinical scores.

METHODS

We repeated quantitative 3 T MR analysis of thigh muscles and clinical testing one year after baseline in 10 treatment-naïve patients with SMA, 5 with Type 2 (21.6 ± 7.0 years) and 5 with Type 3 (33.4 ± 11.9 years). MR protocol consisted of Dixon, T₂ mapping and diffusion tensor imaging (DTI). The temporal relation of parameters was examined with a mixed model.

RESULTS

We detected a significant increase in fat fraction (baseline, 38.2% SE 0.6; follow-up, 39.5% SE 0.6; +1.3%, p = 0.001) in all muscles. Muscles with moderate to high fat infiltration at baseline show a larger increase over time (+1.6%, p < 0.001). We did not find any changes in DTI parameters except for low fat-infiltration muscles (m. adductor longus and m. biceps femoris (short head)). The T₂ of muscles decreased from 28.2 ms to 28.0 ms (p = 0.07). Muscle strength and motor function scores were not significantly different between follow-up and baseline.

CONCLUSION

Longitudinal imaging data show slow disease progression in skeletal muscle of the thigh of (young-) adult patients with SMA despite stable strength and motor function scores. Quantitative muscle imaging demonstrates potential as a biomarker for disease activity and monitoring of therapy response.

Exercise-induced muscle damage: multi-parametric MRI quantitative assessment

BACKGROUND

To explore the value of magnetic resonance quantitative analysis using diffusion tensor imaging, T2 mapping, and intravoxel incoherent motion in the evaluation of eccentric exercise-induced muscle damage and to compare the effects of various eccentric exercise modes at different time points in rats.

METHODS

A total of 174 Sprague-Dawley male rats were randomly divided into five groups: control, once-only exercise, continuous exercise, intermittent exercise, and once-fatigue exercise groups. Each experimental group was divided into seven time-subgroups: 0.5 h, 24 h, 48 h, 72 h, 96 h, 120 h and 168 h after exercise. The quadriceps femoris muscles were then scanned using magnetic resonance imaging. The apparent diffusion coefficient and fractional anisotropy values of diffusion tensor imaging, T2 values of T2 mapping, D and D* values of intravoxel incoherent motion and optical density values of desmin were measured. Associations among different eccentric exercise programmes, magnetic resonance imaging findings, and histopathological results were evaluated. Dunnett's test, two-way repeated measures analysis of variance, and Pearson correlation analysis were used for statistical analysis.

RESULTS

Diffusion tensor imaging showed that the number of muscle fibre bundles decreased to varying degrees with different time points and eccentric exercises. Apparent diffusion coefficient values of the exercise groups showed a trend that first increased and then decreased, the opposite of fractional anisotropy. The specimens in all eccentric exercise programmes showed high signal T2 values after exercise, the highest among which was in the once-fatigue exercise group. D and D* in the experimental groups were significantly higher than those in the control group at 0.5-48 h after exercise. The apparent diffusion coefficient, fractional anisotropy, T2, D and D* values correlated with the optical density values of desmin.

CONCLUSIONS

Diffusion tensor imaging, T2 mapping, and intravoxel incoherent motion technology accurately reflect the degree of skeletal muscle damage and recovery associated with eccentric exercise. The degree of muscle damage was the lowest in the continuous exercise group and the highest in the once-fatigue exercise group, which may provide more information and guidance for the formulation of physical and athletic training programmes.

Evaluation of interrater reliability of different muscle segmentation techniques in diffusion tensor imaging

INTRODUCTION

Muscle diffusion tensor imaging (mDTI) is a quantitative MRI technique that can provide information about muscular microstructure and integrity. Ultrasound and DTI studies have shown intramuscular differences, and therefore separation of different muscles for analysis is essential. The commonly used methods to assess DTI metrics in muscles are manual segmentation and tract-based analysis. Recently methods such as volume-based tractography have been applied to optimize muscle architecture estimation, but can also be used to assess DTI metrics.

PURPOSE

To evaluate diffusion metrics obtained using three different methods-volume-based tractography, manual segmentation-based analysis and tract-based analysis-with respect to their interrater reliability and their ability to detect intramuscular variance.

MATERIALS AND METHODS

30 volunteers underwent an MRI examination in a 3 T scanner using a 16-channel Torso XL coil. Diffusion-weighted images were acquired to obtain DTI metrics. These metrics were evaluated in six thigh muscles using volume-based tractography, manual segmentation and standard tractography. All three methods were performed by two independent raters to assess interrater reliability by ICC analysis and Bland-Altman plots. Ability to assess intramuscular variance was compared using an ANOVA with muscle as a between-subjects factor.

RESULTS

Interrater reliability for all methods was found to be excellent. The highest interrater reliability was found for volume-based tractography (ICC ≥ 0.967). Significant differences for the factor muscle in all examined diffusion parameters were shown in muscles using all methods (main effect p < 0.001).

CONCLUSIONS

Diffusion data can be assessed by volume tractography, standard tractography and manual segmentation with high interrater reliability. Each method produces different results for the investigated DTI parameters. Volume-based tractography was superior to conventional manual segmentation and tractography regarding interrater reliability and detection of intramuscular variance, while tract-based analysis showed the lowest coefficients of variation.

The effect of ageing on skeletal muscle as assessed by quantitative MR imaging: an association with frailty and muscle strength

Background

Skeletal muscles undergo changes with ageing which can cause sarcopenia that can result in frailty. Quantitative MRI may detect the muscle-deficit component of frailty which could help improve the understanding of ageing muscles.

Aims

To investigate whether quantitative MRI measures of T2, fat fraction (FF), diffusion tensor imaging and muscle volume can detect differences within the muscles between three age groups, and to assess how these measures compare with frailty index, gait speed and muscle power.

Methods

18 ‘young’ (18–30 years), 18 ‘middle-aged’ (31–68 years) and 18 ‘older’ (> 69 years) healthy participants were recruited. Participants had an MRI of their dominant thigh. Knee extension and flexion power and handgrip strength were measured. Frailty (English Longitudinal Study of Ageing frailty index) and gait speed were measured in the older participants.

Results

Young participants had a lower muscle MRI T2, FF and mean diffusivity than middle-aged and older participants; middle-aged participants had lower values than older participants. Young participants had greater muscle flexion and extension power, muscle volume and stronger hand grip than middle-aged and older participants; middle-aged participants had greater values than the older participants. Quantitative MRI measurements correlated with frailty index, gait speed, grip strength and muscle power.

Discussion

Quantitative MRI and strength measurements can detect muscle differences due to ageing. Older participants had raised T2, FF and mean diffusivity and lower muscle volume, grip strength and muscle power.

Conclusions

Quantitative MRI measurements correlate with frailty and muscle function and could be used for identifying differences across age groups within muscle.

Effect of b Value on Imaging Quality for Diffusion Tensor Imaging of the Spinal Cord at Ultrahigh Field Strength

Objective

To explore the optimal b value setting for diffusion tensor imaging of rats' spinal cord at ultrahigh field strength (7 T).

Methods

Spinal cord diffusion tensor imaging data were collected from 14 rats (5 healthy, 9 spinal cord injured) with a series of b values (200, 300, 400, 500, 600, 700, 800, 900, and 1000 s/mm²) under the condition that other scanning parameters were consistent. The image quality (including image signal-to-noise ratio and image distortion degree) and data quality (i.e., the stability and consistency of the DTI-derived parameters, referred to as data stability and data consistency) were quantitatively evaluated. The min-max normalization method was used to process the calculation results of the four indicators. Finally, the image and data quality under each b value were synthesized to determine the optimal b value.

Results

b = 200 s/mm² and b = 900 s/mm² ranked in the top two of the comprehensive evaluation, with the best image quality at b = 200 s/mm² and the best data quality at b = 900 s/mm².

Conclusion

Considering the shortcomings of the ability of low b values to reflect the microstructure, b = 900 s/mm² can be used as the optimal b value for 7 T spinal cord diffusion tensor scanning.

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 T₂ 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 ³¹ 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 water T₂ 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 T₂ SD was found in the majority of FR muscles compared with NFR and healthy control muscles. No differences in average water T₂ 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 T₂ heterogeneity and changes in some metabolic indices. No clear differences were found for the DTI indices between groups.

Quantitative MRI in myositis patients: comparison with healthy volunteers and radiological visual assessment

AIM

To assess whether magnetic resonance imaging (MRI)-based measurements of T2, fat fraction, diffusion tensor imaging, and muscle volume can detect differences between the muscles of myositis patients and healthy controls, and to identify how they compare with semi-quantitative MRI diagnosis.

MATERIALS AND METHODS

Sixteen myositis patients and 16 age- and gender-matched healthy controls underwent MRI of their thigh. Quantitative MRI measurements and radiologists' semi-quantitative scores were assessed. Strength was assessed using an isokinetic dynamometer.

RESULTS

Fat fraction and T2 values were higher in myositis patients whereas muscle volume was lower compared to healthy controls. There was no difference in diffusion. Muscle strength was lower in myositis patients compared to healthy controls. In a subgroup of eight patients, scored as unaffected by radiologists, T2 values were still significantly higher in myositis patients.

CONCLUSIONS

Quantitative MRI measurements can detect differences between myositis patients and healthy controls. Changes in the muscles of myositis patients, undetected by visual, semi-quantitative scoring, can be detected using quantitative T2 measurements. This suggests that MRI T2 values may be useful for the management of myositis patients.

Clinical trial readiness study of distal myopathy and dysphagia in nephropathic cystinosis

BACKGROUND

Nephropathic cystinosis is a lysosomal storage disorder with late-onset systemic complications, such as myopathy and dysphagia. Currently employed outcome measures lack sensitivity and responsiveness for dysphagia and myopathy, a limitation to clinical trial readiness.

METHODS

We evaluated 20 patients with nephropathic cystinosis in two visits over the course of a year to identify outcomes sensitive to detect changes over time. Patients also underwent an expiratory muscle strength training program to assess any effects on aspiration and dysphagia.

RESULTS

There were significant differences in the Timed Up and Go Test (TUG) and Timed 25-Foot Walk (25-FW) between baseline and 1-y follow-up (P < .05). Maximum expiratory pressure (MEP) and peak cough flow (PCF) significantly improved following respiratory training (P < .05).

CONCLUSIONS

Improved respiratory outcomes may enhance patients ability to expel aspirated material from the airway, stave off pulmonary sequelae associated with chronic aspiration, and yield an overall improvement in physical health and well-being.

Diffusion tensor imaging of the anterior cruciate ligament graft following reconstruction: a longitudinal study

OBJECTIVE

To longitudinally monitor remodeling of human autograft following anterior cruciate ligament (ACL) reconstruction with DTI.

METHODS

Twenty-eight patients underwent DTI follow-up at 3, 8, and 14 months after clinically successful ACL reconstruction with tendon autograft. Among these, 18 patients had a concomitant lateral extra-articular procedure (LET). DTI data from 7 healthy volunteers was also obtained. Diffusion parameters (fractional anisotropy, FA; mean diffusivity, MD; axial diffusivity, AD; and radial diffusivity, RD) were evaluated within the fiber tractography volumes of the ACL graft and posterior cruciate ligament (PCL) in all patients. Data were analyzed using a linear mixed-effects model with post hoc testing using Bonferroni-Holm correction for multiple testing. The effect of additional LET was studied.

RESULTS

The ACL graft showed a significant decrease of FA over time (F = 4.00, p = 0.025), while the diffusivities did not significantly change over time. For PCL there were no significant DTI changes over time. A different evolution over time between patients with and without LET was noted for all diffusivity values of the ACL graft with reduced AD values in patients with LET at 8 months postoperatively (p = 0.048; adjusted p = 0.387). DTI metrics of the ACL graft differed largely from both native ACL and tendon at 14 months postoperatively.

CONCLUSION

Our study has shown the potential of DTI to longitudinally monitor the remodeling process in human ACL reconstruction. DTI analysis indicates that graft remodeling is incomplete at 14 months postoperatively.

KEY POINTS

• DTI can be used to longitudinally monitor the remodeling process in human ACL reconstruction. • DTI analysis indicates that autograft remodeling is incomplete at 14 months postoperatively. • DTI may be helpful for evaluating new ACL treatments.

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

The focus of this review is the application of advanced MRI to study the effect of aging and disuse related remodeling of the extracellular matrix (ECM) on force transmission in the human musculoskeletal system. Structural MRI includes (i) ultra-low echo times (UTE) maps to visualize and quantify the connective tissue, (ii) diffusion tensor imaging (DTI) modeling to estimate changes in muscle and ECM microstructure, and (iii) magnetization transfer contrast imaging to quantify the macromolecular fraction in muscle. Functional MRI includes dynamic acquisitions during contraction cycles enabling computation of the strain tensor to monitor muscle deformation. Further, shear strain extracted from the strain tensor may be a potential surrogate marker of lateral transmission of force. Biochemical and histological analysis of muscle biopsy samples can provide "gold-standard" validation of some of the MR findings. The review summarizes biochemical studies of ECM adaptations with age and with disuse. A brief summary of animal models is included as they provide experimental confirmation of longitudinal and lateral force transmission pathways. Computational muscle models enable exploration of force generation and force pathways and elucidate the link between structural adaptations and functional consequences. MR image findings integrated in a computational model can explain and predict subject specific functional changes to structural adaptations. Future work includes development and validation of MRI biomarkers using biochemical analysis of muscle tissue as a reference standard and potential translation of the imaging markers to the clinic to noninvasively monitor musculoskeletal disease conditions and changes consequent to rehabilitative interventions.

Quantitative MRI Reveals Microstructural Changes in the Upper Leg Muscles After Running a Marathon

Background

The majority of sports‐related injuries involve skeletal muscle. Unlike acute trauma, which is often caused by a single traumatic event leading to acute symptoms, exercise‐induced microtrauma may remain subclinical and difficult to detect. Therefore, novel methods to detect and localize subclinical exercise‐induced muscle microtrauma are desirable.

Purpose

To assess acute and delayed microstructural changes in upper leg muscles with multiparametric quantitative MRI after running a marathon.

Study Type

Longitudinal; 1‐week prior, 24–48 hours postmarathon and 2‐week follow‐up

Population

Eleven men participants (age: 47–68 years).

Field Strength/Sequence

Spin‐echo echo planar imaging (SE‐EPI) with diffusion weighting, multispin echo, Dixon, and fat‐suppressed turbo spin‐echo (TSE) sequences at 3T. MR datasets and creatine kinase (CK) concentrations were obtained at three timepoints.

Assessment

Diffusion parameters, perfusion fractions, and quantitative (q)T₂ values were determined for hamstring and quadriceps muscles, TSE images were scored for acute injury. The vastus medialis and biceps femoris long head muscles were divided and analyzed in five segments to assess local damage.

Statistical Tests

Differences between timepoints in MR parameters were assessed with a multilevel linear mixed model and in CK concentrations with a Friedman test. Mean diffusivity (MD) and qT₂ for whole muscle and muscle segments were compared using a multivariate analysis of covariance (MANCOVA).

Results

CK concentrations were elevated (1194 U/L [166–3906], P < 0.001) at 24–48 hours postmarathon and returned to premarathon values (323 U/L [56–2216]) at 2‐week follow‐up. Most of the MRI diffusion indices in muscles without acute injury changed at 24–48 hours postmarathon and returned to premarathon values at follow‐up (MD, RD, and λ3; P < 0.006). qT₂ values (P = 0.003) and perfusion fractions (P = 0.003) were higher at baseline compared to follow‐up. Local assessments of MD and qT₂ revealed more pronounced changes than whole muscle assessment (2–3‐fold; P < 0.01).

Data Conclusion

Marathon running‐induced microtrauma was detected with MRI in individual whole upper leg muscles and even more pronounced on local segments.

Level of Evidence

2

Technical Efficacy Stage

3 J. Magn. Reson. Imaging 2020;52:407–417.

Diffusion tensor imaging combined with T2 mapping to quantify changes in the skeletal muscle associated with training and endurance exercise in competitive triathletes

OBJECTIVES

Diffusion tensor magnetic resonance imaging (DTI) and T2 mapping enable the detection of exercise-induced changes in the skeletal muscle microenvironment. This study prospectively quantified DTI metrics and T2 relaxation times of thigh muscles in competitive triathletes at rest and following a triathlon race in comparison with sedentary controls.

METHODS

Twenty-two triathletes (males N = 16, females N = 6) and twenty-three controls (males N = 16, females N = 7) underwent magnetic resonance imaging (MRI) on a 3 T system at baseline (time point 1; 72 h at rest). Twelve triathletes (males N = 8, females N = 4) underwent a second scan (time point 2; 3 h of completing a triathlon race). The tensor eigenvalues (λ₁, λ₂, λ₃), mean diffusivity (MD), fractional anisotropy (FA), and T2 times were compared between controls and triathletes at time point 1 and triathletes at time points 1 and 2 using independent and paired t tests.

RESULTS

In comparison with the controls at time point 1, the T2 times of rectus femoris (RF, p < 0.02), adductor magnus (AM, p = 0.02), biceps femoris (BF, p < 0.001), semitendinosus (ST, p = 0.005), and semimembranosus (SM, p = 0.003) muscles were significantly increased in triathletes. At time point 2 in triathletes, the average tensor metrics (MD, λ₃/ λ₁) of BF, ST, and SM muscles increased (p < 0.05) and FA values in ST and SM muscles decreased (p < 0.03). T2 times were not significantly changed between both time points in triathletes.

CONCLUSION

Our results indicate that this multiparametric MRI protocol allows detection and quantification of changes in the skeletal muscle microenvironment caused by endurance training and acute strenuous exercise.

KEY POINTS

• Endurance training results in changes to the skeletal microstructure, which can be quantified using MRI-based diffusion tensor imaging. • The combined application of MRI diffusion tensor imaging and T2 mapping allows the differentiation of microstructural changes caused by active exercise or endurance training. • Environmental adaptations of the skeletal muscle caused by physical training are influenced by gender.

A 12-channel flexible receiver coil for accelerated tongue imaging

Objective

MRI of the tongue requires acceleration to minimise motion artefacts and to facilitate real-time imaging of swallowing. To accelerate tongue MRI, we designed a dedicated flexible receiver coil.

Materials and methods

We designed a flexible 12-channel receiver coil for tongue MRI at 3T and compared it to a conventional head-and-neck coil regarding SNR and g-factor. Furthermore, two accelerated imaging techniques were evaluated using both coils: multiband (MB) diffusion-tensor imaging (DTI) and real-time MRI of swallowing.

Results

The flexible coil had significantly higher SNR in the anterior (2.1 times higher, P = 0.002) and posterior (2.0 times higher, P < 0.001) parts of the tongue, while the g-factor was lower at higher acceleration. Unlike for the flexible coil, the apparent diffusion coefficient (P = 0.001) and fractional anisotropy (P = 0.008) deteriorated significantly while using the conventional coil after accelerating DTI with MB. The image quality of real-time MRI of swallowing was significantly better for hyoid elevation (P = 0.029) using the flexible coil.

Conclusion

Facilitated by higher SNR and lower g-factor values, our flexible tongue coil allows faster imaging, which was successfully demonstrated in MB DTI and real-time MRI of swallowing.

Electronic supplementary material

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

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

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

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.

Consensus-based technical recommendations for clinical translation of renal diffusion-weighted MRI

Objectives

Standardization is an important milestone in the validation of DWI-based parameters as imaging biomarkers for renal disease. Here, we propose technical recommendations on three variants of renal DWI, monoexponential DWI, IVIM and DTI, as well as associated MRI biomarkers (ADC, D, D*, f, FA and MD) to aid ongoing international efforts on methodological harmonization.

Materials and methods

Reported DWI biomarkers from 194 prior renal DWI studies were extracted and Pearson correlations between diffusion biomarkers and protocol parameters were computed. Based on the literature review, surveys were designed for the consensus building. Survey data were collected via Delphi consensus process on renal DWI preparation, acquisition, analysis, and reporting. Consensus was defined as ≥ 75% agreement.

Results

Correlations were observed between reported diffusion biomarkers and protocol parameters. Out of 87 survey questions, 57 achieved consensus resolution, while many of the remaining questions were resolved by preference (65–74% agreement). Summary of the literature and survey data as well as recommendations for the preparation, acquisition, processing and reporting of renal DWI were provided.

Discussion

The consensus-based technical recommendations for renal DWI aim to facilitate inter-site harmonization and increase clinical impact of the technique on a larger scale by setting a framework for acquisition protocols for future renal DWI studies. We anticipate an iterative process with continuous updating of the recommendations according to progress in the field.

Electronic supplementary material

The online version of this article (10.1007/s10334-019-00790-y) contains supplementary material, which is available to authorized users.

Denoising of diffusion MRI improves peripheral nerve conspicuity and reproducibility

BACKGROUND

Quantitative diffusion MRI is a promising technique for evaluating peripheral nerve integrity but low signal-to-noise ratio (SNR) can impede measurement accuracy.

PURPOSE

To evaluate principal component analysis (PCA) and generalized spherical deconvolution (genSD) denoising techniques to improve within-subject reproducibility and peripheral nerve conspicuity.

STUDY TYPE

Prospective.

SUBJECTS

Seven healthy volunteers and three peripheral neuropathy patients.

FIELD STRENGTH/SEQUENCE

3T/multiband single-shot echo planar diffusion sequence using multishell 55-direction scheme.

ASSESSMENT

Images were processed using four methods: "original" (no denoising), "average" (10 repetitions), "PCA-only," and "PCA + genSD." Tibial and common peroneal nerve segmentations and masks were generated from volunteer diffusion data. Quantitative (SNR and contrast-to-noise ratio [CNR]) values were calculated. Three radiologists qualitatively evaluated nerve conspicuity for each method. The two denoising methods were also performed in three patients with peripheral neuropathies.

STATISTICAL TESTS

For healthy volunteers, calculations included SNR and CNR_FA (computed using FA values). Coefficient of variation (CV%) of CNR_FA quantified within-subject reproducibility. Groups were compared with two-sample t-tests (significance P < 0.05; two-tailed, Bonferroni-corrected). Odds ratios (ORs) quantified the relative rates of each of three radiologists confidently identifying a nerve, per slice, for the four methods.

RESULTS

"PCA + genSD" yielded the highest SNR (mean_overall = 14.83 ± 1.99) and tibial and common peroneal nerve CNR_FA (mean_tibial = 3.45, mean_peroneal = 2.34) compared to "original" (P _SNR < 0.001; P _CNR = 0.011) and "PCA-only" (P _SNR < 0.001, P _CNR < 0.001). "PCA + genSD" had higher within-subject reproducibility (low CV%) for tibial (6.04 ± 1.98) and common peroneal nerves (8.27 ± 2.75) compared to "original" and "PCA-only." The mean FA was higher for "original" than "average" (P < 0.001), but did not differ significantly between "average" and "PCA + genSD" (P = 0.14). "PCA + genSD" had higher tibial and common peroneal nerve conspicuity than "PCA-only" (OR_tibial = 2.50, P < 0.001; OR_peroneal = 1.86, P < 0.001) and "original" (OR_tibial = 2.73, P < 0.001; OR_peroneal = 2.43, P < 0.001).

DATA CONCLUSION

PCA + genSD denoising method improved SNR, CNR_FA , and within-subject reproducibility (CV%) without biasing FA and nerve conspicuity. This technique holds promise for facilitating more reliable, unbiased diffusion measurements of peripheral nerves.

LEVEL OF EVIDENCE

2 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2020;51:1128-1137.

Clinical myopathy in patients with nephropathic cystinosis

BACKGROUND

Nephropathic cystinosis is a lysosomal storage disorder. Patient survival years after renal transplantation has revealed systemic complications including distal myopathy and dysphagia.

METHODS

We evaluated 20 adult patients with nephropathic cystinosis using patient-reported and clinical outcome measures. Standard motor measures, video fluoroscopy swallow studies, and tests of respiratory function were performed. We also used Rasch analysis of an initial survey to design a 16-item survey focused on upper and lower extremity function, which was completed by 31 additional patients.

RESULTS

Distal myopathy and dysphagia were common in patients with nephropathic cystinosis. Muscle weakness ranges from mild involvement of intrinsic hand muscles to prominent distal greater than proximal weakness and contractures.

CONCLUSIONS

In addition to further characterization of underlying dysphagia and muscle weakness, we propose a new psychometrically devised, disease specific, functional outcome measures for distal myopathy in patients with nephropathic cystinosis.

Magnetic resonance imaging of the cervical spinal cord in spinal muscular atrophy

Objective

In this study we investigated the potential value of magnetic resonance imaging (MRI) and diffusion tensor imaging (DTI) in characterizing changes in the cervical spinal cord and peripheral nerve roots in vivo in patients with spinal muscular atrophy (SMA).

Methods

We developed an MRI protocol with 4 sequences to investigate the cervical spinal cord and nerve roots on a 3 Tesla MRI system. We used 2 anatomical MRI sequences to investigate cross-sectional area (CSA) at each spinal segment and the diameter of ventral and dorsal nerve roots, and two diffusion tensor imaging (DTI) techniques to estimate the fractional anisotropy (FA), mean (MD), axial (AD) and radial diffusivity (RD) in 10 SMA patients and 20 healthy controls.

Results

There were no significant differences in CSA (p > .1), although an 8.5% reduction of CSA in patients compared to healthy controls was apparent at segment C7. DTI data showed a higher AD in grey matter of patients compared to healthy controls (p = .033). Significantly lower MD, AD and RD values were found in rostral nerve roots (C3-C5) in patients (p < .045).

Conclusions

We showed feasibility of an advanced 3 T MRI protocol that allowed differences to be determined between patients and healthy controls, confirming the potential of this technique to assess pathological mechanisms in SMA. After further development and confirmation of findings in a larger sample, these techniques may be used to study disease course of SMA in vivo and evaluate response to survival motor neuron (SMN) augmenting therapy.

•

The developed MRI sequences measure (micro)structural spinal cord changes in SMA.

•

cervical spinal cross-sectional area is overall (non-significantly) smaller in SMA.

•

In nerve roots C3-C8 all DTI parameters were lower in patients compared to controls.

•

Largest differences in DTI parameters were located at the rostral cervical segments.

Crossing muscle fibers of the human tongue resolved in vivo using constrained spherical deconvolution

BACKGROUND

Surgical resection of tongue cancer may impair swallowing and speech. Knowledge of tongue muscle architecture affected by the resection could aid in patient counseling. Diffusion tensor imaging (DTI) enables reconstructions of muscle architecture in vivo. Reconstructing crossing fibers in the tongue requires a higher-order diffusion model.

PURPOSE

To develop a clinically feasible diffusion imaging protocol, which facilitates both DTI and constrained spherical deconvolution (CSD) reconstructions of tongue muscle architecture in vivo.

STUDY TYPE

Cross-sectional study.

SUBJECTS/SPECIMEN

One ex vivo bovine tongue resected en bloc from mandible to hyoid bone. Ten healthy volunteers (mean age 25.5 years; range 21-34 years; four female).

FIELD STRENGTH/SEQUENCE

Diffusion-weighted echo planar imaging at 3 T using a high-angular resolution diffusion imaging scheme acquired twice with opposing phase-encoding for B₀ -field inhomogeneity correction. The scan of the healthy volunteers was divided into four parts, in between which the volunteers were allowed to swallow, resulting in a total acquisition time of 10 minutes.

ASSESSMENT

The ability of resolving crossing muscle fibers using CSD was determined on the bovine tongue specimen. A reproducible response function was estimated and the optimal peak threshold was determined for the in vivo tongue. The quality of tractography of the in vivo tongue was graded by three experts.

STATISTICAL TESTS

The within-subject coefficient of variance was calculated for the response function. The qualitative results of the grading of DTI and CSD tractography were analyzed using a multilevel proportional odds model.

RESULTS

Fiber orientation distributions in the bovine tongue specimen showed that CSD was able to resolve crossing muscle fibers. The response function could be determined reproducibly in vivo. CSD tractography displayed significantly improved tractography compared with DTI tractography (P = 0.015).

DATA CONCLUSION

The 10-minute diffusion imaging protocol facilitates CSD fiber tracking with improved reconstructions of crossing tongue muscle fibers compared with DTI.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:96-105.