Diffusion-Tensor MRI Based Skeletal Muscle Fiber Tracking

Is a 3D representation of muscle architecture needed to model craniomaxillofacial skeletal mechanics?

The craniomaxillofacial skeleton (CMFS) is sensitive to the direction and magnitude of loading, making experimental testing challenging due to the complex network of musculature and thin bone. Finite element (FE) models have been used to characterize their mechanical behavior, often utilizing link elements to simulate muscles. Utilizing multimodal imaging, a specimen-specific CMFS FE model incorporating 3D masseter geometry and fiber directions was developed. 3D representation of the masseters resulted in lower peak intensity and smoother strain distribution in the zygomatic region, suggesting link muscle modeling may not sufficiently capture complex load transfer from muscle to bone in CMFS FE models.

Imaging and Image Processing Techniques for High-Resolution Visualization of Connective Tissue with MRI: Application to Fascia, Aponeurosis, and Tendon

Recent interest in musculoskeletal connective tissues like tendons, aponeurosis, and deep fascia has led to a greater focus on in vivo medical imaging, particularly MRI. Given the rapid T2* decay of collagenous tissues, advanced ultra-short echo time (UTE) MRI sequences have proven useful in generating high-signal images of these tissues. To further these advances, we discuss the integration of UTE with Diffusion Tensor Imaging (DTI) and explore image processing techniques to enhance the localization, labeling, and modeling of connective tissues. These techniques are especially valuable for extracting features from thin tissues that may be difficult to distinguish. We present data from lower leg scans of 30 healthy subjects using a non-Cartesian MRI sequence to acquire axial 2D images to segment skeletal muscle and connective tissue. DTI helped differentiate aponeurosis from deep fascia by analyzing muscle fiber orientations. The dual echo imaging methods yielded high-resolution images of deep fascia, where in-plane spatial resolutions were between 0.3 × 0.3 mm to 0.5 × 0.5 mm with a slice thickness of 3-5 mm. Techniques such as K-Means clustering, FFT edge detection, and region-specific scaling were most effective in enhancing images of deep fascia, aponeurosis, and tendon to enable high-fidelity modeling of these tissues.

Development of a three-dimensional muscle-driven lower limb model developed using an improved CFD-FE method

Analysis of the musculoskeletal movements (gait analysis) is needed in many scenarios. The in vivo method has some difficulties. For example, recruiting human subjects for the gait analysis is challenging due to many issues. In addition, when plenty of subjects are required, the follow-up experiments take a long period and the dropout of subjects always occurs. An efficient and reliable in silico simulation platform for gait analysis has been desired for a long time. Therefore, a technique using three-dimensional (3D) muscle modeling to drive the 3D musculoskeletal model was developed and the application of the technique in the simulation of lower limb movements was demonstrated. A finite element model of the lower limb with anatomically high fidelity was developed from the MRI data, where the main muscles, the bones, the subcutaneous tissues, and the skin were reconstructed. To simulate the active behavior of 3D muscles, an active, fiber-reinforced hyperelastic muscle model was developed using the user-defined material (VUMAT) model. Two typical movements, that is, hip abduction and knee lifting, were simulated by activating the responsible muscles. The results show that it is reasonable to use the improved CFD-FE method proposed in the present study to simulate the active contraction of the muscle, and it is feasible to simulate the movements by activating the relevant muscles. The results from the present technique closely match the physiological scenario and thus the technique developed has a great potential to be used in the in silico human simulation platform for many purposes.

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 6-week randomized-controlled field study: effect of isokinetic eccentric resistance training on strength, flexibility and muscle structure of the shoulder external rotators in male junior handball players

