Skeletal muscle diffusion tensor-MRI fiber tracking: rationale, data acquisition and analysis methods, applications and future directions.

The mechanical functions of muscles involve the generation of force and the actuation of movement by shortening or lengthening under load. These functions are influenced, in part, by the internal arrangement of muscle fibers with respect to the muscle's mechanical line of action. This property is known as muscle architecture. In this review, we describe the use of diffusion tensor (DT)-MRI muscle fiber tracking for the study of muscle architecture. In the first section, the importance of skeletal muscle architecture to function is discussed. In addition, traditional and complementary methods for the assessment of muscle architecture (brightness-mode ultrasound imaging and cadaver analysis) are presented. Next, DT-MRI is introduced and the structural basis for the reduced and anisotropic diffusion of water in muscle is discussed. The third section discusses issues related to the acquisition of skeletal muscle DT-MRI data and presents recommendations for optimal strategies. The fourth section discusses methods for the pre-processing of DT-MRI data, the available approaches for the calculation of the diffusion tensor and the seeding and propagating of fiber tracts, and the analysis of the tracking results to measure structural properties pertinent to muscle biomechanics. Lastly, examples are presented of how DT-MRI fiber tracking has been used to provide new insights into how muscles function, and important future research directions are highlighted. Copyright © 2016 John Wiley & Sons, Ltd.

Multi-parametric MRI characterization of healthy human thigh muscles at 3.0 T - relaxation, magnetization transfer, fat/water, and diffusion tensor imaging.

Muscle diseases commonly have clinical presentations of inflammation, fat infiltration, fibrosis, and atrophy. However, the results of existing laboratory tests and clinical presentations are not well correlated. Advanced quantitative MRI techniques may allow the assessment of myo-pathological changes in a sensitive and objective manner. To progress towards this goal, an array of quantitative MRI protocols was implemented for human thigh muscles; their reproducibility was assessed; and the statistical relationships among parameters were determined. These quantitative methods included fat/water imaging, multiple spin-echo T2 imaging (with and without fat signal suppression, FS), selective inversion recovery for T1 and quantitative magnetization transfer (qMT) imaging (with and without FS), and diffusion tensor imaging. Data were acquired at 3.0 T from nine healthy subjects. To assess the repeatability of each method, the subjects were re-imaged an average of 35 days later. Pre-testing lifestyle restrictions were applied to standardize physiological conditions across scans. Strong between-day intra-class correlations were observed in all quantitative indices except for the macromolecular-to-free water pool size ratio (PSR) with FS, a metric derived from qMT data. Two-way analysis of variance revealed no significant between-day differences in the mean values for any parameter estimate. The repeatability was further assessed with Bland-Altman plots, and low repeatability coefficients were obtained for all parameters. Among-muscle differences in the quantitative MRI indices and inter-class correlations among the parameters were identified. There were inverse relationships between fractional anisotropy (FA) and the second eigenvalue, the third eigenvalue, and the standard deviation of the first eigenvector. The FA was positively related to the PSR, while the other diffusion indices were inversely related to the PSR. These findings support the use of these T1 , T2 , fat/water, and DTI protocols for characterizing skeletal muscle using MRI. Moreover, the data support the existence of a common biophysical mechanism, water content, as a source of variation in these parameters.

Preservation of biological function despite oxidative modification of the apolipoprotein A-I mimetic peptide 4F.

