Improved reconstruction of crossing fibers in the mouse optic pathways with orientation distribution function fingerprinting.

PURPOSE: The accuracy of diffusion MRI tractography reconstruction decreases in the white matter regions with crossing fibers. The optic pathways in rodents provide a challenging structure to test new diffusion tractography approaches because of the small crossing volume within the optic chiasm and the unbalanced 9:1 proportion between the contra- and ipsilateral neural projections from the retina to the lateral geniculate nucleus, respectively. METHODS: Common approaches based on Orientation Distribution Function (ODF) peak finding or statistical inference were compared qualitatively and quantitatively to ODF Fingerprinting (ODF-FP) for reconstruction of crossing fibers within the optic chiasm using in vivo diffusion MRI ( n = 18 $$ n=18 $$ healthy C57BL/6 mice). Manganese-Enhanced MRI (MEMRI) was obtained after intravitreal injection of manganese chloride and used as a reference standard for the optic pathway anatomy. RESULTS: ODF-FP outperformed by over 100% all the tested methods in terms of the ratios between the contra- and ipsilateral segments of the reconstructed optic pathways as well as the spatial overlap between tractography and MEMRI. CONCLUSION: In this challenging model system, ODF-Fingerprinting reduced uncertainty of diffusion tractography for complex structural formations of fiber bundles.

Tractography passes the test: Results from the diffusion-simulated connectivity (disco) challenge.

Estimating structural connectivity from diffusion-weighted magnetic resonance imaging is a challenging task, partly due to the presence of false-positive connections and the misestimation of connection weights. Building on previous efforts, the MICCAI-CDMRI Diffusion-Simulated Connectivity (DiSCo) challenge was carried out to evaluate state-of-the-art connectivity methods using novel large-scale numerical phantoms. The diffusion signal for the phantoms was obtained from Monte Carlo simulations. The results of the challenge suggest that methods selected by the 14 teams participating in the challenge can provide high correlations between estimated and ground-truth connectivity weights, in complex numerical environments. Additionally, the methods used by the participating teams were able to accurately identify the binary connectivity of the numerical dataset. However, specific false positive and false negative connections were consistently estimated across all methods. Although the challenge dataset doesn't capture the complexity of a real brain, it provided unique data with known macrostructure and microstructure ground-truth properties to facilitate the development of connectivity estimation methods.

Time-dependent diffusivity and kurtosis in phantoms and patients with head and neck cancer.

PURPOSE: To assess the reliability of measuring diffusivity, diffusional kurtosis, and cellular-interstitial water exchange time with long diffusion times (100-800 ms) using stimulated-echo DWI. METHODS: Time-dependent diffusion MRI was tested on two well-established diffusion phantoms and in 5 patients with head and neck cancer. Measurements were conducted using an in-house diffusion-weighted STEAM-EPI pulse sequence with multiple diffusion times at a fixed TE on three scanners. We used the weighted linear least-squares fit method to estimate time-dependent diffusivity, D ( t ) $$ D(t) $$ , and diffusional kurtosis, K ( t ) $$ K(t) $$ . Additionally, the Kärger model was used to estimate cellular-interstitial water exchange time ( τ ex $$ {\tau}_{ex} $$ ) from K ( t ) $$ K(t) $$ . RESULTS: Diffusivity measured by time-dependent STEAM-EPI measurements and commercial SE-EPI showed comparable results with R2 above 0.98 and overall 5.4 ± 3.0% deviation across diffusion times. Diffusional kurtosis phantom data showed expected patterns: constant D $$ D $$ and K $$ K $$  = 0 for negative controls and slow varying D $$ D $$ and K $$ K $$ for samples made of nanoscopic vesicles. Time-dependent diffusion MRI in patients with head and neck cancer found that the Kärger model could be considered valid in 72% ± 23% of the voxels in the metastatic lymph nodes. The median cellular-interstitial water exchange time estimated for lesions was between 58.5 ms and 70.6 ms. CONCLUSIONS: Based on two well-established diffusion phantoms, we found that time-dependent diffusion MRI measurements can provide stable diffusion and kurtosis values over a wide range of diffusion times and across multiple MRI systems. Moreover, estimation of cellular-interstitial water exchange time can be achieved using the Kärger model for the metastatic lymph nodes in patients with head and neck cancer.

Reproducibility of the Standard Model of diffusion in white matter on clinical MRI systems.