Background: Team handball involves a tremendous amount of shoulder motion with high forces during repeated extended external range of motion. This causes shoulder complaints and overuse injuries. While eccentric training for the lower extremity shows preventive effects by improving strength, range of motion and fascicle length, there is a research gap for the shoulder joint and for advanced tissue characterization using diffusion tensor imaging. Objectives: To investigate the effects of 6-week eccentric isokinetic resistance training on strength, flexibility, and fiber architecture characteristics of the external rotators compared to an active control group in junior male handball players. Methods: 15 subjects were randomly assigned to the eccentric training group and 14 subjects to the active control group (conventional preventive training). Primary outcome measures were eccentric and concentric isokinetic strength of the external rotators, range of motion, and muscle fascicle length and fascicle volume. Results: The intervention group, showed significant changes in eccentric strength (+15%). The supraspinatus and infraspinatus muscles showed significant increases in fascicle length (+13% and +8%), and in fractional anisotropy (+9% and +6%), which were significantly different from the control group. Conclusion: Eccentric isokinetic training has a significant effect on the function and macroscopic structure of the shoulder external rotators in male junior handball players. While strength parameters and muscle structure improved, range of motion did not change. This research helps understanding the physiology of muscle and the role of eccentric training on shoulder function and muscle structure. Furthermore, DTI was found to be a promising tool for advanced tissue characterization, and the in vivo derived data can also serve as model input variables and as a possibility to extend existing ex-vivo muscle models. Future research is needed for functional and structural changes following convenient eccentric field exercises.

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.

A systematic review of the relationship between muscle oxygen dynamics and energy rich phosphates. Can NIRS help?

BACKGROUND

Phosphocreatine dynamics provide the gold standard evaluation of in-vivo mitochondrial function and is tightly coupled with oxygen availability. Low mitochondrial oxidative capacity has been associated with health issues and low exercise performance.

METHODS

To evaluate the relationship between near-infrared spectroscopy-based muscle oxygen dynamics and magnetic resonance spectroscopy-based energy-rich phosphates, a systematic review of the literature related to muscle oxygen dynamics and energy-rich phosphates was conducted. PRISMA guidelines were followed to perform a comprehensive and systematic search of four databases on 02-11-2021 (PubMed, MEDLINE, Scopus and Web of Science). Beforehand pre-registration with the Open Science Framework was performed. Studies had to include healthy humans aged 18-55, measures related to NIRS-based muscle oxygen measures in combination with energy-rich phosphates. Exclusion criteria were clinical populations, laboratory animals, acutely injured subjects, data that only assessed oxygen dynamics or energy-rich phosphates, or grey literature. The Effective Public Health Practice Project Quality Assessment Tool was used to assess methodological quality, and data extraction was presented in a table.

RESULTS

Out of 1483 records, 28 were eligible. All included studies were rated moderate. The studies suggest muscle oxygen dynamics could indicate energy-rich phosphates under appropriate protocol settings.

CONCLUSION

Arterial occlusion and exercise intensity might be important factors to control if NIRS application should be used to examine energetics. However, more research needs to be conducted without arterial occlusion and with high-intensity exercises to support the applicability of NIRS and provide an agreement level in the concurrent course of muscle oxygen kinetics and muscle energetics.

TRIAL REGISTRATION

https://osf.io/py32n/ .

KEY POINTS

1. NIRS derived measures of muscle oxygenation agree with gold-standard measures of high energy phosphates when assessed in an appropriate protocol setting. 2. At rest when applying the AO protocol, in the absence of muscle activity, an initial disjunction between the NIRS signal and high energy phosphates can been seen, suggesting a cascading relationship. 3. During exercise and recovery a disruption of oxygen delivery is required to provide the appropriate setting for evaluation through either an AO protocol or high intensity contractions.

Effect of isokinetic eccentric training on the human shoulder strength, flexibility, and muscle architecture in physically active men: A preliminary study

Strengthening the rotator cuff muscles is important for injury prevention and rehabilitation. Since muscle fascicle length improves motor performance and is suggested to reduce the risk of injury for the hamstring, it may be an important variable to promote multidirectional changes in the function and macroscopic structure for the shoulder. Recent literature reviews overwhelmingly suggest that eccentric exercises improve fascicle length and functional measures for the lower limb. However, there is a research gap for the shoulder. Since ultrasound imaging is the most commonly used imaging technique to quantify muscle structure, but has yielded heterogeneous results in different studies, there is another issue and a research gap for the imaging method. Based on the research gaps, the purpose of this study was to evaluate the effects of standardized eccentric strength training on the function and structure of the external rotator cuff muscles using an isokinetic dynamometer and MRI. Therefore, a preliminary pre-post intervention study was conducted and 16 physically active men were recruited in October 2021. For the right shoulder, an eccentric isokinetic training was performed twice a week for almost six weeks. The primary outcome measures (external rotators) were active and passive range of motion, eccentric and concentric torque at 30, 60, and 180°/s isokinetic speed, and fascicle length and fascicle volume for the supraspinatus and infraspinatus muscles. The findings show a training effect for the absolute mean values of eccentric strength (+24%, p = .008). The torque-angle relationship increased, especially in the final phase of range of motion, although a 4% (p = .002) decrease in passive range of motion was found in the stretch test. Positive changes in muscle structure were shown for the supraspinatus muscle fascicle length (+16%, p = .003) and fascicle volume (+19%, p = .002). Based on the study results, we can conclude that eccentric isokinetic training has a significant positive effect on the shoulder. To our knowledge, this is the first eccentric training study using both isokinetic dynamometer and muscle diffusion tensor imaging to access functional and structural changes in the human shoulder rotator cuff muscles. The methods were shown to be applicable for interventional studies. Based on these results, populations such as high-performance handball players with highly trained shoulders should be included in future studies.

