In vivo measurement of membrane permeability and myofiber size in human muscle using time-dependent diffusion tensor imaging and the random permeable barrier model

Influence of Magnetic Field Strength on Intravoxel Incoherent Motion Parameters in Diffusion MRI of the Calf

Background: The purpose of this study was to investigate the dependence of Intravoxel Incoherent Motion (IVIM) parameters measured in the human calf on B0. Methods: Diffusion-weighted image data of eight healthy volunteers were acquired using five b-values (0-600 s/mm2) at rest and after muscle activation at 0.55 and 7 T. The musculus gastrocnemius mediale (GM, activated) was assessed. The perfusion fraction f and diffusion coefficient D were determined using segmented fits. The dependence on field strength was assessed using Student's t-test for paired samples and the Wilcoxon signed-rank test. A biophysical model built on the three non-exchanging compartments of muscle, venous blood, and arterial blood was used to interpret the data using literature relaxation times. Results: The measured perfusion fraction of the GM was significantly lower at 7 T, both for the baseline measurement and after muscle activation. For 0.55 and 7 T, the mean f values were 7.59% and 3.63% at rest, and 14.03% and 6.92% after activation, respectively. The biophysical model estimations for the mean proton-density-weighted perfusion fraction were 3.37% and 6.50% for the non-activated and activated states, respectively. Conclusions: B0 may have a significant effect on the measured IVIM parameters. The blood relaxation times suggest that 7 T IVIM may be arterial-weighted whereas 0.55 T IVIM may exhibit an approximately equal weighting of arterial and venous blood.

Clinical physiology: the crucial role of MRI in evaluation of peripheral artery disease

Peripheral artery disease (PAD) is a common vascular disease that primarily affects the lower limbs and is defined by the constriction or blockage of peripheral arteries and may involve microvascular dysfunction and tissue injury. Patients with diabetes have more prominent disease of microcirculation and develop peripheral neuropathy, autonomic dysfunction, and medial vascular calcification. Early and accurate diagnosis of PAD and disease characterization are essential for personalized management and therapy planning. Magnetic resonance imaging (MRI) provides excellent soft tissue contrast and multiplanar imaging capabilities and is useful as a noninvasive imaging tool in the comprehensive physiological assessment of PAD. This review provides an overview of the current state of the art of MRI in the evaluation and characterization of PAD, including an analysis of the many applicable MR imaging techniques, describing the advantages and disadvantages of each approach. We also present recent developments, future clinical applications, and future MRI directions in assessing PAD. The development of new MR imaging technologies and applications in preclinical models with translation to clinical research holds considerable potential for improving the understanding of the pathophysiology of PAD and clinical applications for improving diagnostic precision, risk stratification, and treatment outcomes in patients with PAD.

Glymphatic Imaging in Pediatrics

The glymphatic system, which facilitates cerebrospinal fluid (CSF) flow through the brain parenchyma, is important for brain development and waste clearance. Advances in imaging techniques, particularly magnetic resonance imaging, have make it possible to evaluate glymphatic structures and functions in vivo. Recently, several studies have focused on the development and alterations of the glymphatic system in pediatric disorders. This review discusses the development of the glymphatic system, advances of imaging techniques and their applications in pediatric disorders. First, the results of the reviewed studies indicate that the development of the glymphatic system is a long-lasting process that continues into adulthood. Second, there is a need for improved glymphatic imaging techniques that are non-invasive and fast to improve suitability for pediatric applications, as some of existing methods use contrast injection and are susceptible to motion artifacts from long scanning times. Several novel techniques are potentially feasible for pediatric patients and may be used in the future. Third, the glymphatic dysfunction is associated with a large number of pediatric disorders, although only a few have recently been investigated. In conclusion, research on the pediatric glymphatic system remains an emerging field. The preliminary applications of glymphatic imaging techniques have provided unique insight into the pathological mechanism of pediatric diseases, but mainly limited in visualization of enlarged perivascular spaces and morphological measurements on CSF volumes. More in-depth studies on glymphatic functions are required to improve our understanding of the mechanisms underlying brain development and pediatric diseases. LEVEL OF EVIDENCE: 5 TECHNICAL EFFICACY: Stage 3.

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.

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.

Monte Carlo simulation of water diffusion through cardiac tissue models

Monte Carlo diffusion simulations are commonly used to establish a reliable ground truth of tissue microstructure, including for the validation of diffusion-weighted MRI. However, selecting simulation parameters is challenging and affects validity and reproducibility. We conducted experiments to investigate critical conditions in Monte Carlo simulations, such as tissue representation complexity, simulated molecules, update duration, and compartment size. Results show significant changes in microstructure characteristics when parameters are altered, emphasizing the importance of careful control for a reliable ground truth.

In vivo Mapping of Cellular Resolution Neuropathology in Brain Ischemia by Diffusion MRI

Non-invasive mapping of cellular pathology can provide critical diagnostic and prognostic information. Recent developments in diffusion MRI have produced new tools for examining tissue microstructure at a level well below the imaging resolution. Here, we report the use of diffusion time ( t )-dependent diffusion kurtosis imaging ( t DKI) to simultaneously assess the morphology and transmembrane permeability of cells and their processes in the context of pathological changes in hypoxic-ischemic brain (HI) injury. Through Monte Carlo simulations and cell culture organoid imaging, we demonstrate feasibility in measuring effective size and permeability changes based on the peak and tail of t DKI curves. In a mouse model of HI, in vivo imaging at 11.7T detects a marked shift of the t DKI peak to longer t in brain edema, suggesting swelling and beading associated with the astrocytic processes and neuronal neurites. Furthermore, we observed a faster decrease of the t DKI tail in injured brain regions, reflecting increased membrane permeability that was associated with upregulated water exchange upon astrocyte activation at acute stage as well as necrosis with disrupted membrane integrity at subacute stage. Such information, unavailable with conventional diffusion MRI at a single t, can predict salvageable tissues. For a proof-of-concept, t DKI at 3T on an ischemic stroke patient suggested increased membrane permeability in the stroke region. This work therefore demonstrates the potential of t DKI for in vivo detection of the pathological changes in microstructural morphology and transmembrane permeability after ischemic injury using a clinically translatable protocol.

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