Myeloperoxidase (MPO)-derived hypochlorous acid induces changes in HDL function via redox modifications at the level of apolipoprotein A-I (apoA-I). As 4F and apoA-I share structural and functional properties, we tested the hypothesis that 4F acts as a reactive substrate for hypochlorous acid (HOCl). 4F reduced the HOCl-mediated oxidation of the fluorescent substrate APF in a concentration-dependent manner (ED(50) ∼ 56 ± 3 µM). This reaction induced changes in the physical properties of 4F. Addition of HOCl to 4F at molar ratios ranging from 1:1 to 3:1 reduced 4F band intensity on SDS-PAGE gels and was accompanied by the formation of a higher molecular weight species. Chromatographic studies showed a reduction in 4F peak area with increasing HOCl and the formation of new products. Mass spectral analyses of collected fractions revealed oxidation of the sole tryptophan (Trp) residue in 4F. 4F was equally susceptible to oxidation in the lipid-free and lipid-bound states. To determine whether Trp oxidation influenced its apoA-I mimetic properties, we monitored effects of HOCl on 4F-mediated lipid binding and ABCA1-dependent cholesterol efflux. Neither property was altered by HOCl. These results suggest that 4F serves as a reactive substrate for HOCl, an antioxidant response that does not influence the lipid binding and cholesterol effluxing capacities of the peptide.

Characterizing Intracranial Hemodynamics in Sickle Cell Anemia: Impact of Patient-Specific Viscosity.

PURPOSE: Pediatric and adult patients with sickle cell anemia (SCA) are at increased risk of stroke and cerebrovascular accident. In the general adult population, there is a relationship between arterial hemodynamics and pathology; however, this relationship in SCA patients remains to be elucidated. The aim of this work was to characterize circle of Willis hemodynamics in patients with SCA and quantify the impact of viscosity choice on pathophysiologically-relevant hemodynamics measures. METHODS: Based on measured vascular geometries, time-varying flow rates, and blood parameters, detailed patient-specific simulations of the circle of Willis were conducted for SCA patients (n = 6). Simulations quantified the impact of patient-specific and standard blood viscosities on wall shear stress (WSS). RESULTS: These results demonstrated that use of a standard blood viscosity introduces large errors into the estimation of pathophysiologically-relevant hemodynamic parameters. Standard viscosity models overpredicted peak WSS by 55% and 49% for steady and pulsatile flow, respectively. Moreover, these results demonstrated non-uniform, spatial patterns of positive and negative WSS errors related to viscosity, and standard viscosity simulations overpredicted the time-averaged WSS by 32% (standard deviation = 7.1%). Finally, differences in shear rate demonstrated that the viscosity choice alters the simulated near-wall flow field, impacting hemodynamics measures. CONCLUSIONS: This work presents simulations of circle of Willis arterial flow in SCA patients and demonstrates the importance and feasibility of using a patient-specific viscosity in these simulations. Accurately characterizing cerebrovascular hemodynamics in SCA populations has potential for elucidating the pathophysiology of large-vessel occlusion, aneurysms, and tissue damage in these patients.

Combined In Silico and In Vitro Approach Predicts Low Wall Shear Stress Regions in a Hemofilter that Correlate with Thrombus Formation In Vivo.

A major challenge in developing blood-contacting medical devices is mitigating thrombogenicity of an intravascular device. Thrombi may interfere with device function or embolize from the device to occlude distant vascular beds with catastrophic consequences. Chemical interactions between plasma proteins and bioengineered surface occur at the nanometer scale; however, continuum models of blood predict local shear stresses that lead to platelet activation or aggregation and thrombosis. Here, an iterative approach to blood flow path design incorporating in silico, in vitro, and in vivo experiments predicted the occurrence and location of thrombi in an implantable hemofilter. Low wall shear stress (WSS) regions identified by computational fluid dynamics (CFD) predicted clot formation in vivo. Revised designs based on CFD demonstrated superior performance, illustrating the importance of a multipronged approach for a successful design process.

Original article submission: Platelet stress accumulation analysis to predict thrombogenicity of an artificial kidney.