Time-dependent diffusion tensor imaging and diffusion modeling of age-related differences in the medial gastrocnemius and feasibility study of correlations to histopathology

PURPOSE

Implement STEAM-DTI to model time-dependent diffusion eigenvalues using the random permeable barrier model (RPBM) to study age-related differences in the medial gastrocnemius (MG) muscle. Validate diffusion model-extracted fiber diameter for histological assessment.

METHODS

Diffusion imaging at different diffusion times (Δ) was performed on seven young and six senior participants. Time-dependent diffusion eigenvalues (λ2 (t), λ3 (t), and D⊥ (t); average of λ2 (t) and λ3 (t)) were fit to the RPBM to extract tissue microstructure parameters. Biopsy of the MG tissue for histological assessment was performed on a subset of participants (four young, six senior).

RESULTS

λ3 (t) was significantly higher in the senior cohort for the range of diffusion times. RPBM fits to λ2 (t) yielded fiber diameters in agreement to those from histology for both cohorts. The senior cohort had lower values of volume fraction of membranes, ζ, in fits to λ2 (t), λ3 (t), and D⊥ (t) (significant for fit to λ3 (t)). Fits of fiber diameter from RPBM to that from histology had the highest correlation for the fit to λ2 (t).

CONCLUSION

The age-related patterns in λ2 (t) and λ3 (t) could tentatively be explained from RPBM fits; these patterns may potentially arise from a decrease in fiber asymmetry and an increase in permeability with age.

Diffusion tensor imaging: Influence of segmentation on fiber tracking in the supraspinatus muscle-An inter-operator reliability analysis

The ability of muscle to generate force depends on its architecture and health condition. MR-based diffusion tensor imaging of muscle (mDTI) is an innovative approach for showing the fiber arrangement for the whole muscle volume. For accurate calculations of fiber metrics, muscle segmentation prior to tractography is regarded as necessary. Since segmentation is known to be operator dependent, it is important to understand how segmentation affects tractography. The aim of this study was to compare the results of deterministic fiber tracking based on muscle models generated by two independent operators. In addition, this study compares the results with a segmentation-free approach. Fifteen subjects underwent mDTI of the right shoulder. The results showed that mDTI can be successfully applied to complex joints such as the human shoulder. Furthermore, operator segmentation did not influence the results of fiber tracking and fascicle length (FL), fiber volume (FV), fractional anisotropy (FA), axial diffusivity (AD), radial diffusivity (RD), and mean diffusivity (MD) showed excellent intraclass correlation estimates (≥ 0.975). As an exploratory approach, the segmentation-free fiber tracking showed significant differences in terms of mean fascicle length. Based on these findings, we conclude that tractography is not sensitive to small deviations in muscle segmentation. Furthermore, it implies that mDTI and automatic segmentation approaches or even a segmentation-free analysis can be considered for evaluation of muscle architecture.

Imaging biomarkers in the idiopathic inflammatory myopathies

Idiopathic inflammatory myopathies (IIMs) are a group of acquired muscle diseases with muscle inflammation, weakness, and other extra-muscular manifestations. IIMs can significantly impact the quality of life, and management of IIMs often requires a multi-disciplinary approach. Imaging biomarkers have become an integral part of the management of IIMs. Magnetic resonance imaging (MRI), muscle ultrasound, electrical impedance myography (EIM), and positron emission tomography (PET) are the most widely used imaging technologies in IIMs. They can help make the diagnosis and assess the burden of muscle damage and treatment response. MRI is the most widely used imaging biomarker of IIMs and can assess a large volume of muscle tissue but is limited by availability and cost. Muscle ultrasound and EIM are easy to administer and can even be performed in the clinical setting, but they need further validation. These technologies may complement muscle strength testing and laboratory studies and provide an objective assessment of muscle health in IIMs. Furthermore, this is a rapidly progressing field, and new advances are going to equip care providers with a better objective assessment of IIMS and eventually improve patient management. This review discusses the current state and future direction of imaging biomarkers in IIMs.