Estimating intra- and extra-axonal microstructure parameters, such as volume fractions and diffusivities, has been one of the major efforts in brain microstructure imaging with MRI. The Standard Model (SM) of diffusion in white matter has unified various modeling approaches based on impermeable narrow cylinders embedded in locally anisotropic extra-axonal space. However, estimating the SM parameters from a set of conventional diffusion MRI (dMRI) measurements is ill-conditioned. Multidimensional dMRI helps resolve the estimation degeneracies, but there remains a need for clinically feasible acquisitions that yield robust parameter maps. Here we find optimal multidimensional protocols by minimizing the mean-squared error of machine learning-based SM parameter estimates for two 3T scanners with corresponding gradient strengths of 40and80mT/m. We assess intra-scanner and inter-scanner repeatability for 15-minute optimal protocols by scanning 20 healthy volunteers twice on both scanners. The coefficients of variation all SM parameters except free water fraction are ≲10% voxelwise and 1-4% for their region-averaged values. As the achieved SM reproducibility outcomes are similar to those of conventional diffusion tensor imaging, our results enable robust in vivo mapping of white matter microstructure in neuroscience research and in the clinic.

Insights from the IronTract challenge: Optimal methods for mapping brain pathways from multi-shell diffusion MRI.

Limitations in the accuracy of brain pathways reconstructed by diffusion MRI (dMRI) tractography have received considerable attention. While the technical advances spearheaded by the Human Connectome Project (HCP) led to significant improvements in dMRI data quality, it remains unclear how these data should be analyzed to maximize tractography accuracy. Over a period of two years, we have engaged the dMRI community in the IronTract Challenge, which aims to answer this question by leveraging a unique dataset. Macaque brains that have received both tracer injections and ex vivo dMRI at high spatial and angular resolution allow a comprehensive, quantitative assessment of tractography accuracy on state-of-the-art dMRI acquisition schemes. We find that, when analysis methods are carefully optimized, the HCP scheme can achieve similar accuracy as a more time-consuming, Cartesian-grid scheme. Importantly, we show that simple pre- and post-processing strategies can improve the accuracy and robustness of many tractography methods. Finally, we find that fiber configurations that go beyond crossing (e.g., fanning, branching) are the most challenging for tractography. The IronTract Challenge remains open and we hope that it can serve as a valuable validation tool for both users and developers of dMRI analysis methods.

Performance of orientation distribution function-fingerprinting with a biophysical multicompartment diffusion model.

PURPOSE: Orientation Distribution Function (ODF) peak finding methods typically fail to reconstruct fibers crossing at shallow angles below 40°, leading to errors in tractography. ODF-Fingerprinting (ODF-FP) with the biophysical multicompartment diffusion model allows for breaking this barrier. METHODS: A randomized mechanism to generate a multidimensional ODF-dictionary that covers biologically plausible ranges of intra- and extra-axonal diffusivities and fraction volumes is introduced. This enables ODF-FP to address the high variability of brain tissue. The performance of the proposed approach is evaluated on both numerical simulations and a reconstruction of major fascicles from high- and low-resolution in vivo diffusion images. RESULTS: ODF-FP with the suggested modifications correctly identifies fibers crossing at angles as shallow as 10 degrees in the simulated data. In vivo, our approach reaches 56% of true positives in determining fiber directions, resulting in visibly more accurate reconstruction of pyramidal tracts, arcuate fasciculus, and optic radiations than the state-of-the-art techniques. Moreover, the estimated diffusivity values and fraction volumes in corpus callosum conform with the values reported in the literature. CONCLUSION: The modified ODF-FP outperforms commonly used fiber reconstruction methods at shallow angles, which improves deterministic tractography outcomes of major fascicles. In addition, the proposed approach allows for linearization of the microstructure parameters fitting problem.

Three-Dimensional Printed Anatomic Models Derived From Magnetic Resonance Imaging Data: Current State and Image Acquisition Recommendations for Appropriate Clinical Scenarios.

Three-dimensional (3D) printing technologies have been increasingly utilized in medicine over the past several years and can greatly facilitate surgical planning thereby improving patient outcomes. Although still much less utilized compared to computed tomography (CT), magnetic resonance imaging (MRI) is gaining traction in medical 3D printing. The purpose of this study was two-fold: 1) to determine the prevalence in the existing literature of using MRI to create 3D printed anatomic models for surgical planning and 2) to provide image acquisition recommendations for appropriate clinical scenarios where MRI is the most suitable imaging modality. The workflow for creating 3D printed anatomic models from medical imaging data is complex and involves image segmentation of the regions of interest and conversion of that data into 3D surface meshes, which are compatible with printing technologies. CT is most commonly used to create 3D printed anatomic models due to the high image quality and relative ease of performing image segmentation from CT data. As compared to CT datasets, 3D printing using MRI data offers advantages since it provides exquisite soft tissue contrast needed for accurate organ segmentation and it does not expose patients to unnecessary ionizing radiation. MRI, however, often requires complicated imaging techniques and time-consuming postprocessing procedures to generate high-resolution 3D anatomic models needed for 3D printing. Despite these challenges, 3D modeling and printing from MRI data holds great clinical promises thanks to emerging innovations in both advanced MRI imaging and postprocessing techniques. EVIDENCE LEVEL: 2 TECHNICAL EFFICATCY: 5.