BACKGROUND

Becker muscular dystrophy (BMD) is an X-linked disorder characterized by slow, progressive muscle damage and muscle weakness. Hallmarks include fibre-size variation and replacement of skeletal muscle with fibrous and adipose tissues, after repeated cycles of regeneration. Muscle histology can detect these features, but the required biopsies are invasive, are difficult to repeat and capture only small muscle volumes. Diffusion-tensor magnetic resonance imaging (DT-MRI) is a potential non-invasive alternative that can calculate muscle fibre diameters when applied with the novel random permeable barrier model (RPBM). In this study, we assessed muscle fibre diameters using DT-MRI in BMD patients and healthy controls and compared these with histology.

METHODS

We included 13 BMD patients and 9 age-matched controls, who underwent water-fat MRI and DT-MRI at multiple diffusion times, allowing RPBM parameter estimation in the lower leg muscles. Tibialis anterior muscle biopsies were taken from the contralateral leg in 6 BMD patients who underwent DT-MRI and from an additional 32 BMD patients and 15 healthy controls. Laminin and Sirius-red stainings were performed to evaluate muscle fibre morphology and fibrosis. Twelve ambulant patients from the MRI cohort underwent the North Star ambulatory assessment, and 6-min walk, rise-from-floor and 10-m run/walk functional tests.

RESULTS

RPBM fibre diameter was significantly larger in BMD patients (P = 0.015): mean (SD) = 68.0 (25.3) μm versus 59.4 (19.2) μm in controls. Inter-muscle differences were also observed (P ≤ 0.002). Both inter- and intra-individual RPBM fibre diameter variability were similar between groups. Laminin staining agreed with the RPBM, showing larger median fibre diameters in patients than in controls: 72.5 (7.9) versus 63.2 (6.9) μm, P = 0.006. However, despite showing similar inter-individual variation, patients showed more intra-individual fibre diameter variability than controls-mean variance (SD) = 34.2 (7.9) versus 21.4 (6.9) μm, P < 0.001-and larger fibrosis areas: median (interquartile range) = 21.7 (5.6)% versus 14.9 (3.4)%, P < 0.001. Despite good overall agreement of RPBM and laminin fibre diameters, they were not associated in patients who underwent DT-MRI and muscle biopsy, perhaps due to lack of colocalization of DT-MRI with biopsy samples.

CONCLUSIONS

DT-MRI RPBM metrics agree with histology and can quantify changes in muscle fibre size that are associated with regeneration without the need for biopsies. They therefore show promise as imaging biomarkers for muscular dystrophies.

To Explore the Changes and Differences of Microstructure of Vocal Fold in Vocal Fold Paralysis and Cricoarytenoid Joint Dislocation by Diffusion Tensor Imaging

OBJECTIVE

The diffusion characteristics of water molecules were measured in the vocal folds of canines exhibiting unilateral vocal fold paralysis and unilateral cricoarytenoid joint dislocation. These characteristics were used in conjunction with a histological examination of the microstructural changes of vocal fold muscle fibers to explore the feasibility of diffusion tensor imaging (DTI) in distinguishing unilateral vocal fold paralysis and unilateral cricoarytenoid joint dislocation as well as evaluating microstructural changes.

METHODS

Ten beagles were randomly divided into three groups: four in the unilateral vocal fold paralysis group, four in the unilateral cricoarytenoid joint dislocation group, and two in the normal group. Unilateral recurrent laryngeal nerve resection was performed in the vocal fold paralysis group. Unilateral cricoarytenoid joint dislocation surgery was performed in the dislocation group. No intervention was performed in the normal group. Four months postintervention, the larynges were excised and put into a magnetic resonance imaging (MRI) system (9.4T BioSpec MRI, Bruker, German) for scanning, followed by an analysis of diffusion parameters among the different groups for statistical significance. After MRI scanning, the vocal folds were cut into sections, stained with hematoxylin and eosin, and scanned digitally. The mean cross-sectional area of muscle fibers, and the mean diameter of muscle fibers in the vocal folds were calculated by target detection and extraction technology. Mean values of each measurement were used to compare the differences among the three groups. Pearson correlation analysis was performed on the DTI parameters and the results from histological section extraction.

RESULTS

The paralysis group had significantly higher Fractional Anisotropy (FA) compared to the dislocation group and normal group (P = 0.004). The paralysis group also had a significantly lower Tensor Trace value compared to the dislocation group and normal group (P = 0.000). The average cross-sectional area of vocal fold muscle fibers in the paralysis group was significantly smaller than the dislocation group and normal group (P = 0.000). Pearson correlation analysis yielded values of, r = -0.785, P = 0.01 between the average cross-sectional area of vocal muscle fibers and FA, and values of r = 0.881, P = 0.00 between Tensor Trace and the average cross-sectional area of vocal muscle.

CONCLUSION

FA and Tensor Trace can be used as effective parameters to reflect the changes of microstructure in vocal fold paralysis and cricoarytenoid joint dislocation. DTI is an objective and quantitative method to effectively evaluate unilateral vocal fold paralysis and unilateral cricoarytenoid joint dislocation, also capable of noninvasively evaluating vocal fold muscle fiber microstructure.