A clustering-based method for estimating pennation angle from B-mode ultrasound images

B-mode ultrasound (US) is often used to noninvasively measure skeletal muscle architecture, which contains human intent information. Extracted features from B-mode images can help improve closed-loop human-robotic interaction control when using rehabilitation/assistive devices. The traditional manual approach to inferring the muscle structural features from US images is laborious, time-consuming, and subjective among different investigators. This paper proposes a clustering-based detection method that can mimic a well-trained human expert in identifying fascicle and aponeurosis and, therefore, compute the pennation angle. The clustering-based architecture assumes that muscle fibers have tubular characteristics. It is robust for low-frequency image streams. We compared the proposed algorithm to two mature benchmark techniques: UltraTrack and ImageJ. The performance of the proposed approach showed higher accuracy in our dataset (frame frequency is 20 Hz), that is, similar to the human expert. The proposed method shows promising potential in automatic muscle fascicle orientation detection to facilitate implementations in biomechanics modeling, rehabilitation robot control design, and neuromuscular disease diagnosis with low-frequency data stream.

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.

Automated Method for Tracking Human Muscle Architecture on Ultrasound Scans during Dynamic Tasks

Existing approaches for automated tracking of fascicle length (FL) and pennation angle (PA) rely on the presence of a single, user-defined fascicle (feature tracking) or on the presence of a specific intensity pattern (feature detection) across all the recorded ultrasound images. These prerequisites are seldom met during large dynamic muscle movements or for deeper muscles that are difficult to image. Deep-learning approaches are not affected by these issues, but their applicability is restricted by their need for large, manually analyzed training data sets. To address these limitations, the present study proposes a novel approach that tracks changes in FL and PA based on the distortion pattern within the fascicle band. The results indicated a satisfactory level of agreement between manual and automated measurements made with the proposed method. When compared against feature tracking and feature detection methods, the proposed method achieved the lowest average root mean squared error for FL and the second lowest for PA. The strength of the proposed approach is that the quantification process does not require a training data set and it can take place even when it is not possible to track a single fascicle or observe a specific intensity pattern on the ultrasound recording.

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

BACKGROUND

The purpose of this study was to develop a DTI-based method to quantitatively assess fiber angles and changes therein in leg muscles in order to facilitate longitudinal studies on muscle fiber architectural adaptations in healthy subjects.

METHODS

The upper legs of five volunteers were scanned twice on the same day. The right lower legs of five volunteers were scanned twice with the ankle in three positions, i.e. -15° dorsiflexion, 0° neutral position, and 30° plantarflexion. The MRI protocols consisted of a noise scan, a 3-point mDixon scan and a DTI scan. Fiber-angle color maps were generated for four muscles in the upper legs and two muscles in the lower leg. Voxel-wise fiber angles (θ) were calculated from the angle between the principal eigenvector of the diffusion tensor and a reference line defined between the origo and insertion points of each muscle. Bland-Altman analysis, intraclass correlation coefficient (ICC), coefficient of variation (CV%), minimal detectable change (MDC), standard error (SE) and Friedman test were used for assessing the feasibility of this method and in order to have an indication of the repeatability and the sensitivity.

RESULTS

Bland-Altman analysis showed good repeatability (CV%<10 and 0.7≤ICC≤0.9) with exception of the Tibialis Anterior (TA) muscle in dorsiflexion position(CV%: 12.2) and the Semitendinosus (ST) muscle (left leg) (CV%: 11.4). The best repeatability metrics were found for the SOL muscle in neutral position (CV%: 2.6). Changes in average θ in TA and SOL with ankle positions were observed in accordance with expected agonist and antagonist functions of both muscles. For example, for the anterior left compartment the change in fiber angle Δθ with respect to the neutral position Δθ = -1.6° ± 0.8° and 2.2° ± 2.8° (p = 0.008), for dorsiflexion and plantarflexion, respectively.

