Faculty Bibliography

OBJECTIVES/OBJECTIVE:We describe measurement of skeletal muscle kinetics with multiple echo diffusion tensor imaging (MEDITI). This approach allows characterization of the microstructural dynamics in healthy and pathologic muscle. MATERIALS AND METHODS/METHODS:In a Siemens 3-T Skyra scanner, MEDITI was used to collect dynamic DTI with a combination of rapid diffusion encoding, radial imaging, and compressed sensing reconstruction in a multi-compartment agarose gel rotation phantom and within in vivo calf muscle. An MR-compatible ergometer (Ergospect Trispect) was employed to enable in-scanner plantar flexion exercise. In a HIPAA-compliant study with written informed consent, post-exercise recovery of DTI metrics was quantified in eight volunteers. Exercise response of DTI metrics was compared with that of T2-weighted imaging and characterized by a gamma variate model. RESULTS: = 0.303 ± 0.185). Diffusion and T2-weighted response magnitudes were correlated (e.g., r = 0.792, p = 0.019 for nMD vs. nT2w). CONCLUSION/CONCLUSIONS:We have demonstrated the feasibility of MEDITI for capturing spatially resolved diffusion tensor data in dynamic systems including post-exercise skeletal muscle recovery following in-scanner plantar flexion.

A novel approach is presented for group statistical analysis of diffusion weighted MRI datasets through voxelwise Orientation Distribution Functions (ODF). Recent advances in MRI acquisition make it possible to use high quality diffusion weighted protocols (multi-shell, large number of gradient directions) for routine in vivo study of white matter architecture. The dimensionality of these data sets is however often reduced to simplify statistical analysis. While these approaches may detect large group differences, they do not fully capitalize on all acquired image volumes. Incorporation of all available diffusion information in the analysis however risks biasing the outcome by outliers. Here we propose a statistical analysis method operating on the ODF, either the diffusion ODF or fiber ODF. To avoid outlier bias and reliably detect voxelwise group differences and correlations with demographic or behavioral variables, we apply the Low-Rank plus Sparse (L+S) matrix decomposition on the voxelwise ODFs which separates the sparse individual variability in the sparse matrix S whilst recovering the essential ODF features in the low-rank matrix L. We demonstrate the performance of this ODF L+S approach by replicating the established negative association between global white matter integrity and physical obesity in the Human Connectome dataset. The volume of positive findings (p<0.01, 227 cm³) agrees with and expands on the volume found by TBSS (17 cm³), Connectivity based fixel enhancement (15 cm³) and Connectometry (212 cm³). In the same dataset we further localize the correlations of brain structure with neurocognitive measures such as fluid intelligence and episodic memory. The presented ODF L+S approach will aid in the full utilization of all acquired diffusion weightings leading to the detection of smaller group differences in clinically relevant settings as well as in neuroscience applications.

PURPOSE: Diffusion spectrum imaging (DSI) provides us non-invasively and robustly with anatomical details of brain microstructure. To achieve sufficient angular resolution, DSI requires a large number of q-space samples, leading to long acquisition times. This need is mitigated here by combining the beneficial properties of Radial q-space sampling for DSI with a Multi-Echo Stimulated Echo Sequence (MESTIM). METHODS: Full 2D k-spaces for each of several q-space samples, along the same radially outward line in q-space, are acquired in one readout train with one spin and three stimulated echoes. RF flip angles are carefully chosen to distribute spin magnetization over the echoes and the DSI reconstruction is adapted to account for differences in diffusion time among echoes. RESULTS: Individual datasets and bootstrapped reproducibility analysis demonstrate image quality and SNR of the more-than-twofold-accelerated RDSI MESTIM sequence. Orientation distribution functions (ODF) and tractography results benefit from the longer diffusion times of the latter echoes in the echo train. CONCLUSION: A MESTIM sequence can be used to shorten RDSI acquisition times significantly without loss of image or ODF quality. Further acceleration is possible by combination with simultaneous multi-slice techniques. Magn Reson Med, 2017. (c) 2017 International Society for Magnetic Resonance in Medicine.