Application of diffusion microstructure imaging in musculoskeletal radiology - translation from head to shoulders

OBJECTIVES

Quantitative MRI techniques, such as diffusion microstructure imaging (DMI), are increasingly applied for advanced tissue characterization. We determined its value in rotator cuff (RC) muscle imaging by studying the association of DMI parameters to isometric strength and fat fraction (FF).

METHODS

Healthy individuals prospectively underwent 3T-MRI of the shoulder using DMI and chemical shift encoding-based water-fat imaging. RC muscles were segmented and quantitative MRI metrics (V-ISO, free fluid; V-intra, compartment inside of muscle fibers; V-extra, compartment outside of muscle fibers, and FF) were extracted. Isometric shoulder strength was quantified using specific clinical tests. Sex-related differences were assessed with Student's t. Association of DMI-metrics, FF, and strength was tested. A factorial two-way ANOVA was performed to compare the main effects of sex and external/internal strength-ratio and their interaction effects on quantitative imaging parameters ratios of infraspinatus/subscapularis.

RESULTS

Among 22 participants (mean age: 26.7 ± 3.1 years, 50% female, mean BMI: 22.6 ± 1.9 kg/m²), FF of the individual RC muscles did not correlate with strength or DMI parameters (all p > 0.05). Subjects with higher V-intra (r = 0.57 to 0.87, p < 0.01) and lower V-ISO (r = -0.6 to -0.88, p < 0.01) had higher internal and external rotation strength. Moreover, V-intra was higher and V-ISO was lower in all RC muscles in males compared to female subjects (all p < 0.01). There was a sex-independent association of external/internal strength-ratio with the ratio of V-extra of infraspinatus/subscapularis (p = 0.02).

CONCLUSIONS

Quantitative DMI parameters may provide incremental information about muscular function and microstructure in young athletes and may serve as a potential biomarker.

KEY POINTS

• Diffusion microstructure imaging was successfully applied to non-invasively assess the microstructure of rotator cuff muscles in healthy volunteers. • Sex-related differences in the microstructural composition of the rotator cuff were observed. • Muscular microstructural metrics correlated with rotator cuff strength and may serve as an imaging biomarker of muscular integrity and function.

Reduction of the cardiac pulsation artifact and improvement of lesion conspicuity in flow‐compensated diffusion images in the liver—A quantitative evaluation of postprocessing algorithms

PURPOSE

To enhance image quality of flow-compensated diffusion-weighted liver MRI data by increasing the lesion conspicuity and reducing the cardiac pulsation artifact using postprocessing algorithms.

METHODS

Diffusion-weighted image data of 40 patients with liver lesions had been acquired at 1.5 T. These data were postprocessed with 5 different algorithms (weighted averaging, p-mean, percentile, outlier exclusion, and exception set). Four image properties of the postprocessed data were evaluated for optimizing the algorithm parameters. These properties were the lesion to tissue contrast-to-noise ratio (CNR), the reduction of the cardiac pulsation artifact, the data consistency, and the vessel darkness. They were combined into a total quality score ( Q total , $$ {Q}_{\mathrm{total}}, $$ set to 1 for the trace-weighted reference image), which was used to rate the image quality objectively.

RESULTS

The weighted averaging algorithm performed best according to the total quality score ( Q total = 1.111 ± 0.067 $$ {Q}_{\mathrm{total}}=1.111\pm 0.067 $$ ). The further ranking was outlier exclusion algorithm ( Q total = 1.086 ± 0.061 $$ {Q}_{\mathrm{total}}=1.086\pm 0.061 $$ ), p-mean algorithm ( Q total = 1.045 ± 0.049 $$ {Q}_{\mathrm{total}}=1.045\pm 0.049 $$ ), percentile algorithm ( Q total = 1.012 ± 0.049 $$ {Q}_{\mathrm{total}}=1.012\pm 0.049 $$ ), and exception set algorithm ( Q total = 0.957 ± 0.027 $$ {Q}_{\mathrm{total}}=0.957\pm 0.027 $$ ). All optimized algorithms except for the exception set algorithm corrected the pulsation artifact and increased the lesion CNR. Changes in Q total $$ {Q}_{\mathrm{total}} $$ were significant for all optimized algorithms except for the percentile algorithm. Liver ADC was significantly reduced (except for the exception set algorithm), particularly in the left lobe.

CONCLUSION

Postprocessing algorithms should be used for flow-compensated liver DWI. The proposed weighted averaging algorithm seems to be suited best to increase the image quality of artifact-corrupted flow-compensated diffusion-weighted liver data.

Time-Dependent Diffusion MRI for Quantitative Microstructural Mapping of Prostate Cancer