T1 and T2 quantification using magnetic resonance fingerprinting in mild traumatic brain injury.

OBJECTIVES: To assess whether MR fingerprinting (MRF)-based relaxation properties exhibit cross-sectional and prospective correlations with patient outcome and compare the results with those from DTI. METHODS: Clinical imaging, MRF, and DTI were acquired in patients (24 ± 10 days after injury (timepoint 1) and 90 ± 17 days after injury (timepoint 2)) and once in controls. Patient outcome was assessed with global functioning, symptom profile, and neuropsychological testing. ADC and fractional anisotropy (FA) from DTI and T1 and T2 from MRF were compared in 12 gray and white matter regions with Mann-Whitney tests. Bivariate associations between MR measures and outcome were assessed using the Spearman correlation and logistic regression. RESULTS: Data from 22 patients (38 ± 12 years; 17 women) and 18 controls (32 ± 8 years; 12 women) were analyzed. Fourteen patients (37 ± 12 years; 11 women) returned for timepoint 2, while two patients provided only timepoint 2 clinical outcome data. At timepoint 1, there were no differences between patients and controls in T1, T2, and ADC, while FA was lower in mTBI frontal white matter. T1 at timepoint 1 and the change in T1 exhibited more (n = 18) moderate to strong correlations (|r|= 0.6-0.85) with clinical outcome at timepoint 2 than T2 (n = 3), FA (n = 7), and ADC (n = 2). High T1 at timepoint 1, and serially increasing T1, accounted for five of the six MR measures with the highest utility for identification of non-recovered patients at timepoint 2 (AUC > 0.80). CONCLUSION: T1 derived from MRF was found to have higher utility than T2, FA, and ADC for predicting 3-month outcome after mTBI. KEY POINTS: • In a region-of-interest approach, FA, ADC, and T1 and T2 all showed limited utility in differentiating patients from controls at an average of 24 and 90 days post-mild traumatic brain injury. • T1 at 24 days, and the serial change in T1, revealed more and stronger predictive correlations with clinical outcome at 90 days than did T2, ADC, or FA. • T1 showed better prospective identification of non-recovered patients at 90 days than ADC, T2, and FA.

CG-SENSE revisited: Results from the first ISMRM reproducibility challenge.

PURPOSE: The aim of this work is to shed light on the issue of reproducibility in MR image reconstruction in the context of a challenge. Participants had to recreate the results of "Advances in sensitivity encoding with arbitrary k-space trajectories" by Pruessmann et al. METHODS: The task of the challenge was to reconstruct radially acquired multicoil k-space data (brain/heart) following the method in the original paper, reproducing its key figures. Results were compared to consolidated reference implementations created after the challenge, accounting for the two most common programming languages used in the submissions (Matlab/Python). RESULTS: Visually, differences between submissions were small. Pixel-wise differences originated from image orientation, assumed field-of-view, or resolution. The reference implementations were in good agreement, both visually and in terms of image similarity metrics. DISCUSSION AND CONCLUSION: While the description level of the published algorithm enabled participants to reproduce CG-SENSE in general, details of the implementation varied, for example, density compensation or Tikhonov regularization. Implicit assumptions about the data lead to further differences, emphasizing the importance of sufficient metadata accompanying open datasets. Defining reproducibility quantitatively turned out to be nontrivial for this image reconstruction challenge, in the absence of ground-truth results. Typical similarity measures like NMSE of SSIM were misled by image intensity scaling and outlier pixels. Thus, to facilitate reproducibility, researchers are encouraged to publish code and data alongside the original paper. Future methodological papers on MR image reconstruction might benefit from the consolidated reference implementations of CG-SENSE presented here, as a benchmark for methods comparison.

Fingerprinting Orientation Distribution Functions in diffusion MRI detects smaller crossing angles.

Diffusion tractography is routinely used to study white matter architecture and brain connectivity in vivo. A key step for successful tractography of neuronal tracts is the correct identification of tract directions in each voxel. Here we propose a fingerprinting-based methodology to identify these fiber directions in Orientation Distribution Functions, dubbed ODF-Fingerprinting (ODF-FP). In ODF-FP, fiber configurations are selected based on the similarity between measured ODFs and elements in a pre-computed library. In noisy ODFs, the library matching algorithm penalizes the more complex fiber configurations. ODF simulations and analysis of bootstrapped partial and whole-brain in vivo datasets show that the ODF-FP approach improves the detection of fiber pairs with small crossing angles while maintaining fiber direction precision, which leads to better tractography results. Rather than focusing on the ODF maxima, the ODF-FP approach uses the whole ODF shape to infer fiber directions to improve the detection of fiber bundles with small crossing angle. The resulting fiber directions aid tractography algorithms in accurately displaying neuronal tracts and calculating brain connectivity.