An implantable artificial kidney using a hemofilter constructed from an array of silicon membranes to provide ultrafiltration requires a suitable blood flow path to ensure stable operation in vivo. Two types of flow paths distributing blood to the array of membranes were evaluated: parallel and serpentine. Computational fluid dynamics (CFD) simulations were used to guide the development of the blood flow paths. Pressure data from animal tests were used to obtain pulsatile flow conditions imposed in the transient simulations. A key consideration for stable operation in vivo is limiting platelet stress accumulation to avoid platelet activation and thrombus formation. Platelet stress exposure was evaluated by CFD particle tracking methods through the devices to provide distributions of platelet stress accumulation. The distributions of stress accumulation over the duration of a platelet lifetime for each device revealed that stress accumulation for the serpentine flow path exceeded levels expected to cause platelet activation while the accumulated stress for the parallel flow path was below expected activation levels.

Quantitative Magnetic Resonance Imaging of Skeletal Muscle Disease.

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

Anisotropic Smoothing Improves DT-MRI-Based Muscle Fiber Tractography.

PURPOSE: To assess the effect of anisotropic smoothing on fiber tracking measures, including pennation angle, fiber tract length, and fiber tract number in the medial gastrocnemius (MG) muscle in healthy subjects using diffusion-weighted magnetic resonance imaging (DW-MRI). MATERIALS AND METHODS: 3T DW-MRI data were used for muscle fiber tractography in the MG of healthy subjects. Anisotropic smoothing was applied at three levels (5%, 10%, 15%), and pennation angle, tract length, fiber tract number, fractional anisotropy, and principal eigenvector orientation were quantified for each smoothing level. RESULTS: Fiber tract length increased with pre-fiber tracking smoothing, and local heterogeneities in fiber direction were reduced. However, pennation angle was not affected by smoothing. CONCLUSION: Modest anisotropic smoothing (10%) improved fiber-tracking results, while preserving structural features.

Matching of postcontraction perfusion to oxygen consumption across submaximal contraction intensities in exercising humans.

Studying the magnitude and kinetics of blood flow, oxygen extraction, and oxygen consumption at exercise onset and during the recovery from exercise can lead to insights into both the normal control of metabolism and blood flow and the disturbances to these processes in metabolic and cardiovascular diseases. The purpose of this study was to examine the on- and off-kinetics for oxygen delivery, extraction, and consumption as functions of submaximal contraction intensity. Eight healthy subjects performed four 1-min isometric dorsiflexion contractions, with two at 20% MVC and two at 40% MVC. During one contraction at each intensity, relative perfusion changes were measured by using arterial spin labeling, and the deoxyhemoglobin percentage (%HHb) was estimated using the spin- and gradient-echo sequence and a previously published empirical calibration. For the whole group, the mean perfusion did not increase during contraction. The %HHb increased from ∼28 to 38% during contractions of each intensity, with kinetics well described by an exponential function and mean response times (MRTs) of 22.7 and 21.6 s for 20 and 40% MVC, respectively. Following contraction, perfusion increased ∼2.5-fold. The %HHb, oxygen consumption, and perfusion returned to precontraction levels with MRTs of 27.5, 46.4, and 50.0 s, respectively (20% MVC), and 29.2, 75.3, and 86.0 s, respectively (40% MVC). These data demonstrate in human subjects the varied recovery rates of perfusion and oxygen consumption, along with the similar rates of %HHb recovery, across these exercise intensities.

Diffusion-Tensor MRI Based Skeletal Muscle Fiber Tracking.

A skeletal muscle's function is strongly influenced by the internal organization and geometric properties of its fibers, a property known as muscle architecture. Diffusion-tensor magnetic resonance imaging-based fiber tracking provides a powerful tool for non-invasive muscle architecture studies, has three-dimensional sensitivity, and uses a fixed frame of reference. Significant advances have been made in muscle fiber tracking technology, including defining seed points for fiber tracking, quantitatively characterizing muscle architecture, implementing denoising procedures, and testing validity and repeatability. Some examples exist of how these data can be integrated with those from other advanced MRI and computational methods to provide novel insights into muscle function. Perspectives are offered regarding future directions in muscle diffusion-tensor imaging, including needs to develop an improved understanding for the microstructural basis for reduced and anisotropic diffusion, establish the best practices for data acquisition and analysis, and integrate fiber tracking with other physiological data.