Background Recently developed time-dependent diffusion MRI has potential in characterizing cellular tissue microstructures; however, its value in imaging prostate cancer (PCa) remains unknown. Purpose To investigate the feasibility of time-dependent diffusion MRI-based microstructural mapping for noninvasively characterizing cellular properties of PCa and for discriminating between clinically significant PCa and clinically insignificant disease. Materials and Methods Men with a clinical suspicion of PCa were enrolled prospectively between October 2019 and August 2020. Time-dependent diffusion MRI data were acquired with pulsed and oscillating gradient diffusion MRI sequences at an equivalent diffusion time of 7.5-30 msec on a 3.0-T scanner. Time-dependent diffusion MRI-based microstructural parameters, including cell diameter, intracellular volume fraction, cellularity, and diffusivities, were estimated with a two-compartment model. These were compared for different International Society of Urological Pathology grade groups (GGs), and their performance in discriminating clinically significant PCa (GG >1) from clinically insignificant disease (benign and GG 1) was determined with a linear discriminant analysis. The fitted microstructural parameters were validated by means of correlation with histopathologic measurements. Results In the 48 enrolled men, the time-dependent diffusion MRI measurements showed that higher GG was correlated with higher intracellular volume fraction and higher cellularity (intracellular volume fraction = 0.22, 0.36, 0.34, 0.37, and 0.40 in GGs 1-5, respectively; P < .001 at one-way analysis of variance), while lower cell diameter was found at higher GGs (diameter = 23.4, 18.3, 19.2, 17.9, and 18.5 μm in GGs 1-5, respectively; P = .002). Among all measurements derived from time-dependent diffusion MRI, cellularity achieved the highest diagnostic performance, with an accuracy of 92% (44 of 48 participants) and area under the receiver operating characteristic curve of 0.96 (95% CI: 0.87, 0.99) in discriminating clinically significant PCa from clinically insignificant disease. Microstructural mapping was supported by positive correlations between time-dependent diffusion MRI-based and pathologic examination-based intracellular volume fraction (r = 0.83; P < .001). Conclusion Time-dependent diffusion MRI-based microstructural mapping correlates with pathologic findings and demonstrates promise for characterizing prostate cancer. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Chatterjee and Oto in this issue.

Diffusion MRI fiber diameter for muscle denervation assessment

BACKGROUND

To develop and evaluate a diffusion MRI-based apparent muscle fiber diameter (AFD) method in patients with muscle denervation. It was hypothesized that AFD differences between denervated, non-denervated and control muscles would be greater than those from standard diffusion metrics.

METHODS

A spin-echo diffusion acquisition with multi-b-valued diffusion sampling was used. An orientation-invariant dictionary approach utilized a cylinder-based forward model and multi-compartment model for obtaining restricted and free fractions. Simulations were performed to determine precision, bias, and optimize dictionary parameters. In all, 18 exams of patients with muscle denervation and 8 exams of healthy subjects were performed at 3T. Six regions of interests (ROIs) within separate shoulder muscles were selected, yielding three groups consisting 47 control (healthy), 36 non-denervated (patients), and 68 denervated (patients) muscle ROIs. Two-sample t-tests (α=0.05) between groups were performed with Holm-Bonferroni correction. T2- and fat fraction (FF)-mapping were acquired for comparison.

RESULTS

Mean AFD was 89.7±13.6 µm in control, 71.6±15.3 µm in non-denervated, and 60.7±15.9 µm in denervated muscles and were significantly different (P<0.001) in paired comparisons and in 10/12 individual muscle region comparisons. Correlation between AFD and FF (-0.331, P<0.001) was low, but correlation between FA and FF was negligible (0.197, P=0.016). Correlation was low between AFD and T2 (-0.395, P<0.001) and between FA and T2 (0.359, P<0.001).

CONCLUSIONS

Diffusion MRI-based AFD complements T2- and FF-mapping techniques to non-invasively assess muscle denervation.

Combined diffusion-relaxometry microstructure imaging: Current status and future prospects

Microstructure imaging seeks to noninvasively measure and map microscopic tissue features by pairing mathematical modeling with tailored MRI protocols. This article reviews an emerging paradigm that has the potential to provide a more detailed assessment of tissue microstructure-combined diffusion-relaxometry imaging. Combined diffusion-relaxometry acquisitions vary multiple MR contrast encodings-such as b-value, gradient direction, inversion time, and echo time-in a multidimensional acquisition space. When paired with suitable analysis techniques, this enables quantification of correlations and coupling between multiple MR parameters-such as diffusivity, T 1 , T 2 , and T 2 ∗ . This opens the possibility of disentangling multiple tissue compartments (within voxels) that are indistinguishable with single-contrast scans, enabling a new generation of microstructural maps with improved biological sensitivity and specificity.

High Inter-Rater Reliability of Manual Segmentation and Volume-Based Tractography in Healthy and Dystrophic Human Calf Muscle

BACKGROUND

Muscle diffusion tensor imaging (mDTI) is a promising surrogate biomarker in the evaluation of muscular injuries and neuromuscular diseases. Since mDTI metrics are known to vary between different muscles, separation of different muscles is essential to achieve muscle-specific diffusion parameters. The commonly used technique to assess DTI metrics is parameter maps based on manual segmentation (MSB). Other techniques comprise tract-based approaches, which can be performed in a previously defined volume. This so-called volume-based tractography (VBT) may offer a more robust assessment of diffusion metrics and additional information about muscle architecture through tract properties. The purpose of this study was to assess DTI metrics of human calf muscles calculated with two segmentation techniques-MSB and VBT-regarding their inter-rater reliability in healthy and dystrophic calf muscles.

METHODS

20 healthy controls and 18 individuals with different neuromuscular diseases underwent an MRI examination in a 3T scanner using a 16-channel Torso XL coil. DTI metrics were assessed in seven calf muscles using MSB and VBT. Coefficients of variation (CV) were calculated for both techniques. MSB and VBT were performed by two independent raters to assess inter-rater reliability by ICC analysis and Bland-Altman plots. Next to analysis of DTI metrics, the same assessments were also performed for tract properties extracted with VBT.

RESULTS

For both techniques, low CV were found for healthy controls (≤13%) and neuromuscular diseases (≤17%). Significant differences between methods were found for all diffusion metrics except for λ₁. High inter-rater reliability was found for both MSB and VBT (ICC ≥ 0.972). Assessment of tract properties revealed high inter-rater reliability (ICC ≥ 0.974).

CONCLUSIONS