CONCLUSION

Our method facilitates fast inspection and quantification of muscle fiber angles in the lower and upper leg muscles in rest and detection of changes in lower-leg muscle fiber angles with varying ankle angles.

Discovering fiber type architecture over the entire muscle using data-driven analysis

Skeletal muscle function is inferred from the spatial arrangement of muscle fiber architecture, which corresponds to myofiber molecular and metabolic features. Myofiber features are often determined using immunofluorescence on a local sampling, typically obtained from a median region. This median region is assumed to represent the entire muscle. However, it remains largely unknown to what extent this local sampling represents the entire muscle. We present a pipeline to study the architecture of muscle fiber features over the entire muscle, including sectioning, staining, imaging to image quantification and data‐driven analysis with Myofiber type were identified by the expression of myosin heavy chain (MyHC) isoforms, representing contraction properties. We reconstructed muscle architecture from consecutive cross‐sections stained for laminin and MyHC isoforms. Examining the entire muscle using consecutive cross‐sections is extremely laborious, we provide consideration to reduce the dataset without loosing spatial information. Data‐driven analysis with over 150,000 myofibers showed spatial variations in myofiber geometric features, myofiber type, and the distribution of neuromuscular junctions over the entire muscle. We present a workflow to study histological changes over the entire muscle using high‐throughput imaging, image quantification, and data‐driven analysis. Our results suggest that asymmetric spatial distribution of these features over the entire muscle could impact muscle function. Therefore, instead of a single sampling from a median region, representative regions covering the entire muscle should be investigated in future studies.

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.

Established paths and new avenues: a review of the main radiological techniques for investigating sarcopenia

Sarcopenia is a clinical condition mainly affecting the elderly that can be associated in a long run with severe consequences like malnutrition and frailty. Considering the progressive ageing of the world population and the socio-economic impact of this disease, much effort is devoted and has to be further focused on an early and accurate diagnostic assessment of muscle loss. Currently, several radiological techniques can be applied for evaluating sarcopenia. If dual-energy X-ray absorptiometry (DXA) is still considered the main tool and it is even recommended as reference by the most current guidelines of the European working group on sarcopenia in older people (EWGSOP), the role of ultrasound (US), computed tomography (CT), peripheral quantitative CT (pQCT), and magnetic resonance imaging (MRI) should not be overlooked. Indeed, such techniques can provide robust qualitative and quantitative information. In particular, regarding MRI, the use of sequences like diffusion-weighted imaging (DWI), diffusion tensor imaging (DTI), magnetic resonance spectroscopy (MRS) and mapping that could provide further insights into the physiopathological features of sarcopenia, should be fostered. In an era pointing to the quantification and automatic evaluation of diseases, we call for future research extending the application of organ tailored protocols, taking advantage of the most recent technical developments.

In vivo muscle morphology comparison in post-stroke survivors using ultrasonography and diffusion tensor imaging

Skeletal muscle architecture significantly influences the performance capacity of a muscle. A DTI-based method has been recently considered as a new reference standard to validate measurement of muscle structure in vivo. This study sought to quantify muscle architecture parameters such as fascicle length (FL), pennation angle (PA) and muscle thickness (tm) in post-stroke patients using diffusion tensor imaging (DTI) and to quantitatively compare the differences with 2D ultrasonography (US) and DTI. Muscle fascicles were reconstructed to examine the anatomy of the medial gastrocnemius, posterior soleus and tibialis anterior in seven stroke survivors using US- and DTI-based techniques, respectively. By aligning the US and DTI coordinate system, DTI reconstructed muscle fascicles at the same scanning plane of the US data can be identified. The architecture parameters estimated based on two imaging modalities were further compared. Significant differences were observed for PA and tm between two methods. Although mean FL was not significantly different, there were considerable intra-individual differences in FL and PA. On the individual level, parameters measured by US agreed poorly with those from DTI in both deep and superficial muscles. The significant differences in muscle parameters we observed suggested that the DTI-based method seems to be a better method to quantify muscle architecture parameters which can provide important information for treatment planning and to personalize a computational muscle model.

Multimodal imaging assessment and histologic correlation of the female rat pelvic floor muscles' anatomy