A diffusion measurement in the short-time surface-to-volume ratio (S/V) limit (Mitra et al., Phys Rev Lett. 1992;68:3555) can disentangle the free diffusion coefficient from geometric restrictions to diffusion. Biophysical parameters, such as the S/V of tissue membranes, can be used to estimate microscopic length scales non-invasively. However, due to gradient strength limitations on clinical MRI scanners, pulsed gradient spin echo (PGSE) measurements are impractical for probing the S/V limit. To achieve this limit on clinical systems, an oscillating gradient spin echo (OGSE) sequence was developed. Two phantoms containing 10 fiber bundles, each consisting of impermeable aligned fibers with different packing densities, were constructed to achieve a range of S/V values. The frequency-dependent diffusion coefficient, D(omega), was measured in each fiber bundle using OGSE with different gradient waveforms (cosine, stretched cosine, and trapezoidal), while D(t) was measured from PGSE and stimulated-echo measurements. The S/V values derived from the universal high-frequency behavior of D(omega) were compared against those derived from quantitative proton density measurements using single spin echo (SE) with varying echo times, and from magnetic resonance fingerprinting (MRF). S/V estimates derived from different OGSE waveforms were similar and demonstrated excellent correlation with both SE- and MRF-derived S/V measures (rho >/= 0.99). Furthermore, there was a smoother transition between OGSE frequency f and PGSE diffusion time when using teffS/V=9/64f, rather than the commonly used teff = 1/(4f), validating the specific frequency/diffusion time conversion for this regime. Our well-characterized fiber phantom can be used for the calibration of OGSE and diffusion modeling techniques, as the S/V ratio can be measured independently using other MR modalities. Moreover, our calibration experiment offers an exciting perspective of mapping tissue S/V on clinical systems.

PURPOSE: Diffusion spectrum imaging (DSI) has been shown to be an effective tool for noninvasively depicting the anatomical details of brain microstructure. Existing implementations of DSI sample the diffusion encoding space using a rectangular grid. Here we present a different implementation of DSI whereby a radially symmetric q-space sampling scheme for DSI is used to improve the angular resolution and accuracy of the reconstructed orientation distribution functions. METHODS: Q-space is sampled by acquiring several q-space samples along a number of radial lines. Each of these radial lines in q-space is analytically connected to a value of the orientation distribution functions at the same angular location by the Fourier slice theorem. RESULTS: Computer simulations and in vivo brain results demonstrate that radial diffusion spectrum imaging correctly estimates the orientation distribution functions when moderately high b-values (4000 s/mm2) and number of q-space samples (236) are used. CONCLUSION: The nominal angular resolution of radial diffusion spectrum imaging depends on the number of radial lines used in the sampling scheme, and only weakly on the maximum b-value. In addition, the radial analytical reconstruction reduces truncation artifacts which affect Cartesian reconstructions. Hence, a radial acquisition of q-space can be favorable for DSI. Magn Reson Med, 2015. (c) 2015 Wiley Periodicals, Inc.

PURPOSE: To examine heterogeneous breast cancer through intravoxel incoherent motion (IVIM) histogram analysis. MATERIALS AND METHODS: This HIPAA-compliant, IRB-approved retrospective study included 62 patients (age 48.44 +/- 11.14 years, 50 malignant lesions and 12 benign) who underwent contrast-enhanced 3 T breast MRI and diffusion-weighted imaging. Apparent diffusion coefficient (ADC) and IVIM biomarkers of tissue diffusivity (Dt), perfusion fraction (fp), and pseudo-diffusivity (Dp) were calculated using voxel-based analysis for the whole lesion volume. Histogram analysis was performed to quantify tumour heterogeneity. Comparisons were made using Mann-Whitney tests between benign/malignant status, histological subtype, and molecular prognostic factor status while Spearman's rank correlation was used to characterize the association between imaging biomarkers and prognostic factor expression. RESULTS: The average values of the ADC and IVIM biomarkers, Dt and fp, showed significant differences between benign and malignant lesions. Additional significant differences were found in the histogram parameters among tumour subtypes and molecular prognostic factor status. IVIM histogram metrics, particularly fp and Dp, showed significant correlation with hormonal factor expression. CONCLUSION: Advanced diffusion imaging biomarkers show relationships with molecular prognostic factors and breast cancer malignancy. This analysis reveals novel diagnostic metrics that may explain some of the observed variability in treatment response among breast cancer patients. KEY POINTS: * Novel IVIM biomarkers characterize heterogeneous breast cancer. * Histogram analysis enables quantification of tumour heterogeneity. * IVIM biomarkers show relationships with breast cancer malignancy and molecular prognostic factors.