Both segmentation techniques can be used in the evaluation of DTI metrics in healthy controls and different NMD with low rater dependency and high precision but differ significantly from each other. Our findings underline that the same segmentation protocol must be used to ensure comparability of mDTI data.

Intravoxel Incoherent Motion Magnetic Resonance Imaging in Skeletal Muscle: Review and Future Directions

Throughout the body, muscle structure and function can be interrogated using a variety of noninvasive magnetic resonance imaging (MRI) methods. Recently, intravoxel incoherent motion (IVIM) MRI has gained momentum as a method to evaluate components of blood flow and tissue diffusion simultaneously. Much of the prior research has focused on highly vascularized organs, including the brain, kidney, and liver. Unique aspects of skeletal muscle, including the relatively low perfusion at rest and large dynamic range of perfusion between resting and maximal hyperemic states, may influence the acquisition, postprocessing, and interpretation of IVIM data. Here, we introduce several of those unique features of skeletal muscle; review existing studies of IVIM in skeletal muscle at rest, in response to exercise, and in disease states; and consider possible confounds that should be addressed for muscle-specific evaluations. Most studies used segmented nonlinear least squares fitting with a b-value threshold of 200 sec/mm² to obtain IVIM parameters of perfusion fraction (f), pseudo-diffusion coefficient (D*), and diffusion coefficient (D). In healthy individuals, across all muscles, the average ± standard deviation of D was 1.46 ± 0.30 × 10^-3  mm² /sec, D* was 29.7 ± 38.1 × 10^-3  mm² /sec, and f was 11.1 ± 6.7%. Comparisons of reported IVIM parameters in muscles of the back, thigh, and leg of healthy individuals showed no significant difference between anatomic locations. Throughout the body, exercise elicited a positive change of all IVIM parameters. Future directions including advanced postprocessing models and potential sequence modifications are discussed. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 2.

Clinically feasible diffusion MRI in muscle: Time dependence and initial findings in Duchenne muscular dystrophy

PURPOSE

To characterize the diffusion time-dependence in muscle in healthy adult volunteers, boys with Duchenne's muscular dystrophy (DMD), and age-matched controls in a clinically feasible acquisition time for pediatric applications.

METHODS

Diffusion data were acquired using a pulsed gradient stimulated echo diffusion preparation at 5 different diffusion times (70, 130, 190, 250, and 330 ms), at 4 different b-values (0, 200, 400, 600, and 800 s/mm² ) and 6 directions (orthogonal x, y, and z and diagonal xy, xz, and yz) and processed to obtain standard diffusion indices (mean diffusivity [MD] and fractional anisotropy [FA]) at each diffusion time.

RESULTS

Time-dependent diffusion was seen in muscle in healthy adult volunteers, boys with DMD, and age-matched controls. Boys with DMD showed reduced MD and increased FA values in comparison to age matched controls across a range of diffusion times. A diffusion time of Δ = 190 ms had the largest effect size.

CONCLUSIONS

These results could be used to optimize diffusion imaging in this disease further and imply that these diffusion indices may become an important biomarker in monitoring progression in DMD in the future.

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

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

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.

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

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

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 present and the future of microstructure MRI: From a paradigm shift to normal science

The aspiration of imaging tissue microstructure with MRI is to uncover micrometer-scale tissue features within millimeter-scale imaging voxels, in vivo. This kind of super-resolution has fueled a paradigm shift within the biomedical imaging community. However, what feels like an ongoing revolution in MRI, has been conceptually experienced in physics decades ago; from this point of view, our current developments can be seen as Thomas Kuhn's "normal science" stage of progress. While the concept of model-based quantification below the nominal imaging resolution is not new, its possibilities in neuroscience and neuroradiology are only beginning to be widely appreciated. This disconnect calls for communicating the progress of tissue microstructure MR imaging to its potential users. Here, a number of recent research developments are outlined in terms of the overarching concept of coarse-graining the tissue structure over an increasing diffusion length. A variety of diffusion models and phenomena are summarized on the phase diagram of diffusion MRI, with the unresolved problems and future directions corresponding to its unexplored domains.

Lattice Boltzmann method for simulation of diffusion magnetic resonance imaging physics in multiphase tissue models

We report an implementation of the lattice Boltzmann method (LBM) to integrate the Bloch-Torrey equation, which describes the evolution of the transverse magnetization vector and the fate of the signal of diffusion magnetic resonance imaging (dMRI). Motivated by the need to interpret dMRI experiments in biological tissues, and to offset the small time-step limitation of classical LBM, a hybrid LBM scheme is introduced and implemented to solve the Bloch-Torrey equation. A membrane boundary condition is presented which is able to accurately represent the effects of thin curvilinear membranes typically found in biological tissues. As implemented, the hybrid LBM scheme accommodates piece-wise uniform transport, dMRI parameters, periodic and mirroring outer boundary conditions, and finite membrane permeabilities on non-boundary-conforming inner boundaries. By comparing with analytical solutions of limiting cases, we demonstrate that the hybrid LBM scheme is more accurate than the classical LBM scheme. The proposed explicit LBM scheme maintains second-order spatial accuracy, stability, and first-order temporal accuracy for a wide range of parameters. The parallel implementation of the hybrid LBM code in a multi-CPU computer system, as well as on GPUs, is straightforward and efficient. Along with offering certain advantages over finite element or Monte Carlo schemes, the proposed hybrid LBM constitutes a flexible scheme that can by easily adapted to model more complex interfacial conditions and physics in heterogeneous multiphase tissue models and to accommodate sophisticated dMRI sequences.

Novel Muscle Imaging in Inflammatory Rheumatic Diseases-A Focus on Ultrasound Shear Wave Elastography and Quantitative MRI