Pelvic floor disorders negatively impact millions of women worldwide. Although there is a strong epidemiological association with childbirth, the mechanisms leading to the dysfunction of the integral constituents of the female pelvic floor, including pelvic floor skeletal muscles, are not well understood. This is in part due to the constraints associated with directly probing these muscles, which are located deep in the pelvis. Thus, experimental models and non-invasive techniques are essential for advancing knowledge of various phenotypes of pelvic floor muscle injury and pathogenesis of muscle dysfunction, as well as developing minimally invasive approaches for the delivery of novel therapeutics. The most widely used animal model for pelvic floor disorders is the rat. However, the radiological anatomy of rat pelvic floor muscles has not been described. To remedy this gap, the current study provides the first detailed description of the female rat pelvic floor muscles' radiological appearance on MR and ultrasound images, validated by correlation with gross anatomy and histology. We also demonstrate that ultrasound guidance can be used to target rat pelvic floor muscles for possible interventional therapies.

Diffusion tensor imaging and diffusion modeling: Application to monitoring changes in the medial gastrocnemius in disuse atrophy induced by unilateral limb suspension

BACKGROUND

Diffusion tensor imaging (DTI) assesses underlying tissue microstructure, and has been applied to studying skeletal muscle. Unloading of the lower leg causes decreases in muscle force, mass, and muscle protein synthesis as well as changes in muscle architecture.

PURPOSE

To monitor the change in DTI indices in the medial gastrocnemius (MG) after 4-week unilateral limb suspension (ULLS) and to explore the feasibility of extracting tissue microstructural parameters based on a two-compartment diffusion model.

STUDY TYPE

Prospective cohort study.

SUBJECTS

Seven moderately active subjects (29.1 ± 5.7 years).

FIELD STRENGTH/SEQUENCE

3T, single-shot fat-suppressed echo planar spin echo sequence.

ASSESSMENT

Suspension-related changes in the DTI indices (eigenvalues: λ₁ , λ₂ , λ₃ , fractional anisotropy; coefficient of planarity) were statistically analyzed. Changes in model-derived tissue parameters (muscle fiber circularity and diameter, intracellular volume fraction, and residence time) after suspension are qualitatively discussed.

STATISTICAL TESTS

Changes in the DTI indices of the MG between pre- and postsuspension were assessed using repeated-measures two-way analysis of variance (ANOVA).

RESULTS

All the eigenvalues (λ₁ : P = 0.025, λ₂ : P = 0.035, λ₃ : P = 0.049) as well as anisotropic diffusion coefficient (P = 0.029) were significantly smaller post-ULLS. Diffusion modeling revealed that fibers were more circular (circularity index increased from 0.55 to 0.95) with a smaller diameter (diameter decreased from 82-60 μm) postsuspension.

DATA CONCLUSION

We have shown that DTI indices change with disuse and modeling can relate these voxel level changes to changes in the tissue microarchitecture.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2018.

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.

Weighted Mean of Signal Intensity for Unbiased Fiber Tracking of Skeletal Muscles: Development of a New Method and Comparison With Other Correction Techniques

OBJECTIVES

The aim of this study was to investigate the origin of random image artifacts in stimulated echo acquisition mode diffusion tensor imaging (STEAM-DTI), assess the role of averaging, develop an automated artifact postprocessing correction method using weighted mean of signal intensities (WMSIs), and compare it with other correction techniques.

MATERIALS AND METHODS

Institutional review board approval and written informed consent were obtained. The right calf and thigh of 10 volunteers were scanned on a 3 T magnetic resonance imaging scanner using a STEAM-DTI sequence.Artifacts (ie, signal loss) in STEAM-based DTI, presumably caused by involuntary muscle contractions, were investigated in volunteers and ex vivo (ie, human cadaver calf and turkey leg using the same DTI parameters as for the volunteers). An automated postprocessing artifact correction method based on the WMSI was developed and compared with previous approaches (ie, iteratively reweighted linear least squares and informed robust estimation of tensors by outlier rejection [iRESTORE]). Diffusion tensor imaging and fiber tracking metrics, using different averages and artifact corrections, were compared for region of interest- and mask-based analyses. One-way repeated measures analysis of variance with Greenhouse-Geisser correction and Bonferroni post hoc tests were used to evaluate differences among all tested conditions. Qualitative assessment (ie, images quality) for native and corrected images was performed using the paired t test.

