Search

PURPOSE: This study implements a compressed sensing (CS) 3-directional velocity encoded phase contrast (VE-PC) imaging for studying skeletal muscle kinematics within 40 s. METHODS: Independent variable density random sampling in the phase encoding direction for each temporal frame was implemented for various combinations of CS-factors and views per segment. CS reconstruction was performed for the combined multicoil, temporal datasets using temporal Fourier transform followed by temporal principal component analysis sparsifying transformations. The method was tested on a flow phantom and in vivo, on velocity and strain rate of the medial gastrocnemius muscle of 11 subjects performing isometric contractions. RESULTS: For the flow phantom, velocity from 8 undersampled sequences matched very well with the flowmeter values over a range of velocities spanning in vivo muscle velocities. Bland-Altman plots of the peak strain rate eigenvalues comparing 7 undersampled sequences was in good agreement with the reference (full k-space) scan. CS-factor of 4 combined with views per segment of 4 (scan times reduced by 4) yielded images with no visual artifacts allowing and yielded velocities and strain rate maps in the lower leg muscle in 40 s. CONCLUSION: This study shows that a reduction in scan time of velocity encoded phase contrast imaging up to a factor of 4 is possible using the proposed CS reconstruction.

3D strain or strain rate tensor mapping comprehensively captures regional muscle deformation. While compressive strain along the muscle fiber is a potential measure of the force generated, radial strains in the fiber cross-section may provide information on the material properties of the extracellular matrix. Additionally, shear strain may potentially inform on the shearing of the extracellular matrix; the latter has been hypothesized as the mechanism of lateral transmission of force. Here, we implement a novel fast MR method for velocity mapping to acquire multi-slice images at different % maximum voluntary contraction (MVC) for 3D strain mapping to explore deformation in the plantar-flexors under isometric contraction in a cohort of young and senior subjects. 3D strain rate and strain tensors were computed and eigenvalues and two invariants (maximum shear and volumetric strain) were extracted. Strain and strain rate indices (contractile and in-plane strain/strain rate, shear strain/strain rate) changed significantly with %MVC (30 and 60% MVC) and contractile and shear strain with age in the medial gastrocnemius. In the soleus, significant differences with age in contractile and shear strain were seen. Univariate regression revealed weak but significant correlation ofin-planeandshear strainandshear strain rateindices to %MVC and correlation ofcontractileandshear strainindices to force. The ability to map strain tensor components provides unique insights into muscle physiology: with contractile strain providing an index of the force generated by the muscle fibers while the shear strain could potentially be a marker of lateral transmission of force.

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.

Age- and disuse- related loss of muscle force is disproportionately larger than the loss of muscle mass. Earlier studies reported that comparing concentric and eccentric contractions, there is a significant age-related decrease in force only in concentric contractions. Magnetic Resonance Imaging enables mapping of muscle deformation and has been used to study isometric but not eccentric contractions. We report MRI based strain rate mapping of the medial gastrocnemius in subjects pre- and post-unloading induced by Unilateral Limb Suspension. In contrast to isometric contraction, no difference in strain rate indices were observed post-unloading, in conformance with preserved force during eccentric contractions.

BACKGROUND Strain rate (SR) is a measure of the rate of regional deformation that can be computed by analyzing velocity-encoded phase-contrast 2D images. Recent studies have explored the changes in normal components of the strain tensor in aging muscle, while shear strain may also provide valuable information. PURPOSE To compute the shear SR from velocity-encoded MRI of the lower leg and to study the correlation of SR parameters measured in the medial gastrocnemius (MG) to muscle force in a cohort of young and senior subjects. STUDY TYPE Prospective cohort study. SUBJECTS Six young (26.1 ± 2.3 years) and six senior (76.7 ± 8.3 years) healthy females; two other subjects were scanned on three separate occasions for repeatability studies. FIELD STRENGTH/SEQUENCE 1.5T using a single oblique sagittal slice with velocity-encoding in three directions (velocity-encoded phase contrast gradient echo sequence). ASSESSMENT Age-related and regional differences in the SR eigenvalues (SRfiber , SRin-plane ), normal SRs (SRff , SRcc ), and shear SRs (SRfc , SRfc_max ) were statistically analyzed. STATISTICAL TESTS Difference between young and senior cohorts were assessed using two-way analysis of variance (ANOVAs). The coefficient of variation and repeatability coefficient were calculated from repeat studies. Univariate and stepwise multivariable linear regression was performed to identify predictors of force. RESULTS During isometric plantarflexion contraction, SRs in the principal basis (SRfiber , SRin-plane ) and maximum shear SR (SRfc_max ) was significantly lower in the senior cohort (P < 0.05). On multiple variable regression, maximum shear SR (SRfc_max ) and normal SR in the fiber cross-section (SRcc ) were significantly associated with force (R = 0.681, F = 14.034, P < 0.001). DATA CONCLUSION This study establishes that computation of shear strain is feasible and is a significant predictor of force variability with age. LEVEL OF EVIDENCE 2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2018;47:1351-1357.