In recent years, imaging has played an increasing role in the clinical management of patients with rheumatic diseases with respect to aiding diagnosis, guiding therapy and monitoring disease progression. These roles have been underpinned by research which has enhanced our understanding of disease pathogenesis and pathophysiology of rheumatology conditions, in addition to their key role in outcome measurement in clinical trials. However, compared to joints, imaging research of muscles is less established, despite the fact that muscle symptoms are very common and debilitating in many rheumatic diseases. Recently, it has been shown that even though patients with rheumatoid arthritis may achieve clinical remission, defined by asymptomatic joints, many remain affected by lingering constitutional systemic symptoms like fatigue, tiredness, weakness and myalgia, which may be attributed to changes in the muscles. Recent improvements in imaging technology, coupled with an increasing clinical interest, has started to ignite new interest in the area. This perspective discusses the rationale for using imaging, particularly ultrasound and MRI, for investigating muscle pathology involved in common inflammatory rheumatic diseases. The muscles associated with rheumatic diseases can be affected in many ways, including myositis-an inflammatory muscle condition, and myopathy secondary to medications, such as glucocorticoids. In addition to non-invasive visual assessment of muscles in these conditions, novel imaging techniques like shear wave elastography and quantitative MRI can provide further useful information regarding the physiological and biomechanical status of the muscle.

Real-time measurement of diffusion exchange rate in biological tissue

Diffusion exchange spectroscopy (DEXSY) provides a means to isolate the signal attenuation associated with exchange from other sources of signal loss. With the total diffusion weighting b₁+b₂=b_s held constant, DEXSY signals acquired with b₁=0 or b₂=0 have no exchange weighting, while a DEXSY signal acquired with b₁=b₂ has maximal exchange weighting. The exchange rate can be estimated by fitting a diffusion exchange model to signals acquired with variable mixing times. Conventionally, acquired signals are normalized by a signal with b₁=0 and b₂=0 to remove the decay due to spin-lattice relaxation. Instead, division by a signal with equal b_s but b₁=0 or b₂=0 reduces spin-lattice relaxation weighting of the apparent exchange rate (AXR). Furthermore, apparent diffusion-weighted R₁ relaxation rates can be estimated from non-exchange-weighted DEXSY signals. Estimated R₁ values are utilized to remove signal decay due to spin-lattice relaxation from exchange-weighted signals, permitting a more precise estimate of AXR with less data. Data reduction methods are proposed and tested with regards to statistical accuracy and precision of AXR estimates on simulated and experimental data. Simulations show that the methods are capable of accurately measuring the ground-truth exchange rate. The methods remain accurate even when the assumption that DEXSY signals attenuate with b is violated, as occurs for restricted diffusion. Experimental data was collected from fixed neonatal mouse spinal cord samples at 25 and 7°C using the strong static magnetic field gradient produced by a single-sided permanent magnet (i.e., an NMR MOUSE). The most rapid method for exchange measurements requires only five data points (an 80 s experiment as implemented) and achieves a similar level of accuracy and precision to the baseline method using 44 data points. This represents a significant improvement in acquisition speed, overcoming a barrier which has limited the use of DEXSY on living specimen.

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.

Time-dependent diffusion MRI as a probe of microstructural changes in a mouse model of Duchenne muscular dystrophy

Dystrophic muscles show a high variability of fibre sizes and altered sarcolemmal integrity, which are typically assessed by histology. Time-dependent diffusion MRI is sensitive to tissue microstructure and its investigation through age-related changes in dystrophic and healthy muscles may help the understanding of the onset and progression of Duchenne muscular dystrophy (DMD). We investigated the capability of time-dependent diffusion MRI to quantify age and disease-related changes in hind-limb muscle microstructure between dystrophic (mdx) and wild-type (WT) mice of three age groups (7.5, 22 and 44 weeks). Diffusion time-dependent apparent diffusion coefficients (ADCs) of the gastrocnemius and tibialis anterior muscles were determined versus age and diffusion-gradient orientation at six diffusion times (Δ; range: 25-350 ms). Mean muscle ADCs were compared between groups and ages, and correlated with T₂ , using Student's t test, one-way analysis of variance and Pearson correlation, respectively. Muscle fibre sizes and sarcolemmal integrity were evaluated by histology and compared with diffusion measurements. Hind-limb muscle ADC showed characteristic restricted diffusion behaviour in both mdx and WT animals with decreasing ADC values at longer Δ. Significant differences in ADC were observed at long Δ values (≥ 250 ms; p < 0.05, comparison between groups; p < 0.01, comparison between ages) with ADC increased by 5-15% in dystrophic muscles, indicative of reduced diffusion restriction. No significant correlation was found between T₂ and ADC. Additionally, muscle fibre size distributions showed higher variability and lower mean fibre size in mdx than WT animals (p < 0.001). The extensive Evans Blue Dye uptake shown in dystrophic muscles revealed substantial sarcolemmal damage, suggesting diffusion measurements as more consistent with altered permeability rather than changes in muscle fibre sizes. This study shows the potential of diffusion MRI to non-invasively discriminate between dystrophic and healthy muscles with enhanced sensitivity when using long Δ.

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.

Stimulated echo double diffusion encoded imaging of closed pores: Influence and removal of unbalanced terms