RESULTS

Randomly localized and shaped artifacts affected all volunteer data sets. Artifact burden during voluntary muscle contractions increased on average from 23.1% to 77.5% but were absent ex vivo. Diffusion tensor imaging metrics (mean diffusivity, fractional anisotropy, radial diffusivity, and axial diffusivity) had a heterogeneous behavior, but in the range reported by literature. Fiber track metrics (number, length, and volume) significantly improved in both calves and thighs after artifact correction in region of interest- and mask-based analyses (P < 0.05 each). Iteratively reweighted linear least squares and iRESTORE showed equivalent results, but WMSI was faster than iRESTORE. Muscle delineation and artifact load significantly improved after correction (P < 0.05 each).

CONCLUSIONS

Weighted mean of signal intensity correction significantly improved STEAM-based quantitative DTI analyses and fiber tracking of lower-limb muscles, providing a robust tool for musculoskeletal applications.

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.

Quantitative Magnetic Resonance Imaging of Skeletal Muscle Disease

Quantitative magnetic resonance imaging (qMRI) describes the development and use of MRI to quantify physical, chemical, and/or biological properties of living systems. Neuromuscular diseases often exhibit a temporally varying, spatially heterogeneous, and multi-faceted pathology. The goal of this protocol is to characterize this pathology using qMRI methods. The MRI acquisition protocol begins with localizer images (used to locate the position of the body and tissue of interest within the MRI system), quality control measurements of relevant magnetic field distributions, and structural imaging for general anatomical characterization. The qMRI portion of the protocol includes measurements of the longitudinal and transverse relaxation time constants (T1 and T2, respectively). Also acquired are diffusion-tensor MRI data, in which water diffusivity is measured and used to infer pathological processes such as edema. Quantitative magnetization transfer imaging is used to characterize the relative tissue content of macromolecular and free water protons. Lastly, fat-water MRI methods are used to characterize fibro-adipose tissue replacement of muscle. In addition to describing the data acquisition and analysis procedures, this paper also discusses the potential problems associated with these methods, the analysis and interpretation of the data, MRI safety, and strategies for artifact reduction and protocol optimization.

Loss of electrical anisotropy is an unrecognized feature of dystrophic muscle that may serve as a convenient index of disease status

OBJECTIVE

We sought to understand the alteration in the anisotropic, or direction dependent, character of muscle as measured by electrical impedance myography (EIM) in subjects with Duchenne muscular dystrophy (DMD) and its potential to serve as a biomarker of disease status.

METHODS

Thirty-six boys with DMD and 27 healthy controls were measured with EIM, with electrical current applied both parallel and perpendicular to the major muscle fiber direction. In addition, muscle extracted from 10 mdx and 10 wild-type mice were measured analogously.

RESULTS

Normalized reactance anisotropy, a direction-dependent measure of membrane charge storage capability, was significantly lower in the four muscles of DMD subjects as compared to controls (p<0.01). Normalized reactance anisotropy also decreased with increasing age in DMD subjects (r=-0.36, p=0.031), but not in healthy boys. Analogous changes were observed in mdx mouse gastrocnemius as compared to wild type (p=0.019).

CONCLUSION

These results support that loss of electrical anisotropy is a previously unrecognized feature of dystrophic muscle.

SIGNIFICANCE

Anisotropic alterations may offer novel indices to assist in neuromuscular disease diagnosis and to serve as easy-to-obtain biomarkers in clinical therapeutic trials.

Human skeletal muscle behavior in vivo: Finite element implementation, experiment, and passive mechanical characterization

In this study, useful methods for active human skeletal muscle material parameter determination are provided. First, a straightforward approach to the implementation of a transversely isotropic hyperelastic continuum mechanical material model in an invariant formulation is presented. This procedure is found to be feasible even if the strain energy is formulated in terms of invariants other than those predetermined by the software's requirements. Next, an appropriate experimental setup for the observation of activation-dependent material behavior, corresponding data acquisition, and evaluation is given. Geometry reconstruction based on magnetic resonance imaging of different deformation states is used to generate realistic, subject-specific finite element models of the upper arm. Using the deterministic SIMPLEX optimization strategy, a convenient quasi-static passive-elastic material characterization is pursued; the results of this approach used to characterize the behavior of human biceps in vivo indicate the feasibility of the illustrated methods to identify active material parameters comprising multiple loading modes. A comparison of a contact simulation incorporating the optimized parameters to a reconstructed deformed geometry of an indented upper arm shows the validity of the obtained results regarding deformation scenarios perpendicular to the effective direction of the nonactivated biceps. However, for a valid, activatable, general-purpose material characterization, the material model needs some modifications as well as a multicriteria optimization of the force-displacement data for different loading modes.

