The effect of noise and lipid signals on determination of Gaussian and non-Gaussian diffusion parameters in skeletal muscle

Feasibility of Multiparametric Perfusion Assessment in Diabetic Foot Ulcer Using Intravoxel Incoherent Motion and Blood Oxygenation-Level Dependent MRI

BACKGROUND

Patients with type-2 diabetes (T2DM) are at increased risk of developing diabetic foot ulcers (DFU) and experiencing impaired wound healing related to underlying microvascular disease.

PURPOSE

To evaluate the sensitivity of intra-voxel incoherent motion (IVIM) and blood oxygen level dependent (BOLD) MRI to microvascular changes in patients with DFUs.

STUDY TYPE

Case-control.

POPULATION

20 volunteers who were age and body mass index matched, including T2DM patients with DFUs (N = 10, mean age = 57.5 years), T2DM patients with controlled glycemia and without DFUs (DC, N = 5, mean age = 57.4 years) and healthy controls (HC, N = 5, mean age = 52.8 years).

FIELD STRENGTH/SEQUENCE

3T/multi-b-value IVIM and dynamic BOLD.

ASSESSMENT

Resting IVIM parameters were obtained using a multi-b-value diffusion-weighted imaging sequence and two IVIM models were fit to obtain diffusion coefficient (D), pseudo-diffusion coefficient (D*), perfusion fraction (f) and microvascular volume fraction (MVF) parameters. Microvascular reactivity was evaluated by inducing an ischemic state in the foot with a blood pressure cuff during dynamic BOLD imaging. Perfusion indices were assessed in two regions of the foot: the medial plantar (MP) and lateral plantar (LP) regions.

STATISTICAL TESTS

Effect sizes of group mean differences were assessed using Hedge's g adjusted for small sample sizes.

RESULTS

DFU participants exhibited elevated D*, f, and MVF values in both regions (g ≥ 1.10) and increased D (g = 1.07) in the MP region compared to DC participants. DC participants showed reduced f and MVF compared to HC participants in the MP region (g ≥ 1.06). Finally, the DFU group showed reduced tolerance for ischemia in the LP region (g = -1.51) and blunted reperfusion response in both regions (g < -2.32) compared to the DC group during the cuff-occlusion challenge.

DATA CONCLUSION

The combined use of IVIM and BOLD MRI shows promise in differentiating perfusion abnormalities in the feet of diabetic patients and suggests hyperperfusion in DFU patients.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY: Stage 1.

Assessment of fitting methods and variability of IVIM parameters in muscles of the lumbar spine at rest

Intravoxel incoherent motion (IVIM) MRI provides insight into tissue diffusion and perfusion. Here, estimates of perfusion fraction ( f ), pseudo-diffusion coefficient (D * ), and diffusion coefficient ( D ) obtained via different fitting methods are compared to ascertain (1) the optimal analysis strategy for muscles of the lumbar spine and (2) repeatability of IVIM parameters in skeletal muscle at rest. Diffusion-weighted images were acquired in the lumbar spine at rest in 15 healthy participants. Data were fit to the bi-exponential IVIM model to estimate f , D * and D using three variably segmented approaches based on non-linear least squares fitting, and a Bayesian fitting method. Assuming that perfusion and diffusion are temporally stable in skeletal muscle at rest, and spatially uniform within a spinal segment, the optimal analysis strategy was determined as the approach with the lowest temporal or spatial variation and smallest residual between measured and fit data. Inter-session repeatability of IVIM parameters was evaluated in a subset of 11 people. Finally, simulated IVIM signal at varying signal to noise ratio were evaluated to understand precision and bias. Experimental results showed that IVIM parameter values differed depending on the fitting method. A three-step non-linear least squares fitting approach, where D , f , andD * were estimated sequentially, generally yielded the lowest spatial and temporal variation. Solving all parameters simultaneously yielded the lowest residual between measured and fit data, however there was substantial spatial and temporal variability. Results obtained by Bayesian fitting had high spatial and temporal variability in addition to a large residual between measured and fit data. Simulations showed that all fitting methods can fit the IVIM data at signal to noise ratios >35, and thatD * was the most challenging to accurately obtain. Overall, this study motivates use of a three-step non-linear least squares fitting strategy to quantify IVIM parameters in skeletal muscle.

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.

Effect of compression garments on muscle perfusion in delayed-onset muscle soreness: A quantitative analysis using intravoxel incoherent motion MR perfusion imaging