Nuclear magnetic resonance (NMR) diffusion pore imaging has been proposed to study the shape of arbitrary closed pores filled with an NMR-detectable medium by use of nonclassical diffusion encoding schemes. Potential applications can be found in biomedical imaging and porous media research. When studying non-point-symmetric pores, NMR signals with nonvanishing imaginary parts arise containing the pore shape information, which is lost for classical diffusion encoding schemes. Key limitations are the required high magnetic field gradient amplitudes and T2 relaxation while approaching the diffusion long-time limit. To benefit from the slower T1 decay, we demonstrate the feasibility of diffusion pore imaging with stimulated echoes using Monte Carlo simulations and experiments with hyperpolarized xenon-129 gas in well-defined geometries and show that the necessary complex-valued signals can be acquired. Analytical derivation of the stimulated echo double diffusion encoded signal was performed to investigate the effect of the additionally arising undesired terms on the complex phase information. These terms correspond to signals arising for spin-echo sequences with unbalanced gradients. For most possible applications, the unbalanced terms can be neglected. If non-negligible, selection of the appropriate signal component using a phase cycling scheme was demonstrated experimentally. Using stimulated echoes may be a step towards application of diffusion pore imaging to larger pores with gradient amplitudes available today in preclinical systems.

Quantifying brain microstructure with diffusion MRI: Theory and parameter estimation

We review, systematize and discuss models of diffusion in neuronal tissue, by putting them into an overarching physical context of coarse-graining over an increasing diffusion length scale. From this perspective, we view research on quantifying brain microstructure as occurring along three major avenues. The first avenue focusses on transient, or time-dependent, effects in diffusion. These effects signify the gradual coarse-graining of tissue structure, which occurs qualitatively differently in different brain tissue compartments. We show that transient effects contain information about the relevant length scales for neuronal tissue, such as the packing correlation length for neuronal fibers, as well as the degree of structural disorder along the neurites. The second avenue corresponds to the long-time limit, when the observed signal can be approximated as a sum of multiple nonexchanging anisotropic Gaussian components. Here, the challenge lies in parameter estimation and in resolving its hidden degeneracies. The third avenue employs multiple diffusion encoding techniques, able to access information not contained in the conventional diffusion propagator. We conclude with our outlook on future directions that could open exciting possibilities for designing quantitative markers of tissue physiology and pathology, based on methods of studying mesoscopic transport in disordered systems.

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.

Comparison of diffusion-weighted MRI and anti-Stokes Raman scattering (CARS) measurements of the inter-compartmental exchange-time of water in expression-controlled aquaporin-4 cells

We performed multi-b and multi-diffusion-time diffusion-weighted magnetic resonance imaging on aquaporin-4-expressing (AQ) and -non-expressing (noAQ) cells, and demonstrated a clear difference between the signals from the two cell types. The data were interpreted using a two-compartment (intra and extracellular spaces) model including inter-compartmental exchange. It was also assumed that restricted diffusion of water molecules inside the cells leads to the intracellular diffusion coefficient being inversely proportional to the diffusion-time. Estimates of the water-exchange-times obtained with this model are compared to those measured using an independent optical imaging technique (coherent anti-Stokes Raman scattering imaging, CARS). For both techniques it was found that the exchange-time estimated for the noAQ cells was significantly longer than that for the AQ cells.

Quantifying myofiber integrity using diffusion MRI and random permeable barrier modeling in skeletal muscle growth and Duchenne muscular dystrophy model in mice

PURPOSE

To measure the microstructural changes during skeletal muscle growth and progressive pathologies using the random permeable model with diffusion MRI, and compare findings to conventional imaging modalities such as three-point Dixon and T2 imaging.

METHODS

In vivo and ex vivo DTI experiments with multiple diffusion times (20-700 ms) were completed on wild-type (n = 22) and muscle-dystrophic mdx mice (n = 8) at various developmental time points. The DTI data were analyzed with the random permeable model framework that provides estimates of the unrestricted diffusion coefficient (D0 ), membrane surface-to-volume ratio (S/V), and membrane permeability (κ). In addition, the MRI experiments included conventional measures, such as tissue fat fractions and T2 relaxation.

RESULTS

During normal muscle growth between week 4 and week 13, the in vivo S/V, fractional anisotropy, and fat fraction correlated positively with age (ρ = 0.638, 0.664, and 0.686, respectively), whereas T2 correlated negatively (ρ = -0.847). In mdx mice, all DTI random permeable model parameters and fat fraction had significant positive correlation with age, whereas fractional anisotropy and T2 did not have significant correlation with age. Histological measurements of the perimeter-to-area ratio served as a proxy for the model-derived S/V in the cylindrical myofiber geometry, and had a significant correlation with the ex vivo S/V (r = 0.71) as well as the in vivo S/V (r = 0.56).

CONCLUSION

The present study demonstrates that DTI at multiple diffusion times with the random permeable model analysis allows for noninvasively quantifying muscle fiber microstructural changes during both normal muscle growth and disease progression. Future studies can apply our technique to evaluate current and potential treatments to muscle myopathies.

Characterization of prostate microstructure using water diffusion and NMR relaxation