Establishing the biomechanical properties of the pelvic soft tissues through an inverse finite element analysis using magnetic resonance imaging

The mechanical characteristics of the female pelvic floor are relevant when explaining pelvic dysfunction. The decreased elasticity of the tissue often causes inability to maintain urethral position, also leading to vaginal and rectal descend when coughing or defecating as a response to an increase in the internal abdominal pressure. These conditions can be associated with changes in the mechanical properties of the supportive structures—namely, the pelvic floor muscles—including impairment. In this work, we used an inverse finite element analysis to calculate the material constants for the passive mechanical behavior of the pelvic floor muscles. The numerical model of the pelvic floor muscles and bones was built from magnetic resonance axial images acquired at rest. Muscle deformation, simulating the Valsalva maneuver with a pressure of 4 KPa, was compared with the muscle displacement obtained through additional dynamic magnetic resonance imaging. The difference in displacement was of 0.15 mm in the antero-posterior direction and 3.69 mm in the supero-inferior direction, equating to a percentage error of 7.0% and 16.9%, respectively. We obtained the shortest difference in the displacements using an iterative process that reached the material constants for the Mooney–Rivlin constitutive model ( c10=11.8 KPa and c20=5.53  E−02 KPa). For each iteration, the orthogonal distance between each node from the group of nodes which defined the puborectal muscle in the numerical model versus dynamic magnetic resonance imaging was computed. With the methodology used in this work, it was possible to obtain in vivo biomechanical properties of the pelvic floor muscles for a specific subject using input information acquired non-invasively.

Multi-parametric MRI characterization of inflammation in murine skeletal muscle

Myopathies often display a common set of complex pathologies that include muscle weakness, inflammation, compromised membrane integrity, fat deposition, and fibrosis. Multi-parametric, quantitative, non-invasive imaging approaches may be able to resolve these individual pathological components. The goal of this study was to use multi-parametric MRI to investigate inflammation as an isolated pathological feature. Proton relaxation, diffusion tensor imaging (DTI), quantitative magnetization transfer (qMT-MRI), and dynamic contrast enhanced (DCE-MRI) parameters were calculated from data acquired in a single imaging session conducted 6-8 hours following the injection of λ-carrageenan, a local inflammatory agent. T2 increased in the inflamed muscle and transitioned to bi-exponential behavior. In diffusion measurements, all three eigenvalues and the apparent diffusion coefficient increased, but λ3 had the largest relative change. Analysis of the qMT data revealed that the T1 of the free pool and the observed T1 both increased in the inflamed tissue, while the ratio of exchanging spins in the solid pool to those in the free water pool (the pool size ratio) significantly decreased. DCE-MRI data also supported observations of an increase in extracellular volume. These findings enriched the understanding of the relation between multiple quantitative MRI parameters and an isolated inflammatory pathology, and may potentially be employed for other single or complex myopathy models.

Quantitative effects of inclusion of fat on muscle diffusion tensor MRI measurements

PURPOSE

To determine the minimum water percentage in a muscle region of interest that would allow diffusion tensor (DT-) MRI data to reflect the diffusion properties of pure muscle accurately.

MATERIALS AND METHODS

Proton density-weighted images with and without fat saturation were obtained at the mid-thigh in four subjects. Co-registered DT-MR images were used to calculate the diffusion tensor's eigenvalues and fractional anisotropy.

RESULTS

The eigenvalues transitioned monotonically as a function of water signal percentage from values near to those expected for pure fat to those for pure muscle. Also, the fractional anisotropy transitioned monotonically from 0.50 (fat) to 0.20 (muscle). For water signal percentages >55%, none of the diffusion indices differed significantly from those for regions of >90% muscle.

CONCLUSION

Accounting for the T1 and T2 values of muscle and fat and the pulse sequence properties, it is concluded that, as a conservative estimate, regions must contain at least 76% muscle tissue to reflect the diffusion properties of pure muscle accurately.