BACKGROUND Urinary incontinence is a major clinical problem arising primarily from age-related degenerative changes to the sphincter muscles. However, the precise anatomy of the normal male sphincter muscles has yet to be established. Diffusion tensor imaging (DTI) may offer a unique insight into muscle microstructure and fiber architecture. PURPOSE To explore the anatomy of the urethral sphincter muscles pertinent to urinary continence function using DT-MRI. STUDY TYPE Prospective cohort study. SUBJECTS Eleven normal male subjects (mean age: 25.4 years); two subjects were scanned in three separate sessions to assess reproducibility. FIELD STRENGTH/SEQUENCE 3T; using a diffusion-weighted spin echo planar sequence. ASSESSMENT DT parameters including fractional anisotropy (FA), primary (λ1 ), secondary (λ2 ), and tertiary (λ3 ) eigenvalues, Apparent diffusion coefficient and radial diffusivity were analyzed statistically, while tracked muscle fibers were assessed visually. STATISTICAL TESTS Regional differences (sphincters and longitudinal muscle of the urethra) in the DTI indices were assessed by one-way analysis of variance. A Tukey post-hoc test was used to identify significant differences between muscle regions. RESULTS Two sphincter muscles, one proximal near the base of the bladder, corresponding to the lisso-sphincter, and the other distal to the end of the prostate corresponding to the rhabdo-sphincter, surrounding a central urethral muscle fiber bundle, were clearly identified. FA was higher and λ3 lower in the proximal sphincter muscle compared to the central urethral muscle and the distal sphincter (P

Purpose This study explores changes in strain rate (SR) (rate of regional deformation) parameters extracted from velocity-encoded MRI and their relationship to muscle force loss following 4-week unilateral lower limb suspension in healthy humans. Methods Two-dimensional SR maps were derived from three directional velocity-encoded MR phase-contrast images of the medial gastrocnemius in seven subjects. Atrophy-related and regional differences in the SR eigenvalues, angle between the SR and muscle fiber (SR-fiber angle), and strain rates in the fiber basis were statistically analyzed using analysis of variance and linear regression. Results During isometric contraction, SR in the fiber cross section (SRin-plane) was significantly lower, and the SR-fiber angle was significantly higher postsuspension (P

Purpose To develop and evaluate an automated algorithm to segment intramuscular adipose (IMAT) and connective (IMCT) tissue from musculoskeletal MRI images acquired with a dual echo Ultrashort TE (UTE) sequence. Theory and methods The dual echo images and calculated structure tensor images are the inputs to the multichannel fuzzy cluster mean (MCFCM) algorithm. Modifications to the basic multichannel fuzzy cluster mean include an adaptive spatial term and bias shading correction. The algorithm was tested on digital phantoms simulating IMAT/IMCT tissue under varying conditions of image noise and bias and on ten subjects with varying amounts of IMAT/IMCT. Results The MCFCM including the adaptive spatial term and bias shading correction performed better than the original MCFCM and adaptive spatial MCFCM algorithms. IMAT/IMCT was segmented from the unsmoothed simulated phantom data with a mean Dice coefficient of 0.933 ±0.001 when contrast-to-noise (CNR) was 140 and bias was varied between 30% and 65%. The algorithm yielded accurate in vivo segmentations of IMAT/IMCT with a mean Dice coefficient of 0.977 ±0.066. Conclusion The proposed algorithm is completely automated and yielded accurate segmentation of intramuscular adipose and connective tissue in the digital phantom and in human calf data. Magn Reson Med 77:870-883, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

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: λ1 , λ2 , λ3 , 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 (λ1 : P = 0.025, λ2 : P = 0.035, λ3 : 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.