The aim of this prospective cohort study was to evaluate the effect of compression garments under resting conditions and after the induction of delayed-onset muscle soreness (DOMS) by MR perfusion imaging using intravoxel incoherent motion (IVIM). Magnetic resonance imaging of both lower legs of 16 volunteers was performed before and after standardized eccentric exercises that induced DOMS. A compression garment (21-22 mmHg) was worn during and for 6 h after exercise on one randomly selected leg. IVIM MR imaging, represented as total muscle perfusion D*f, perfusion fraction f and tissue diffusivity D, were compared between baseline and directly, 30 min, 6 h and 48 h after exhausting exercise with and without compression. Creatine kinase levels and T2-weighted images were acquired at baseline and after 48 h. DOMS was induced in the medial head of the gastrocnemius muscle (MGM) in all volunteers. Compression garments did not show any significant effect on IVIM perfusion parameters at any time point in the MGM or the tibialis anterior muscle (p > 0.05). Microvascular perfusion in the MGM increased significantly in both the compressed and noncompressed leg between baseline measurements and those taken directly after and 30 min after the exercise: the relative median f increased by 31.5% and 24.7% in the compressed and noncompressed leg, respectively, directly after the exercise compared with the baseline value. No significant change in tissue perfusion occurred 48 h after the induction of DOMS compared with baseline. It was concluded that compression garments (21-22 mmHg) do not alter microvascular muscle perfusion at rest, nor do they have any significant effect during the regeneration phase of DOMS. In DOMS, only a short-term effect of increased muscle perfusion (30 min after exercise) was observed, with normalization occurring during regeneration after 6-48 h. The normalization of perfusion independently of compression after 6 h may have implications for diagnostic and therapeutic strategies and for the better understanding of pathophysiological pathways in DOMS.

Parsimonious modeling of skeletal muscle perfusion: Connecting the stretched exponential and fractional Fickian diffusion

PURPOSE

To develop an anomalous (non-Gaussian) diffusion model for characterizing skeletal muscle perfusion using multi-b-value DWI.

THEORY AND METHODS

Fick's first law was extended for describing tissue perfusion as anomalous superdiffusion, which is non-Gaussian diffusion exhibiting greater particle spread than that of the Gaussian case. This was accomplished using a space-fractional derivative that gives rise to a power-law relationship between mean squared displacement and time, and produces a stretched exponential signal decay as a function of b-value. Numerical simulations were used to estimate parameter errors under in vivo conditions, and examine the effect of limited SNR and residual fat signal. Stretched exponential DWI parameters, α and D , were measured in thigh muscles of 4 healthy volunteers at rest and following in-magnet exercise. These parameters were related to a stable distribution of jump-length probabilities and used to estimate microvascular volume fractions.

RESULTS

Numerical simulations showed low dispersion in parameter estimates within 1.5% and 1%, and bias errors within 3% and 10%, for α and D , respectively. Superdiffusion was observed in resting muscle, and to a greater degree following exercise. Resting microvascular volume fraction was between 0.0067 and 0.0139 and increased between 2.2-fold and 4.7-fold following exercise.

CONCLUSIONS

This model captures superdiffusive molecular motions consistent with perfusion, using a parsimonious representation of the DWI signal, providing approximations of microvascular volume fraction comparable with histological estimates. This signal model demonstrates low parameter-estimation errors, and therefore holds potential for a wide range of applications in skeletal muscle and elsewhere in the body.

Diffusion-weighted MRI with intravoxel incoherent motion modeling for assessment of muscle perfusion in the thigh during post-exercise hyperemia in younger and older adults

Aging is associated with impaired endothelium-dependent vasodilation that leads to muscle perfusion impairment and contributes to organ dysfunction. Impaired muscle perfusion may result in inadequate delivery of oxygen and nutrients during and after muscle contraction, leading to muscle damage. The ability to study the relationship between perfusion and muscle damage has been limited using traditional muscle perfusion measures, which are invasive and risky. To overcome this limitation, we optimized a diffusion-weighted MRI sequence and validated an intravoxel incoherent motion (IVIM) analysis based on Monte Carlo simulation to study muscle perfusion impairment with aging during post-exercise hyperemia. Simulation results demonstrated that the bias of IVIM-derived perfusion fraction (f_p ) and diffusion of water molecules in extra-vascular tissue (D) ranged from -3.3% to 14% and from -16.5% to 0.002%, respectively, in the optimized experimental condition. The dispersion in f_p and D ranged from 3.2% to 9.5% and from 0.9% to 1.1%, respectively. The mid-thigh of the left leg of four younger (21-30 year old) and four older (60-90 year old) healthy females was studied using the optimized protocol at baseline and at seven time increments occurring every 3.25 min following in-magnet dynamic knee extension exercise performed using a MR-compatible ergometer with a workload of 0.4 bar for 2.5 min. After exercise, both f_p and D significantly increased in the rectus femoris (active muscle during exercise) but not in adductor magnus (inactive muscle), reflecting the fact that the local increase in perfusion with both groups showed a maximum value in the second post-exercise time-point. A significantly greater increase in perfusion from the baseline (p < 0.05) was observed in the younger group (37 ± 12.05%) compared with the older group (17.57 ± 15.92%) at the first post-exercise measurement. This work establishes a reliable non-invasive method that can be used to study the effects of aging on dynamic changes in muscle perfusion as they relate to important measures of physical function.

The Role of Muscle Perfusion in the Age-Associated Decline of Mitochondrial Function in Healthy Individuals