For many pathologies, early structural tissue changes occur at the cellular level, on the scale of micrometers or tens of micrometers. Magnetic resonance imaging (MRI) is a powerful non-invasive imaging tool used for medical diagnosis, but its clinical hardware is incapable of reaching the cellular length scale directly. In spite of this limitation, microscopic tissue changes in pathology can potentially be captured indirectly, from macroscopic imaging characteristics, by studying water diffusion. Here we focus on water diffusion and NMR relaxation in the human prostate, a highly heterogeneous organ at the cellular level. We present a physical picture of water diffusion and NMR relaxation in the prostate tissue, that is comprised of a densely-packed cellular compartment (composed of stroma and epithelium), and a luminal compartment with almost unrestricted water diffusion. Transverse NMR relaxation is used to identify fast and slow T ₂ components, corresponding to these tissue compartments, and to disentangle the luminal and cellular compartment contributions to the temporal evolution of the overall water diffusion coefficient. Diffusion in the luminal compartment falls into the short-time surface-to-volume (S/V) limit, indicating that only a small fraction of water molecules has time to encounter the luminal walls of healthy tissue; from the S/V ratio, the average lumen diameter averaged over three young healthy subjects is measured to be 217.7±188.7 μm. Conversely, the diffusion in the cellular compartment is highly restricted and anisotropic, consistent with the fibrous character of the stromal tissue. Diffusion transverse to these fibers is well described by the random permeable barrier model (RPBM), as confirmed by the dynamical exponent ϑ = 1/2 for approaching the long-time limit of diffusion, and the corresponding structural exponent p = -1 in histology. The RPBM-derived fiber diameter and membrane permeability were 19.8±8.1 μm and 0.044±0.045 μm/ms, respectively, in agreement with known values from tissue histology and membrane biophysics. Lastly, we revisited 38 prostate cancer cases from a recently published study, and found the same dynamical exponent ϑ = 1/2 of diffusion in tumors and benign regions. Our results suggest that a multi-parametric MRI acquisition combined with biophysical modeling may be a powerful non-invasive complement to prostate cancer grading, potentially foregoing biopsies.

[Diffusion Tensor Imaging of Rotator Cuff: Influence of Arm Position]

Diffusion tensor imaging (DTI) of skeletal muscles has been reported as capable to characterize physiological properties, tissue microstructure and architectural organization. However, the DTI indices may vary with the contractile state of the muscles, and in the rotator cuff muscles, a change in forearm position can result in variation of the DTI indices. The purpose of this study was to examine the influence of forearm position on the major DTI indices of the rotator cuff muscles. The DTI of right rotator cuff was acquired under the neutral position and external and internal rotation of the forearm in nine healthy volunteers. Fractional anisotropy (FA) and mean diffusivity (MD) of each muscle were calculated and compared among the three forearm positions. FA and MD were significantly different between external and internal rotation in infraspinatus, teres minor and subscapularis (p<0.05). We considered that this difference was due to the change in cross-sectional area of muscle fibers based on their contractile state. That is, when the muscle is contracted, its cross-sectional area is increased and the muscle fiber density in the short axis direction becomes less. This causes a change in FA and MD due to increase in λ₂ and λ₃ through increased diffusion of intercellular water in the short axis direction. In conclusion, the DTI indices of the rotator cuff muscles are affected by the forearm position.

The effects of ageing on mouse muscle microstructure: a comparative study of time-dependent diffusion MRI and histological assessment

The investigation of age-related changes in muscle microstructure between developmental and healthy adult mice may help us to understand the clinical features of early-onset muscle diseases, such as Duchenne muscular dystrophy. We investigated the evolution of mouse hind-limb muscle microstructure using diffusion imaging of in vivo and in vitro samples from both actively growing and mature mice. Mean apparent diffusion coefficients (ADCs) of the gastrocnemius and tibialis anterior muscles were determined as a function of diffusion time (Δ), age (7.5, 22 and 44 weeks) and diffusion gradient direction, applied parallel or transverse to the principal axis of the muscle fibres. We investigated a wide range of diffusion times with the goal of probing a range of diffusion lengths characteristic of muscle microstructure. We compared the diffusion time-dependent ADC of hind-limb muscles with histology. ADC was found to vary as a function of diffusion time in muscles at all stages of maturation. Muscle water diffusivity was higher in younger (7.5 weeks) than in adult (22 and 44 weeks) mice, whereas no differences were observed between the older ages. In vitro data showed the same diffusivity pattern as in vivo data. The highlighted differences in diffusion properties between young and mature muscles suggested differences in underlying muscle microstructure, which were confirmed by histological assessment. In particular, although diffusion was more restricted in older muscle, muscle fibre size increased significantly from young to adult age. The extracellular space decreased with age by only ~1%. This suggests that the observed diffusivity differences between young and adult muscles may be caused by increased membrane permeability in younger muscle associated with properties of the sarcolemma.

Advanced MRI Techniques for Muscle Imaging

Advanced magnetic resonance imaging (MRI) techniques can evaluate a wide array of muscle pathologies including acute or chronic muscle injury, musculotendinous response to injury, intramuscular collections and soft tissue masses, and others. In recent years, MRI has played a more important role in muscle disease diagnosis and monitoring. MRI provides excellent spatial and contrast resolution and helps direct optimal sites for muscle biopsy. Whole-body MRI now helps identify signature patterns of muscular involvement in large anatomical regions with relative ease. Quantitative MRI has advanced the evaluation and disease tracking of muscle atrophy and fatty infiltration in entities such as muscular dystrophies. Multivoxel magnetic resonance spectroscopy (MRS) now allows a more thorough, complete evaluation of a muscle of interest without the inherent sampling bias of single-voxel MRS or biopsy. Diffusion MRI allows quantification of muscle inflammation and capillary perfusion as well as muscle fiber tracking.

Fundamentals of diffusion MRI physics

Diffusion MRI is commonly considered the "engine" for probing the cellular structure of living biological tissues. The difficulty of this task is threefold. First, in structurally heterogeneous media, diffusion is related to structure in quite a complicated way. The challenge of finding diffusion metrics for a given structure is equivalent to other problems in physics that have been known for over a century. Second, in most cases the MRI signal is related to diffusion in an indirect way dependent on the measurement technique used. Third, finding the cellular structure given the MRI signal is an ill-posed inverse problem. This paper reviews well-established knowledge that forms the basis for responding to the first two challenges. The inverse problem is briefly discussed and the reader is warned about a number of pitfalls on the way.