PURPOSE: To compare fitting methods and sampling strategies, including the implementation of an optimized b-value selection for improved estimation of intravoxel incoherent motion (IVIM) parameters in breast cancer. METHODS: Fourteen patients (age, 48.4 +/- 14.27 years) with cancerous lesions underwent 3 Tesla breast MRI examination for a HIPAA-compliant, institutional review board approved diffusion MR study. IVIM biomarkers were calculated using "free" versus "segmented" fitting for conventional or optimized (repetitions of key b-values) b-value selection. Monte Carlo simulations were performed over a range of IVIM parameters to evaluate methods of analysis. Relative bias values, relative error, and coefficients of variation (CV) were obtained for assessment of methods. Statistical paired t-tests were used for comparison of experimental mean values and errors from each fitting and sampling method. RESULTS: Comparison of the different analysis/sampling methods in simulations and experiments showed that the "segmented" analysis and the optimized method have higher precision and accuracy, in general, compared with "free" fitting of conventional sampling when considering all parameters. Regarding relative bias, IVIM parameters fp and Dt differed significantly between "segmented" and "free" fitting methods. CONCLUSION: IVIM analysis may improve using optimized selection and "segmented" analysis, potentially enabling better differentiation of breast cancer subtypes and monitoring of treatment. Magn Reson Med, 2014. (c) 2014 Wiley Periodicals, Inc.

OBJECTIVE: To measure background parenchymal enhancement (BPE) and compare with other contrast enhancement values and diffusion-weighted MRI parameters in healthy and cancerous breast tissue at the clinical level. MATERIALS AND METHODS: This HIPAA-compliant, IRB approved retrospective study enrolled 77 patients (38 patients with breast cancer - mean age 51.8+/-10.0 years; 39 high-risk patients for screening evaluation - mean age 46.3+/-11.7 years), who underwent contrast-enhanced 3T breast MRI. Contrast enhanced MRI and diffusion-weighted imaging were performed to quantify BPE, lesion contrast enhancement, and apparent diffusion coefficient (ADC) metrics in fibroglandular tissue (FGT) and lesions. RESULTS: BPE did not correlate with ADC values. Mean BPE for the lesion-bearing patients was higher (43.9%) compared to that of the high-risk screening patients (28.3%, p=0.004). Significant correlation (r=0.37, p<0.05) was found between BPE and lesion contrast enhancement. CONCLUSION: No significant association was observed between parenchymal or lesion enhancement with conventional apparent diffusion metrics, suggesting that proliferative processes are not co-regulated in cancerous and parenchymal tissue.

When diffusion biomarkers display transient changes, i.e. in muscle following exercise, traditional diffusion-tensor imaging (DTI) methods lack the temporal resolution to resolve the dynamics. This article presents an MRI method for dynamic diffusion-tensor acquisitions on a clinical 3T scanner. This method, the Single-Line Multiple-Echo Diffusion-Tensor Acquisition Technique (SL-MEDITATE), achieves a high temporal resolution (4 s) by rapid diffusion encoding through the acquisition of multiple echoes with unique diffusion sensitization and limiting the readout to a single line volume. The method is demonstrated in a rotating anisotropic phantom, a flow phantom with adjustable flow speed and in vivo skeletal calf muscle of healthy volunteers following a plantar flexion exercise. The rotating and flow-varying phantom experiments show that SL-MEDITATE correctly identifies the rotation of the first diffusion eigenvector and the changes in diffusion-tensor parameter magnitudes, respectively. Immediately following exercise, the in vivo mean diffusivity (MD) time courses show, before the well-known increase, an initial decrease that is not typically observed in traditional DTI. In conclusion, SL-MEDITATE can be used to capture transient changes in tissue anisotropy in a single line. Future progress might allow for dynamic DTI when combined with appropriate k-space trajectories and compressed sensing reconstruction

Radial spin-echo diffusion imaging allows motion-robust imaging of tissues with very low T2 values like articular cartilage with high spatial resolution and signal-to-noise ratio (SNR). However, in vivo measurements are challenging, due to the significantly slower data acquisition speed of spin-echo sequences and the less efficient k-space coverage of radial sampling, which raises the demand for accelerated protocols by means of undersampling. This work introduces a new reconstruction approach for undersampled diffusion-tensor imaging (DTI). A model-based reconstruction implicitly exploits redundancies in the diffusion-weighted images by reducing the number of unknowns in the optimization problem and compressed sensing is performed directly in the target quantitative domain by imposing a total variation (TV) constraint on the elements of the diffusion tensor. Experiments were performed for an anisotropic phantom and the knee and brain of healthy volunteers (three and two volunteers, respectively). Evaluation of the new approach was conducted by comparing the results with reconstructions performed with gridding, combined parallel imaging and compressed sensing and a recently proposed model-based approach. The experiments demonstrated improvements in terms of reduction of noise and streaking artifacts in the quantitative parameter maps, as well as a reduction of angular dispersion of the primary eigenvector when using the proposed method, without introducing systematic errors into the maps. This may enable an essential reduction of the acquisition time in radial spin-echo diffusion-tensor imaging without degrading parameter quantification and/or SNR