Maximum oxidative capacity of skeletal muscle measured by in vivo phosphorus magnetic resonance spectroscopy (³¹P-MRS) declines with age, and negatively affects whole-body aerobic capacity. However, it remains unclear whether the loss of oxidative capacity is caused by reduced volume and function of mitochondria or limited substrate availability secondary to impaired muscle perfusion. Therefore, we sought to elucidate the role of muscle perfusion on the age-related decline of muscle oxidative capacity and ultimately whole-body aerobic capacity. Muscle oxidative capacity was assessed by ³¹P-MRS post-exercise phosphocreatine recovery time (τ_PCr), with higher τ_PCr reflecting lower oxidative capacity, in 75 healthy participants (48 men, 22–89 years) of the Genetic and Epigenetic Signatures of Translational Aging Laboratory Testing study. Muscle perfusion was characterized as an index of blood volume at rest using a customized diffusion-weighted MRI technique and analysis method developed in our laboratory. Aerobic capacity (peak-VO₂) was also measured during a graded treadmill exercise test in the same visit. Muscle oxidative capacity, peak-VO₂, and resting muscle perfusion were significantly lower at older ages independent of sex, race, and body mass index (BMI). τ_PCr was significantly associated with resting muscle perfusion independent of age, sex, race, and BMI (p-value = 0.004, β = −0.34). τ_PCr was also a significant independent predictor of peak-VO₂ and, in a mediation analysis, significantly attenuated the association between muscle perfusion and peak-VO₂ (34% reduction for β in perfusion). These findings suggest that the age-associated decline in muscle oxidative capacity is partly due to impaired muscle perfusion and not mitochondrial dysfunction alone. Furthermore, our findings show that part of the decline in whole-body aerobic capacity observed with aging is also due to reduced microvascular blood volume at rest, representing a basal capacity of the microvascular system, which is mediated by muscle oxidative capacity. This finding suggests potential benefit of interventions that target an overall increase in muscle perfusion for the restoration of energetic capacity and mitochondrial function with aging.

A simple and fast adaptive nonlocal multispectral filtering algorithm for efficient noise reduction in magnetic resonance imaging

PURPOSE

We recently introduced a multispectral (MS) nonlocal (NL) filter based on maximum likelihood estimation (MLE) of voxel intensities, termed MS-NLML. While MS-NLML provides excellent noise reduction and improved image feature preservation as compared to other NL or MS filters, it requires considerable processing time, limiting its application in routine analyses. In this work, we introduced a fast, simple, and robust filter, termed nonlocal estimation of multispectral magnitudes (NESMA), for noise reduction in multispectral (MS) magnetic resonance imaging (MRI).

METHODS

Through extensive simulation and in-vivo analyses, we compared the performance of NESMA and MS-NLML in terms of noise reduction and processing efficiency. Further, we introduce two simple adaptive methods that permit spatial variation of similar voxels, R, used in the filtering. The first method is semi-adaptive and permits variation of R across the image by using a relative Euclidean distance (RED) similarity threshold. The second method is fully adaptive and filters the raw data with several RED similarity thresholds to spatially determine the optimal threshold value using an unbiased criterion.

RESULTS

NESMA shows very similar filtering performance as compared to MS-NLML, however, with much simple implementation and very fast processing time. Further, for both filters, the adaptive methods were shown to further reduce noise in comparison with the conventional non-adaptive method in which R is set to a constant value throughout the image.

CONCLUSIONS

NESMA is fast, robust, and straightforward to implement filter. These features render it suitable for routine clinical use and analysis of large MRI datasets.

Use of the NESMA Filter to Improve Myelin Water Fraction Mapping with Brain MRI

BACKGROUND AND PURPOSE

Myelin water fraction (MWF) mapping permits direct visualization of myelination patterns in the developing brain and in pathology. MWF is conventionally measured through multiexponential T₂ analysis which is very sensitive to noise, leading to inaccuracies in derived MWF estimates. Although noise reduction filters may be applied during postprocessing, conventional filtering can introduce bias and obscure small structures and edges. Advanced nonblurring filters, while effective, exhibit a high level of complexity and the requirement for supervised implementation for optimal performance. The purpose of this paper is to demonstrate the ability of the recently introduced nonlocal estimation of multispectral magnitudes (NESMA) filter to greatly improve the determination of MWF parameter estimates from gradient and spin echo (GRASE) imaging data.

METHODS

We evaluated the performance of the NESMA filter for MWF mapping from clinical GRASE imaging data of the human brain, and compared the results to those calculated from unfiltered images. Numerical and in vivo analyses of the brains of three subjects, representing different ages, were conducted.

RESULTS

Our results demonstrated the potential of the NESMA filter to permit high-quality in vivo MWF mapping. Indeed, NESMA permits substantial reduction of random variation in derived MWF estimates while preserving accuracy and detail.

CONCLUSIONS

In vivo estimation of MWF in the human brain from GRASE imaging data was markedly improved through use of the NESMA filter. The use of NESMA may contribute to the goal of high-quality MWF mapping in clinically feasible imaging times.