Low Rank plus Sparse decomposition of ODFs for improved detection of group-level differences and variable correlations in white matter.

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,227cm3, agrees with and expands on the volume found by TBSS (17 cm3), Connectivity based fixel enhancement (15 cm3) and Connectometry (212 cm3). 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.

Accelerated radial diffusion spectrum imaging using a multi-echo stimulated echo diffusion sequence.

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 79:306-316, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Validation of surface-to-volume ratio measurements derived from oscillating gradient spin echo on a clinical scanner using anisotropic fiber phantoms.

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(ω), 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(ω) 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 (ρ ≥ 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.

Radial q-space sampling for DSI.

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 76:769-780, 2016. © 2015 Wiley Periodicals, Inc.

Evaluation of breast cancer using intravoxel incoherent motion (IVIM) histogram analysis: comparison with malignant status, histological subtype, and molecular prognostic factors.

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.

Comparison of fitting methods and b-value sampling strategies for intravoxel incoherent motion in breast cancer.

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.

Dynamic diffusion-tensor measurements in muscle tissue using the single-line multiple-echo diffusion-tensor acquisition technique at 3T.

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.

A model-based reconstruction for undersampled radial spin-echo DTI with variational penalties on the diffusion tensor.

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.

Time-dependent diffusion in skeletal muscle with the random permeable barrier model (RPBM): application to normal controls and chronic exertional compartment syndrome patients.

The purpose of this work was to carry out diffusion tensor imaging (DTI) at multiple diffusion times Td in skeletal muscle in normal subjects and chronic exertional compartment syndrome (CECS) patients and analyze the data with the random permeable barrier model (RPBM) for biophysical specificity. Using an institutional review board approved HIPAA-compliant protocol, seven patients with clinical suspicion of CECS and eight healthy volunteers underwent DTI of the calf muscle in a Siemens MAGNETOM Verio 3 T scanner at rest and after treadmill exertion at four different T(d) values. Radial diffusion values λ(rad) were computed for each of seven different muscle compartments and analyzed with RPBM to produce estimates of free diffusivity D(0), fiber diameter a, and permeability κ. Fiber diameter estimates were compared with measurements from literature autopsy reference for several compartments. Response factors (post/pre-exercise ratios) were computed and compared between normal controls and CECS patients using a mixed-model two-way analysis of variance. All subjects and muscle compartments showed nearly time-independent diffusion along and strongly time-dependent diffusion transverse to the muscle fibers. RPBM estimates of fiber diameter correlated well with corresponding autopsy reference. D(0) showed significant (p < 0.05) increases with exercise for volunteers, and a increased significantly (p < 0.05) in volunteers. At the group level, response factors of all three parameters showed trends differentiating controls from CECS patients, with patients showing smaller diameter changes (p = 0.07), and larger permeability increases (p = 0.07) than controls. Time-dependent diffusion measurements combined with appropriate tissue modeling can provide enhanced microstructural specificity for in vivo tissue characterization. In CECS patients, our results suggest that high-pressure interfiber edema elevates free diffusion and restricts exercise-induced fiber dilation. Such specificity may be useful in differentiating CECS from other disorders or in predicting its response to either physical therapy or fasciotomy.

FastMRI Prostate: A public, biparametric MRI dataset to advance machine learning for prostate cancer imaging.

Magnetic resonance imaging (MRI) has experienced remarkable advancements in the integration of artificial intelligence (AI) for image acquisition and reconstruction. The availability of raw k-space data is crucial for training AI models in such tasks, but public MRI datasets are mostly restricted to DICOM images only. To address this limitation, the fastMRI initiative released brain and knee k-space datasets, which have since seen vigorous use. In May 2023, fastMRI was expanded to include biparametric (T2- and diffusion-weighted) prostate MRI data from a clinical population. Biparametric MRI plays a vital role in the diagnosis and management of prostate cancer. Advances in imaging methods, such as reconstructing under-sampled data from accelerated acquisitions, can improve cost-effectiveness and accessibility of prostate MRI. Raw k-space data, reconstructed images and slice, volume and exam level annotations for likelihood of prostate cancer are provided in this dataset for 47468 slices corresponding to 1560 volumes from 312 patients. This dataset facilitates AI and algorithm development for prostate image reconstruction, with the ultimate goal of enhancing prostate cancer diagnosis.